631. [NOIP2011] 聪明的质监员

★★   输入文件：qc.in   输出文件：qc.out   简单对比
时间限制：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。

　　这道题一开始看到Σ还以为是数学题，果断跳过，然后发现自己错了……

　　然后又以为是平衡树，然而昨天虽然刚打了平衡树却因为模板错误而没能学习区间操作，果断心虚。于是打了最暴力的暴力，强行水了5分。

　　这道题最重要的有两点，第一，前缀和，第二，二分。

　　先解释前缀和，这个很容易理解，就是sumv[i]为在i之前所有high不小于w的value的前缀和sumw[i]为在i之前高度不小于w的个数。

　　然后就是二分，考试的时候以为是三分，之前没打过，打挂了。如果以abs(Y-S)为y轴，w为x轴，那么图像是一个绝对值函数，类似二次函数，很容易让人想到三分，然而值得注意的是Y本身是随w增大而递减的，而Y-S也是如此，我们就可以利用这个性质不断地二分w，使Y无限的去接近S,Y大于0就右移，反之左移，记得开long long。

　　

 #include<iostream>

 #include<cstdlib>

 #include<cstdio>

 #include<cstring>

 #include<queue>

 #include<algorithm>

 #include<cmath>

 #include<map>

 using namespace std;

 int n,m;

 long long s;

 struct no{

     int va;

     int w;

 }node[];

 struct qu{

     int li,ri;

 }que[];

 long long sumw[],sumv[];

 long long check(int w){

     memset(sumv,,sizeof(sumv));

     memset(sumw,,sizeof(sumw));

     for(int i=;i<=n;i++)

     {

         sumw[i]=sumw[i-];

         sumv[i]=sumv[i-];

         if(node[i].w>=w)

         {

             sumw[i]++;

             sumv[i]+=node[i].va;

         }

     }

     long long ans=;

     for(int i=;i<=m;i++)

     {

         ans+=(sumv[que[i].ri]-sumv[que[i].li-])*(sumw[que[i].ri]-sumw[que[i].li-]);

     }

     return ans;

 }

 long long minn(long long a,long long b){

     if(a>b)return b;

     return a;

 }

 long long llab(long long a){

     if(a<)

     return -a;

     return a;

 }

 int main(){

     scanf("%d%d%lld",&n,&m,&s);

     int mx=;

     for(int i=;i<=n;i++)

     {

         scanf("%d%d",&node[i].w,&node[i].va);

         mx=max(mx,node[i].w);

     }

     for(int i=;i<=m;i++)

         scanf("%d%d",&que[i].li,&que[i].ri);

     int li=,ri=mx;

     long long ans=-;

     while(li<=ri)

     {

         int mid=(li+ri)/;

         long long x=check(mid);

         if(ans==-)ans=llab(x-s);

         ans=minn(ans,llab(x-s));

         if(x>s)

             li=mid+;

         else

             ri=mid-;

     }

     printf("%lld\n",ans);

     //while(1);

     return ;

 }

NOIP2011聪明的质监员题解的相关教程结束。