poj1160 Post Office（区间dp+四边形不等式优化）

题目链接

分析：
一看这道题，就有一个疑问进入我的脑海：
我们要怎么确定一个村庄到底去哪个邮局呢？

实际上邮局的存在就把村庄分成了若干区间
一个区间内的村中都去一个邮局
如果我们确定了区间，那么要找到区间最优解是不难的（枚举暴力计算）

也就是说整个序列被我们划分成了P个区间，每个区间内的最优解我们可以预处理出
所以设计状态：f[i][j]表示前i个村庄被划分成了j个区间（建了j个邮局）
转移方程：f[i][j]=min{f[k][j-1]+c[k+1][i]}
c[i][j]表示在区间(i,j)内建一个邮局的最小代价

c的预处理要注意一下（推柿子就好了）

//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int INF=0x33333333;
int dis[305],sum[305];
int f[305][35],c[305][305];
int n,m;

void pre()
{
    int i,j,k;
    for (i=1;i<=n;i++)
    {
        c[i][i]=0;
        for (j=i+1;j<=n;j++)
        {
            int s=sum[j]-sum[i]-dis[i]*(j-i);
            c[i][j]=s;
            for (k=i+1;k<=j;k++)
            {
                s=s+(2*k-i-j-1)*(dis[k]-dis[k-1]);
                c[i][j]=min(c[i][j],s);
            }
        }
    }
}

void doit()
{
    int i,j,k;
    memset(f,0,sizeof(f));
    for (int i=1;i<=n;i++) f[i][1]=c[1][i];
    for (i=2;i<=n;i++)
        for (j=2;j<=min(i,m);j++)
        {
            f[i][j]=INF;
            for (k=1;k<i;k++)
                f[i][j]=min(f[i][j],f[k][j-1]+c[k+1][i]);
        }
    printf("%d\n",f[n][m]);
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
        scanf("%d",&dis[i]);
    int a=dis[1];
    for (int i=1;i<=n;i++) 
    {
        dis[i]-=a;                //每个点到第一个点的距离
        sum[i]=sum[i-1]+dis[i];   //dis的前缀和
    }
    pre();
    doit();
    return 0;
}

---

之后我又发现c数组的计算完全不用这么麻烦，
就像UVa11300中说的，取最优值得位置一定是区间的中点
所以我们实际上可以O(1)出解

但是这道题还没完，因为我们现在的效率是O(n^3)的
但是我们可以利用四边形不等式优化到O(n^2)

//这里写代码片
#include<cstdio>
#include<cstring>
#include<iostream>

using namespace std;

const int INF=0x33333333;
int dis[305],sum[305];
int f[305][35],c[305][305],s[305][305];
int n,m;

int get(int i,int j)
{
    if ((i+j)%2==0) return sum[j]-sum[(i+j)/2]-sum[(i+j)/2-1]+sum[i-1];
    else return sum[j]-sum[(i+j)/2]-sum[(i+j)/2-1]+sum[i-1]-dis[(i+j)/2];
}

void pre()
{
    int i,j,k;
    for (i=1;i<=n;i++)
    {
        c[i][i]=0;
        for (j=i+1;j<=n;j++)
            c[i][j]=get(i,j);
    }
}

void doit()
{
    int i,j,k;
    memset(f,0x33,sizeof(f));
    for (i=1;i<=n;i++) f[i][1]=c[1][i];

    for (j=2;j<=m;j++)
    {
        s[n+1][j]=n-1;
        for (i=n;i>=j;i--)
            for (k=s[i][j-1];k<=s[i+1][j];k++)
            if (f[i][j]>f[k][j-1]+c[k+1][i])
            {
                f[i][j]=f[k][j-1]+c[k+1][i];
                s[i][j]=k;
            }
    }

    printf("%d\n",f[n][m]);
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&dis[i]);
        sum[i]=sum[i-1]+dis[i];
    }   
    pre();
    doit();
    return 0;
}