Luogu P1436 棋盘分割 暴力DP

我的天，，，，，n=8,k<=15，，，这怕不是暴力DP+高维数组。。。。

---

开一个五维数组f[k][i][j][p][q]表示从(i,j)到(p,q)中分成k个矩形最小的平方和。

然后初始化时用上二维前缀和

DP时有些像区间DP，枚举分割线，分成两部分，取min就好

#include<cstdio>
#include<iostream>
#include<cstring>
#define R register int
#define For(i,a,b) for(R i=a;i<=b;++i)
#define ffor(i,a,b) for(R i=b;i>=a;--i)
const int N=;
using namespace std;
inline int g() {
    R ret=,fix=; register char ch; while(!isdigit(ch=getchar())) fix=ch=='-'?-:fix;
    do ret=ret*+(ch^); while(isdigit(ch=getchar())); return ret*fix;
}
inline void ckmn(int& a,int b) {b<a?a=b:a=a;}
int n;
int f[N<<][N][N][N][N],a[N][N],s[N][N];
signed main() { memset(f,0x3f,sizeof(f)); n=g();
    For(i,,) For(j,,) a[i][j]=g();
    For(i,,) For(j,,) s[i][j]=s[i-][j]+s[i][j-]-s[i-][j-]+a[i][j];
    For(i,,) For(j,,) For(p,i,) For(q,j,) {
        R tmp=s[p][q]-s[p][j-]-s[i-][q]+s[i-][j-]; f[][i][j][p][q]=tmp*tmp;
    } For(tk,,n) ffor(i,,) ffor(j,,) For(p,i,) For(q,i,) {
        For(t,i+,p) ckmn(f[tk][i][j][p][q],min(f[tk-][i][j][t-][q]+f[][t][j][p][q],f[][i][j][t-][q]+f[tk-][t][j][p][q]));
        For(t,j+,q) ckmn(f[tk][i][j][p][q],min(f[tk-][i][j][p][t-]+f[][i][t][p][q],f[][i][j][p][t-]+f[tk-][i][t][p][q]));
    } printf("%d\n",f[n][][][][]);
}

第一次define for循环，太多了。。QAQ。。。

---

2019.04.28