《第章高精度计算C版课件》由会员分享，可在线阅读，更多相关《第章高精度计算C版课件（31页珍藏版）》请在金锄头文库上搜索。
1、第一章第一章 高精度计算高精度计算利用计算机进行数值计算，有时会遇到这样的问题：有些计算要求精度高，希望计算的数的位数可达几十位甚至几百位，虽然计算机的计算精度也算较高了，但因受到硬件的限制，往往达不到实际问题所要求的精度。我们可以利用程序设计的方法去实现这样的高精度计算。介绍常用的几种高精度计算的方法。高精度计算中需要处理好以下几个问题：(1)数据的接收方法和存贮方法数据的接收和存贮：当输入的数很长时，可采用字符串方式输入，这样可输入数字很长的数，利用字符串函数和操作运算，将每一位数取出，存入数组中。另一种方法是直接用循环加数组方法输入数据。voidinit(inta)/传入一个数组stri
2、ngs;cins;/读入字符串sa0=s.length();/用a0计算字符串s的位数for(i=1;i=10)ci%=10;+ci+1;减法借位：if(aibi)-ai+1;ai+=10;ci=ai-bi;乘法进位：ci+j-1=ai*bj+x+ci+j-1;x=ci+j-1/10;ci+j-1%=10;(4)商和余数的求法商和余数处理：视被除数和除数的位数情况进行处理.3PPT学习交流【例1】高精度加法。输入两个正整数，求它们的和。【分析】输入两个数到两个变量中，然后用赋值语句求它们的和，输出。但是，我们知道，在C+语言中任何数据类型都有一定的表示范围。而当两个被加数很大时，上述算法显然不
3、能求出精确解，因此我们需要寻求另外一种方法。在读小学时，我们做加法都采用竖式方法，如图1。这样，我们方便写出两个整数相加的算法。 8 5 6 + 2 5 5 1 1 1 1 图1 A3 A2 A1+ B3 B2 B1 C4 C3 C2 C1 图2如果我们用数组A、B分别存储加数和被加数，用数组C存储结果。则上例有A1=6，A2=5，A3=8，B1=5，B2=5，B3=2，C4=1，C3=1，C2=1，C1=1，两数相加如图2所示。4PPT学习交流因此，算法描述如下：intc100;voidadd(inta,intb)/a,b,c都为数组，分别存储被加数、加数、结果inti=1,x=0;/x是进
4、位while(i=a数组长度)|(i=b数组的长度)ci=ai+bi+x;/第i位相加并加上次的进位x=ci/10;/向高位进位ci%=10;/存储第i位的值i+;/位置下标变量5PPT学习交流通常，读入的两个整数用可用字符串来存储，程序设计如下：#include#include#includeusingnamespacestd;intmain()chara1100,b1100;inta100,b100,c100,lena,lenb,lenc,i,x;memset(a,0,sizeof(a);memset(b,0,sizeof(b);memset(c,0,sizeof(c);gets(a1);
5、gets(b1);/输入加数与被加数lena=strlen(a1);lenb=strlen(b1);for(i=0;i=lena-1;i+)alena-i=a1i-48;/加数放入a数组for(i=0;i=lenb-1;i+)blenb-i=b1i-48;/加数放入b数组lenc=1;x=0;6PPT学习交流while(lenc=lena|lenc=1;i-)coutci;/输出结果coutendl;return0;7PPT学习交流【例2】高精度减法。输入两个正整数，求它们的差。【算法分析】类似加法，可以用竖式求减法。在做减法运算时，需要注意的是：被减数必须比减数大，同时需要处理借位。高精度减
6、法的参考程序：#include#include#includeusingnamespacestd;intmain()inta256,b256,c256,lena,lenb,lenc,i;charn256,n1256,n2256;memset(a,0,sizeof(a);memset(b,0,sizeof(b);memset(c,0,sizeof(c);8PPT学习交流printf(Inputminuend:);gets(n1);/输入被减数printf(Inputsubtrahend:);gets(n2);/输入减数if(strlen(n1)strlen(n2)|(strlen(n1)=str
7、len(n2)&strcmp(n1,n2)n2时，返回正整数；n1n2时，返回负整数/处理被减数和减数，交换被减数和减数strcpy(n,n1);/将n1数组的值完全赋值给n数组strcpy(n1,n2);strcpy(n2,n);cout-;/交换了减数和被减数，结果为负数lena=strlen(n1);lenb=strlen(n2);for(i=0;i=lena-1;i+)alena-i=int(n1i-0);/被减数放入a数组for(i=0;i=lenb-1;i+)blenb-i=int(n2i-0);/减数放入b数组9PPT学习交流 i=1;while(i=lena|i=lenb)if
8、(ai1)lenc-;/最高位的0不输出for(i=lenc;i=1;i-)coutci;/输出结果coutendl;return0;10PPT学习交流【例3】高精度乘法。输入两个正整数，求它们的积。【算法分析】类似加法，可以用竖式求乘法。在做乘法运算时，同样也有进位，同时对每一位进行乘法运算时，必须进行错位相加，如图3、图4。分析c数组下标的变化规律，可以写出如下关系式：ci=ci+c”i+由此可见，ci跟ai*bj乘积有关，跟上次的进位有关，还跟原ci的值有关，分析下标规律，有ci+j-1=ai*bj+x+ci+j-1;x=ci+j-1/10;ci+j-1%=10; 8 5 6 2 5 4
9、 2 8 0 1 7 1 2 2 1 4 0 0 图3A 3 A 2 A 1 B 2 B 1 C4C3 C2 C1 C”5C”4C”3C”2 C 6 C 5 C 4 C 3 C 2 C 1 图411PPT学习交流高精度乘法的参考程序：#include#include#includeusingnamespacestd;intmain()chara1100,b1100;inta100,b100,c100,lena,lenb,lenc,i,j,x;memset(a,0,sizeof(a);memset(b,0,sizeof(b);memset(c,0,sizeof(c);gets(a1);gets(
10、b1);lena=strlen(a1);lenb=strlen(b1);for(i=0;i=lena-1;i+)alena-i=a1i-48;for(i=0;i=lenb-1;i+)blenb-i=b1i-48;12PPT学习交流for(i=1;i=lena;i+)x=0;/用于存放进位for(j=1;j1)/删除前导0lenc-;for(i=lenc;i=1;i-)coutci;coutendl;return0;13PPT学习交流【例4】高精度除法。输入两个正整数，求它们的商（做整除）。【算法分析】做除法时，每一次上商的值都在，每次求得的余数连接以后的若干位得到新的被除数，继续做除法。因此，
11、在做高精度除法时，要涉及到乘法运算和减法运算，还有移位处理。当然，为了程序简洁，可以避免高精度除法，用09次循环减法取代得到商的值。这里，我们讨论一下高精度数除以单精度数的结果，采取的方法是按位相除法。14PPT学习交流#include#include#includeusingnamespacestd;intmain()chara1100,c1100;inta100,c100,lena,i,x=0,lenc,b;memset(a,0,sizeof(a);memset(c,0,sizeof(c);gets(a1);cinb;lena=strlen(a1);for(i=0;i=lena-1;i+)
12、ai+1=a1i-48;15PPT学习交流for(i=1;i=lena;i+)/按位相除ci=(x*10+ai)/b;x=(x*10+ai)%b;lenc=1;while(clenc=0&lenclena)lenc+;/删除前导0for(i=lenc;i=lena;i+)coutci;coutendl;return0;16PPT学习交流实质上，在做两个高精度数运算时候，存储高精度数的数组元素可以不仅仅只保留一个数字，而采取保留多位数（例如一个整型或长整型数据等），这样，在做运算（特别是乘法运算）时，可以减少很多操作次数。例如图5就是采用4位保存的除法运算，其他运算也类似。具体程序可以修改上述例
13、题予以解决，程序请读者完成。示例：12345678945=12345678945=27434841/45=0，1%45=1取12345/45=27412345%45=15取156789/45=3484答案为2743484,余数为156789%45=9图517PPT学习交流【例5】高精除以高精，求它们的商和余数。【算法分析】高精除以低精是对被除数的每一位（这里的“一位”包含前面的余数，以下都是如此）都除以除数，而高精除以高精则是用减法模拟除法，对被除数的每一位都减去除数，一直减到当前位置的数字（包含前面的余数）小于除数（由于每一位的数字小于10，所以对于每一位最多进行10次计算）具体实现程序如下
14、：18PPT学习交流#include#includeusingnamespacestd;inta101,b101,c101,d,i;voidinit(inta)strings;cins;/读入字符串sa0=s.length();/用a0计算字符串s的位数for(i=1;i=a0;i+)ai=sa0-i-0;/将数串s转换为数组a，并倒序存储voidprint(inta)/打印输出if(a0=0)cout00;i-)coutai;coutb则为1，ab0)return1;/a的位数大于b则a比b大if(a00;i-)/从高位到低位比较if(aibi)return1;if(aibi)return-
15、1;return0;/各位都相等则两数相等。voidnumcpy(intp,intq,intdet)/复制p数组到q数组从det开始的地方for(inti=1;i=p0;i+)qi+det-1=pi;q0=p0+det-1;20PPT学习交流voidjian(inta,intb)/计算a=a-bintflag,i;flag=compare(a,b);/调用比较函数判断大小if(flag=0)a0=0;return;/相等if(flag=1)/大于for(i=1;i=a0;i+)if(ai0&aa0=0)a0-;/修正a的位数return;21PPT学习交流voidchugao(inta,int
16、b,intc)inttmp101;c0=a0-b0+1;for(inti=c0;i0;i-)memset(tmp,0,sizeof(tmp);/数组清零numcpy(b,tmp,i);while(compare(a,tmp)=0)ci+;jian(a,tmp);/用减法来模拟while(c00&cc0=0)c0-;return;22PPT学习交流intmain() memset(a,0,sizeof(a); memset(b,0,sizeof(b); memset(c,0,sizeof(c); init(a);init(b); chugao(a,b,c); print(c); print(a)
17、; return0;23PPT学习交流【例6】回文数【问题描述】若一个数（首位不为零）从左向右读与从右向左读都是一样，我们就将其称之为回文数。例如：给定一个10进制数56，将56加65（即把56从右向左读），得到121是一个回文数。又如，对于10进制数87，STEPl：8778=165STEP2：165561=726STEP3：7266271353STEP4：1353+3531=4884在这里的一步是指进行了一次N进制的加法，上例最少用了4步得到回文数4884。写一个程序，给定一个N（2N10或N=16）进制数M求最少经过几步可以得到回文数。如果在30步以内（包含30步）不可能得到回文数，则输
18、出“Impossible”【输入样例】：987【输出样例】：6【算法分析】N进制运算1、当前位规范由%10改为%n2、进位处理由/10改为/n3、其他运算规则不变24PPT学习交流【参考程序】#include#includeusingnamespacestd;intn,a101,b101,ans,i;voidinit(inta)/将数串s转化为整数数组astrings; cinns;/读入字符串s memset(a,0,sizeof(a);/数组a清0 a0=s.length();/用a0计算字符串s的位数 for(i=1;i=0&sa0-i=9)ai=sa0-i-0;elseai=sa0-i
19、-A+10;boolcheck(inta)/判别整数数组a是否为回文数for(i=1;i=a0;i+)if(ai!=aa0-i+1)returnfalse; returntrue;25PPT学习交流voidjia(inta)/整数数组a与其反序数b进行n进制加法运算for(inti=1;i=a0;i+)bi=aa0-i+1;/反序数bfor(inti=1;i=a0;i+)ai+=bi;/逐位相加for(inti=1;i0)a0+;/修正新的a的位数（a+b最多只能的一个进位）intmain() init(a);if(check(a)cout0endl;return0;ans=0;/步数初始化为
20、0while(ans+=30)jia(a);if(check(a)coutansendl;return0;coutImpossible;/输出无解信息return0;26PPT学习交流【上机练习】1、求N！的值【问题描述】用高精度方法，求N！的精确值(N以一般整数输入)。【输入样例】ni.in10【输出样例】ni.out36288002、求A/B高精度值【问题描述】计算A/B的精确值，设A，B是以一般整数输入，计算结果精确小数后20位。【输入样例】ab.in43【输出样例】ab.out4/3=1.33333333333333333333【输入样例】ab.in65【输出样例】ab.out6/5=
21、1.227PPT学习交流3、求n累加和【问题描述】用高精度方法，求s=1+2+3+n的精确值(n以一般整数输入)。【输入样例】ja.in10【输出样例】ja.out554、阶乘和(sum.pas)【问题描述】已知正整数N（N=100）,设S=1!+2!+3!+.N!。其中!表示阶乘,即N!=1*2*3*(N-1)*N，如:3!=1*2*3=6。请编程实现：输入正整数N,输出计算结果S的值。【输入样例】sum.in4【输出样例】sum.out335、高精度求积(MULTIPLY.PAS)【问题描述】输入两个高精度正整数M和N（M和N均小于100位）。【问题求解】求这两个高精度数的积。【输入样例】
22、MULTIPLY.IN363【输出样例】MULTIPLY.OUT10828PPT学习交流6、天使的起誓（YUBIKILI.pas）【问题描述】TENSHI非常幸运的被选为掌管智慧之匙的天使。在正式任职之前，她必须和其他新当选的天使一样，要宣誓。宣誓仪式是每位天使各自表述自己的使命，她们的发言稿被放在N个呈圆形排列的宝盒中。这些宝盒按顺时针方向被编上号码、N-1、N。一开始天使们站在编号为N的宝盒旁。她们各自手上都有一个数字，代表她们自己的发言稿所在的盒子是从号盒子开始按顺时针方向的第几个。例如：有个盒子，那么如果TENSHI手上的数字为，那么她的发言稿所在盒子就是第个。现在天使们开始按照自己手
23、上的数字来找发言稿，先找到的就可以先发言。TENSHI一下子就找到了，于是她最先上台宣誓：“我将带领大家开启NOI之门”TENSHI宣誓结束以后，陆续有天使上台宣誓。可是有一位天使找了好久都找不到她的发言稿，原来她手上的数字M非常大，她转了好久都找不到她想找的宝盒。【问题求解】请帮助这位天使找到她想找的宝盒的编号。【输入格式】从文件YUBIKILI.IN的第一、二行分别读入正整数N和M，其中N、M满足2N108，2M101000【输出格式】把所求宝盒的编号输出到文件YUBIKILI.OUT，文件只有一行（包括换行符）。样例一YUBIKILI.IN79YUBIKILI.OUT2样例二YUBIKI
24、LI.IN11108YUBIKILI.OUT929PPT学习交流7、Hanoi双塔问题(Noip2007)【问题描述】给定A、B、C三根足够长的细柱，在A柱上放有2n个中间有孔的圆盘，共有n个不同的尺寸，每个尺寸都有两个相同的圆盘，注意这两个圆盘是不加区分的（下图为n=3的情形）。现要将这些圆盘移到C柱上，在移动过程中可放在B柱上暂存。要求：（1）每次只能移动一个圆盘；（2）A、B、C三根细柱上的圆盘都要保持上小下大的顺序；任务：设An为2n个圆盘完成上述任务所需的最少移动次数，对于输入的n，输出An。【输入格式】输入文件hanoi.in为一个正整数n，表示在A柱上放有2n个圆盘。【输出格式】输出文件hanoi.out仅一行，包含一个正整数,为完成上述任务所需的最少移动次数An。样例一hanoi.in1hanoi.out2样例二hanoi.inhanoi.out26【限制】对于50%的数据，1=n=25对于100%的数据，1=n=200【提示】设法建立An与An-1的递推关系式。30PPT学习交流此课件下载可自行编辑修改，供参考！感谢您的支持，我们努力做得更好！31PPT学习交流
