原题链接：https://www.luogu.com.cn/problem/P1045

求2^p-1的最后500位数 （P<3100000）

以前碰到高精都是python大法好
想写一下c的高精，第一次写，踩了好多坑
这个题是高精乘法+快速幂

第一小问是求位数，可以用一个求位数的公式
k=log10(N)+1
2^p-1和2^p的位数是一样的，所以直接求2^p的位数即可
log10（2^p）=p*log10(2)

代码：

#include <bits/stdc++.h>
#define ll long long
using namespace std;
ll a[1010],s[1010],t[1010];
void result1()
{
    memset(t,0,sizeof(t));
    for(int i=1;i<=500;i++)
    {
        for(int j=1;j<=500;j++)
        {
            t[i+j-1]+=s[i]*a[j];
        }
    }
    for(int i=1;i<=500;i++)
    {
        t[i+1]+=t[i]/10;
        t[i]%=10;
        s[i]=t[i];
    }
}
void result2()
{
    memset(t,0,sizeof(t));
    for(int i=1;i<=500;i++)
    {
        for(int j=1;j<=500;j++)
        {
            t[i+j-1]+=a[j]*a[i];
        }
    }
    for(int i=1;i<=500;i++)
    {
        t[i+1]+=t[i]/10;
        t[i]%=10;
        a[i]=t[i];
    }
}
int main()
{
    ll p;
    cin>>p;
    ll kk=p*log10(2)+1;
    s[1]=1;a[1]=2;
    while(p)
    {
        if(p&1)
        {result1();}
        result2();
        p/=2;
    }
    s[1]--;
    cout<<kk<<endl;
    for(int i=1;i<=500;i++)
    {
        int c=500-i+1;
        if(i%50==0)
        {cout<<s[c]<<endl;}
        else
        {cout<<s[c];}
    }
    return 0;
}