BZOJ 1002 轮状病毒（生成树计数+高精度）

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

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

现给定n(N<=100)，编程计算有多少个不同的n轮状病毒
Input
第一行有1个正整数n
Output
计算出的不同的n轮状病毒数输出
Sample Input
3
Sample Output
16
Solution
生成树计数，需要用到Kirchhoff Matrix Tree定理：设G为一无向图，称其度数矩阵与邻接矩阵之差R为该图的基尔霍夫矩阵，R的任一n-1阶主子式的行列式的值为G的生成树个数
令圆心为0点，其他点顺次编号1~n，那么由以上定理


知道G[n]=3G[n-1]-G[n-2]-2后用高精度即可得到答案
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 1111
struct BigInt
{
    const static int mod=10000;
    const static int LEN=4;
    int a[maxn],len;
    BigInt() 
    {
        memset(a,0,sizeof(a));
        len=1;
    }
    void init(int x)
    {
        memset(a,0,sizeof(a));
        len=0;
        do
        {
            a[len++]=x%mod;
            x/=mod;
        }while(x);
    }
    int Compare(const BigInt &b)
    {
        if(len<b.len)return -1;
        if(len>b.len)return 1;
        for(int i=len-1;i>=0;i--)
            if(a[i]<b.a[i])return -1;
            else if(a[i]>b.a[i])return 1;
        return 0;
    }
    BigInt operator +(const BigInt &b)const
    {
        BigInt ans;
        ans.len=max(len,b.len);
        for(int i=0;i<=ans.len;i++)ans.a[i]=0;
        for(int i=0;i<ans.len;i++)
        {
            ans.a[i]+=((i<len)?a[i]:0)+((i<b.len)?b.a[i]:0);
            ans.a[i+1]+=ans.a[i]/mod;
            ans.a[i]%=mod;
        }
        if(ans.a[ans.len]>0)ans.len++;
        return ans;
    }
    BigInt operator -(const BigInt &b)const
    {
        BigInt ans;
        ans.len=len;
        int k=0;
        for(int i=0;i<ans.len;i++)
        {
            ans.a[i]=a[i]+k-b.a[i];
            if(ans.a[i]<0)ans.a[i]+=mod,k=-1;
            else k=0;           
        }
        while(ans.a[ans.len-1]==0&&ans.len>1)ans.len--;
        return ans;
    }
    BigInt operator *(const BigInt &b)const
    {
        BigInt ans;
        for(int i=0;i<len;i++)
        {
            int k=0;
            for(int j=0;j<b.len;j++)
            {
                int temp=a[i]*b.a[j]+ans.a[i+j]+k;
                ans.a[i+j]=temp%mod;
                k=temp/mod;
            }
            if(k!=0)ans.a[i+b.len]=k;
        }
        ans.len=len+b.len;
        while(ans.a[ans.len-1]==0&&ans.len>1)ans.len--;
        return ans;
    }
    void output()
    {
        printf("%d",a[len-1]);
        for(int i=len-2;i>=0;i--)
            printf("%04d",a[i]);
        printf("\n");
    }
}ans[111];
int main()
{
    BigInt a,b;
    a.init(3),b.init(2);
    ans[1].init(1),ans[2].init(5);
    for(int i=3;i<=100;i++)
        ans[i]=ans[i-1]*a-ans[i-2]+b;
    int n;
    while(~scanf("%d",&n))ans[n].output();
    return 0;
}