poj 1320 Street Numbers （佩尔方程）

题目链接：哆啦A梦传送门

佩尔方程：参考链接：维基

最小解，记作(x1,y1)，则所有的解(xi,yi)由以下的递回关系式得到：

 

各个变量的通项公式：

题意：求满足1+2+3+……+(n-1)=(n+1)+(n+2)+……+m的前10项n和m。

我们化简可得：

那么就可以得结果了。x1=3，y1=1。

故

 

代码：

///佩尔方程应用

#include<cstdio>
#include<algorithm>
#include<cstring>

using namespace std;
typedef long long LL;

int main()
{
    LL px=3,py=1,x1=3,y1=1,n=8;

    for(int i=1;i<=10;i++)
    {
        LL x=3*px+8*py;
        LL y=3*py+px;
        printf("%10lld%10lld\n",y,(x-1)/2);
        px=x;
        py=y;
    }
    return 0;
}