Pascal语言高精度计算资料

1、高精度计算2.1 高精度加法2.2 高精度减法2.3 高精度乘法2.4 高精度除法高精度加法高精度加法程序如下：programHighPrecision1_Plus;constfn_inp=hp1.inp;fn_out=hp1.out;maxlen=100;maxlengthofthenumbertypehp=recordlen:integer;lengthofthenumbers:array1.maxlenofintegers1isthelowestpositionslenisthehighestpositionend;varx:array1.2ofhp;y:hp;x:input;y:out

2、putprocedurePrintHP(constp:hp);vari:integer;beginfori:=p.lendownto1dowrite(p.si);end;procedureinit;varst:string;j,i:integer;beginassign(input,fn_inp);reset(input);forj:=1to2dobeginreadln(st);xj.len:=length(st);fori:=1toxj.lendochangestringtoHPxj.si:=ord(stxj.len+1-i)-ord(0);end;close(input);end;proc

3、edurePlus(a,b:hp;varc:hp);c:=a+bvari,len:integer;beginfillchar(c,sizeof(c),0);ifa.lenb.lenthenlen:=a.lengetthebiggerlengthofa,belselen:=b.len;fori:=1tolendoplusfromlowtohighbegininc(c.si,a.si+b.si);ifc.si=10thenbegindec(c.si,10);inc(c.si+1);add1toahigherpositionend;end;ifc.slen+10theninc(len);c.len:

4、=len;end;proceduremain;beginPlus(x1,x2,y);end;procedureout;beginassign(output,fn_out);rewrite(output);PrintHP(y);writeln;close(output);end;begininit;main;out;end.2.2 高精度减法高精度减法程序如下:programHighPrecision2_Subtract;constfn_inp=hp2.inp;fn_out=hp2.out;maxlen=100;maxlengthofthenumbertypehp=recordlen:integ

5、er;lengthofthenumbers:array1.maxlenofintegers1isthelowestpositionslenisthehighestpositionend;varx:array1.2ofhp;y:hp;x:input;y:outputpositive:boolean;procedurePrintHP(constp:hp);vari:integer;beginfori:=p.lendownto1dowrite(p.si);end;procedureinit;varst:string;j,i:integer;beginassign(input,fn_inp);rese

6、t(input);forj:=1to2dobeginreadln(st);xj.len:=length(st);fori:=1toxj.lendochangestringtoHPxj.si:=ord(stxj.len+1-i)-ord(0);end;close(input);end;procedureSubtract(a,b:hp;varc:hp);c:=a-b,supposea=bvari,len:integer;beginfillchar(c,sizeof(c),0);ifa.lenb.lenthenlen:=a.lengetthebiggerlengthofa,belselen:=b.l

7、en;fori:=1tolendosubtractfromlowtohighbegininc(c.si,a.si-b.si);ifc.si1)and(c.slen=0)dodec(len);c.len:=len;end;functionCompare(consta,b:hp):integer;ifab0ifa=b-1ifab.lenthenlen:=a.lengetthebiggerlengthofa,belselen:=b.len;while(len0)and(a.slen=b.slen)dodec(len);findapositionwhichhaveadifferentdigitifle

8、n=0thencompare:=0nodifferenceelsecompare:=a.slen-b.slen;end;proceduremain;beginifCompare(x1,x2)=10)dobegininc(c.slen+1,c.slendiv10);c.slen=c.slenmod10;inc(len);end;while(len1)and(c.slen=0)dodec(len);c.len:=len;end;proceduremain;beginMultiply(x,z,y);end;procedureout;beginassign(output,fn_out);rewrite

9、(output);PrintHP(y);writeln;close(output);end;begininit;main;out;end.高精度乘一个整型数据(integer)只需要将上述程序的 hp 类型定义如下即可：typehp=recordlen:integerlengthofthenumbers:array1.maxlenoflongints1isthelowestpositionslenisthehighestpositionend;.高精度乘高精度程序如下：programHighPrecision4_Multiply2;constfn_inp=hp4.inp;fn_out=hp4.

10、out;maxlen=100;maxlengthofthenumbertypehp=recordlen:integer;lengthofthenumbers:array1.maxlenofintegers1isthelowestpositionslenisthehighestpositionend;varx:array1.2ofhp;y:hp;x:input;y:outputprocedurePrintHP(constp:hp);vari:integer;beginfori:=p.lendownto1dowrite(p.si);end;procedureinit;varst:string;j,

