1002: [FJOI2007]轮状病毒

题目链接

题目大意：求有多少个n轮状病毒

题解：用奇怪的东西推出f[i]=f[i-1]*3-f[i-2]+2，无限Orz

我的收获：练习了下高精度模板……

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;

#define rad 10000
#define B 4

int n;

struct Bi{
    int s[2005],l;
    Bi(){memset(s,0,sizeof(s));l=1;}
}f[105];

Bi operator +(Bi a,int x){
    a.s[1]+=x;
    for(int i=1;i<=a.l;i++)
        if(a.s[i]>=rad) a.s[i]-=rad,a.s[i+1]++;
    if(a.s[a.l+1]) a.l++;
    return a;
}

Bi operator -(Bi a,Bi b){
    for(int i=1;i<=a.l;i++){
        a.s[i]-=b.s[i];
        if(a.s[i]<0) a.s[i]+=rad,a.s[i+1]--;
    }
    while(a.l>1&&!a.s[a.l]) a.l--;
    return a;
}

Bi operator *(Bi a,int x){
    for(int i=1;i<=a.l;i++) a.s[i]*=x;
    for(int i=1;i<=a.l;i++){
        a.s[i+1]+=a.s[i]/rad;
        a.s[i]%=rad;
    }
    if(a.s[a.l+1]) a.l++;
    return a;
}

void print(Bi a){printf("%d",a.s[a.l]);for(int i=a.l-1;i;i--) printf("%04d",a.s[i]);}

void init()
{
    cin>>n;
    f[1].s[1]=1,f[2].s[1]=5;
    for(int i=3;i<=n;i++)
        f[i]=f[i-1]*3-f[i-2]+2;
    print(f[n]);
}

int main()
{
    init();
    return 0;
}