BZOJ1002轮状病毒

1002: [FJOI2007]轮状病毒
Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 3852 Solved: 2106
Description
给定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
高精乘，高精减
附上本蒟蒻的代码：

//f[i]=(f[i-1]*3-f[i-2]+2)
#include<cstdio>
#include<iostream>
using namespace std;
struct kx
{
    int a[101],len;
};
int n;
kx mul(kx 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;
}
kx sub(kx a,kx 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 read()
{
    int w=0,c=1; char ch=getchar();
    while (ch<'0' || ch>'9')
      {
        if (ch=='-') c=-1;
        ch=getchar();
      }
    while (ch>='0' && ch<='9')
      w=w*10+ch-'0',ch=getchar();
    return w*c;
}

int main()
{
    kx f[101];
    f[1].a[1]=1,f[2].a[1]=5,f[1].len=f[2].len=1;
    n=read();
    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;
}