松江1843路

松江1843路

（来源：Luogu.org）

题目描述

涞坊路是一条长L米的道路，道路上的坐标范围从0到L，路上有N座房子，第i座房子建在坐标为x[i]的地方，其中住了r[i]人。
松江1843路公交车要在这条路上建一个公交站，市政府希望让最多的人得到方便，因此希望所有的每一个的居民，从家到车站的距离的总和最短。
公交站应该建在哪里呢？

输入输出格式

输入格式：
第一行输入L、N。

接下来N行，每行两个整数x[i]和r[i]。

输出格式：
一个整数，最小的每个人从家到车站的距离的总和。

输入输出样例

输入样例#1：
100 3
20 3
50 2
70 1
输出样例#1：
110

输入样例#2：
100 2
0 1
100 10
输出样例#2：
100
输入样例#3：
10000000000 5
3282894320 391
4394338332 929
6932893249 181
7823822843 440
9322388365 623
输出样例#3：
5473201404068

说明

样例解释1

当建在坐标40的时候，所有人距离车站的距离总和为 |20−40|×3+|50−40|×2+|70−40|×1=110。

数据范围和约定

对于10%的数据，1≤N≤50，R[i]=1。

对于30%的数据，1≤N≤100，R[i]≤10，1≤L≤1000。

对于70%的数据，1≤N≤1000，R[i]≤100，1≤L≤10^6。

对于全部数据，1≤L≤10^10，1≤N≤10^5，0≤x[i]≤L，1≤r[i]≤1000

思路

（来自蒟蒻的心路历程）

二分思想

看到这个题，我率先想到了二分法去解决这个问题，然后满心欢喜的打出来了二分代码。做完之后发现只得了10分。其余全部Wrong Answer。

#include<iostream>
#include<cmath>
using namespace std;
long long i,j,m,n,l,sum,mid,suml,sumr,temp;

struct data
{
    long long x,v;
}a[100001];

long long r()
{
    char ch=getchar();
    long long ans=0;
    while(ch<'0'||ch>'9')
    {
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        ans*=10;
        ans+=ch-'0';
        ch=getchar();
    }
    return ans;
}

long long qiuhe(long long xx)
{
    temp=0;
    for(int i=1;i<=n;i++)
    {
        temp+=abs(a[i].x-xx)*a[i].v;
    }
    return temp;
}

void erfen(long long le,long long ri)
{
    if(ri-le<=1) {sum=min(suml,sumr);return;}
    mid=(le+ri)/2;
    suml=qiuhe(le),sumr=qiuhe(ri);
    if(suml<=sumr)
    {erfen(le,mid);return;}
    else
    {erfen(mid,ri);return;}
}

int main()
{
    l=r(),n=r();
    for(i=1;i<=n;i++)
    a[i].x=r(),a[i].v=r();
    erfen(0,l);
    cout<<sum;
}

分段二分思想

这时，旁边的dalao告诉我，题述花费f(x)的图像应该是不单调的，但对于每一座房子之间的花费函数图像应该是单调的。所以我们只需要分段将它二分求最值，然后再比较即可。结果发现有超时和WA。

#include<iostream>
#include<cmath>
using namespace std;
long long i,j,m,n,l,sum=0x7fffffffffff,mid,suml,sumr,temp;

struct data
{
    long long x,v;
}a[100001];

long long r()
{
    char ch=getchar();
    long long ans=0;
    while(ch<'0'||ch>'9')
    {
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        ans*=10;
        ans+=ch-'0';
        ch=getchar();
    }
    return ans;
}

long long qiuhe(long long xx)
{
    temp=0;
    for(int i=1;i<=n;i++)
    {
        temp+=abs(a[i].x-xx)*a[i].v;
    }
    return temp;
}

void erfen(long long le,long long ri)
{
    if(ri-le<=1) {sum=min(suml,sum);sum=min(sumr,sum);return;}
    mid=(le+ri)/2;
    suml=qiuhe(le),sumr=qiuhe(ri);
    if(suml<=sumr)
    {erfen(le,mid);return;}
    else
    {erfen(mid,ri);return;}
}

int main()
{
    l=r(),n=r();
    for(i=1;i<=n;i++)
    a[i].x=r(),a[i].v=r();
    for(i=1;i<n;i++)
    erfen(a[i].x,a[i+1].x);
    cout<<sum;
}

简单的思想

既然是求最大值，而且函数在每一个小区间内单调，所以我们就把每一个房子对应的函数值求出来比较即可。但是这样有可能超时，因为N比较大。因而要在每一次操作的时候将0点对应的函数值初始化出来，然后向后递推操作。然后就很开心的AC了。

#include<iostream>
#include<cmath>
using namespace std;
long long i,j,m,n,l,sum,mid,le,ri,temp;

struct data
{
    long long x,v;
}a[100001];

long long qiuhe(long long xx)
{
    temp=0;
    for(int i=1;i<=n;i++)
    {
        temp+=abs(a[i].x-xx)*a[i].v;
    }
    return temp;
}

long long r()
{
    char ch=getchar();
    long long ans=0;
    while(ch<'0'||ch>'9')
    {
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        ans*=10;
        ans+=ch-'0';
        ch=getchar();
    }
    return ans;
}

int main()
{
    l=r(),n=r();
    for(i=1;i<=n;i++)
    {
    a[i].x=r(),a[i].v=r();
    ri+=a[i].v;
    sum+=a[i].x*a[i].v;
    }
    long long ans=sum;
    for(i=1;i<=n;i++)
    {
        ans+=le*(a[i].x-a[i-1].x);
        ans-=ri*(a[i].x-a[i-1].x);
        ri-=a[i].v;
        le+=a[i].v;
        sum=min(sum,ans);
    }
    cout<<sum;
}