首页 > 代码库 > Vijos P1740聪明的质检员
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 151 52 53 54 55 51 52 43 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 }
Vijos P1740聪明的质检员