poj 1191 棋盘分割 记忆化搜索

纠结了一天，奇葩的代码，dp初始化的时候，用memset(dp, 0x3f, sizeof(dp)),记忆化搜索判断用 if( dp[tim][x1][y1][x2][y2] == 0x3f3f3f3f),TLE了整整一晚上，后来随便换成用-1做判断条件就A掉了，还0ms……又去试了下是不是poj里面dp[tim][x1][y1][x2][y2] == 0x3f3f3f3f会出问题，用a+b代码试了下，完全没问题啊……

还是不知道到底哪挫了…………………………

白天写的时候，没有化简方差公式，用一个area存每一块的面积，用一个th数组存储dfs累加的 (xi - avg)*(xi - avg) 的最小值，不知道哪有问题，写了好多遍也不对，后来看黑书上，按化简的方差公式写，过样例，然后就开始上面的 endless TLE …………

很简单的记忆化搜索

for( k= x1; k<x2; k++){

mm =min( mm, dfs( tim-1, x1, y1, k, y2) + get_var( k+1, y1, x2, y2));

mm =min( mm, dfs( tim-1, k+1, y1, x2, y2) + get_var( x1, y1, k, y2));

}

for( k= y1; k<y2; k++){

mm =min( mm, dfs( tim-1, x1, y1, x2, k) + get_var( x1, k+1, x2, y2));

mm =min( mm, dfs( tim-1, x1, k+1, x2, y2) + get_var( x1, y1, x2, k));

}

dp[tim][x1][y1][x2][y2] = mm;

把目前的这块（范围 x1, y1, x2, y2 ) 分为两部分， 一部分算分出去的，计算其大小，供最后处理方差，另一部分dfs，继续进行分割

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring>
#include <cmath>
using namespace std;
const int INF= 0x3f3f3f3f;
int dp[16][9][9][9][9], n, mit[9][9], ans[9][9], area[9][9][9][9];
int get_var(int x1, int y1, int x2, int y2){
if( area[x1][y1][x2][y2] != -1) return area[x1][y1][x2][y2];
int tmp = ans[x2][y2] - ans[x1-1][y2] - ans[x2][y1-1] + ans[x1-1][y1-1];
return area[x1][y1][x2][y2] = tmp*tmp;
}
int dfs(int tim, int x1, int y1, int x2, int y2){
if( dp[tim][x1][y1][x2][y2] != -1) return dp[tim][x1][y1][x2][y2] ;
if( tim == 0){
dp[tim][x1][y1][x2][y2] = get_var( x1, y1, x2, y2);
return dp[tim][x1][y1][x2][y2] ;
}
int k, s, s1, s2, mm= INF;
for( k= x1; k<x2; k++){
mm =min( mm, dfs( tim-1, x1, y1, k, y2) + get_var( k+1, y1, x2, y2));
mm =min( mm, dfs( tim-1, k+1, y1, x2, y2) + get_var( x1, y1, k, y2));
}
for( k= y1; k<y2; k++){
mm =min( mm, dfs( tim-1, x1, y1, x2, k) + get_var( x1, k+1, x2, y2));
mm =min( mm, dfs( tim-1, x1, k+1, x2, y2) + get_var( x1, y1, x2, k));
}
return ( dp[tim][x1][y1][x2][y2] = mm);
}
int main(){
//  freopen("1.txt", "r", stdin);
int i, j, k;
int aa[2][2][2][2][2];
double avg, ss;
while( scanf("%d", &n) != EOF){
memset( ans, 0, sizeof(ans));
memset( area, -1, sizeof( area));
for( i= 1; i< 9; i++){
for( j= 1; j< 9; j++){
scanf("%d", &mit[i][j]);
ans[i][j]= ans[i-1][j] + ans[i][j-1] + mit[i][j] - ans[i-1][j-1];
}
}
avg= ans[8][8]*1.0/n;
memset(dp, -1, sizeof( dp));
dfs(n-1, 1, 1, 8, 8);
memset(aa, INF, sizeof( INF));
ss= sqrt(fabs(dp[n-1][1][1][8][8]*1.0/n - avg*avg));
printf("%.3f\n", ss);
}
return 0;
}