bzoj 1002: [FJOI2007]轮状病毒

#include<cstdio>
#include<iostream>
using namespace std;
int f[102][102],n;
void mul(int a1)
{
f[a1][0]=f[a1-1][0];
for(int i=1;i<=f[a1-1][0];i++)
{
f[a1][i]+=f[a1-1][i]*3;
f[a1][i+1]+=f[a1][i]/10;
f[a1][i]%=10;
}
if(f[a1][f[a1][0]+1])
f[a1][0]++;
return;
}
void jia(int a1)
{
f[a1][1]+=2;
for(int i=1;i<=f[a1][0];i++)
{
f[a1][i]-=f[a1-2][i];
if(f[a1][i]<0)
{
f[a1][i+1]--;
f[a1][i]+=10;
}
}
for(;f[a1][f[a1][0]]==0;)
f[a1][0]--;
return;
}
int main()
{
scanf("%d",&n);
f[1][0]=1;
f[1][1]=1;
f[2][0]=1;
f[2][1]=5;
for(int i=3;i<=n;i++)
{
mul(i);
jia(i);
}
for(int i=f
[0];i;i--)
printf("%d",f
[i]);
return 0;
}