原题链接:https://www.luogu.com.cn/problem/P1045
求2p-1的最后500位数 (P<3100000)
以前碰到高精都是python大法好
想写一下c的高精,第一次写,踩了好多坑
这个题是高精乘法+快速幂
第一小问是求位数,可以用一个求位数的公式
k=log10(N)+1
2p-1和2p的位数是一样的,所以直接求2p的位数即可
log10(2p)=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;
}