poj 1191(棋盘分割)

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<cmath>
using namespace std;
double col[10][10],s[10][10][10][10],var[10][10][10][10][20];
int main()
{
    double n,a;
    double sum=0;
    memset(col,0,sizeof(col));
    memset(var,0,sizeof(var));
    scanf("%lf",&n);
    for(int i=1;i<=8;i++)
    for(int j=1;j<=8;j++)
    {
        scanf("%lf",&a);
        col[i][j]=col[i-1][j]+a;
        //cout<<a<<endl;
        sum+=a;
    }
    double avg=sum/n;
    //cout<<avg<<endl;
    //for(int i=1;i<=8;i++,cout<<endl)
    //for(int j=1;j<=8;j++)
    //{
       // printf("%.0lf ",col[i][j]);
    //}
    for(int i=1;i<=8;i++)
    for(int j=1;j<=8;j++)
    for(int k=i;k<=8;k++)
    for(int h=j;h<=8;h++)
    {
        if(h>1)
        s[i][j][k][h]=s[i][j][k][h-1]+col[k][h]-col[i-1][h];
        else s[i][j][k][h]=col[k][h]-col[i-1][h];
        var[i][j][k][h][1]=((s[i][j][k][h]-avg)*(s[i][j][k][h]-avg))/n;
        //printf("%.0lf ",s[i][j][k][h]);
    }
    //int x,y,z,p;
    //while(cin>>x>>y>>z>>p)
    //cout<<s[x][y][z][p]<<endl;
    for(int k=2;k<=n;k++)
        for(int i=1;i<=8;i++)
            for(int j=1;j<=8;j++)
                for(int f=i;f<=8;f++)
                    for(int h=j;h<=8;h++)
                    {
                        double mi=9999999999;
                        for(int i0=i;i0<f;i0++)
                        {
                            mi=min(var[i][j][i0][h][k-1]+var[i0+1][j][f][h][1],mi);
                            mi=min(var[i][j][i0][h][1]+var[i0+1][j][f][h][k-1],mi);
                        }
                        for(int j0=j;j0<h;j0++)
                        {
                            mi=min(var[i][j][f][j0][k-1]+var[i][j0+1][f][h][1],mi);
                            mi=min(var[i][j][f][j0][1]+var[i][j0+1][f][h][k-1],mi);
                        }
                        var[i][j][f][h][k]=mi;
                    }
    printf("%.3lf\n",sqrt(var[1][1][8][8][(int)n]));
    return 0;
}