【高精】wikioi 1087 麦森数 （noip2003）

题目链接：http://wikioi.com/problem/1087/

分析：

高精，很明显。

第一问是很简单的，位数可以用公式算出来，只需要求一个对数而已。

数学原理：十进制正整数n的位数为int(log10(n))+1。所以2^P-1的位数int(log10(2)*p)+1 。

对于第二问，算到前500位就可以了，后面不用管它。

不加优化的高精肯定超时，因此我们可以压位处理，因为题目要求保留500位，因此我们压5位，便于最后输出。

压位后仍超时。

快速幂可以，但是我不会。。。

所以换个思路。

本题是求2的幂，每次乘2要乘很多次，所以我们预处理一下，把乘数变成2的k次方，减少相乘次数，从n次变为n/k次，这样会快很多。

对于k的取值，尽量大，别爆掉就好。

先把乘数变成2的k次方，然后用n对k取余得到s，把初值赋成2的s次方，最后再和扩大后的乘数相乘n/k次。

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#define mo 100000
using namespace std;
int main()
{
	long long n,a[101];
	memset(a,0,sizeof(a));
	cin>>n;
	cout<<int(n*log10(2)+1)<<endl;
	a[1]=1;
	long long jj=n%30,ji=1073741824;
	for (long long i=1;i<=jj;i++)	a[1]*=2;
	n-=jj;
	n/=30;
	for (long long i=1;i<=n;i++)
	 {
	 	for (long long j=1;j<=100;j++) 
	 	 {
		  a[j]*=ji;
	    }
		 for (long long j=1;j<=100;j++)
	     if (a[j]>mo)
	      {
	      	a[j+1]+=a[j]/mo;
	      	a[j]%=mo;
	      }
	 }
	 a[1]-=1;
	 for (long long i=1;;i++)
	  {
	  	if (a[i]>0)break;
	  	a[i]+=mo;
	  	a[i-1]--;
	  }
  printf("%05lld",a[100]);
  for (long long i=99;i>=1;i--) 
  {
  printf("%05lld",a[i]);	
  if ((i-1)%10==0) cout<<'\n';
  }
}