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;}