题目描述
小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 。
关键点:
(x为检查参数,y为各个区间的和) 、递减函数、 时间复杂度(nlog(w))
ac代码:
#include<iostream>
#include<bits/stdc++.h>
using namespace std;
using ll = long long;
const int N = 2 * 1e5 + 10;
int n, m;
int w[N], v[N], l[N], r[N];
ll ans = 1e18;
ll s;
ll sw[N], sv[N];
bool check(int x) {
memset(sw, 0, sizeof(sw));
memset(sv, 0, sizeof(sv));
for (int i = 1; i <= n; i++) {
if (w[i] >= x)sw[i] = sw[i - 1] + 1, sv[i] = sv[i - 1] + v[i];
else sw[i] = sw[i - 1], sv[i] = sv[i - 1];
}
ll y = 0;
for (int i = 1; i <= m; i++) {
y += (sw[r[i]] - sw[l[i] - 1]) * (sv[r[i]] - sv[l[i] - 1]);
}
ans = min(ans, llabs(y - s));
return y>=s;
}
void find() {
int l = 0, r = 1e6 + 1;
while (l + 1 < r) {
int mid = l + r >> 1;
if (check(mid))l=mid;
else r=mid;
}
}
int main() {
cin >> n >> m >> s;
for (int i = 1; i <= n; i++)cin >> w[i] >> v[i];
for (int i = 1; i <= m; i++)cin >> l[i] >> r[i];
find();
cout << ans;
}