【bzoj1002】[FJOI2007]轮状病毒

Description

　　轮状病毒有很多变种，所有轮状病毒的变种都是从一个轮状基产生的。一个N轮状基由圆环上N个不同的基原子
和圆心处一个核原子构成的，2个原子之间的边表示这2个原子之间的信息通道。如下图所示

　　N轮状病毒的产生规律是在一个N轮状基中删去若干条边，使得各原子之间有唯一的信息通道，例如共有16个不
同的3轮状病毒，如下图所示

　　现给定n(N<=100)，编程计算有多少个不同的n轮状病毒
Input

　　第一行有1个正整数n

Output

　　计算出的不同的n轮状病毒数输出

Sample Input

3
Sample Output

16

代码
f[i]=(f[i-1]*3-f[i-2]+2)

#include<cstring>
#include<iostream>
#include<cstdio>
#include<cmath>
int a[10000],b[10000],c[10000];
using namespace std;
int main()
{
    int n;
    scanf("%d",&n);
    if (n==1){cout<<"1";return 0;}
    if (n==2){cout<<"5";return 0;}
    a[1]=1;b[1]=5;
    int lena=1,lenb=1,lenc=0;
    for (int l=3;l<=n;l++)
    {
        int z=0;
        for (int i=1;i<=lenb||z!=0;i++)
        {
            c[i]=b[i]*3+z;
            z=c[i]/10;
            c[i]%=10;
            lenc=max(lenc,i);
        }
        z=0;
        for (int i=1;i<=lena||z!=0;i++)
        {
            c[i]=c[i]-z-a[i];
            if (c[i]<0) z=1;else z=0;
            c[i]=(c[i]+10)%10;
        }
        while (c[lenc]==0) lenc--;
        c[1]+=2;z=c[1]/10;c[1]%=10;
        for (int i=2;z!=0;i++)
        {
            c[i]+=z;
            z=c[i]/10;
            c[i]%=10;
        }
        for (int i=1;i<=lenb;i++)
        {
            a[i]=b[i];
        }
        for (int i=1;i<=lenc;i++)
        {
            b[i]=c[i];
        }
        lena=lenb;lenb=lenc;
    }
    for (int i=lenb;i>=1;i--)
    {
        printf("%d",b[i]);
    }
    return 0;
}