注册 登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

那一日,泪水打湿雪花冰冷的心

曾经的滋味,回忆在一次次的离合中,回眸时,那一刻,如涟漪般在一刹那融化

 
 
 

日志

 
 

spoj EQU2  

2010-10-19 22:28:12|  分类: spoj |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
pell方程,x^2-d*y^2=1,用求解连分数
import java.util.*;
import java.math.*;
public class Main {
static BigInteger p,q,p1,p2,p3,q1,q2,q3,a1,a2,a0,h1,h2,g1,g2,n0;
static int n,t;
static void solve(){
p2=BigInteger.ONE; p1=BigInteger.ZERO;
q2=BigInteger.ZERO; q1=BigInteger.ONE;
a0=a1=BigInteger.valueOf((long)Math.sqrt(n));
g1=BigInteger.ZERO;h1=BigInteger.ONE;
n0=BigInteger.valueOf(n);
while(true){
g2=a1.multiply(h1).subtract(g1);
h2=(n0.subtract(g2.multiply(g2))).divide(h1);
a2=(g2.add(a0)).divide(h2);
p=p2.multiply(a1).add(p1);
q=q2.multiply(a1).add(q1);
if(p.multiply(p).subtract(n0.multiply(q.multiply(q))).equals(BigInteger.ONE))
return ;
a1=a2;
g1=g2;
h1=h2;
p1=p2;p2=p;
q1=q2;q2=q;
}
}
public static void main(String[] args) {
Scanner in=new Scanner(System.in);
t=in.nextInt();
for(int i=0;i<t;++i){
n=in.nextInt();
solve();
System.out.println(p+" "+q);
}
}

}
  评论这张
 
阅读(359)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018