[NOIP2003]麦森数

  题目描述

形如2^p-1的素数称为麦森数，这时p一定也是个素数。但反过来不一定，即如果p是个素数，2^p-1不一定也是素数。 到1998年底，人们已找到了37个麦森数。最大的一个是p=3021377，它有909526位。 麦森数有许多重要应用，它与完全数密切相关。 现要求输入p(1000＜ P ＜ 3100000)，计算2^p-1的位数和最后500位数字(用十进制高精度数表示)。

   输入

第1行：1个整数p(1000

输出

第1行：十进制高精度数2^p-1的位数； 第2..11行：十进制高精度数2^p-1的最后500位数字(每行输出50位，共输出10行，不足500位时高位补0)； 注意：不必验证2^p-1与p是否为素数。

样例输入

Copy (如果复制到控制台无换行，可以先粘贴到文本编辑器，再复制)

1279

样例输出

386

00000000000000000000000000000000000000000000000000

00000000000000000000000000000000000000000000000000

38615262247266704805319112350403608059673360298012

00000000000000104079321946643990819252403273640855

23944173232418484242161395428100779138356624832346

50138443826029173234888531116082853841658502825560

49081399066056773207629241295093892203457731833496 61583550472959420547689811211693677147548478866962 46662248318909188018470682222031405210266984354887

32958028878050869736186900714720710555703168729087

    分析： 重定义运算符+快速幂

#include
#include
#include
#include
#include
#define LL long long
using namespace std;
//#define BIG 100000LL
const LL BIG=100000LL; 
struct INT{
    LL data[500];
    int len;
    INT(){memset(data,0,sizeof(data));len=0;}
    void operator=(int x){
        while(x){
            data[len++]=x%BIG;
            x/=BIG;
        }
    }
}ans,tmp;
int p;
INT operator*(const INT &a,const INT &b){
    INT ret; 
    int i, j, add; 
    for(i=0;i>=1;
        tmp=tmp*tmp;
    }
}
int main(){
    int i;
    scanf("%d",&p);
    printf("%d\n",int(log10(2)*p)+1);
    work();
    ans.data[0]--;
    for(i=99;i>=0;i--){
        printf("%05lld",ans.data[i]);
        if(i%10==0)
            printf("\n");
    }
}