小T 是一名质量监督员,最近负责检验一批矿产的质量。这批矿产共有 n 个矿石,从 1到n 逐一编号,每个矿石都有自己的重量 wi 以及价值vi 。检验矿产的流程是:
1 、给定m 个区间[Li,Ri];
2 、选出一个参数 W;
3 、对于一个区间[Li,Ri],计算矿石在这个区间上的检验值Yi:
这批矿产的检验结果Y 为各个区间的检验值之和。即:Y1+Y2...+Ym
若这批矿产的检验结果与所给标准值S 相差太多
不想费时间去检验另一批矿产,所以他想通过调整参数W 的值,让检验结果尽可能的靠近
标准值S,即使得S-Y 的绝对值最小。请你帮忙求出这个最小值。
输入
第一行包含三个整数n,m,S,分别表示矿石的个数、区间的个数和标准值。
接下来的n 行,每行2个整数,中间用空格隔开,第i+1 行表示 i 号矿石的重量 wi 和价值vi。
接下来的m 行,表示区间,每行2 个整数,中间用空格隔开,第i+n+1 行表示区间[Li,Ri]的两个端点Li 和Ri。注意:不同区间可能重合或相互重叠。
#include <iostream>
#include <algorithm>
using namespace std;
struct interval{
int l,r;
interval(){}
interval(int ll,int rr):l(ll),r(rr){}
};
interval range[200020];
struct mineNode{
int w,v,index;
mineNode(){}
mineNode(int ww,int vv,int ii):w(ww),v(vv),index(ii){} //结构体
};
bool operator < (mineNode a,mineNode b){
return a.w < b.w;
}
mineNode mine[200020];
int n,m;
long long s;
long long prefixSum[200020];// i 前面满足条件的物品的个数
long long prefixSumValue[200020]; // i 前面满足条件的物品的价值总和
long long calculatorY(int w){
//预处理出前缀和
int idx = 0;
long long tmpSum = 0;
for (int j = 1; j <= n; ++j) {
if(mine[j].w >= w){
idx++;
tmpSum += mine[j].v;
}
prefixSum[j] = idx;
prefixSumValue[j] = tmpSum;
}
long long sigma = 0,sigmav = 0,sumY = 0;
for (int i = 0; i < m; ++i) {//遍历所有区间
sigma = prefixSum[range[i].r] - prefixSum[range[i].l - 1];//注意此处的-1
sigmav = prefixSumValue[range[i].r] - prefixSumValue[range[i].l - 1];
sumY += sigma * sigmav;
}
return sumY;
}
int main() {
cin >> n >> m >> s;
int minW = 1 << 30,maxW = -1;
int w,v;
for (int i = 1; i <= n; ++i) {
cin >> w >> v;
mine[i] = mineNode(w,v,i);
minW = min(minW,w); //预处理
maxW = max(maxW,w);
}
int t1,t2;
for (int j = 0; j < m; ++j) {
cin >> t1 >> t2;
range[j] = interval(t1,t2);
}
int l = minW,r = maxW + 1;
int mid;
long long y,ans = 100000000000000;
while (l <= r){
mid = (l + r) / 2;
y = calculatorY(mid);
//注意此处的处理
ans = min(ans,abs(y - s));
if(y == s) break;
if(y > s){
l = mid + 1; //二分走起
} else {
r = mid - 1;
}
}
cout << ans << endl; //输出
return 0;
}