【IOI2000第五题】邮局问题

描述

一些村庄建在一条笔直的高速公路边上，我们用一条坐标轴来描述这条公路，每个村庄的坐标都是整数，没有两个村庄的坐标相同。两个村庄的距离定义为坐标之差的绝对值。我们需要在某些村庄建立邮局。使每个村庄使用与它距离最近的邮局，建立邮局的原则是：所有村庄到各自使用的邮局的距离总和最小。
数据规模：1<=村庄数<=300, 1<=邮局数<=30, 1<=村庄坐标<=10000
格式

输入格式

2行
第一行:n m {表示有n个村庄，建立m个邮局}
第二行:a1 a2 a3 .. an {表示n个村庄的坐标}
输出格式

1行
第一行:l {l表示最小距离总和}
样例1

样例输入1

10 5
1 2 3 6 7 9 11 22 44 50

样例输出1

9

前言：
前天晚上loli讲过一个问题，有几头牛站成一排，所有的牛都去一个地方集合，问在哪个地方集合所有的牛跑的距离最短。
答案是选中位数——都去中间标号的牛那里去集合（注意不是中间位置），可以证明：如果不选在中间，超过一半的牛的跑的距离在增加，小于一半的牛的数量在减少，肯定不是最优。
如果每个位置有多头牛，那就是加权中位数（数量中间的那头牛），证明略。

所以这道题选邮局就一定选在中间的村庄。
分析：
本题是一道线性动归。首先我们考虑可以参与划状态的变量：（1）所服务的村庄数（2）所设立的邮局数。于是我们就考虑到子状态为在没有放置这个邮局之前的一个子区间的最小距离数值，如果加上在来一个邮局所需要的区间距离和就是一个方案。其中我们把以前几个邮局所服务的区间看做是一个独立的区间，而新加区间又是独立的，所以满足无后效性，处理后每一个最小的数值都可以看做上述划分，且上一个扩充的来的子区间已是最优，所以满足最优子结构。

方程：
f[i][j]表示前i个村庄中放j个邮局的最小路程 。data（i，j）表示第i个村庄到第j个村庄之间建一个邮局耗费的最小距离总和（中位数）。

f[i][j]=min(f[i][j],f[k][j-1]+data(k+1,i));

这一类的线性dp就是把最后一个拿出来讨论。
初始化：极大值
边界：f[i][1]=data(1,i); （1<=i<=n）
输出：f[n][m]

code:

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
int n,m,a[301];
int f[301][31];

int Abs(int x) {return x<0?-x:x;}//原来是abs函数写错了！ 

int data(int x,int y)
{
    if (x==y) return 0;
    int mid=(x+y)>>1; int sum=0;
    for (int i=x; i<=y; i++)
        sum+=Abs(a[i]-a[mid]);
    return sum;
}

int main()
{
    scanf("%d%d",&n,&m);
    for (int i=1; i<=n; i++)
        scanf("%d",&a[i]);
    memset(f,0x3f,sizeof(f));
    for (int i=1; i<=n; i++)
        f[i][1]=data(1,i);
    for (int i=2; i<=n; i++)
        for (int j=1; j<=m; j++)
        {
            for (int k=j; k<=i-1; k++)//枚举断点 
            f[i][j]=min(f[i][j],f[k][j-1]+data(k+1,i));//这类方程就是看把最后一个分给谁的问题 
        }
/*  for (int i=1; i<=n; i++)
    {
        for (int j=1; j<=m; j++) printf("%d ",f[i][j]);
        printf("\n");
    }*/
    printf("%d",f[n][m]);
    return  0;
}
/*
10 5
1 2 3 6 7 9 11 22 44 50

9
*/

PS 好像可以用四边形不等式优化？