poj 1320 Street Numbers 解pell方程

题意：

解（2*n+1）^2-8*a^2=1.

分析：

令x=2*n+1,y=a,则转化为x^2-8*y^2=1的pell方程。pell方程详解

代码：

//poj 1320
//sep9
#include <iostream> 
using namespace std;int main()
{int x1=3,y1=1,t=10;while(t--){int x2,y2;x2=3*x1+8*y1;y2=3*y1+x1;printf("%10d%10d\n",y2,(x2-1)/2);x1=x2;y1=y2;}return 0;	
}