POJ 1320 Street Number(佩尔方程)

形如x^2-d*y^2=1(d>1且d不为完全平方数)的不定方程称为佩尔方程。

若佩尔方程的最小特解为(x1,y1)

可以有迭代公式

Xn=X n-1*x1+d*Y n-1*y1

Yn=X n-1*y1+Y n-1*x1

HIT的《数论及应用》上对于此题的描述存在问题,应该是

求解不相等的n,m使得1+2+3+.....+(n-1)=(n+1)+...+m

化简得n(n-1)/2=(n-m)(n+m+1)/2.==>(2*m+1)^2-8*y^2=1(佩尔方程)

这里x1=3,y1=1显然是最小正整数解。

#include <cstdio>
using namespace std;
int main()
{
    int px=3,py=1,d=8,x,y;
    for(int i=1;i<=10;i++)
    {
        x=px*x1+d*y1*py;
        y=py*x1+px*y1;
        printf("%10d%10d\n",y,(x-1)/2);
        px=x;
        py=y;
    }
    return 0;
}



  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值