#NOIP2011#聪明的质检员(二分)

8 篇文章 0 订阅

【NOIP2011】聪明的质监员

时间限制: 1 Sec   内存限制: 128 MB

题目描述

       小T 是一名质量监督员,最近负责检验一批矿产的质量。这批矿产共有n 个矿石,从1到n 逐一编号,每个矿石都有自己的重量wi 以及价值vi。检验矿产的流程是:
       1、给定m 个区间[Li,Ri];
       2、选出一个参数W;
       3、对于一个区间[Li,Ri],计算矿石在这个区间上的检验值Yi :


j是矿石编号。这批矿产的检验结果为各个区间的校验值之和。即


若这批矿产的检验结果与所给标准值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。
输出只有一行,包含一个整数,表示所求的最小值。

样例输入

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,每次check,统计值满足单调。


Code:

#include<iostream> 
#include<cstdio> 
#include<cstdlib> 
#include<cstring> 
#include<algorithm> 
using namespace std; 
  
typedef long long LL; 
  
const int Max = 200000; 
const int INF = 0x3f3f3f3f; 
  
int N, M, Maxw, Minw; 
LL S; 
int num[Max + 5]; 
LL sum[Max + 5]; 
int W[Max + 5], V[Max + 5]; 
int L[Max + 5], R[Max + 5]; 
  
bool getint(int & num){ 
    char c; int flg = 1;    num = 0; 
    while((c = getchar()) < '0' || c > '9'){ 
        if(c == '-')   flg = -1; 
        if(c == -1) return 0; 
    } 
    while(c >= '0' && c <= '9'){ 
        num = num * 10 + c - 48; 
        if((c = getchar()) == -1)   return 0; 
    } 
    num *= flg; 
    return 1; 
} 
  
LL Cnt(int w){ 
    num[0] = sum[0] = 0; 
    for(int i = 1; i <= N; ++ i){ 
        num[i] = num[i - 1] + (W[i] >= w); 
        sum[i] = sum[i - 1] + (W[i] >= w ? V[i] : 0); 
    } 
    LL rt = 0; 
    for(int i = 1; i <= M; ++ i) 
        rt += (num[R[i]] - num[L[i] - 1]) * (sum[R[i]] - sum[L[i] - 1]); 
    return rt; 
} 
  
int main(){ 
    getint(N), getint(M); 
    scanf("%lld", & S); 
    Maxw = -INF, Minw = INF; 
    for(int i = 1; i <= N; ++ i){ 
        getint(W[i]), getint(V[i]); 
        Maxw = max(Maxw, W[i]); 
        Minw = min(Minw, W[i]); 
    } 
    for(int i = 1; i <= M; ++ i) 
        getint(L[i]), getint(R[i]); 
    LL a = Cnt(Minw), b = Cnt(Maxw); 
    if(a <= S){ 
        printf("%lld\n", S - a); 
        return 0; 
    } 
    if(b >= S){ 
        printf("%lld\n", b - S); 
        return 0; 
    } 
    int l = Minw, r = Maxw, mid; 
    LL rt = 0x3f3f3f3f3f3f3f3fll; 
    while(l < r){ 
        mid = (l + r) >> 1; 
        LL now = Cnt(mid); 
        rt = min(rt, abs(now - S)); 
        if(now < S) r = mid - 1; 
        else l = mid + 1; 
    } 
    rt = min(rt, min(abs(Cnt(l) - S), abs(Cnt(r) - S))); 
    printf("%lld\n", rt); 
    return 0; 
} 








  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值