BZOJ-1002: [FJOI2007]轮状病毒(打表找规律or递推 + 高精度)

1002: [FJOI2007]轮状病毒

Time Limit: 1 Sec  Memory Limit: 162 MB
Submit: 5803  Solved: 3154
[Submit][Status][Discuss]

Description

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

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

　　现给定n(N<=100)，编程计算有多少个不同的n轮状病毒

 

Input

　　第一行有1个正整数n

Output

　　计算出的不同的n轮状病毒数输出

Sample Input

3

Sample Output

16

HINT

 

Source

我承认我承认……这题是在路边沾花惹草遇到的……因为看到那个递推的式子极其的简单所以就写了 _(:зゝ∠)_ 对了这貌似是我第一次写高精度减法 *▽*

http://blog.csdn.net/Clove_unique/article/details/51438878

#include "bits/stdc++.h"
using namespace std;
typedef long long LL;
const int MAX=105;
int n;
struct Node{
    int a[505];
    Node (){memset(a,0,sizeof(a));}
    Node operator * (const Node tt) {
        int i,j,k;
        Node an;an.a[0]=a[0]+tt.a[0];
        for (i=1;i<=a[0];i++){
            for (j=1;j<=tt.a[0];j++){
                an.a[i+j-1]+=a[i]*tt.a[j];
                an.a[i+j]+=an.a[i+j-1]/10;
                an.a[i+j-1]%=10;
            }
        }
        while (an.a[an.a[0]]==0) an.a[0]--;
        return an;
    }
    Node operator + (const Node &tt) {
        int i,j;
        Node an;an.a[0]=max(a[0],tt.a[0])+1;
        for (i=1;i<=max(a[0],tt.a[0]);i++){
            an.a[i]+=a[i]+tt.a[i];
            an.a[i+1]+=an.a[i]/10;
            an.a[i]%=10;
        }
        while (an.a[an.a[0]]==0) an.a[0]--;
        return an;
    }
    Node operator - (const Node &tt) {
        int i,j;
        Node an;an.a[0]=a[0];
        for (i=1;i<=a[0];i++){
            if (a[i]<tt.a[i]) a[i+1]--,a[i]+=10;
            an.a[i]=a[i]-tt.a[i];
        }
        while (an.a[an.a[0]]==0) an.a[0]--;
        return an;
    }
}f[MAX];
int main(){
    freopen ("virus.in","r",stdin);
    freopen ("virus.out","w",stdout);
    int i,j;
    scanf("%d",&n);
    f[1].a[0]=f[1].a[1]=f[2].a[0]=1;f[2].a[1]=5;
    Node x;x.a[0]=1,x.a[1]=2;
    Node y;y.a[0]=1,y.a[1]=3;
    for (i=3;i<=n;i++)
        f[i]=f[i-1]*y-f[i-2]+x;
    for (i=f[n].a[0];i>0;i--)
        printf("%d",f[n].a[i]);
    return 0;
}