HDU3480：Division——题解

http://acm.hdu.edu.cn/showproblem.php?pid=3480

将一列数划分成几个集合，这些集合的并集为该数列，求每个数列的（最大值-最小值）^2的和的最小值。

简单的dp都会写，就不讲了。

然后就是四边形优化了，参考：https://blog.csdn.net/noiau/article/details/72514812

事实上四边形优化的条件一般是靠打表打出来的。

于是简单记录下吧：

先排序。

设dp[i][j]为前j个数划分成i个集合的最小值，cost[i][j]为i～j的集合价值。

显然有dp[i][j]=min{dp[i][j],dp[i][k]+cost[k+1][j]}

接着打表得出（就是打一个矩阵，观察矩阵每行每列都是递增的）：

s[i-1][j]<=s[i][j]<=s[i][j+1]

然后就可以利用第三条结论来优化了。

（此外能否用四边形不等式优化还和你如何定义dp也是有关系的……我就是被坑了把dp两个状态倒换一下才行。）

还有一些注意事项看一下https://www.cnblogs.com/mlystdcall/p/6525962.html吧。

#include<cstdio>
#include<iostream>
#include<vector>
#include<queue>
#include<cstring>
#include<algorithm>
#include<map>
using namespace std;
const int N=;
const int M=;
const int INF=1e9;
inline int read(){
    int X=,w=;char ch=;
    while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
    while(isdigit(ch))X=(X<<)+(X<<)+(ch^),ch=getchar();
    return w?-X:X;
}
inline int sqr(int k){return k*k;}
int a[N],dp[M][N],s[M][N];
int main(){
    int t=read();
    for(int cas=;cas<=t;cas++){
    printf("Case %d: ",cas);
    int n=read(),m=read();
    for(int i=;i<=n;i++)a[i]=read();
    sort(a+,a+n+);
    for(int i=;i<=n;i++){
        dp[][i]=sqr(a[i]-a[]);
        s[][i]=;
    }
    for(int i=;i<=m;i++){
        s[i][n+]=n-;
        for(int j=n;j>=i;j--){
        dp[i][j]=INF;
        for(int k=s[i-][j];k<=s[i][j+];k++){
            if(dp[i][j]>dp[i-][k]+sqr(a[j]-a[k+])){
            dp[i][j]=dp[i-][k]+sqr(a[j]-a[k+]);
            s[i][j]=k;
            }
        }
        }
    }
    printf("%d\n",dp[m][n]);
    }
    return ;
}

+++++++++++++++++++++++++++++++++++++++++++

+本文作者：luyouqi233。 　　　　　　　　　　　　　　+

+欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++