1 //形如2n-1的素数称为麦森数，这时n一定也是个素数。但反过来不一定，即如果n是个素数。2n-1不一定也是素数。
2
3 #include<iostream>
4 #include<cmath>
5 #include<cstdio>
6 #include<cstring>
7 #define N 126
8 using namespace std;
9 int ans[N],anspow[N];
10 void mult(int ans[],int anspow[])
11 {
12 int i,j;
13 int c[N];
14 memset(c,0,sizeof(c));
15 for(i=0;i<N;i++)
16 {
17 for(j=0;j<N;j++)
18 {
19 if(i+j<N)//超出500位的部分不计算
20 {
21 c[i+j]+=ans[j]*anspow[i];
22 }
23 }
24 for(j=0;j<N-1;j++)
25 {
26 if(c[j]>=10000)
27 {
28 c[j+1]+=c[j]/10000;//压4位
29 c[j]%=10000;
30 }
31 }
32 }
33 memcpy(ans,c,N*sizeof(int)); //复制函数
34 }
35 int main()
36 {
37 int P,i,j;
38 while(cin>>P)
39 {
40 memset(ans,0,sizeof(ans));
41 memset(anspow,0,sizeof(anspow));
42 printf("%d\n",(int)(P*log10(2)+1));
43 ans[0]=1;
44 anspow[0]=2;
45 /************关键部分计算2^P*******
46 2^p=（2^1）*（2^2）*（2^3）*（2^4）*（2^5）………………
47 简单说下：P=5 -----101（二进制）
48 p & 1 =1（最右边一位） -->>>ans=2 anspow=2
49 p>>1=110 anspow=2^2
50 p & 1 =0 此时表明2^3不存在 ans=2 anspow=2^4
51 p>>1=1
52 p & 1 =1 ------>>>>> ans=2^5
53 p>>1=0 ------结束
54 ************************************/
55 while(P)
56 {
57 if( P & 1)
58 mult(ans,anspow);
59 P>>=1;
60 mult(anspow,anspow);
61 }
62 ans[0]--;//2^P的个位为2，4，6，8，故可以-1
63 /****************输出格式的控制************************/
64 for(i=124;i>=0;i--)
65 {
66 if(i%25==12)
67 {
68 printf("%02d\n%02d",ans[i]/100,ans[i]%100);
69 }
70 else
71 {
72 printf("%04d",ans[i]);
73 if(i%25==0) printf("\n");
74 }
75 }
76 /***************************************************/
77 }
78 return 0;
79 }
1 /*************麦森数****************/
2 #include<iostream>
3 #include<cstdio>
4 #include<cmath>
5 #define N 100 //压5位
6 using namespace std;
7 int ans[N];
8 void mult(int t)
9 {
10 int i,temp,last=0;
11 for ( i=N-1; i>=0; i--)
12 {
13 temp=(ans[i]<<t)+last; //乘2^t，加进位
14 last=temp/100000;
15 ans[i]=temp-last*100000; //temp%100000
16 }
17 }
18 void output()
19 {
20 int i;
21 for(i=1;i<=N;i++)
22 { printf("%05d",ans[i-1]);
23 if (i%10==0)cout<<endl;
24 }
25 }
26 int main()
27 {
28 int P,times;
29 while(cin>>P)
30 {
31 memset(ans,0,sizeof(ans));
32 ans[N-1]=1;
33 cout<<(int)(P*log10(2)+1)<<endl;
34 /***********************关键部分*****************
35 2^P=2^(14*times)*2^(P%14) 用移位
36 之所以取14的原因 2^14*(99999)=1.638382*10^9 在（int）
37 范围，15以后都会超int，主要体现在mult()中
38 **********************************************/
39 times=P/14; //只能到14 15以后压缩时会超范围
40 while(times--)
41 mult(14);
42 mult(P%14);
43 --ans[99];
44 output();
45 }
46 return 0;
47 }
