HDU2281(Pell方程的妙用)
来源:互联网 发布:js 字符串加上单引号 编辑:程序博客网 时间:2024/05/08 07:23
HDU2281(Pell方程的妙用)
分类: 数论2013-07-29 10:21 145人阅读 评论(0) 收藏 举报
题目:Square Number
题意:给出一个数N,求1到N的范围内,找到一个最大的n,满足,注意这里的N很大啊,10^18
解析:由于2n+1是奇数,所以3|(2n+1),这样的话2|(n+1),得到n=6k+1,其中k=0,1,2。。。
则n=6k+1<=N,所以:
我们可以得到3k+1与4k+1互素,证明:
现在就是有:
令:
消去k得到:,继续令x=2a,y=b
得到,哈哈,Pell方程,在10^18范围内由于x,y都是10^9级别的,所以实际上处理到Pell方程的第17个
解就可以了。另外这里的x是偶数才行。
然后先求最小解,然后递推式,最后程序中再判断一下就行了。
- #include <iostream>
- #include <string.h>
- #include <stdio.h>
- using namespace std;
- typedef long long LL;
- const int N=25;
- LL x[N],y[N];
- LL tmpx,tmpy,n,k;
- void Init()
- {
- x[1]=2;y[1]=1;
- for(int i=2;i<20;i++)
- {
- x[i]=2*x[i-1]+3*y[i-1];
- y[i]=x[i-1]+2*y[i-1];
- }
- }
- int main()
- {
- Init();
- while(cin>>n)
- {
- if(n==0) break;
- for(int i=17;i>=0;i--)
- {
- if(x[i]%2==0&&x[i]*x[i]<=2*n+2)
- {
- tmpx=x[i]/2;
- tmpy=y[i];
- break;
- }
- }
- k=(tmpx*tmpx-1)/3;
- cout<<6*k+1<<" "<<tmpx*tmpy<<endl;
- }
- return 0;
- }
- HDU2281(Pell方程的妙用)
- 连分数求解Pell方程
- 解Pell方程(java)
- HDU 3292 pell方程
- POJ 2427 pell方程
- HDU 2281 pell方程
- poj1320(解pell方程)
- 连分数求解Pell方程
- 求Pell方程解
- 佩尔(Pell)方程
- MapleSim在模型方程的妙用
- Pell方程解法+连分数自己整理
- Pell Equation (佩尔方程)
- UVa 138 Street Numbers(数论&Pell方程)
- HDU 2281 Square Number (Pell方程&数论)
- poj 1320 Street Numbers(解pell方程)
- poj 1320 Street Numbers 解pell方程
- PE 66 Diophantine equation(Pell方程)
- HDU 1151 二分图
- 面试题33:把字符串转换为整数
- Meta标签详解
- 【Android】两步搞定AndroidSDKManager在线更新SDK版本失败问题
- Ubuntu开启ssh服务
- HDU2281(Pell方程的妙用)
- HDU 1281 二分图
- AFNetworking 学习笔记
- 调试PHP扩展
- 如何遵循函数内部实现规则,提高函数质量?
- 关于RM中的X3014错误,以及mul() 、天空盒
- hdu 1236 排名
- 相似度计算概述
- UVA 10341 Solve It