POJ1160 Post Office 四边形不等式优化

题目大意
在一条直线上有n个村庄，选出m个村庄，在其中每个村庄建立一个邮局，要求每个村庄到最近邮局的距离和最小。

f[m,n]为问题的解
f[i,j]表示在前j个村庄建立i个邮局的最小代价
状态转移方程和边界条件f[1,j]=w[1,j]
m[i,j]=min(m[i-1,k]+w[k+1,j]) i<=j

其中w[i,j]表示在d[i..j]之间建立一个邮局的最小距离和，
可以证明，当仅建立一个邮局时，最优解出现在中间，
即设建立邮局的村庄为k，则k=(i+j)/2

证明函数w满足四边形不等式，即
w[i,j]+w[i’,j’]<=w[i’,j]+w[i,j’] , i<=i’<=j<=j’

算法降次优化 s[i-1,j]<=s[i,j]<=s[i,j+1]

#include <iostream>  
#include <cstdio>
#include <algorithm>
#include <cstring>
using namespace std; 

int f[35][305], dis[305], w[305][305], s[35][305];
int n, m;

int main(){
    while(scanf("%d%d", &n, &m) == 2){
        for(int i=1; i<=n; i++){
            scanf("%d", &dis[i]);
            w[i][i]=0;//从第i个村庄到第j个村庄中只建一个邮局的代价 
        }
        for(int i=1; i<=n; i++)
            for(int j=i; j<=n; j++)
                w[i][j] = w[i][j-1] + dis[j] - dis[(i+j)>>1];//邮局一定建在区间中点 
        memset(f,0x7f,sizeof(f));
        for(int i=1; i<=n; i++){
            f[1][i] = w[1][i];//f[i][j]表示在1~j中建i个邮局的代价 
            s[1][i] = 0;//边界条件 
        }
        for(int i=2; i<=m; i++){  
            s[i][n+1]=n;//边界条件 
            for(int j=n; j>i; j--){//因为范围是s[i-1][j]~s[i][j+1]，所以j要反向遍历 
                for(int k=s[i-1][j]; k<=s[i][j+1]; k++)  {  
                    if(f[i-1][k] + w[k+1][j] < f[i][j]){
                        s[i][j] = k;
                        f[i][j] = f[i-1][k] + w[k+1][j]; 
                    }  
                }  
            }  
        } 
        printf("%d\n",f[m][n]);
    }
    return 0;
}