P1314 聪明的质监员
难度 提高+/省选-
时空限制 1000ms / 128MB
题目描述
小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 。
——————————————————————————————————————————-
注意
看数据范围,有些变量要开long long
。
思路1
这是一个水思路:
在矿石重量Wi中找到一个最大值和一个最小值,作为二分的起点和终点,然后进行二分,因为这道题终答案具有递增的性质,所以当算出来的值大于标准时,r = mid
,反之l = mid
。然后写一个函数,用来计算每一次的答案,不断更新
bool check(ll x) {
ll sum = 0;
ll he = 0;
ll ans = 0;
for(ll i = 1; i <= m; ++i) {
sum = 0, he = 0;
for(ll j = b[i].L; j <= b[i].R; ++j) {
if(a[j].w >= x) {
sum ++;
he += a[j].v;
}
}
ans += sum * he;
}
if(abs(ans - S) < ANS){
ANS = abs(ans - S);
}
if(ans > S) return true;
else return false;
}
-------------------------------
ll l = v[1], r = v[n];
ll mid;
while(l + 1 < r) {
mid = (l + r) / 2;
if(check(mid)){
l = mid;
}
else r = mid;
}
思路二
二分的基本思路不变,只是在计算的时候每次预处理一个前缀和,这样可以仅剩不少时间
for(int i=1;i<=n;++i){
if(W[i]>=w) v[i]=v[i-1]+V[i],num[i]=num[i-1]+1;
else v[i]=v[i-1],num[i]=num[i-1];
}
for(int i=1;i<=m;++i){
Y+=(v[R[i]]-v[L[i]-1])*(num[R[i]]-num[L[i]-1]);
}
最终代码
#include<bits/stdc++.h>
using namespace std;
const int N = 2e6 + 9;
long long n,m;
long long L[N],R[N];
long long W[N],V[N],num[N];
long long minx=1e7+1,maxx=-1;
long long ans=1e13+1,s;
long long Y=0;
long long v[N];
inline long long check(long long w){
Y=0;
memset(v,0,sizeof(v));
memset(num,0,sizeof(num));
for(int i=1;i<=n;++i){
if(W[i]>=w) v[i]=v[i-1]+V[i],num[i]=num[i-1]+1;
else v[i]=v[i-1],num[i]=num[i-1];
}
for(int i=1;i<=m;++i){
Y+=(v[R[i]]-v[L[i]-1])*(num[R[i]]-num[L[i]-1]);
}
return Y ;
}
int main (){
scanf("%lld%lld%lld",&n,&m,&s);
for(int i=1;i<=n;++i){
scanf("%lld %lld",&W[i],&V[i]);
minx=min(W[i],minx);
maxx=max(W[i],maxx);
}
for(int i=1;i<=m;++i){
scanf("%lld %lld",&L[i],&R[i]);
}
long long l=minx,r=maxx+1;
while(l+1<r){
long long mid=l+r>>1;
ans=min(ans,llabs(check(mid)-s));
if(Y>s) l=mid;
else r=mid;
}
long long ls=min(llabs(check(l)-s),llabs(check(r)-s));
printf("%lld\n",ls>ans?ans:ls);
}