摘要:本文主要向大家介绍了 C/C++知识点之佩尔方程,通过具体的内容向大家展示,希望对大家学习C/C++知识点有所帮助。
本文主要向大家介绍了 C/C++知识点之佩尔方程,通过具体的内容向大家展示,希望对大家学习C/C++知识点有所帮助。
1.连分数
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
using namespace std;
typedef long long ll;
ll a[20000];
bool pell_minimum_solution(ll n,ll &x0,ll &y0){
ll m=(ll)sqrt((double)n);
double sq=sqrt(n);
int i=0;
if(m*m==n)return false;//当n是完全平方数则佩尔方程无解
a[i++]=m;
ll b=m,c=1;
double tmp;
do{
c=(n-b*b)/c;
tmp=(sq+b)/c;
a[i++]=(ll)(floor(tmp));
b=a[i-1]*c-b;
//printf("%lld %lld %lld\n",a[i-1],b,c);
}while(a[i-1]!=2*a[0]);
ll p=1,q=0;
for(int j=i-2;j>=0;j--){
ll t=p;
p=q+p*a[j];
q=t;
//printf("a[%d]=%lld %lld %lld\n",j,a[j],p,q);
}
if((i-1)%2==0){x0=p;y0=q;}
else{x0=2*p*p+1;y0=2*p*q;}
return true;
}
int main(){
ll n,x,y;
while(~scanf("%lld",&n)){
if(pell_minimum_solution(n,x,y)){
printf("%lld^2-%lld*%lld^2=1\t",x,n,y);
printf("%lld-%lld=1\n",x*x,n*y*y);
}
}
2.暴力
递推式
x[n]=x[n-1]*x[1]+d*y[n-1]*y[1];
y[n]=x[n-1]*y[1]+y[n-1]*x[1];
本文由职坐标整理并发布,希望对同学们有所帮助。了解更多详情请关注职坐标编程语言C/C+频道!
您输入的评论内容中包含违禁敏感词
我知道了
请输入正确的手机号码
请输入正确的验证码
您今天的短信下发次数太多了,明天再试试吧!
我们会在第一时间安排职业规划师联系您!
您也可以联系我们的职业规划师咨询:
版权所有 职坐标-一站式IT培训就业服务领导者 沪ICP备13042190号-4
上海海同信息科技有限公司 Copyright ©2015 www.zhizuobiao.com,All Rights Reserved.
沪公网安备 31011502005948号