COGS41. [NOIP2003] 麦森数

【题目描述】
形如$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)
【输出格式】
第1行：十进制高精度数$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

对于第一问，利用对数$log$的性质可以直接计算位数$n=log_{10}{(2^P)}+1=P*log_{10}{2}+1$
（注，本公式其实把$2^P-1$中的$-1$忽略了，因为它并不影响位数）
（不太熟悉$log$性质的同学请参见百度百科基本性质第二条）
对于第二问，用快速幂+高精乘+半个高精减。

    #include<cstdio>
    #include<math.h>
    #include<string.h>
    #include<iostream>
    using namespace std;
    int P;double n;
    int a[110000],ans[110000];
    void read(int &x)
    {
        int f=1;char t=getchar();x=0;
        while((t<48)or(t>57)){if(t=='-')f=-1;t=getchar();}
        while((t>=48)and(t<=57)){x=x*10+t-48;t=getchar();}
        x*=f;
    }
    void jian(int a[])
    {
        a[1]-=1;
        for (int i=1;i<=a[0];++i)
        {
            if (a[i]<0) {a[i]+=10;--a[i+1];}
            else break;
        }
    }
    void cheng(int a[])
    {
        int len=a[0]+1;
        int c[510];memset(c,0,sizeof(c));
        for (int i=1;i<=a[0];++i)
        {
            int t=0;
            c[i]=c[i]+t+a[i]*2;
            t=c[i]/10;
            c[i]%=10;
            c[i+1]=t;
        }
        while ((c[len]==0)and(len>1)) --len;
        for (int i=1;i<=len;++i) a[i]=c[i];
        a[0]=len;
    }
    void cheng1(int a[],int b[])
    {
        int c[110000];memset(c,0,sizeof(c));
        int n=min(a[0],500),m=min(b[0],500),len=n+m;
        for (int i=1;i<=n;++i)
        {
            int t=0;
            for (int j=1;j<=m;++j)
            {
                c[i+j-1]=c[i+j-1]+t+a[i]*b[j];
                t=c[i+j-1]/10;
                c[i+j-1]%=10;
            }
            c[i+m]=t;
        }
        while ((c[len]==0)and(len>1)) --len;
        for (int i=1;i<=len;++i) a[i]=c[i];
        a[0]=len;
    }
    int main()
    {
        freopen("mason.in","r",stdin);
        freopen("mason.out","w",stdout);
        read(P);
        n=(double)P*(double)(log(2)/log(10));
        n+=1;printf("%d",(int)n);
        a[1]=2;a[0]=1;
        ans[1]=1;ans[0]=1;
        for (;P;P>>=1,cheng1(a,a))
        if (P&1) cheng1(ans,a);
        jian(ans);
        for (int i=500;i>=1;--i)
        {
            if (i%50==0) printf("\n");
            printf("%d",ans[i]);
        }
        return 0;
    }