luogu 1045 麦森数

题目大意：

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

思路：

一道高精度练习题

其中位数是一个结论 位数=[P*log2]+1

然后就是高精度，因为作死的压位打了好久

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#define inf 2139062143
#define ll long long
#define MAXN 6010
#define MOD 10000
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=getchar();}
    while(isdigit(ch)) {x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int n;
struct bign
{
    int len,num[125];
    bign(){len=0;memset(num,0,sizeof(num));}
    void Print()
    {
        for(int i=124;i>0;i--)
            if(i%25==12) printf("%02d
%02d",num[i]/100,num[i]%100);
            else if(i%25==0) printf("%04d
",num[i]);
            else printf("%04d",num[i]);
        printf("%04d",num[0]-1);
    }
};
bign mul(bign a,bign b)
{
    bign res;int t;
    //cout<<"a: ";a.Print();cout<<"b: ";b.Print();
    for(int i=0;i<=a.len;i++)
        for(int j=0;j<=b.len;j++)
        {
            t=a.num[i]*b.num[j];
            if(i+j<=124) 
                if(i+j<124) res.num[i+j]+=t,res.num[i+j+1]+=res.num[i+j]/MOD,res.num[i+j]%=MOD;
                else res.num[i+j]+=t,res.num[i+j]%=MOD;
        }
    for(int i=0;i<=124;i++)
        if(res.num[i]) res.len=i;
    return  res;
}
void q_pow()
{
    bign ans,t;ans.num[0]=1,t.num[0]=2;
    while(n)
    {
        if(n&1) ans=mul(ans,t);
        t=mul(t,t);
        n>>=1;
    }
    ans.Print();
}
int main()
{
    n=read();
    printf("%d
",(int)(n*0.30103)+1);
    q_pow();
}