P1045メイソン数（数論＆高精度）

P1045メイソン数（数論＆高精度）

タイトルポータル

トピック：2 ^ p-1桁、最後の500桁を出力
アイデア：2 ^ p-1は2 ^ pと同じです（2 ^ p桁は0ではないため）。したがって、2 ^ p 10 ^ n形式に変換するだけ

ACコード：

#include<bits/stdc++.h>
using namespace std;
const int N=1e3+5;
int ans[N],tmp[N],a[N],p;
void muti(int *x,int *y){ //高精x*y 
	memset(tmp,0,sizeof tmp);
	for(int i=1;i<=500;i++)
		for(int j=1;j<=500;j++)
			tmp[i+j-1]+=x[i]*y[j];
	for(int i=1;i<=500;i++)
	{
		tmp[i+1]+=tmp[i]/10;
		tmp[i]%=10;
	}
	memcpy(ans,tmp,sizeof a);
}
void fun(int *f){ //高精平方 
	memset(tmp,0,sizeof tmp);
	for(int i=1;i<=500;i++)
		for(int j=1;j<=500;j++)
			tmp[i+j-1]+=f[i]*f[j];
	for(int i=1;i<=500;i++)
	{
		tmp[i+1]+=tmp[i]/10;
		tmp[i]%=10;
	}
	memcpy(a,tmp,sizeof a);
}  
void ksm(int *x,int n){ //快速幂 
	while(n){
		if(n&1) muti(ans,x);
		fun(x);
		n>>=1;
	}
}
int main(){
	cin>>p;
	cout<<(int)(log10(2)*p+1)<<endl;//求位数 
	a[1]=2,ans[1]=1;//初始化 
	ksm(a,p);
	ans[1]-=1;
	for(int i=500;i>=1;i--){
		if(i<500&&i%50==0) printf("\n%d",ans[i]);//控制格式 
		else printf("%d",ans[i]);
	}
	return 0;
}