[NOIP2011 提高组] 聪明的质监员
题目描述
小T
是一名质量监督员,最近负责检验一批矿产的质量。这批矿产共有
n
n
n 个矿石,从
1
1
1 到
n
n
n 逐一编号,每个矿石都有自己的重量
w
i
w_i
wi 以及价值
v
i
v_i
vi 。检验矿产的流程是:
- 给定$ m$ 个区间 [ l i , r i ] [l_i,r_i] [li,ri];
- 选出一个参数 W W W;
- 对于一个区间 [ l i , r i ] [l_i,r_i] [li,ri],计算矿石在这个区间上的检验值 y i y_i yi:
y i = ∑ j = l i r i [ w j ≥ W ] × ∑ j = l i r i [ w j ≥ W ] v j y_i=\sum\limits_{j=l_i}^{r_i}[w_j \ge W] \times \sum\limits_{j=l_i}^{r_i}[w_j \ge W]v_j yi=j=li∑ri[wj≥W]×j=li∑ri[wj≥W]vj
其中 j j j 为矿石编号。
这批矿产的检验结果 y y y 为各个区间的检验值之和。即: ∑ i = 1 m y i \sum\limits_{i=1}^m y_i i=1∑myi
若这批矿产的检验结果与所给标准值
s
s
s 相差太多,就需要再去检验另一批矿产。小T
不想费时间去检验另一批矿产,所以他想通过调整参数
W
W
W 的值,让检验结果尽可能的靠近标准值
s
s
s,即使得
∣
s
−
y
∣
|s-y|
∣s−y∣ 最小。请你帮忙求出这个最小值。
输入格式
第一行包含三个整数 n , m , s n,m,s n,m,s,分别表示矿石的个数、区间的个数和标准值。
接下来的 n n n 行,每行两个整数,中间用空格隔开,第 i + 1 i+1 i+1 行表示 i i i 号矿石的重量 w i w_i wi 和价值 v i v_i vi。
接下来的 m m m 行,表示区间,每行两个整数,中间用空格隔开,第 i + n + 1 i+n+1 i+n+1 行表示区间 [ l i , r i ] [l_i,r_i] [li,ri] 的两个端点 l i l_i li 和 r i r_i ri。注意:不同区间可能重合或相互重叠。
输出格式
一个整数,表示所求的最小值。
样例 #1
样例输入 #1
5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3
样例输出 #1
10
提示
【输入输出样例说明】
当 W W W 选 4 4 4 的时候,三个区间上检验值分别为 20 , 5 , 0 20,5 ,0 20,5,0 ,这批矿产的检验结果为 25 25 25,此时与标准值 S S S 相差最小为 10 10 10。
【数据范围】
对于 $10% $ 的数据,有 1 ≤ n , m ≤ 10 1 ≤n ,m≤10 1≤n,m≤10;
对于 $30% $的数据,有 1 ≤ n , m ≤ 500 1 ≤n ,m≤500 1≤n,m≤500 ;
对于 $50% $ 的数据,有 $ 1 ≤n ,m≤5,000$;
对于 70 % 70\% 70% 的数据,有 1 ≤ n , m ≤ 10 , 000 1 ≤n ,m≤10,000 1≤n,m≤10,000 ;
对于 100 % 100\% 100% 的数据,有 $ 1 ≤n ,m≤200,000$, 0 < w i , v i ≤ 1 0 6 0 < w_i,v_i≤10^6 0<wi,vi≤106, 0 < s ≤ 1 0 12 0 < s≤10^{12} 0<s≤1012, 1 ≤ l i ≤ r i ≤ n 1 ≤l_i ≤r_i ≤n 1≤li≤ri≤n 。
详解
二分题,找单调性:w↑,y↓
check基准值W
需要求(|y-s|)min
如何处理[l,r]
用前缀和处理大于基准值的个数,然后同时处理另一项的前缀和
最后求和
记得记录最小值
代码实现
#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=2e5+5;
int n,m;
int s;
int l=0,r=maxn*10;
int ans=1e14;
int sx[maxn],sv[maxn];
struct rock{ int w,v; }a[maxn];
struct q{ int l,r; }b[maxn];
bool check(int d){
int as=0;
for(int i=1;i<=n;i++)
if(a[i].w>=d) sx[i]=sx[i-1]+1,sv[i]=sv[i-1]+a[i].v;
else sx[i]=sx[i-1],sv[i]=sv[i-1];
for(int i=1;i<=m;i++)
as+=(sx[b[i].r]-sx[b[i].l-1])*(sv[b[i].r]-sv[b[i].l-1]);
if(abs(s-as)<=ans) ans=abs(s-as);
if(s>as) return false;
else return true;
}
signed main()
{
cin>>n>>m>>s;
for(int i=1;i<=n;i++) cin>>a[i].w>>a[i].v;
for(int i=1;i<=m;i++) cin>>b[i].l>>b[i].r;
while(l<=r){
int mid=l+r>>1;
if(check(mid)) l=mid+1;
else r=mid-1;
}
cout<<ans;
return 0;
}