http://poj.org/problem?id=2427
//* @author
import java.util.*;
import java.math.*;
public class Main {
static BigInteger p1,p2,p3,q1,q2,q3,a1,a2,a0,h1,h2,g1,g2,p,q;
static int nn;
static void solve(){
p2=BigInteger.ONE; p1=BigInteger.ZERO;
q2=BigInteger.ZERO; q1=BigInteger.ONE;
a1=BigInteger.valueOf((long)Math.sqrt(nn));
a0=BigInteger.valueOf((long)Math.sqrt(nn));
g1=BigInteger.ZERO;h1=BigInteger.ONE;
BigInteger n0=BigInteger.valueOf(nn);
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))).compareTo(BigInteger.ONE)==0)
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);
while(in.hasNext()){
nn=in.nextInt();
int t=(int)Math.sqrt(nn);
if(t*t==nn)
System.out.println("No solution!");
else{
solve();
System.out.println(p+" "+q);
}
}
}
}
Meta
-
Recent Posts
Recent Comments
Archives
- May 2024
- April 2023
- February 2023
- January 2023
- December 2022
- November 2022
- September 2022
- June 2022
- July 2021
- January 2021
- February 2020
- September 2019
- March 2018
- February 2018
- August 2016
- July 2016
- June 2016
- May 2016
- April 2016
- March 2016
- February 2016
- January 2016
- December 2015
- November 2015
- October 2015
- September 2015
- August 2015
- July 2015
- June 2015
- May 2015
- April 2015
- March 2015
- February 2015
- January 2015
- December 2014
- November 2014
- October 2014
- September 2014
- August 2014
- July 2014
- June 2014
- May 2014
- April 2014
- March 2014
- February 2014
- January 2014
- December 2013
- November 2013
- October 2013
- September 2013
- August 2013
- July 2013
- June 2013
- May 2013
- April 2013
- March 2013
- February 2013
- January 2013
- December 2012
- November 2012
- October 2012
- September 2012
- August 2012
- July 2012
- June 2012
- May 2012
- April 2012
- March 2012
- February 2012
- January 2012
- December 2011
- November 2011
- October 2011
- September 2011
- August 2011
- July 2011
- June 2011
- May 2011
- April 2011
- March 2011
- February 2011
- January 2011
- December 2010
- November 2010
- October 2010
- September 2010
- August 2010
- July 2010
- June 2010
- May 2010
- April 2010
- March 2010
- February 2010
- January 2010
- December 2009
- November 2009
- October 2009
- September 2009
- August 2009
- July 2009
- June 2009
- May 2009
- April 2009
- March 2009
- February 2009
Categories
