bzoj1505: [NOI2004]小H的小屋

传送门
用$g[i][j]$表示南方向上，x轴方向上长度为i的区间，划分成j个矩形的最小面积。
用$f[k][i][j]$表示x轴方向上长度为k的区间，北方向上放i个矩形，南方向上放j个矩形的最小面积。
可以得到下面的dp方程
$g[i][j]=min(g[i′][j−1]+(i−i′)^2K2)$
$f[k][i][j]=min(f[k′][i−1][j′]+(k−k′^)2K1+g[k−k′][j−j′])$
求$g[][]$的过程复杂度是$O(n^3)$
，但求$f[][][]$的复杂度是$O(n^5)$，太慢了。
我们考虑优化$k',j'$枚举范围
可以发现，假设$f[k][i][j]$是从$f[k1][i−1][j1]$递推而来，$f[k][i−1][j]$是从$f[k2][i−2][j2]$递推而来，那么有$k1>k2,j1>j2。$
这是显然的，同样长度的区间，划分的话，一定是划分得尽量均匀，让每一段矩形的大小尽量接近。那么显然划分成i段，北面最后一段矩形的长度比划分成$i−1$段，北面最后一段矩形的长度短。
利用这个特点，每次DP枚举$k′,j′$时，我们只需要从$k2,j2$开始枚举即可。这样复杂度就降到了$O(n^3)$

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cstdlib>
#define N 105
#define inf 1e100
using namespace std;
double g[N][N],f[N][N][N],k1,k2;
int lastj[N][N][N],lastk[N][N][N],n,m;
int main(){
    scanf("%lf%lf%d%d",&k1,&k2,&m,&n);
    for (int i=0;i<N;i++)
        for (int j=0;j<N;j++) g[i][j]=inf;
    for (int i=0;i<N;i++)
        for (int j=0;j<N;j++)
            for (int k=0;k<N;k++) f[i][j][k]=inf;
    g[0][0]=0;
    for (int i=1;i<=100;i++)
        for (int j=1;j<=i;j++)
            for (int k=j-1;k<i;k++)
                g[i][j]=min(g[i][j],g[k][j-1]+(i-k)*(i-k)*k2);
    f[0][0][0]=0;
    for (int k=1;k<=100;k++)
        for (int i=1;i<=m;i++)
            for (int j=1;j<=n;j++)
                for (int lk=lastk[k][i-1][j];lk<k;lk++)
                    for (int lj=lastj[k][i-1][j];lj<j;lj++)
                        if (f[k][i][j]>f[lk][i-1][lj]+(k-lk)*(k-lk)*k1+g[k-lk][j-lj]){
                            f[k][i][j]=f[lk][i-1][lj]+(k-lk)*(k-lk)*k1+g[k-lk][j-lj];
                            lastk[k][i][j]=lk;
                            lastj[k][i][j]=lj;
                        }
    printf("%.1lf",f[100][m][n]);
}