[NOIP2011]聪明的质监员

ThearticlediscussesanITexpertsproblemwheretheyusebinarysearchtominimizethedifferencebetweentheinspectionresultandthegivenstandardvalueforabatchofminerals,optimizingtheinspectionprocess.
摘要由CSDN通过智能技术生成

小T是一名质量监督员,最近负责检验一批矿产的质量。这批矿产共有n个矿石,从1到 n 逐一编号,每个矿石都有自己的重量wi以及价值vi。检验矿产的流程是:
1、给定m个区间[Li,Ri];
2、选出一个参数W;

3、对于一个区间[Li,Ri],计算矿石在这个区间上的检验值 Yi :

Yi=∑j1∗∑jvjY_{i}=\sum_{j}{1}*\sum_{j}{v_{j}}Yi​=∑j​1∗∑j​vj​且wj≥W,j是矿石编号

这批矿产的检验结果Y为各个区间的检验值之和。即:Y=∑i=1mYiY=\sum_{i=1}^{m}{Y_{i}}Y=∑i=1m​Yi​
若这批矿产的检验结果与所给标准值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

输入
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 ≤ 106,0 < S ≤ 1012,1 ≤ Li ≤ Ri ≤ n。

思路:这里的区间检验值是指在编号为Li,Ri的范围内,(质量达到参数w的个数)与(质量达到w的矿石的价值之和)的乘积,检验结果是各个检验值之和。显然在这里,w的值并不能直接得出,它与多个量都有关。那么我们尝试去枚举每一个w,直接枚举效率很低,所以我们可以运用二分将w的值确定,也就是通过二分w的值,使得s-y的绝对值最小。那么二分可不可行呢,这里我们要求的是最小的绝对值,因此二分左右边界改变的条件是s与y的大小关系,而w越大,符合要求的越少,相应的检验值y也应该越少,这样我们就可以进行二分了。

接下来我们要做的是对于每个w,我们期望能反馈出一个检验结果y,为此,我定义了一个函数pd。它的返回值就是检验结果,输入量是w。因为有多个区间的关系,如果每一次都枚举的话容易爆,因此做一下预处理,求前缀和。

#include<bits/stdc++.h>
using namespace std;
int n,m;
long long s;//s比较大,要用ll,不然会爆
int num[200010];//num[i]代表前i个中符合条件的矿石数目
long long sum[200010];//代表前i个矿石的价值和
int L[200010],R[200010];//存每个检验区间
struct ks
{
	int v;
	int w;	
}a[200010];//结构体
long long pd(int w)
{
	memset(num,0,sizeof(num));//初始化
	memset(sum,0,sizeof(sum));
	long long k=0;
	for(int i=1;i<=n;i++)
	{
		if(a[i].w>=w)//这里ifelse可以并成两个式子,不过看起来难受,用了比较一般的方法
		{
			num[i]=num[i-1]+1;
			sum[i]=sum[i-1]+a[i].v;
		}
		else
		{
			num[i]=num[i-1];
			sum[i]=sum[i-1];
		}
	}
	for(int i=0;i<m;i++)
	{
		k+=(num[R[i]]-num[L[i]-1])*(sum[R[i]]-sum[L[i]-1]);//右边是每个区间的检验值,-1不要漏
	}
	return k;
}
int main()
{
    cin>>n>>m>>s;
    for(int i=1;i<=n;i++)cin>>a[i].w>>a[i].v;//读入数据
    for(int i=0;i<m;i++)cin>>L[i]>>R[i];
    int l=0,r=1000010,mid;
    long long ss,ans=1000000000010;//ans初值要大
    while(l<=r)
    {
    	mid=(l+r)/2;
    	ss=pd(mid);
    	ans=min(ans,abs(ss-s));
    	if(ss<s)r=mid-1;
    	else    l=mid+1;
	}
	cout<<ans<<endl;
    return 0;   
}

  • 5
    点赞
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值