Vijos P1740聪明的质检员

题目

描述

小 T 是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有n个矿石，从1到n逐一编号，每个矿石都有自己的重量wi以及价值vi。检验矿产的流程是：
1、给定m个区间[Li，Ri]；
2、选出一个参数W；
3、对于一个区间[Li，Ri]，计算矿石在这个区间上的检验值Yi：
Yi = ∑1*∑vj，j∈[Li, Ri]且wj ≥ W，j是矿石编号
这批矿产的检验结果Y 为各个区间的检验值之和。即：Y = ∑Yi，i ∈[1, m]
若这批矿产的检验结果与所给标准值S相差太多，就需要再去检验另一批矿产。小T不想费时间去检验另一批矿产，所以他想通过调整参数W的值，让检验结果尽可能的靠近标准值S，即使得S-Y的绝对值最小。请你帮忙求出这个最小值。

格式

输入格式

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

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

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

输出格式

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

样例1

样例输入1[复制]

 

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

样例输出1[复制]

 

10

限制

1s

提示

样例说明：当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。

来源

NOIp2011提高组Day2第二题

题解

这道题目我们可以发现随着w的增大，最后的Y也是增大的，所以我们可以二分出最接近S的几个w试一下。算法复杂度是O(nlogw)的。是可以做的。

代码

 1 /*Author:WNJXYK*/
 2 #include<cstdio>
 3 #include<iostream>
 4 using namespace std;
 5 const int Maxn=200000;
 6 long long w[Maxn+10],v[Maxn+10];
 7 long long num[Maxn+10],vs[Maxn+10];
 8 int l[Maxn+10],r[Maxn+10];
 9 long long S;
10 int n,m;
11 long long maxw;
12 inline long long abs(long long x){
13     if (x<0) return -x;
14     return x;
15 }
16 inline long long getAns(int x){
17     num[0]=vs[0]=0;
18     for (int i=1;i<=n;i++){
19         if (w[i]>=x){
20             num[i]=1;
21             vs[i]=v[i];
22         }else{
23             num[i]=0;
24             vs[i]=0;
25         }
26         num[i]+=num[i-1];
27         vs[i]+=vs[i-1];
28     }
29     long long Ans=0;
30     for (int i=1;i<=m;i++){
31         Ans+=(num[r[i]]-num[l[i]-1])*(vs[r[i]]-vs[l[i]-1]);
32     }
33     return Ans;
34 }
35 int main(){
36     scanf("%d%d%lld",&n,&m,&S);
37     for (int i=1;i<=n;i++) scanf("%lld%lld",&w[i],&v[i]);
38     for (int i=1;i<=n;i++) if (maxw<w[i]) maxw=w[i];
39     for (int i=1;i<=m;i++) scanf("%d%d",&l[i],&r[i]);
40     int left=0,right=maxw;
41     while(left+3<right){
42         int mid=(left+right)/2;
43         if (getAns(mid)>S){
44             left=mid;
45         }else{
46             right=mid;
47         }
48     } 
49     long long Ans=abs(getAns(left)-S);
50     for (int i=left+1;i<=right;i++){
51         long long tmp=abs(getAns(i)-S);
52         if (tmp<Ans)Ans=tmp;
53     }
54     printf("%lld\n",Ans);
55     return 0;
56 }

View Code