LG P1314 聪明的质监员

题目描述

小T 是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有 n 个矿石，从 1到n 逐一编号，每个矿石都有自己的重量 wi 以及价值vi 。检验矿产的流程是：

1 、给定m 个区间[Li,Ri]；

2 、选出一个参数 W；

3 、对于一个区间[Li,Ri]，计算矿石在这个区间上的检验值Yi：

这批矿产的检验结果Y 为各个区间的检验值之和。即：Y1+Y2…+Ym

若这批矿产的检验结果与所给标准值S 相差太多，就需要再去检验另一批矿产。小T

不想费时间去检验另一批矿产，所以他想通过调整参数W 的值，让检验结果尽可能的靠近

标准值S，即使得S-Y 的绝对值最小。请你帮忙求出这个最小值。

输入输出格式

输入格式：

输入文件qc.in 。

第一行包含三个整数n，m，S，分别表示矿石的个数、区间的个数和标准值。

接下来的n 行，每行2个整数，中间用空格隔开，第i+1 行表示 i 号矿石的重量 wi 和价值vi。

接下来的m 行，表示区间，每行2 个整数，中间用空格隔开，第i+n+1 行表示区间[Li,Ri]的两个端点Li 和Ri。注意：不同区间可能重合或相互重叠。

输出格式：

输出文件名为qc.out。

输出只有一行，包含一个整数，表示所求的最小值。

输入输出样例

输入样例#1： 复制

5 3 15 
1 5 
2 5 
3 5 
4 5 
5 5 
1 5 
2 4 
3 3 

输出样例#1： 复制

10

说明

【输入输出样例说明】

当W 选4 的时候，三个区间上检验值分别为 20、5 、0 ，这批矿产的检验结果为 25，此

时与标准值S 相差最小为10。

【数据范围】

对于10% 的数据，有 1 ≤n ，m≤10；

对于30% 的数据，有 1 ≤n ，m≤500 ；

对于50% 的数据，有 1 ≤n ，m≤5,000；

对于70% 的数据，有 1 ≤n ，m≤10,000 ；

对于100%的数据，有 1 ≤n ，m≤200,000，0 < wi, vi≤10^6，0 < S≤10^12，1 ≤Li ≤Ri ≤n 。

---

/*WA了两天，不知道为什么
最后发现居然是绝对值函数的参数和读入优化没有开long long，气人*/

---

一开始的思路居然是二分abs(s-y)然后再验证可行性…只有大概25分左右
然后想到了二分w，因为y是关于w的单调函数
w增大时，y减小（因为w增大会使得w[i]>=w[i]的矿石个数减小）
因此可以二分w来求
然后求区间里的w[i]>=w的矿石的个数和价值之和可以用前缀和来求，加个条件就ok

代码

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>
#define ei else if
#define ll long long int
#define For(i,l,r) for(long long int i=l;i<=r;++i)
using namespace std;
ll read()
{
    char c;
    ll a=0;
    bool flag=0;
    while((c=getchar())==' '||c=='\r'||c=='\n');
    if(c=='-')
    {
        flag=1;
        c=getchar();
    }
    while(isdigit(c))
    {
        a*=10;
        a+=c;
        a-='0';
        c=getchar();
    }
    return a*(flag?-1:1);
}
ll pv[200001],pc[200001],s,ans=-1,tmp;
int n,m,w[200001],v[200001],l[200001],r[200001],maxx,minn=2147483647,maxr,minl=2147483647;
ll ab(ll x)
{
    return x>0?x:-x;
}
int main()
{
    //freopen("wa.out","w",stdout);
    n=read();m=read();s=read();
    For(i,1,n)
    {w[i]=read();v[i]=read();if(w[i]<minn)minn=w[i];if(w[i]>maxx)maxx=w[i];}
    For(i,1,m)
    {l[i]=read();r[i]=read();if(r[i]>maxr)maxr=r[i];if(l[i]<minl)minl=l[i];}
    minn--;maxx+=2;
    ll mid;
    while(minn<=maxx)
    {
        memset(pv,0,sizeof pv);
        memset(pc,0,sizeof pc);
        tmp=0;
        mid=((minn+maxx)>>1);
        For(i,minl,maxr)
        {
            pv[i]=pv[i-1];
            pc[i]=pc[i-1];
            if(w[i]>=mid)
            {
                pv[i]+=v[i];
                pc[i]++;
            }
        }
        For(i,1,m)
        {
            tmp+=((pc[r[i]]-pc[l[i]-1])*(pv[r[i]]-pv[l[i]-1]));
        }
        if(ab(s-tmp)<ans||ans<0)
         ans=ab(s-tmp);
        //cout<<minn<<' '<<mid<<' '<<maxx<<' '<<ans<<' '<<tmp<<'\n';
        if(tmp<s)
         maxx=mid-1;
        ei(tmp==s)
         break;
        else
         minn=mid+1;
    }
    printf("%lld",ans);
    return 0;
}