bzoj1002: [FJOI2007]轮状病毒

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define rep(i,n) for(int i=1;i<=n;i++)
#define clr(x,c) memset(x,c,sizeof(x))
#define REP(i,s,t) for(int i=s;i<=t;i++)
struct node{
	int a[100];int len;
};
node f[105];
node mly(node a,int x){
	rep(i,a.len) a.a[i]*=x;
	rep(i,a.len) a.a[i+1]+=a.a[i]/10,a.a[i]%=10;
	if(a.a[a.len+1]) a.len++;
	return a;
}
node down(node a,node b){
	rep(i,a.len) {
		a.a[i]-=b.a[i];
		if(a.a[i]<0) a.a[i]+=10,a.a[i+1]--;
	}
	a.a[1]+=2;int cur=1;
	while(a.a[cur]>=10) a.a[cur]%=10,a.a[cur+1]++,cur++;
	while(!a.a[a.len]) a.len--;
	return a;
}
int main(){
	int n;scanf("%d",&n);
	f[1].a[1]=1;f[2].a[1]=5;f[1].len=f[2].len=1;
	REP(i,3,n) f[i]=down(mly(f[i-1],3),f[i-2]);
	for(int i=f[n].len;i;i--) printf("%d",f[n].a[i]);
	return 0;
}