高精度 麦森数

【原题链接】

NOIP 2003

POJ2706 http://bailian.openjudge.cn/practice/2706/

洛谷1045 http://www.luogu.org/problem/show?pid=1045

【题目描述】

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

 

        最开始做的时候用高精度，但把每一位都算了出来，题目只要求输出后500位，所以只用算到Long，否则就会超时。还有一个求位数的快速方法：log10（2）*p，即  log10(2)*p   将p提到了前面。

【代码实现】

#include<iostream>

#include<cmath>

#define MAXL 10000

#define Long 505

using namespace std;

int n;

inline string mul(string s1,string s2)//高精乘，后面再快速幂里用到

{

    if(s1.size()

    inta[MAXL]={0},b[MAXL]={0},c[MAXL<<1]={0};//左移一位，相当于乘2,三个数组都必须清零

    intl1=s1.size(),l2=s2.size();

    for(int i=1;i<=l2; ++i)a[i]=s1[l1-i]-48,b[i]=s2[l2-i]-48;

    for(int i=l2+1;i<=l1; ++i)a[i]=s1[l1-i]-48;

    for(int i=1;i<=Long; ++i)//只用算到Long 否则会超时

        for(intj=1; j<=Long; ++j)

            c[i+j-1]+=a[i]*b[j];

    for(int i=1;i<=Long; ++i)

        if(c[i]>=10)

        {

            c[i+1]+=c[i]/10;

            c[i]%=10;

        }

    string s;

    //if(c[l1+l2])s+=c[l1+l2]+48;

    for(int i=Long;i>=1; --i)

        s+=c[i]+48;

    return s;

}

string qpow(string s,int n)//高精快速幂

{

    stringret="1";

    while(n)

    {

        if(n%2)

            ret=mul(ret,s);

        s=mul(s,s);

        n/=2;

    }

    return ret;

}

int main()

{

    cin>>n;

    cout<<ceil(log10(2)*n)<<endl;//快速计算位数调用cmath

    strings=qpow("2",n);

    int len=s.size();

    if(len>500)s.erase(0,len-500);//若大于500位抹掉前面的位数

    {

        for(inti=0;i<=498;i++)

        {

            cout<<s[i];

            if(iP==49)cout<<endl;//每到第50位换行

        }

        cout<<s[len-1]-'1';//最后一位-1后输出

    }

    return 0;

}