【题目描述】
形如 2P−1 的素数称为麦森数,这时P一定也是个素数。但反过来不一定,即如果 P 是个素数。2P−1 不一定也是素数。到1998年底,人们已找到了37个麦森数。最大的一个是P=3021377,它有909526位。麦森数有许多重要应用,它与完全数密切相关。
任务:从文件中输入P(1000 < P < 3100000) ,计算 2P−1 的位数和最后500位数字(用十进制高精度数表示)
【输入格式】
文件中只包含一个整数P(1000 < P < 3100000)
【输出格式】
第1行:十进制高精度数 2P−1 的位数。
第2-11行:十进制高精度数 2P−1 的最后500位数字。(每行输出50位,共输出10行,不足500位时高位补0)
不必验证 2P−1 与P是否为素数。
【样例输入】
1279
【样例输出】
386
00000000000000000000000000000000000000000000000000
00000000000000000000000000000000000000000000000000
00000000000000104079321946643990819252403273640855
38615262247266704805319112350403608059673360298012
23944173232418484242161395428100779138356624832346
49081399066056773207629241295093892203457731833496
61583550472959420547689811211693677147548478866962
50138443826029173234888531116082853841658502825560
46662248318909188018470682222031405210266984354887
32958028878050869736186900714720710555703168729087
对于第一问,利用对数
log
的性质可以直接计算位数
n=log10(2P)+1=P∗log102+1
(注,本公式其实把
2P−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;
}