NOI1999 JZYZOJ1289 棋盘分割 dp 方差的数学结论

http://172.20.6.3/Problem_Show.asp?id=1289

除了下标一坨一坨屎一样挺恶心其他都还挺容易的dp，这道题才发现scanf保留小数位是四舍五入的，惊了。

f[k][x1][y1][x2][y2]

嗯写的时候猜错结论了，本来以为是求下属分配方案中平方和与平均数平方*k的差最小的方案赋给f，没想到是直接找最小的。

代码

 #include<cstdio>
 #include<cstring>
 #include<iostream>
 #include<algorithm>
 #include<cmath>
 #include<queue>
 using namespace std;
 int n;
 int a[][]={};
 double f[][][][][]={};
 bool vis[][][][][]={};
 int main(){
     scanf("%d",&n);
     for(int i=;i<=;i++){
         for(int j=;j<=;j++){
             scanf("%d",&a[i][j]);
             a[i][j]+=a[i][j-];
             a[i][j]+=a[i-][j];
             a[i][j]-=a[i-][j-];
         }
     }
     for(int i=;i<=;i++){
         for(int j=;j<=;j++){
             for(int w=i;w<=;w++){
                 for(int q=j;q<=;q++){
                     f[][i][j][w][q]=1.0*(a[w][q]-a[i-][q]-a[w][j-]+a[i-][j-]);
                     f[][i][j][w][q]*=f[][i][j][w][q];vis[][i][j][w][q]=;
                 }
             }
         }
     }
     double ro=1.0*a[][]/(1.0*n);
     ro=ro*ro;
     for(int k=;k<=n;k++){
         for(int i=;i<=;i++){
             for(int j=;j<=;j++){
                 for(int w=i;w<=;w++){
                     for(int q=j;q<=;q++){
                         for(int z=i;z<w;z++){
                             if((vis[k-][i][j][z][q]&&vis[][z+][j][w][q])&&(vis[k][i][j][w][q]==||f[k][i][j][w][q]-(f[k-][i][j][z][q]+f[][z+][j][w][q])>)){
                                 f[k][i][j][w][q]=f[k-][i][j][z][q]+f[][z+][j][w][q];vis[k][i][j][w][q]=;
                             }
                             if((vis[][i][j][z][q]&&vis[k-][z+][j][w][q])&&(vis[k][i][j][w][q]==||f[k][i][j][w][q]-(f[][i][j][z][q]+f[k-][z+][j][w][q])>)){
                                 f[k][i][j][w][q]=f[][i][j][z][q]+f[k-][z+][j][w][q];vis[k][i][j][w][q]=;
                             }
                         }
                         for(int z=j;z<q;z++){
                             if((vis[k-][i][j][w][z]&&vis[][i][z+][w][q])&&(vis[k][i][j][w][q]==||f[k][i][j][w][q]-(f[k-][i][j][w][z]+f[][i][z+][w][q])>)){
                                 f[k][i][j][w][q]=f[k-][i][j][w][z]+f[][i][z+][w][q];vis[k][i][j][w][q]=;
                             }
                             if((vis[][i][j][w][z]&&vis[k-][i][z+][w][q])&&(vis[k][i][j][w][q]==||f[k][i][j][w][q]-(f[][i][j][w][z]+f[k-][i][z+][w][q])>)){
                                 f[k][i][j][w][q]=f[][i][j][w][z]+f[k-][i][z+][w][q];vis[k][i][j][w][q]=;
                             }
                         }
                     }
                 }
             }
         }
     }
     double z=f[n][][][][]/n-ro;
     z=sqrt(z);
     printf("%.3f\n",z);
     return ;
 }