poj 1320 佩尔方程递推

题意：求方程 1+2+3+...+n=(n+1)+(n+2)+...+m 求这个方程的前十个正整数解，分别按格式输出。

思路：化简一下就是 （2m+1）^2 - 8n^2 =1。满足佩尔方程，令x=2m+1 ,y=n,D=8。那么 x(n)=x(n-1)*x1+D*y(n-1)*y1 y(n)=x(n-1)*y1+y(n-1)*x1。最小解是x1=3,y1=1 所以剩下就直接推十个答案出来就行了。

佩尔方程是满足x^2 - a*y^2=1形式的方程。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int res[15][2];
int main(){
int i,j;
res[0][0] = 3;
res[0][1] = 1;
for(i = 1;i<=10;i++){
res[i][0] = res[i-1][0]*res[0][0] + 8*res[i-1][1]*res[0][1];
res[i][1] = res[i-1][0]*res[0][1] + res[i-1][1]*res[0][0];
printf("%10d%10d\n",res[i][1],(res[i][0]-1)>>1);
}
return 0;
}