POJ3154 Graveyard

　　原题传送：http://poj.org/problem?id=3154

　　模型转换一下，把圆形的路从某个雕塑处剪开拉成一条直线，这样就很容易处理出添加雕塑前后的位置数组a[N]和c[N+M]。对于原状态的一个雕塑a[i]，在新位置数组c[N+M]中找出小于等于a[i]的位置c[j]，然后取位置j和其相邻两点j-1、j+1(如果有的话)距离的最小值(a[i] - min(c[j-1], c[j], c[j+1]))为a[i]须移动的最短距离。

　　原来想的时候会顾虑这样一个问题：在原来的位置数组a[N]中会不会出现相邻的两个点a[i],a[i+1]同时和c[j]距离最近，也就是说上述思想会不会导致a[i]和a[i+1]移到同一个位置？其实是不会出现这种情况的，可以画个数轴分析一下，如果雕塑增多了，那么原来雕塑肯定在新位置数组的两点之间，如果同时位于中点，上述分析显然成立，如果不同时位于中点，那么，便会偏向同一侧而不会取这两点中间的c[j]。

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <algorithm>
#define N 1005
using namespace std;

int n, m;
double a[N], c[N * 2];

int main()
{
    int i;
    double j, d, ans;
    while(scanf("%d%d", &n, &m) == 2)
    {
        d = 10000.0 / n;
        for(j = 0.0, i = 0; i < n; i ++, j += d)
            a[i] = j;
        d = 10000.0 / (n + m);
        for(j = 0.0, i = 0; i < n + m; i ++, j += d)
            c[i] = j;

        double k;
        for(ans = 0.0, i = 0; i < n; i ++)
        {
            int loc = lower_bound(c, c + n + m, a[i]) - c;
            k = fabs(a[i] - c[loc]);
            if(loc > 0)
                k = min(k, fabs(a[i] - c[loc - 1]));
            if(loc < n + m - 1)
                k = min(k, fabs(a[i] - c[loc + 1]));
            ans += k;
        }
        printf("%.4f\n", ans);
    }
    return 0;
}