POJ1160 Post Office 四边形不等式优化

　　题目链接：http://poj.org/problem?id=1160

　　f[i][j]表示前 i 个村庄放 j 个邮局的最小花费，w[i][j]表示 i-j 村庄之间放一个邮局的最小花费，则转移方程：f[i][j]= Min{ f[k][j-1] + w[k+1][i] }，把f[i][j]转换一下，表示前 j 个村庄放 i 个邮局的最小花费，则 f[i][j]=Min{ f[i-1][k] + w[k+1][i] }。

　　四边形不等式优化有如下定理：

　　　　1.当决策代价函数w[i][j]满足w[i][j]+w[i’][j’]<=w[I;][j]+w[i][j’](i<=i’<=j<=j’)时,称w满足四边形不等式.当函数w[i][j]满足w[i’][j]<=w[i][j’] i<=i’<=j<=j’)时,称w关于区间包含关系单调.

　　　　2.如果状态转移方程m为且决策代价w满足四边形不等式的单调函数(可以推导出m亦为满足四边形不等式的单调函数),则可利用四边形不等式推出最优决策s的单调函数性,从而减少每个状态的状态数,将算法的时间复杂度由原来的O(n^3)降低为O(n^2).方法是通过记录子区间的最优决策来减少当前的决策量.令:

s[i][j]=max{k | ma[i][j] = m[i][k-1] + m[k][j] + w[i][j]}

由于决策s具有单调性,因此状态转移方程可修改为:

　　因此只要证明w满足凸和包含关系就可以用四边形不等式优化了。

//STATUS:C++_AC_16MS_584KB
#include<stdio.h>
#include<string.h>
#define min(a,b) ( (a)<(b)?(a):(b) )
const int MAX1=310,MAX2=35,INF=10000000;
int dp[MAX1][MAX2],sum[MAX1][MAX1],p[MAX1],P,V;
int main()
{
//    freopen("in.txt","r",stdin);
    int i,j,k,min;
    while(~scanf("%d%d",&V,&P))
    {
        min=0x7fffffff;
        memset(sum,0,sizeof(sum));
        for(i=1;i<=V;i++)
            scanf("%d",&p[i]);

        for(i=1;i<V;i++){
            for(j=i+1;j<=V;j++){
                sum[i][j]=sum[i][j-1]+p[j]-p[(i+j)/2];
            }
        }
        for(i=1;i<=V;i++){
            dp[i][1]=sum[1][i];
            dp[i][i]=0;
        }

        for(j=2;j<=P;j++){
            for(i=j+1;i<=V;i++){
                dp[i][j]=INF;
                for(k=j-1;k<i;k++)
                     dp[i][j]=min(dp[i][j],dp[k][j-1]+sum[k+1][i]);
            }
        }

        printf("%d\n",dp[V][P]);
    }
    return 0;
}