DSP

高精度练习之超大整数开根

2019-07-13 21:00发布

点击就送屠龙宝刀 题目描述 Description 给出一个正整数n,求n开根号后的整数部分的值。n的位数不超过1000位。
输入描述 Input Description 读入一个不超过1000位的正整数n。
输出描述 Output Description 输出所求答案
样例输入 Sample Input 17
样例输出 Sample Output 4
数据范围及提示 Data Size & Hint n为不超过1000位的正整数 满足 1<=n <10^1000
首先我是一个蒟蒻23333所以我不会一些奇怪的东西比如暴力二分做这个题
then
有一个奇怪的东西叫做手动开根:
手动开根号的基本方法:
1. 将被开方数从右向左每隔2位用撇号分开;
2. 从左边第一段求得算数平方根的第一位数字;
3. 从第一段减去这个第一位数字的平方,再把被开方数的第二段写下来,作为第一个余数;
4. 把所得的第一位数字乘以20,去除第一个余数,所得的商的整数部分作为试商(如果这个整数部分大于或等于10,就改用9左试商,如果第一个余数小于第一位数字乘以20的积,则得试商0);
5. 把第一位数字的20倍加上试商的和,乘以这个试商,如果所得的积大于余数时,就要把试商减1再试,直到积小于或等于余数为止,这个试商就是算数平方根的第二位数字;
6. 用同样方法继续求算数平方根的其他各位数字。 example:
这里写图片描述
so 我们就获得了一种神奇的方法来做这道题切拥有更快的速度和更高的复杂度。 #include #include #include using namespace std; int l; int work(int o,char *O,int I) { char c, *D=O ; if(o>0) { for(l=0;D[l];D[l++]-=10) { D[l++]-=120; D[l]-=110; while(!work(0,O,l)) D[l]+=20; putchar((D[l]+1032)/20); } putchar(10); } else { c=o+(D[I]+82)%10-(I>l/2)*(D[I-l+I]+72)/10-9; D[I]+=I<0 ? 0 : !(o=work(c/10,O,I-1))*((c+999)%10-(D[I]+92)%10); } return o; } int main() { char s[1200];s[0]='0'; scanf("%s",s+1); if(strlen(s)%2 == 1) work(2,s+1,0); else work(2,s,0); return 0; }