洛谷1314 聪明的质监员

NOIP2011 Day2 T1

原题地址：

https://www.luogu.org/problem/show?pid=1314

题目描述

小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：

53 15

15

25

35

45

55

15

24

33

输出样例#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 。

 

因为暑假培训的时候就讲过了，所以正解算法被剧透得一干二净=_=，现在打起来明显顺手多了。不过还是有细节没有把握好，调了好久才A过。。。

解题思路

二分+前缀和，看懂公式一切好说，看不懂……也没办法……

就是求该区间符合要求的矿石数num*sum(满足要求的矿石的价值)，再对所有区间内的这个值求和，跟s比较。

参考代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstring>
#include<iomanip>
#define ll longlong
using namespacestd;
llw[200005],v[200005],x[200005],y[200005];
llnum[200005],val[200005]; //数组开大不要钱！
long long n,m,s;
 
ll ef(ll l,ll r)
{
    memset(num,0,sizeof(num));
    memset(val,0,sizeof(val));
    ll mid=(l+r)/2,ans=0,y1,y2;
    for (int i=1;i<=n;i++)
    {
       num[i]=num[i-1];
       val[i]=val[i-1];
       if (w[i]>=mid)
       {
           num[i]++;
           val[i]+=v[i];
       }
    }
    for(int i=1;i<=m;i++)
    {
       y1=num[y[i]]-num[x[i]-1];
       y2=val[y[i]]-val[x[i]-1];
       ans+=y1*y2;
    }
    return ans;
}
 
int main()
{
    ios::sync_with_stdio(false);
    long long mx=0,mi=1000005;
    cin>>n>>m>>s;
    for (int i=1;i<=n;i++)
    {
       cin>>w[i]>>v[i];
       if (w[i]<mi) mi=w[i];
       if (w[i]>mx) mx=w[i];
    }
    for (int i=1;i<=m;i++)
       cin>>x[i]>>y[i];
    ll ans=2100000000000,l=mi,r=mx;  //极大值往大里赋，我会说就因为这个WA了7个点吗！！！
    while (l<=r)
    {
       ll tmp=ef(l,r);
       ll mid=(l+r)/2;
       if (abs(s-tmp)<ans) ans=abs(s-tmp);  //这一步一定要有。
       if (tmp==s) {cout<<0; return 0;}
       if (tmp<s) {r=mid-1;}
       if (tmp>s) {l=mid+1;}
    }
    cout<<ans;
    return 0;
}