BZOJ1002: [FJOI2007]轮状病毒

n<=100的形如的图有多少生成树。不取模。

$f(i)=3*f(i-1)-f(i-2)+2$，VFK的题解

#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
//#include<queue>
#include<math.h>
//#include<time.h>
//#include<iostream>
using namespace std;

int n;
#define maxn 311
const int base=100000000;
struct LLL
{
    int num[maxn],len;
    LLL() {len=0; memset(num,0,sizeof(num));}
    void operator = (int x)
    {
        if (x==0) {num[len=1]=0; return;}
        while (x) {num[++len]=x%base; x/=base;}
    }
    LLL operator + (const LLL &b)
    {
        LLL ans;
        for (int i=1,top=ans.len=max(len,b.len);i<=top;i++)
        {
            ans.num[i]+=(i<=len?num[i]:0)+(i<=b.len?b.num[i]:0);
            if (ans.num[i]>=base)
            {
                ans.num[i]-=base;
                ans.num[i+1]++;
            }
        }
        ans.len+=(ans.num[ans.len+1]>0);
        return ans;
    }
    LLL operator - (const LLL &b)
    {
        LLL ans;
        for (int i=1,top=ans.len=len;i<=top;i++)
        {
            ans.num[i]+=num[i]-(i<=b.len?b.num[i]:0);
            if (ans.num[i]<0)
            {
                ans.num[i]+=base;
                ans.num[i+1]--;
            }
        }
        ans.len-=(ans.num[ans.len]==0);
        return ans;
    }
    LLL operator + (int b)
    {
        LLL tmp; tmp=b;
        return *this+tmp;
    }
    void print()
    {
        printf("%d",num[len]);
        for (int i=len-1;i;i--) printf("%08d",num[i]);
    }
};

LLL a[maxn];
int main()
{
    scanf("%d",&n);
    a[1]=1; a[2]=5;
    for (int i=3;i<=n;i++) a[i]=a[i-1]+a[i-1]+a[i-1]-a[i-2]+2;
    a[n].print();
    return 0;
}