11、i:integer;beginassign(input,fn_inp);reset(input);forj:=1to2dobeginreadln(st);xj.len:=length(st);fori:=1toxj.lendochangestringtoHPxj.si:=ord(stxj.len+1-i)-ord(0);end;close(input);end;procedureMultiply(a,b:hp;varc:hp);c:=a+bvari,j,len:integer;beginfillchar(c,sizeof(c),0);fori:=1toa.lendoforj:=1tob.len

12、dobegininc(c.si+j-1,a.si*b.sj);inc(c.si+j,c.si+j-1div10);c.si+j-1:=c.si+j-1mod10;end;len:=a.len+b.len+1;theproductofanumberwithidigitsandanumberwithjdigitscanonlyhaveatmosti+j+1digitswhile(len1)and(c.slen=0)dodec(len);c.len:=len;end;proceduremain;beginMultiply(x1,x2,y);end;procedureout;beginassign(o

13、utput,fn_out);rewrite(output);PrintHP(y);writeln;close(output);end;begininit;main;out;end.高精度除以整型数据(integer);程序如下：programHighPrecision3_Multiply1;constfn_inp=hp5.inp;fn_out=hp5.out;maxlen=100;maxlengthofthenumbertypehp=recordlen:integer;lengthofthenumbers:array1.maxlenofintegers1isthelowestpositions

14、lenisthehighestpositionend;varx,y:hp;z,w:integer;procedurePrintHP(constp:hp);vari:integer;beginfori:=p.lendownto1dowrite(p.si);end;procedureinit;varst:string;i:integer;beginassign(input,fn_inp);reset(input);readln(st);x.len:=length(st);fori:=1tox.lendochangestringtoHPx.si:=ord(stx.len+1-i)-ord(0);re

15、adln(z);close(input);end;procedureDivide(a:hp;b:integer;varc:hp;vard:integer);c:=adivb;d:=amodbvari,len:integer;beginfillchar(c,sizeof(c),0);len:=a.len;d:=0;fori:=lendownto1dofromhightolowbegind:=d*10+a.si;c.si:=ddivb;d:=dmodb;end;while(len1)and(c.slen=0)dodec(len);c.len:=len;end;proceduremain;begin

16、Divide(x,z,y,w);end;procedureout;beginassign(output,fn_out);rewrite(output);PrintHP(y);writeln;writeln(w);close(output);end;begininit;main;out;end.高精度除以高精度程序如下：programHighPrecision4_Multiply2;constfn_inp=hp6.inp;fn_out=hp6.out;maxlen=100;maxlengthofthenumbertypehp=recordlen:integer;lengthofthenumber

17、s:array1.maxlenofintegers1isthelowestpositionslenisthehighestpositionend;varx:array1.2ofhp;y,w:hp;x:input;y:outputprocedurePrintHP(constp:hp);vari:integer;beginfori:=p.lendownto1dowrite(p.si);end;procedureinit;varst:string;j,i:integer;beginassign(input,fn_inp);reset(input);forj:=1to2dobeginreadln(st

18、);xj.len:=length(st);fori:=1toxj.lendochangestringtoHPxj.si:=ord(stxj.len+1-iD-ord(0);end;close(input);end;procedureSubtract(a,b:hp;varc:hp);c:=a-b,supposea=bvari,len:integer;beginfillchar(c,sizeof(c),0);ifa.lenb.lenthenlen:=a.lengetthebiggerlengthofa,belselen:=b.len;fori:=1tolendosubtractfromlowtoh

19、ighbegininc(c.si,a.si 卜 b.si);ifc.si1)and(c.slen=0)dodec(len);c.len:=len;end;functionCompare(consta,b:hp):integer;(ifab0ifa=b-1ifab.lenthenlen:=a.lengetthebiggerlengthofa,belselen:=b.len;while(len0)and(a.slen=b.slen)dodec(len);findapositionwhichhaveadifferentdigitiflen=0thencompare:=0nodifferenceelsecompare:=a.slen-b.slen;end;procedureMultiply10(vara:hp);a:=a*10vari:Integer;beginfor