POJ 1320Street Numbers(佩尔方程定理)

题目： LINK

题意：最终是这个意思1+2+3+.....+a == a + a+1 + a+2 + ..... + b;求满足这个式子的(a, b) b由大到小的前10个。
 利用求和公式可以得到 2*a*a == b*b + b; 令b = (y-1)/2,带入得到y*y - 8*a*a = 1.(化成佩尔方程的形式)
可以利用佩尔方程定理，如果D是一个正整数且不是完全平方数
        x*x - D*y*y = 1 总是有正整数解.
    如果(x1, y1)是使x1最小的解，则每个解课通过取幂得到：
    xk + yk*sqrt(D) = (x1 + y1*sqrt(D)) ^k , k = 1, 2, 3, ....
可以知道得到的式子最小解(1, 3) ，可以不断带入求得前10个解

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<iostream>
#include<string>
#include<vector>
#include<cmath>
#include<queue>
#include<map>
#include<set>
using namespace std;
#define INF 1000000000
//typedef __int64 LL;
int main()
{
#ifndef ONLINE_JUDGE
	freopen("in.txt", "r", stdin);
#endif // ONLINE_JUDGE
	int x, y;
	int a, b;
	// y*y - 8*a*a = 1;
	x = 3, y = 1;
	a = 3, b = 1;
	for(int i = 1; i <= 10; i++) {
        int aa = a*x + 8 * b * y;
        int bb = x*b + a*y;
        printf("%10d%10d\n", bb, (aa-1)/2);
        a = aa; b = bb;
	}
	return 0;
}