POJ 1320 Street Numbers [佩尔方程]
2013-07-11 17:33
405 查看
题目可转述为:求解两个不相等的正整数n,m(n<m),使得 1+2+……+ n = n +(n+1)+ …… + m;输出前十组满足条件的(n,m)(从小到大),每组占一行,每个数占十个格,向右对齐。
分析:要使1+2+……+ n = n +(n+1)+ …… + m,那么n*(n-1)/2=(m-n)(m+n+1)/2,即(2*m+1)^2-8*n^2=1,令 x =2*m+1,y=n;有x^2-8*y^2=1;这是典型的佩尔方程(形如x^2-d*y^2=1的不定方程称为佩尔方程,其中d>1且d不为完全平方数) 已知最小特解为:x1=3,y1=1,由迭代公式有:xn=xn-1 *x1 + d*yn-1*y1;yn=xn-1*y1+yn-1*x1;
那么 xn+1=3*xn+8*yn; yn+1=xn + 3*yn;
代码如下:
分析:要使1+2+……+ n = n +(n+1)+ …… + m,那么n*(n-1)/2=(m-n)(m+n+1)/2,即(2*m+1)^2-8*n^2=1,令 x =2*m+1,y=n;有x^2-8*y^2=1;这是典型的佩尔方程(形如x^2-d*y^2=1的不定方程称为佩尔方程,其中d>1且d不为完全平方数) 已知最小特解为:x1=3,y1=1,由迭代公式有:xn=xn-1 *x1 + d*yn-1*y1;yn=xn-1*y1+yn-1*x1;
那么 xn+1=3*xn+8*yn; yn+1=xn + 3*yn;
代码如下:
#include<cstdio> #include<cmath> #include<cstring> int main() { int x,y,x1,y1,px,py,d; x1=3; y1=1; px=3; py=1; d=8; for(int i=1;i<=10;i++) { x=px*x1+d*py*y1; y=px*y1+py*x1; printf("%10d%10d\n",y,(x-1)/2); px=x; py=y; } return 0; }
相关文章推荐
- POJ 1320 Street Numbers 佩尔方程
- POJ 1320Street Numbers(佩尔方程定理)
- POJ1320 Street Numbers【佩尔方程】
- POJ 1320 Street Numbers 解佩尔方程
- POJ - 1320 - Street Numbers - (佩尔方程)
- POJ1320 Street Numbers(佩尔方程)
- POJ1320 Street Numbers【佩尔方程】
- POJ 1320:Street Numbers——佩尔方程
- Poj1320佩尔方程见数论书115或者110页(暴力求解)
- POJ1320 UVA138 Street Numbers【Pell方程+数学】
- POJ 1320 Street Number(佩尔方程)
- poj 1320 Street Numbers 解pell方程
- poj-1320 Street Number 佩尔方程
- poj 1320 Street Numbers(解pell方程)
- 毕达哥拉斯三元组:poj 1305+佩尔方程:poj 1320,hdu 3292(特殊不定方程)
- poj 1320 佩尔方程递推
- POJ 1320 Street Numbers 笔记
- POJ 1320-Street Numbers(解佩尔方程)
- poj 1320——Street Numbers
- POJ 1320:Street Numbers 佩尔方程