小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=∑j1∗∑jvj且wj≥W,j是矿石编号
这批矿产的检验结果Y为各个区间的检验值之和。即:Y=∑i=1mYiY=\sum_{i=1}^{m}{Y_{i}}Y=∑i=1mYi
若这批矿产的检验结果与所给标准值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;
}