题目链接: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';
}
}