麦森数

 麦森数（Mason.pas）

【问题描述】形如2^P-1的素数称为麦森数，这时P一定也是个素数。但反过来不一定，即如果P是个素数，2^P-1不一定也是素数。到1998年底，人们已找到了37个麦森数。最大的一个是P=3021377，它有909526位。麦森数有许多重要应用，它与完全数密切相关。

任务：从文件中输入P（1000<P<3100000），计算2^P-1的位数和最后500位数字（用十进制高精度数表示）

 

【输入格式】

文件中只包含一个整数P（1000<P<3100000）

 

【输出格式】

第一行：十进制高精度数2^P-1的位数。

第2-11行：十进制高精度数2^P-1的最后500位数字。（每行输出50位，共输出10行，不足500位时高位补0）

不必验证2^P-1与P是否为素数。

 

【输入样例】

1279

 

【输出样例】

386

00000000000000000000000000000000000000000000000000

00000000000000000000000000000000000000000000000000

00000000000000104079321946643990819252403273640855

38615262247266704805319112350403608059673360298012

23944173232418484242161395428100779138356624832346

49081399066056773207629241295093892203457731833496

61583550472959420547689811211693677147548478866962

50138443826029173234888531116082853841658502825560

46662248318909188018470682222031405210266984354887

32958028878050869736186900714720710555703168729087

==================================================

这道题考察了2个知识点

1.高精度乘高精度

2.二分法..

=============================

var
  p:longint;
  //ans,num:array[0..501]of longint;
  ans,num,x:array[1..1001]of longint;
procedure init;
begin
  assign(input,'p1037.in');
  assign(output,'p1037.out');
  reset(input); rewrite(output);
end;

procedure terminate;
begin
  close(input); close(output);
  halt;
end;

procedure f(p:longint);  //快速幂高精度模块..
var
  i,j:longint;
  t:longint;
begin
  fillchar(num,sizeof(num),0);
  num[1]:=2;
  fillchar(ans,sizeof(ans),0);
  ans[1]:=1;
  
  while p<>0 do
    begin
      t:=p and 1;
      if (p and 1)=1 then
        begin
          //ans*num
          fillchar(x,sizeof(x),0);
          for i:=1 to 500 do
            for j:=1 to 500 do
              x[i+j-1]:=ans[i]*num[j]+x[i+j-1];
          for i:=1 to 500 do
            begin
              x[i+1]:=x[i+1]+x[i] div 10;
              x[i]:=x[i] mod 10;
            end;
          ans:=x;
        end;
      //num*num
      fillchar(x,sizeof(x),0);
      for i:=1 to 500 do
        for j:=1 to 500 do
          x[i+j-1]:=num[i]*num[j]+x[i+j-1];
      for i:=1 to 500 do
        begin
          x[i+1]:=x[i+1]+x[i] div 10;
          x[i]:=x[i] mod 10;
        end;
      num:=x;
      
      p:=p shr 1;
    end;
end;

procedure main;
var
  i:longint;
begin
  readln(p);
  write(trunc(p*ln(2)/ln(10))+1);
  f(p);
  dec(ans[1]);
  i:=1;
  while ans[i]<0 do
    begin
      inc(ans[i],10);
      dec(ans[i+1]);
      inc(i);
    end;
  for i:=500 downto 1 do
    begin
      if i mod 50=0 then writeln;
      write(ans[i]);
    end;
end;

begin
  init;
  main;
  terminate;
end.