时间限制：1 s 内存限制：128 MB
【问题描述】
小 T 是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有n个矿石，从 1 到n逐一编号，每个矿石都有自己的重量wi以及价值vi。检验矿产的流程是：
1. 给定 m个区间[Li,Ri]；
2. 选出一个参数W；
3. 对于一个区间[Li,Ri]，计算矿石在这个区间上的检验值Yi：Yi=∑j1×∑jvj, j∈[Li,Ri]且 wj≥W,j是矿石编号
这批矿产的检验结果Y为各个区间的检验值之和。即：Y=∑i=1mYi
- 若这批矿产的
检验结果与所给标准值 S 相差太多，就需要再去检验另一批矿产。小 T 不想费时间去检验另一批矿产，所以他想通过调整参数 W 的值，让检验结果尽可能的靠近标准值 S，即使得S−Y的绝对值最小。请你帮忙求出这个最小值。
【输入】
输入文件 qc.in。
第一行包含三个整数n，m，S，分别表示矿石的个数、区间的个数和标准值。
接下来的n 行，每行2 个整数，中间用空格隔开，第i+1 行表示i 号矿石的重量wi 和价值vi 。
接下来的m 行，表示区间，每行2 个整数，中间用空格隔开，第i+n+1 行表示区间[Li,Ri]的两个端点Li 和Ri。注意：不同区间可能重合或相互重叠。
【输出】
输出文件名为qc.out。
输出只有一行，包含一个整数，表示所求的最小值。
【输入输出样例】
qc.in
5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3
qc.out
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≤10^6，0 < S≤10^12，1≤Li≤Ri≤n。
【题解】
对于我来说今天的题只有这一道相当可做，但是当时并没有意识到这一点……考试的时候想得太复杂，刚开始画绝对值函数想三分+线段树，打完了才发现线段树行不通，坚信这是一道线段树套Splay或无旋Treap区间求和的题，发现一来不太会打二来数据范围根本承受不了，然后万念俱灰直接弃疗打了个O(n^2)的暴力查询+三分。
虽然这是个单峰函数，但是并不是必须三分。|y-s|与w的绝对值函数是单峰的，但是绝对值函数是根据y与w反相关的单调函数得来的，所以二分也能解决这个问题，只要寻找与最小值最接近的值即可。至于check，甚至不用什么数据结构，只要前缀和就能解决。把[l,r]转化成[1,r]-[1,l-1]，再分别统计到每个位置大于w的矿石个数、价值和，可以在O(n)的时间内计算出取此w的答案。
一道NOIP题，十分清楚十分简单，却被我生生做成了“强行数据结构，自己还不会打”。前缀和、区间 转化这种方法在ad爷的数数题里面已经出现过，当时没想到，现在还是没想到。要充分利用已有的知识，以之为工具解决问题，这里的工具不仅有算法和数据结构，更有许多转化变通的思想。
#include<iostream> #include<cstdio> #include<cstring> #include<cmath> using namespace std; const int sj=200010; int n,m,jd,w[sj],v[sj],l[sj],r[sj],zj,yj,cnt[sj]; long long s,jg,vh[sj]; long long check(int x) { memset(cnt,0,sizeof(cnt)); memset(vh,0,sizeof(vh)); for(int i=1;i<=n;i++) { cnt[i]=cnt[i-1]+((w[i]>=x)?1:0); vh[i]=vh[i-1]+((w[i]>=x)?v[i]:0); } long long ans=0; for(int i=1;i<=m;i++) ans+=(vh[r[i]]-vh[l[i]-1])*(cnt[r[i]]-cnt[l[i]-1]); return ans; } int main() { scanf("%d%d%lld",&n,&m,&s); jg=s; for(int i=1;i<=n;i++) { scanf("%d%d",&w[i],&v[i]); if(w[i]>jd) jd=w[i]; } for(int i=1;i<=m;i++) scanf("%d%d",&l[i],&r[i]); yj=jd; zj=1; while(zj<yj-1) { int mid=(zj+yj)>>1; long long ans; ans=check(mid); if(ans>=s) zj=mid; else yj=mid; if(abs(ans-s)<jg) jg=abs(ans-s); if(jg==0) break; } long long ans=check(zj); if(abs(ans-s)<jg) jg=abs(ans-s); ans=check(yj); if(abs(ans-s)<jg) jg=abs(ans-s); printf("%lld",jg); return 0; }