HDU 3516 Tree Construction (四边形优化DP)

题目大意：给出n个点的坐标，对于 i < j 满足x[i] < x[j]，y[i] > y[j]。问如何将这些点连起来，线只能向上或者向右。求线最短是多大。

定义状态dp[i][j]表示从点i到点j，连起来需要线的最小值。则得到转移方程为：

dp[i][j] = min(dp[i][k]+dp[k+1][j]+cost[i][j]);其中i<k<j，转移量为O(n)，状态为O(n^2)，总复杂度为O(n^3)，TLE。

考虑到该式子跟一般四边形优化的DP是一样的，所以只要证明cost[i][j]是单调递减的即可。

而cost[i][j] = y[k]-y[j] + x[k+1]-x[i]；固定j，有cost[i+1][j]-cost[i][j] = x[i]-x[i+1]；由于对于i<j有x[i]<x[j]，所以x[i]-x[i+1]<0，得证cost[i][j]是单调递减的。故cost[i][j]是凸的。再证明对于[i,j]属于[i',j']，有cost[i][j] <= cost[i'][j']；考虑式子cost[i][j] = y[k]-y[j] + x[k+1]-x[i];则可以直接比较x[i]+y[j]，若x[i]+y[j] > x[i']+[j']则说明了假设是得证的。再看条件i<j有x[i]<x[j],y[i]>y[j]，而i' <= i <= j <= j'，得x[i'] < x[i], y[j]' < y[j]，得证原假设。

然后直接四边形优化了。

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

using namespace std;

const int inf = 0xfffffff;
const int SIZE = 1024;
int n;
int x[SIZE],y[SIZE];
int dp[SIZE][SIZE],s[SIZE][SIZE];

int Abs(int v)
{
    return v>0?v:(-v);
}

int getDis(int s,int e)
{
    return Abs(x[s]-x[e])+Abs(y[s]-y[e]);
}

int main()
{
    while(~scanf("%d",&n))
    {
        for(int i=1; i<=n; i++)
            scanf("%d%d",&x[i],&y[i]);
        for(int i=0; i<SIZE; i++)
            for(int k=0; k<SIZE; k++)
                dp[i][k] = inf,s[i][k] = 0;
        for(int i=1; i<=n; i++)
        {
            dp[i][i] = 0;
            s[i][i+1] = i;
            if(i < n)
                dp[i][i+1] = Abs(x[i]-x[i+1])+Abs(y[i]-y[i+1]);
        }
        for(int p=2; p<n; p++)
        {
            for(int i=1; i+p<=n; i++)
            {
                int r = i+p;
                for(int k=s[i][r-1]; k<=s[i+1][r]; k++)
                {
                    int tem = dp[i][k]+dp[k+1][r]+getDis(i,r)-y[i]+y[k]-x[r]+x[k+1];
                    if(tem < dp[i][r])
                    {
                        dp[i][r] = tem;
                        s[i][r] = k;
                    }
                }
            }
        }
        printf("%d\n",dp[1][n]);
    }
    return 0;
}