百练 2706 麦森数

/*
这道题过得好艰难啊
一直TLE
后来看到大神们都用位运算什么的
敬仰之心顿生。。。 
之前直接模拟大数乘法
导致输入3021377时很久跳不出来
后来把2的指数mod10算，结果还是不行
最后就采用了下面的方法
很巧妙
能大大减少循环次数
输入3021377一瞬间就出来了
提交后status上显示0ms
算法才是王道啊 
*/
#define LOCAL
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#define N 126
using namespace std;
void mutiply(int a[N],int b[N])
{
       int i,j,c[N]={0};
       for(i=0;i<N;i++)
       {
              for(j=0;j<N;j++)                
              {
                 if(i+j<N)   //这里加个判断条件一是不需要算超过500位的二是会提示内存读取错误 
                    c[j+i]+=a[j]*b[i];                           
              }
              for(j=0;j<N-1;j++)
              {
                     if(c[j]>=10000)
                     {
                           c[j+1]+=c[j]/10000;
                           c[j]%=10000;            
                     }                  
              }
       }     
       memcpy(a,c,N*sizeof(int));
}
int main()
{
#ifdef LOCAL
       freopen("input.txt","r",stdin);
       freopen("output.txt","w",stdout);
#endif

    
    int n,p,i,ans[N]={1},two[N]={2};
    cin>>p;n=p;
    while(p)
    {
           if(p&1)        //当p的第一位为1时才乘否则直接跳过 
             mutiply(ans,two);
           mutiply(two,two);//这里是求2^(2^(i+1))
           p=p>>1;        
    }
    ans[0]--;
    cout<<(int)(log10(2.0)*n+1)<<endl;//用公式求 
    for(i=124;i>=0;i--)
    {
            if(i%25==12)
                printf("%02d\n%02d",ans[i]/100,ans[i]%100);
            else
            {printf("%04d",ans[i]);if(i%25==0) printf("\n");}                
    }
    return 0;
}