Problem Statement
小T 是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有 $n$ 个矿石，从 $1$ 到 $n$ 逐一编号，每个矿石都有自己的重量 $w_i$ 以及价值 $v_i$ 。检验矿产的流程是：
- 给定$ m$ 个区间 $[l_i,r_i]$；
- 选出一个参数 $W$；
- 对于一个区间 $[l_i,r_i]$，计算矿石在这个区间上的检验值 $y_i$：
$$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$$其中 $j$ 为矿石编号。
这批矿产的检验结果 $y$ 为各个区间的检验值之和。即：$\sum\limits_{i=1}^m y_i$
若这批矿产的检验结果与所给标准值 $s$ 相差太多，就需要再去检验另一批矿产。小T 不想费时间去检验另一批矿产，所以他想通过调整参数 $W$ 的值，让检验结果尽可能的靠近标准值 $s$，即使得 $|s-y|$ 最小。请你帮忙求出这个最小值。
第一行包含三个整数 $n,m,s$，分别表示矿石的个数、区间的个数和标准值。
接下来的 $n$ 行，每行两个整数，中间用空格隔开，第 $i+1$ 行表示 $i$ 号矿石的重量 $w_i$ 和价值 $v_i$。
接下来的 $m$ 行，表示区间，每行两个整数，中间用空格隔开，第 $i+n+1$ 行表示区间 $[l_i,r_i]$ 的两个端点 $l_i$ 和 $r_i$。注意：不同区间可能重合或相互重叠。
Output
一个整数，表示所求的最小值。
1
2
3
4
5
6
7
8
9
|
5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3
|
Sample Output
Explanation of example
当 $W$ 选 $4$ 的时候，三个区间上检验值分别为 $20,5 ,0$ ，这批矿产的检验结果为 $25$，此时与标准值 $S$ 相差最小为 $10$。
Constraints
对于 $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 < w_i,v_i≤10^6$，$0 < s≤10^{12}$，$1 ≤l_i ≤r_i ≤n$ 。
Solving
$$
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
$$上式的意思是在[l,r]区间内所有满足[w$_{i}\ge$W]的个数总和×其各自对应的矿物价值的区间和。
由于本题数据过于刁钻，部分变量不开long long会直接WA，建议遇到此类题目全程使用LL或者直接：
#define int long long
以下是AC代码：
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
|
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
const int M = 0;
const int N = 200005;
int t;
ll n,m,s;
ll l=0,r=1e6,mid;
ll ans = numeric_limits<ll>::max();
ll minecnt[N];
ll mineval[N];
struct mineral
{
ll w,v;
}mine[N];
struct ad
{
ll l,r;//有时也可以先存储一下l,r后期再进行区间操作
}a[N];
int main()
{
ios::sync_with_stdio(false);cin.tie(nullptr);
cin >> n >> m >> s;
for(int i=1;i<=n;i++)
{
cin >> mine[i].w >> mine[i].v;
}
for(int i=1;i<=m;i++)
{
cin >> a[i].l >> a[i].r;
}
//二分答案
while(l+1!=r)
{
mid = (l+r)>>1;
for(int i=1;i<=n;i++)
{
if(mine[i].w>=mid)
{
//两个前缀和数组，分别记录满足条件的数量和价值。
minecnt[i] = minecnt[i-1] + 1;
mineval[i] = mineval[i-1] + mine[i].v;
}
else
{
//不满足则不做更新
minecnt[i] = minecnt[i-1];
mineval[i] = mineval[i-1];
}
}
ll sum = 0;
for(int i=1;i<=m;i++)
{
//计算总和
sum += (minecnt[a[i].r] - minecnt[a[i].l-1])*(mineval[a[i].r]-mineval[a[i].l-1]);
}
ans = min(abs(s-sum),ans);//不断更新最小值
//已知W（mid）越大 sum越小。
if(sum <= s)
{
r = mid;//小于收缩右边界
}
else
l = mid;//大于收缩左边界
}
cout << ans;
return 0;
}
|