【NOIP2011】聪明的质监员
时间限制: 1 Sec 内存限制: 128 MB题目描述
小T 是一名质量监督员,最近负责检验一批矿产的质量。这批矿产共有n 个矿石,从1到n 逐一编号,每个矿石都有自己的重量wi 以及价值vi。检验矿产的流程是:
1、给定m 个区间[Li,Ri];
2、选出一个参数W;
3、对于一个区间[Li,Ri],计算矿石在这个区间上的检验值Yi :
j是矿石编号。这批矿产的检验结果为各个区间的校验值之和。即
若这批矿产的检验结果与所给标准值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。
输出只有一行,包含一个整数,表示所求的最小值。
样例输入
5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3
样例输出
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。
二分答案的W,每次check,统计值满足单调。
Code:
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long LL;
const int Max = 200000;
const int INF = 0x3f3f3f3f;
int N, M, Maxw, Minw;
LL S;
int num[Max + 5];
LL sum[Max + 5];
int W[Max + 5], V[Max + 5];
int L[Max + 5], R[Max + 5];
bool getint(int & num){
char c; int flg = 1; num = 0;
while((c = getchar()) < '0' || c > '9'){
if(c == '-') flg = -1;
if(c == -1) return 0;
}
while(c >= '0' && c <= '9'){
num = num * 10 + c - 48;
if((c = getchar()) == -1) return 0;
}
num *= flg;
return 1;
}
LL Cnt(int w){
num[0] = sum[0] = 0;
for(int i = 1; i <= N; ++ i){
num[i] = num[i - 1] + (W[i] >= w);
sum[i] = sum[i - 1] + (W[i] >= w ? V[i] : 0);
}
LL rt = 0;
for(int i = 1; i <= M; ++ i)
rt += (num[R[i]] - num[L[i] - 1]) * (sum[R[i]] - sum[L[i] - 1]);
return rt;
}
int main(){
getint(N), getint(M);
scanf("%lld", & S);
Maxw = -INF, Minw = INF;
for(int i = 1; i <= N; ++ i){
getint(W[i]), getint(V[i]);
Maxw = max(Maxw, W[i]);
Minw = min(Minw, W[i]);
}
for(int i = 1; i <= M; ++ i)
getint(L[i]), getint(R[i]);
LL a = Cnt(Minw), b = Cnt(Maxw);
if(a <= S){
printf("%lld\n", S - a);
return 0;
}
if(b >= S){
printf("%lld\n", b - S);
return 0;
}
int l = Minw, r = Maxw, mid;
LL rt = 0x3f3f3f3f3f3f3f3fll;
while(l < r){
mid = (l + r) >> 1;
LL now = Cnt(mid);
rt = min(rt, abs(now - S));
if(now < S) r = mid - 1;
else l = mid + 1;
}
rt = min(rt, min(abs(Cnt(l) - S), abs(Cnt(r) - S)));
printf("%lld\n", rt);
return 0;
}