poj1191 棋盘分割。 dp

连接：http://poj.org/problem?id=1191

思路：额，其实就是直接搞记录一下就可以了。

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <stdlib.h>
#include <vector>
#include <queue>
#include <stack>
#include <cmath>
#define loop(s,i,n) for(i = s;i < n;i++)
#define cl(a,b) memset(a,b,sizeof(a))
int map[10][10];
int sum[10][10][10][10];
int dp[15][10][10][10][10];
double average;
int n;
using namespace std;
void init()
{
    int i,j,ii,jj,iii,jjj;
    loop(1,i,9)
    {
        loop(1,j,9)
        {
            loop(i,ii,9)
            {
                loop(j,jj,9)
                {
                    int ans;
                    ans = 0;
                    loop(i,iii,ii+1)
                        loop(j,jjj,jj+1)
                        ans+=map[iii][jjj];
                    ans = (ans)*(ans);
                    sum[i][j][ii][jj] += ans;
                }

            }
        }
    }
}
int deal(int k,int x,int y,int xx,int yy)
{

    if(dp[k][x][y][xx][yy] >= 0) return dp[k][x][y][xx][yy];
    if(k == n)
    return sum[x][y][xx][yy];

    int i,j,ans;
    dp[k][x][y][xx][yy] = 20000000;
    for(i = x;i <= xx-1;i++)
    {
        ans = min(deal(k+1,x,y,i,yy)+sum[i+1][y][xx][yy],deal(k+1,i+1,y,xx,yy)+sum[x][y][i][yy]);
        dp[k][x][y][xx][yy] = min(ans,dp[k][x][y][xx][yy]);
    }

    for(j = y;j <= yy-1;j++)
    {
        ans = min(deal(k+1,x,y,xx,j)+sum[x][j+1][xx][yy],deal(k+1,x,j+1,xx,yy)+sum[x][y][xx][j]);
        dp[k][x][y][xx][yy] = min(ans,dp[k][x][y][xx][yy]);
    }

    return dp[k][x][y][xx][yy];
}
int main()
{
    //int n;
    //
   // freopen("in.txt","r",stdin);
    while(~scanf("%d",&n))
    {
        int i,j;
        average = 0;
        loop(1,i,9)
        {
            loop(1,j,9)
                scanf("%d",map[i]+j),average += map[i][j];
        }

        average = 1.0*average/n;
        cl(sum,0);
        cl(dp,-1);
        init();
        deal(1,1,1,8,8);
       // puts("yes");

        double ans = sqrt(dp[1][1][1][8][8]*1.0/n-average*average);
        printf("%.3f",ans);
    }
    return 0;
}