dp四边形优化 Hdu 3480 Division 题解

累加器传送门：

http://blog.csdn.net/NOIAu/article/details/71775000

---

题目传送门：

https://cn.vjudge.net/problem/HDU-3480

---

dp状态转移分析博客：

http://blog.csdn.net/NOIAu/article/details/72428364

---

由于在斜率优化的博客里已经讲解了如何dp转移，这里就不在赘述，显然cost[i][j]=（a[j]-a[i]）^2 满足凸性和包含原则（排序后），所以直接四边形优化，如果需要看dp转移，请进入上述的dp转移分析博客

---

代码如下：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>

#define MAXN 10000+10
int T;
int n,m,cn;
int a[MAXN],cnt[MAXN];
int dp[MAXN][MAXN],s[MAXN][MAXN];

int sqr(int i){ return i*i; }

using namespace std;

void init(){
    scanf("%d%d",&n,&m);
    for(register int i=1;i<=n;i++) scanf("%d",&a[i]);
    sort(a+1,a+n+1);
}

void dpp(){
    for(register int i=1;i<=n;i++){
        dp[1][i]=sqr(a[i]-a[1]);
        s[1][i]=1;
    }
    for(register int i=2;i<=m;i++){
        s[i][n+1]=n-1;
        for(register int j=n;j>=i;j--){
            int temp=0x7fffffff;
            int te;
            for(int k=s[i-1][j];k<=s[i][j+1];k++){
                if(temp>dp[i-1][k]+sqr(a[j]-a[k+1])){
                    temp=dp[i-1][k]+sqr(a[j]-a[k+1]);
                    te=k;
                }
            }
            dp[i][j]=temp;
            s[i][j]=te;
        }
    }
    printf("Case %d: %d\n",++cn,dp[m][n]);
}

int main(){
    scanf("%d",&T);
    while(T--){
        init();
        dpp();
    }
}