麦森数解题报告.docVIP

{先用2^24为一次来优化，再用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:array[1..1001]of longint; procedure fb;{高精度乘法 高*高} var i,j:integer; begin for i:=1 to l do for j:=1 to l do begin b[i+j-1]:=b[i+j-1]+a[i]*a[j]; b[i+j]:=b[i+j]+b[i+j-1] div 10; b[i+j-1]:=b[i+j-1] mod 10; end; if i+j-1=500 then l:=500 else if b[i+j]0 then l:=i+j else l:=i+j-1;{l记录最高位} f:=2{f判断，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