COJ 1211 大整数开平方

手写求大整数开根号所得到的值，具体计算过程参考别人的资料，最后利用java的大整数得到答案

 别人博客链接：http://www.cnblogs.com/Rinyo/archive/2012/12/16/2820450.html

1.举例

上式意为65536的开平方为256。手开方过程类似于除法计算。为了方便表述，以下仍称类似位置的数为“被除数”、“除数”、“商”。

以65536为例，其具体计算过程如下：

Step1：将被开方数(为了形象，表述成“被除数”，此例中即为65536)从个位往高位每两位一断写成6,55,35的形式，为了方便表述，以下每一个“,”称为一步

Step2：从高位开始计算开方。例如第一步为6，由于2²=4<6<9=3²，因此只能商2(这就是和除法不同的地方，“除数”和“商”的计算位必须相同)。于是将2写在根号上方，计算开方余项。即高位余项加一步低位，此例中，即为高位余项2和低位一步55，余项即为255。

Step3：将Step2得到的第一步开方得数2乘以20(原理在后面证明)作为第二步除数的高位。即本步除数是4x(四十几)。按照要求，本步的商必须是x。因为45×5=225<255<46×6=276，所以本步商5。

Step4：按照类似方法，继续计算以后的各步。其中，每一步的除数高位都是20×已求出的部分商。例如第三步的除数高位就是25×20=500，所以第三步除数为50x。本例中，506×6=3036恰好能整除，所以256就是最终计算结果。

 1 import java.util.*;  2 import java.math.*;  3  4 public class Main {  5 public static void main(String [] args)  6  {  7 Scanner input = new Scanner(System.in);  8  9 while(input.hasNext()) 10  { 11 String str = input.next(); 12 BigInteger num = new BigInteger(str); 13 int len = str.length(); 14 int len1 = (len+1)/2; 15 BigInteger factor = BigInteger.valueOf(20); 16 BigInteger factor1 = BigInteger.valueOf(10); 17 18 BigInteger div = BigInteger.valueOf(0); 19 BigInteger ans = BigInteger.valueOf(0); 20 BigInteger left = BigInteger.valueOf(0); 21 int t = 0 , index = 0; 22 if(len%2 == 1){ 23 t = str.charAt(0)-'0'; 24 index+=1; 25 }else{ 26 t = str.charAt(0)-'0'; 27 t = t*10 + str.charAt(1)-'0'; 28 index+=2; 29  } 30 int m = (int)Math.sqrt(t*1.0+0.1); 31 ans = BigInteger.valueOf(m); 32 div = BigInteger.valueOf(m).multiply(factor); 33 left = BigInteger.valueOf(t-m*m); 34 // System.out.println(ans); 35 BigInteger tmp = BigInteger.valueOf(0); 36 37 for(int i=1 ; i<len1 ; i++ , index+=2){ 38 left = left.multiply(BigInteger.valueOf(100)); 39 int val = str.charAt(index)-'0'; 40 val = val*10+str.charAt(index+1)-'0'; 41 left = left.add(BigInteger.valueOf(val)); 42 // System.out.println("div: "+i+" "+div+" "+left); 43 int j=9; 44 for(; j>=0 ; j--){ 45 tmp = div.add(BigInteger.valueOf(j)); 46 if(left.compareTo(tmp.multiply(BigInteger.valueOf(j)))>=0) break; 47  } 48 left = left.subtract(tmp.multiply(BigInteger.valueOf(j))); 49 ans = ans.multiply(factor1).add(BigInteger.valueOf(j)); 50 div = ans.multiply(factor); 51  } 52  System.out.println(ans); 53  } 54  } 55 }