BZOJ1002: [FJOI2007]轮状病毒 (DP)

标准做法似乎应该是计算生成树数量的基尔霍夫矩阵之类的..

我看到的做法是一个神奇的高精度dp，当然以后这个blahblahblah矩阵还是要搞一下。。

 

参考(抄袭)网址

 

这个dp的原理就是把环拆成一条含特定点的链和剩下部分(可用dp解决),这样就避免了环具有的一些dp不好解决的奇怪判定.

非常神奇

%想出这个办法的dalao

 

附上非常不走心的非常丑的自己的代码..

 

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
using namespace std;
int n;
int f[110][2][35]={};//1 到i的单独一段都和中间连上了 0 到i的单独一段没有和中间连上
int ans[35]={};
int z[35]={};
int a[35]={};
int b[35]={};
int w=0;
void plu(){//a+b存z
    int e=0;
    for(int i=1;i<=30;i++){
        e=a[i]+b[i]+e;
        z[i]=e%10000;
        e/=10000;
    }
}
void mul(){//b*w存z
    int e=0;
    for(int i=1;i<=30;i++){
        e=w*b[i]+e;
        z[i]=e%10000;
        e/=10000;
    }
}
int main(){
    scanf("%d",&n);
    f[1][1][1]=f[1][0][1]=1;
    f[0][1][1]=1;
    for(int i=2;i<=n;i++){
        for(int j=1;j<=30;j++){
            a[j]=f[i-1][0][j];
            b[j]=f[i-1][1][j];
        }
        plu();
        for(int j=1;j<=30;j++){
            f[i][0][j]=z[j];
            z[j]=0;
        }
        
        w=2;
        mul();
        for(int j=1;j<=30;j++){
            b[j]=z[j];
            z[j]=0;
        }
        plu();
        for(int j=1;j<=30;j++){
            f[i][1][j]=z[j];
            z[j]=0;
        }
    }
    for(int i=1;i<=n;i++){
        w=i*i;
        for(int j=1;j<=30;j++){
            b[j]=f[n-i][1][j];
            a[j]=ans[j];
        }
        mul();
        for(int j=1;j<=30;j++){
            b[j]=z[j];
            z[j]=0;
        }
        plu();
        for(int j=1;j<=30;j++){
            ans[j]=z[j];
            z[j]=0;
        }
    }
    int f=0;
    for(int i=30;i>=1;i--){
        if(ans[i]!=0&&f==0){
            f=1;
            printf("%d",ans[i]);
            continue;
        }
        if(f){
            if(ans[i]>999){
                cout<<ans[i];
            }
            else if(ans[i]>99){
                cout<<0<<ans[i];
            }
            else if(ans[i]>9){
                cout<<0<<0<<ans[i];
            }
            else{
                cout<<0<<0<<0<<ans[i];
            }
        }
    }
    cout<<endl;
    return 0;
}