POJ 1320 Street Numbers [佩尔方程]

题目可转述为：求解两个不相等的正整数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；

代码如下：

#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;}