3028. 【NOIP2011DAY2】聪明的质监员

Description

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

 

 

Input

第一行包含三个整数 n，m，S，分别表示矿石的个数、区间的个数和标准值。
接下来的 n 行，每行2 个整数，中间用空格隔开，第i+1 行表示i 号矿石的重量wi 和价
值vi 。
接下来的 m 行，表示区间，每行2 个整数，中间用空格隔开，第i+n+1 行表示区间[Li,
Ri]的两个端点Li 和Ri。注意：不同区间可能重合或相互重叠。

 

 

Output

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

 

 

Sample Input

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

 

Sample Output

10

 

Hint

【输入输出样例说明】
当 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≤106，0 < S≤1012，1≤Li≤Ri≤n。

 

Solution 

二分/三分。

二分：二分的根据是随着W的变小，Y的值会变大。我们可以二分W的值，求出Y值，再判断如果Y要比S小就将Y变大（即将W变大），若Y的值比S要大就将Y变小（即将W变小），每次求出Y值更新ans。时间复杂度O(log2(n)*n)。

三分：此依据是|S-Y|。因为随着Y的值的变大，|S-Y|是先变小后变大的，它的图像类似于一个“U”，满足三分性，我们所要找的min{|S-Y|}就是“U”的最低点。所以我们三分|S-Y|的值，每次用对应的W求出Y，最后找到min即为ans。时间复杂度O(2log3(n)*n)。

优化：二分或三分时都可以不用二分一个值，可以二分数组的下标，即将W数组抽出来排个序，再在最后多加一个位置为无穷大，再直接用这些W[i]作为W的候选。

 

二分Code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define N 200010
using namespace std;
int n,m,l,r,mid,x[N],y[N],w[N],v[N],a[N],b[N];
ll s,ans=115292150460684697;
int check(int p){
	memset(a,0,sizeof(a));
	memset(b,0,sizeof(b));
	for(int i=1;i<=n;i++){
		a[i]=a[i-1];b[i]=b[i-1];
		if(w[i]>=p){a[i]++;b[i]+=v[i];}
	}
	ll sum=0;
	for(int i=1;i<=m;i++) sum+=(a[y[i]]-a[x[i]-1])*(b[y[i]]-b[x[i]-1]);
	ans=min(ans,abs(sum-s));
	if(sum<s) return 1;else return 0;
}
int main(){
	scanf("%d%d%lld",&n,&m,&s);
	for(int i=1;i<=n;i++) scanf("%d%d",&w[i],&v[i]);
	for(int i=1;i<=m;i++) scanf("%d%d",&x[i],&y[i]);
	l=0;r=1000000;
	while(l<=r){
		mid=(l+r)>>1;
		if(check(mid)) r=mid-1;else l=mid+1;
	}
	printf("%lld",ans);
	return 0;
}

 

三分(优化)Code 

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
#define N 200010
using namespace std;
int n,m,x[N],y[N];
ll s,l,r,mid1,mid2,w[N],v[N],q[N],a[N],b[N],ans=115292150460684697;
ll check(ll p){
	memset(a,0,sizeof(a));
	memset(b,0,sizeof(b));
	for(int i=1;i<=n;i++){
		a[i]=a[i-1];b[i]=b[i-1];
		if(w[i]>=p){a[i]++;b[i]+=v[i];}
	}
	ll sum=0;
	for(int i=1;i<=m;i++) sum+=(a[y[i]]-a[x[i]-1])*(b[y[i]]-b[x[i]-1]);
	ans=min(ans,abs(sum-s));
	return abs(sum-s);
}
int main(){
	scanf("%d%d%lld",&n,&m,&s);
	for(int i=1;i<=n;i++) scanf("%d%d",&w[i],&v[i]),q[i]=w[i];
	sort(q+1,q+1+n);
	q[n+1]=1000001;
	for(int i=1;i<=m;i++) scanf("%d%d",&x[i],&y[i]);
	l=0;r=n+1;
	while(l<=r){
		mid1=(l+r)>>1,mid2=(mid1+r)>>1;
		ll x1=check(q[mid1]),x2=check(q[mid2]);
		if(x1>x2) l=mid1+1;else r=mid2-1;
	}
	printf("%lld",ans);
	return 0;
}

---

---

作者：zsjzliziyang 
QQ:1634151125 
转载及修改请注明 
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/86541687