麦森数解题报告.doc

先用224为一次来优化，再用2(p div 2)*2(p div 2) 高精度 当 p为奇数时 用 2(p div 2)*2(p div 2)*2 var p,q,f,i,j,x,t,l,k,n:longint;s:int64; a,b,c:array1.1001of longint;procedure fb;高精度乘法 高*高 var i,j:integer; begin for i:=1 to l do for j:=1 to l do begin bi+j-1:=bi+j-1+ai*aj; bi+j:=bi+j+bi+j-1 div 10; bi+j-1:=bi+j-1 mod 10; end; if i+j-1=500 then l:=500 else if bi+j0 then l:=i+j else l:=i+j-1;l记录最高位 f:=2f判断，f=2则记录在b数组中 end;procedure fa; var i:integer; begin fb; 2(p div 2)*2(p div 2) for i:=1 to l do begin ci:=ci+bi*2; ci+1:=ci div 10; ci:=ci mod 10; end;高精度*单精度 f:=1;f=1则记录在c数组中 end;procedure f1; begin for i:=1 to (p mod 48) div 2 do24 的余数 begin t:=0; for j:=1 to l do begin aj:=aj*2+t; t:=aj div 10;t：储存aj div 10,并赋值给下一个aj aj:=aj mod 10; end;高精度*单精度 al+1:=t; if (al+10)and(l0 do begin inc(k); x:=ak+t mod 10; t:=t div 10+x div 10; ak:=x mod 10; end;把t的值一位位的放入ak,k=l+1,l+2 if k500 then l:=k else l:=500;l：当l到500时不再往上计算 end; end;beginassign(input,Mason.in);reset(input);assign(output,Mason.out);rewrite(output); read(p); n:=8; n：记总位数 f24; 以224 为一次运算（减少O） f1;以2为一次运算解决 p/24的余数 if odd(p) then fa 做2(p div 2)*2(p div 2)*2else fb; 做2(p div 2)*2(p div 2) 判断p的奇偶 if f=1 then for i:=1 to 500 do bi:=ci;当记录在c数组时，把c赋值给bdec(b1);b1-1i:=1; while (bi100000000 do当s108时，n加一位，s去除末位，以防超出数据结构原先的n=8：108为8位 begin n:=n+1;n：记录总位数 s:=s div 10; end; end; for q:=1 to p mod 11 do一次*2来做剩余的 begin s:=s*2; while s100000000 do begin n:=n+1; s:=s div 10; end; end; write(n);输出总位数 for i:=500 downto 1 do begin if