hdu3516 Tree Construction（四边形不等式优化dp）

题目链接

题目描述：给定一些点$(x_i,y_i)$满足：$i<j,x_i<x_j,y_i>y_j$
用向右和向上的边连起来，使得所有边的长度最小
（从整棵树的最低点开始的所有边的方向只有向左或向右）

分析：
点的分布呈下降趋势
显然第一个用向上的边连，最后一个点用向右的边连
其实每个点的的花费只有两种：纵坐标-1，与前一个点的横坐标之差

---

一开始不成熟的想法：
如果我们知道了向上的边的位置
那么在两个向上边之间的点肯定是连向ta之前的上边

于是就设计了一个状态：
$f[i]$表示第$i$个点用向上的边连接
设$sum[i]=\sum_{j=1}^{i}x_j$
$f[i]=min(f[j]+(sum[i-1]-sum[j])-(i-j-1)*x_j)+y_i-y_n$
（一号点和n号点需要特殊处理一下）

很快就写好了
然而WA了。。。

找到反例了。。。

---

实际上这道题是区间dp

我们可以把整棵树视为两棵子树的合并
$f[i][j]$表示$i-j$的点全部连通的最小花费
$f[i][j]=min(f[i][k]+f[k+1][j]+abs(y[k]-y[j])+abs(x[k+1]-x[i]))$

$abs(y[k]-y[j])+abs(x[k+1]-x[i])$是连接两棵树的代价，
由题目性质决定，得到的树一定是一个直角三角形

式子符合四边形不等式的限制
我们可以用四边形不等式优化

• 对于函数$w[i][j]$，若$w[i][j]+w[i+1][j+1]<=w[i][j+1]+w[i+1][j]$，则$w$满足凸多边形不等式（交叠的小于等于包含）
• 对于函数$w[i][j]$，若$w[a,b]<=w[c,d](c<=a<=b<=d)$，那么我们称$w$关于区间包含状态单调

定理1：$w$满足凸多边形不等式和区间包含状态单调，那么dp也满足四边形不等式
定理2：最优决策$k[i][j]$满足$k[i][j−1]<=k[i][j]<=k[i+1][j]$，这是四边形不等式优化dp的关键，利用这个定理可以每次缩小$k$的枚举范围，可以使时间复杂度从$O(n^3)$降到$O(n^2)$
定理3：$w$为凸当且仅当$w[i][j]+w[i+1][j+1]<=w[i+1][j]+w[i][j+1]$

tip

第一次自己写四边形不等式优化，感天动地

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

using namespace std;

const int N=1002;
int n,x[N],y[N],f[N][N],s[N][N];

int dis(int i,int k1,int k2,int j) {
    return abs(y[k1]-y[j])+abs(x[k2]-x[i]);
}

int main() {
    while (scanf("%d",&n)!=EOF) {
        for (int i=1;i<=n;i++) 
            scanf("%d%d",&x[i],&y[i]);
        memset(f,0x33,sizeof(f));
        for (int i=1;i<=n;i++) {
            f[i][i]=0;
            s[i][i]=i;
        }
        for (int i=n;i>=1;i--) {
            s[i+1][n]=n;
            for (int j=i+1;j<=n;j++) 
                for (int k=s[i][j-1];k<=s[i+1][j];k++) {
                    int t=f[i][k]+f[k+1][j]+dis(i,k,k+1,j);
                    if (f[i][j]>t) {
                        f[i][j]=t;
                        s[i][j]=k;
                    }
                }
        }
        printf("%d\n",f[1][n]);
    }
    return 0;
}