CodeVs 1087 麦森树

1087 麦森树

1. 刚开始的时候想着用log 10加上换底公式 就可以得到位数，500位结果可以用数组模拟，但发现p的次数太大了，会超时，想到了快速幂，但是不知道具体如何实施，因为这里涉及了大数
2. 然后在网上看了一份代码，是用高精度乘法再配合快速幂进行运算的，觉得太漂亮了！这个快速幂主要还是把p当作二进制这样来考虑，可以降低很多的复杂度，减少运算次数。代码转载如下：
   
int a[501],b[501];
int pow(int a[],int b[])//高精度乘法
{ int c[1001];//c要开大一点否则数组越位
memset(c,0,sizeof(c));
for(int i=0;i<500;i++)
{int x=0;
for(int j=0;j<500;j++)
{
c[i+j]=c[i+j]+a[i]*b[j]+x;
x=c[i+j]/10;
c[i+j]%=10;
}
}
memcpy(a,c,500*sizeof(int));//只要管500位之前的就行了
}
int main()
{ long long n;
cin>>n;
cout<<int(n*log10(2))+1<<endl;
a[0]=1;b[0]=2;
while(n)//快速幂
{
if(n%2==1){pow(a,b);}
pow(b,b);
n=n/2;
}
a[0]--;//记得要减一
for(int i=499;i>=0;i--)
{
cout<<a[i];
if(i%50==0){cout<<endl;}
}
return 0;
}