居民迁移 HihoCoder - 1053

公元2411年，人类开始在地球以外的行星建立居住点。在第1326号殖民星上，N个居住点分布在一条直线上。为了方便描述，我们设第i个居住点的位置是Xi，其中居住着Yi位居民。随着冬季的到来，一些人口较多的居住点的生态循环系统已经开始超负荷运转。为了顺利度过严冬，殖民星上的居民一致同意通过转移到人口较少的居住点来减轻人口众多的居住点的负荷。

遗憾的是，1326殖民星的环境非常恶劣。在冬季到来前，每个居民点的居民最远能迁移到距离不超过R的居民点。1326殖民星的居民希望知道，如何安排迁移才能使完成迁移后人口最多的居民点人口最少？

注意有可能存在多个居民点位置相同。

Input

第一行包含一个整数T(1 <= T <= 10)，代表测试数据的组数。

每组数据的第一行包含2个整数N(1 <= N <= 100000)和R(0 <= R <= 10^9)。

以下N行每行包含两个整数，Xi和Yi(0 <= Xi, Yi, <= 10^9)。

Output

对于每组数据输出迁移后人口最多的居民点人口最少可能的数目。

Sample Input

3  
5 1  
10 80  
20 20  
30 100  
40 30  
50 10  
5 10  
10 80  
20 20  
30 100   
40 30  
50 10  
5 20  
10 80  
50 10  
20 20  
30 100  
40 30 

Sample Output

100  
50  
48

思路

这是一道典型，经典的双指针下标题+二分题。
由于题目中有着

对于每组数据输出迁移后人口最多的居民点人口最少可能的数目。

所以易想到二分答案。
而在check()里，我们用i与j下标来分别存人与城市的下标。
b[]来存人，c[]来存城市。
这就如同将每个城市的人都全部迁出来放进b里。
则此时c[]，即每个城市都是空城，即值为0。
而i之前的b[]都成了0，j之前的c[]都成了mid。

代码

#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
const int MAXN=100000;
int K,n,R;
long long l,r,mid;
long long b[MAXN+5],c[MAXN+5];
struct sa{long long x,y;}s[MAXN+5];
bool cmp(sa A,sa B)
{return A.x<B.x||(A.x==B.x&&A.y>B.y);}
bool check()
{
    for(int i=0;i<=n+1;i++)
        b[i]=c[i]=0;
    int i=1,j=1;
    b[i]=s[i].y;
    while(i<=n&&j<=n)
    {
        if(abs(s[i].x-s[j].x)> R)
        {
            if(i>j)j++;
            else return 0;
        }
        else if(c[j]+b[i]> mid)
        {
            b[i]-=mid-c[j];
            c[j++]=mid;
        }
        else
        {
            c[j]+=b[i];
            b[i++]=0;
            b[i]=s[i].y;
        }
    }
    for(int i=1;i<=n;i++)
        if(b[i])return 0;
    return 1;
}
int main()
{
    scanf("%d",&K);
    while(--K>=0)
    {
        scanf("%d%d",&n,&R);
        for(int i=1;i<=n;i++)
        {
            scanf("%lld%lld",&s[i].x,&s[i].y);
            r=max(r,s[i].y);
        }
        sort(s+1,s+n+1,cmp);
        while(l+1<r)
        {
            mid=(r-l)/2+l;
            bool f=check();
            if(f) r=mid;
            else l=mid;
        }
        printf("%lld\n",r);
        for(int i=0;i<=n+1;i++)
            s[i].x=s[i].y=0;
        n=R=r=l=mid=0;
    }
}