計算PI值到一億位的演算法 (轉)

計算PI值到一億位的演算法 (轉)[@more@]

我大體上考慮了一下用計算PI值到一億位的演算法，得到一個大體的演算法，也好用來交流一下。這是一個構造一種新的長四則運算的演算法。所謂長四則運算，是指用的欄位來作一個小數，用一個記錄來作一個小數數位的演算法。先作出longadd(),longsub(),longmulti(),longdiv(),然後再用求PI公式來求即可。下面我大體上作出了longadd()，並尚未經過上機的。先作六個全域性書籤變數：
public
Xieof,Xpeof,Yieof,Ypeof,Zieof,Zpeof:Tbookmarkstr;
procedure longadd(add1,add2,sum:integer);//欄位的值。
var
add1i,add1p,add2i,add2p,sumi,sump:integer;
add1iend,add1pend,add2iend,add2pend,temp,jump:integer;
begin////
add1i:=add1;add1p:=add1+1;add2i:=add2;add2p:=add2+1;
sumi:=sum;sump:=sum+1;temp:=0;jump:=0;//初值，欄位的索引值。
longtable.bookmark:=xpeof;add1pend:=longtable.recno;
longtable.bookmark:=ypeof;add2pend:=longtable.recno;
if add1pendif add1pend>=add2pend then  longtable.bookmark:=xpeof;
while not bof()  do begin//////##
if (longtable.recnoadd2pend)then  begin
fields(sump):=fields(add2p);longtable.prior;continue;
end;
if (longtable.recnoadd1pend) then  begin
fields(sump):=fields(add1p);longtable.prior;continue;
end;
temp:=fields(add1p)+fields(add2p);
if  (jump=1) then begin
temp:=temp+1;
jump:=0;
end;
if (temp>=10)  then  begin
fields(sump):=temp-10;jump:=1;end;
if temp<10 then  fields(sump):=temp;
longtable.prior;
end;//////##
while not longtable.recno>max(add1iend,add2iend) do begin//@@
if (longtable.recno>min(add1iend,add2iend)) and (add1iendfields(sumi):=fields(add2i);longtable.next;continue;
end;
if (longtable.recno>min(add1iend,add2iend)) and(add1iend>add2iend) then  begin
fields(sumi):=fields(add1i);longtable.next;continue;
end;
temp:=fields(add1i)+fields(add2i);
if  (jump=1) then begin
temp:=temp+1;
jump:=0;
end;
if (temp>=10)  then  begin
fields(sumi):=temp-10;jump:=1;end;
if temp<10 then  fields(sumi):=temp;
longtable.next;
end;//@@

end;////