NOIP2011提高组 聪明的质检员（二分答案）

【题目描述】

小T 是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有 n 个矿石，从 1到n 逐一编号，每个矿石都有自己的重量 wi 以及价值vi 。检验矿产的流程是：

1 、给定m 个区间[Li,Ri]；

2 、选出一个参数 W；

3 、对于一个区间[Li,Ri]，计算矿石在这个区间上的检验值Yi：

这批矿产的检验结果Y 为各个区间的检验值之和。即：Y1+Y2…+Ym

若这批矿产的检验结果与所给标准值S 相差太多，就需要再去检验另一批矿产。小T

不想费时间去检验另一批矿产，所以他想通过调整参数W 的值，让检验结果尽可能的靠近

标准值S，即使得S-Y 的绝对值最小。请你帮忙求出这个最小值。

【输入格式】
第一行包含三个整数n，m，S，分别表示矿石的个数、区间的个数和标准值。

接下来的n 行，每行2个整数，中间用空格隔开，第i+1 行表示 i 号矿石的重量 wi 和价值vi。

接下来的m 行，表示区间，每行2 个整数，中间用空格隔开，第i+n+1 行表示区间[Li,Ri]的两个端点Li 和Ri。注意：不同区间可能重合或相互重叠。

【输出格式】
输出只有一行，包含一个整数，表示所求的最小值。

【输入样例】
5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3
【输出样例】
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 。

这道题明显要二分猜答案，并且通过二分W的值来求Y的值。所以在二分W的时候要依据S与Y的关系调整W的值。我之前想的是abs(S-Y)变小时再改W，但是取绝对值就意味着Y有可能比S大，也可能比S小。因此可以二分找一个Y1要大于S的最小值，再二分找一个Y2要小于S的最大值，最后取abs(S-Y1),abs(S-Y2)的较小值。算Y时要用前缀和节省时间。整个时间复杂度为（(n+m)*40）。

#include<stdio.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 2000002;
const long long inf = 40000000000000002ll;
int n, m;
long long s, t = inf, t1[maxn], t2[maxn];
//二分猜答案 

struct data1
{
    int v, w;
};
data1 g[maxn];

struct data2
{
    int l, r;
};
data2 d[maxn];

long long calc(long long W) //求y 
{
    long long t = 0, tt = 0;
    for (int i = 1; i <= n; i++) //计算前缀和
    {
        if (g[i].w >= W)
        {
            t++, tt += g[i].v;
        }
        t1[i] = t;
        t2[i] = tt;
    }

    long long Y = 0;
    for (int i = 1; i <= m; i++)
    {
        Y += (t1[d[i].r] - t1[d[i].l - 1])*(t2[d[i].r] - t2[d[i].l - 1]);
    }
    return Y;
}

int main()
{
    freopen("1.txt", "r", stdin);
    //freopen("qc.out", "w", stdout);

    scanf("%d %d %I64d", &n, &m, &s);
    int maxw = -1;
    for (int i = 1; i <= n; i++)
    {
        scanf("%d %d", &g[i].w, &g[i].v);
        maxw = max(maxw, g[i].w);
    }
    for (int i = 1; i <= m; i++)
    {
        scanf("%d %d", &d[i].l, &d[i].r);
    }

    int a = 0, b = maxw,m1=0,m2=0;
    long long y1 = inf, y2 = -inf, ans1 = inf, ans2 = -inf;
    while (a <= b)
    {
        int mid = a + b >> 1; //s-y<0的最大值 
        y1 = calc(mid);
        //ans1 = min(ans1, abs(s - y1));
        if (y1<s)
        {
            b = mid - 1; //y取大
            m1 = mid;
        }
        else a = mid + 1;
    }
    ans1 =abs(s - calc(m1));
    a = 0, b = maxw; //s-y>0的最小值 
    while (a <= b)
    {
        int mid = a + b >> 1;
        y2 = calc(mid);
        //ans2 = max(ans2, abs(s - y2));
        if (s<y2)
        {
            a = mid + 1;//y取小
            m2 = mid;
        }
        else b = mid - 1;
    }
    ans2 = abs(s - calc(m2));
    printf("%I64d",min(ans1,ans2));
    return 0;
}