[noip2003]麦森数（高精度+ksm）

题目：

我是超链接

题解：

那个后500位很好求啊，只乘后500位就行啊
问题是怎么求位数，因为是$2^p$，尾数一定不是0，所以-1对位数并不会有影响，那么我们就可以求出$k=2^p$的位数就行啦
我们知道$10^p$的位数是p+1问题是怎么把$2^p$转化为$10^p$
$10^{log_{10}2}=2$我们把这个带进去$k=2^p=(10^{log_{10}2})^p=10^{p*log_{10}2}$
那么位数就是$p*log_{10}2+1$咯

代码：

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int k;
struct bignum{int l,a[505];void cl(){memset(a,0,sizeof(a));l=0;}}mul,ans;
bignum gjj(bignum a,bignum b)
{
    int x=0;bignum c;c.cl();
    c.l=min(500,a.l+b.l);
    for (int i=1;i<=a.l && i<=500;i++)
    {
        x=0;
        for (int j=1;j<=b.l && i+j-1<=500;j++)
        {
        c.a[i+j-1]+=a.a[i]*b.a[j]+x;
        x=c.a[i+j-1]/10;
        c.a[i+j-1]%=10;
        }
        if (i+b.l<=500) c.a[i+b.l]+=x;
    }
    return c;
}
void ksm(int k)
{
    int kk=k;
    mul.l=1;mul.a[1]=2;ans.l=1;ans.a[1]=1;
    for (;k;k>>=1,mul=gjj(mul,mul))
      if (k&1) ans=gjj(ans,mul);
    ans.a[1]--;
    printf("%d\n",(int)(log10(2)*kk+1));
    for (int i=500;i>=1;i-=50) 
    {
        for (int j=i;j>i-50;j--)
          printf("%d",ans.a[j]);
        printf("\n");
    }
}
int main()
{
    scanf("%d",&k);
    ksm(k);
}