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[][][]={};//1 到i的单独一段都和中间连上了 0 到i的单独一段没有和中间连上
 int ans[]={};
 int z[]={};
 int a[]={};
 int b[]={};
 int w=;
 void plu(){//a+b存z
     int e=;
     for(int i=;i<=;i++){
         e=a[i]+b[i]+e;
         z[i]=e%;
         e/=;
     }
 }
 void mul(){//b*w存z
     int e=;
     for(int i=;i<=;i++){
         e=w*b[i]+e;
         z[i]=e%;
         e/=;
     }
 }
 int main(){
     scanf("%d",&n);
     f[][][]=f[][][]=;
     f[][][]=;
     for(int i=;i<=n;i++){
         for(int j=;j<=;j++){
             a[j]=f[i-][][j];
             b[j]=f[i-][][j];
         }
         plu();
         for(int j=;j<=;j++){
             f[i][][j]=z[j];
             z[j]=;
         }

         w=;
         mul();
         for(int j=;j<=;j++){
             b[j]=z[j];
             z[j]=;
         }
         plu();
         for(int j=;j<=;j++){
             f[i][][j]=z[j];
             z[j]=;
         }
     }
     for(int i=;i<=n;i++){
         w=i*i;
         for(int j=;j<=;j++){
             b[j]=f[n-i][][j];
             a[j]=ans[j];
         }
         mul();
         for(int j=;j<=;j++){
             b[j]=z[j];
             z[j]=;
         }
         plu();
         for(int j=;j<=;j++){
             ans[j]=z[j];
             z[j]=;
         }
     }
     int f=;
     for(int i=;i>=;i--){
         if(ans[i]!=&&f==){
             f=;
             printf("%d",ans[i]);
             continue;
         }
         if(f){
             if(ans[i]>){
                 cout<<ans[i];
             }
             else if(ans[i]>){
                 cout<<<<ans[i];
             }
             else if(ans[i]>){
                 cout<<<<<<ans[i];
             }
             else{
                 cout<<<<<<<<ans[i];
             }
         }
     }
     cout<<endl;
     return ;
 }