POJ 1191 棋盘分割（区间DP）题解

#include<cmath>
#include<stack>
#include<cstdio>
#include<vector>
#include<cstring>
#include <iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int maxn = 10 + 10;
const int INF = 0x3f3f3f3f;
const int MOD = 1000000007;
int n;
double w[maxn][maxn], dp[maxn][maxn][maxn][maxn][maxn], sum[maxn][maxn];
double get(int x1, int y1, int x2, int y2){
    return sum[x2][y2] - sum[x2][y1 - 1] - sum[x1 - 1][y2] + sum[x1 - 1][y1 - 1];
}
int main(){
    scanf("%d", &n);
    memset(sum, 0, sizeof(sum));
    for(int i = 1; i <= 8; i++){
        for(int j = 1; j <= 8; j++){
            scanf("%lf", &w[i][j]);
            sum[i][j] = sum[i - 1][j] + sum[i][j - 1] - sum[i - 1][j - 1] + w[i][j];
        }
    }
    double per = sum[8][8] / n;
    for(int x1 = 1; x1 <= 8; x1++){
        for(int y1 = 1; y1 <= 8; y1++){
            for(int x2 = x1; x2 <= 8; x2++){
                for(int y2 = y1; y2 <= 8; y2++){
                    double ret = get(x1, y1, x2, y2);
                    dp[x1][y1][x2][y2][1] = ret * ret;
                }
            }
        }
    }
    for(int k = 2; k <= n; k++){
        for(int x1 = 1; x1 <= 8; x1++){
            for(int y1 = 1; y1 <= 8; y1++){
                for(int x2 = x1; x2 <= 8; x2++){
                    for(int y2 = y1; y2 <= 8; y2++){
                        dp[x1][y1][x2][y2][k] = INF;
                        for(int t = x1; t < x2; t++){
                            dp[x1][y1][x2][y2][k] = min(dp[x1][y1][x2][y2][k], dp[x1][y1][t][y2][1] + dp[t + 1][y1][x2][y2][k - 1]);
                            dp[x1][y1][x2][y2][k] = min(dp[x1][y1][x2][y2][k], dp[x1][y1][t][y2][k - 1] + dp[t + 1][y1][x2][y2][1]);
                        }
                        for(int t = y1; t < y2; t++){
                            dp[x1][y1][x2][y2][k] = min(dp[x1][y1][x2][y2][k], dp[x1][y1][x2][t][1] + dp[x1][t + 1][x2][y2][k - 1]);
                            dp[x1][y1][x2][y2][k] = min(dp[x1][y1][x2][y2][k], dp[x1][y1][x2][t][k - 1] + dp[x1][t + 1][x2][y2][1]);
                        }
                    }
                }
            }
        }
    }
    printf("%.3lf\n", sqrt(dp[1][1][8][8][n] / n - per * per));
    return 0;
}