bzoj1002[FJOI2007]轮状病毒 递推+高精度

题意比较简单就不说了。
一开始我用组合数做，发现WA了。。并不知道为什么错，感觉应该是会算重。

主要是用基尔霍夫矩阵推出一个结论，但是我发现n=1,n=2的时候怎么大家的答案都不太一样。结论就是f[i]=(f[i-1]*3-f[i-2]+2),注意要用高精度不然会爆。。

#include
#include
using namespace std;
struct data{
       int a[101],len;
       };
int n;
data mul(data a,int k)
{
    for(int i=1;i<=a.len;i++)
            a.a[i]*=k;
    for(int i=1;i<=a.len;i++)
    {
            a.a[i+1]+=a.a[i]/10;
            a.a[i]%=10;
            }
    if(a.a[a.len+1]!=0)a.len++;
    return a;
}
data sub(data a,data b)
{
    a.a[1]+=2;
    int j=1;
    while(a.a[j]>=10){a.a[j]%=10;a.a[j+1]++;j++;}
    for(int i=1;i<=a.len;i++)
    {
           a.a[i]-=b.a[i];
           if(a.a[i]<0){a.a[i]+=10;a.a[i+1]--;}
    }
    while(a.a[a.len]==0)a.len--;
    return a;
}
int main()
{
    data f[101];f[1].a[1]=1;f[2].a[1]=5;
    f[1].len=f[2].len=1;
    scanf("%d",&n);
    for(int i=3;i<=n;i++)
            f[i]=sub(mul(f[i-1],3),f[i-2]);
    for(int i=f[n].len;i>0;i--)
       printf("%d",f[n].a[i]);
    return 0;
}