题意

题目描述

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

给定(m)个区间([L_i,R_i])
选出一个参数(W)；
对于一个区间([L_i,R_i])，计算矿石在这个区间上的检验值(Y_i)

P1314

这批矿产的检验结果(Y)为各个区间的检验值之和。即：(Y_1+Y_2...+Y_m)

若这批矿产的检验结果与所给标准值(S)相差太多，就需要再去检验另一批矿产。小(T)不想费时间去检验另一批矿产，所以他想通过调整参数(W)的值，让检验结果尽可能的靠近标准值(S)，即使得(S-Y)的绝对值最小。请你帮忙求出这个最小值。

输入输出格式

输入格式：

第一行包含三个整数(n,m,S)，分别表示矿石的个数、区间的个数和标准值。

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

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

输出格式：

一个整数，表示所求的最小值。

输入输出样例

输入样例：

输出样例：

说明

【输入输出样例说明】

当(W)选(4)的时候，三个区间上检验值分别为(20,5,0)，这批矿产的检验结果为 (25)，此时与标准值(S)相差最小为(10)。

【数据范围】

对于(10 \%)的数据，有(1 leq n,m leq 10)；

对于(30 \%)的数据，有(1 leq n,m leq 500)；

对于(50 \%)的数据，有(1 leq n,m leq 5,000)；

对于(70 \%)的数据，有(1 ≤n,m leq 10,000)；

对于(100\%)的数据，有(1 leq n,m leq 200,000, 0<w_i,v_i≤10^6, 0<S leq 10^{12},1 leq L_i leq R_i leq n)。

思路

我点错(ans)文件了。 --alecli

今天下午应alecli大佬之邀来刷(NOIP)的各种蓝题，然后就做到了这一道。

开头的那张公式图片看得我一脸懵逼，不过想一想也还是能像清楚的。一共有(n)个物品，每个物品有两个属性(w)和(v)。我们选取一个(W)，然后对于每一段区间，其中所有(w)属性大于(W)的物品的数量乘上这些物品的(v)值之和，就是这一段区间的(Y)值。把每个区间的(Y)值加起来，就是总(Y)值。题目要求的，就是总(Y)值与给定的数(S)的最小差值。

有两点显然的特性：

选的(W)越大，每个区间所满足(w)属性大于(W)的物品的数量越大；
选的(W)越大，这些物品的(v)值之和也越大。

因为(Y)与这两个属性正相关，所以(Y)关于(W)单调递增，那么我们就可以二分(W)的值，来计算(Y)了。

第二个问题就是如何快速地计算(Y)。因为对于没一段区间我们都只询问区间和，所以很容易想到预处理(O(n))，查询(O(1))的前缀和。具体来说我们这样操作：

tmp=0;//统计Y值
memset(s,0,sizeof s);//s是v的前缀和
memset(cnt,0,sizeof cnt);//cnt是数量的前缀和
for(LL i=1;i<=n;i++)//预处理
{
    if(w[i]>=W) s[i]=s[i-1]+v[i],cnt[i]=cnt[i-1]+1;//满足w>=W，对答案有贡献
    else s[i]=s[i-1],cnt[i]=cnt[i-1];//对答案无贡献
}
for(LL i=1;i<=m;i++) tmp+=(s[r[i]]-s[l[i]-1])*(cnt[r[i]]-cnt[l[i]-1]);//查询每一段区间的Y值

顺利(AC)。

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const LL MAXN=2e5+5;
LL n,m,S,sum,tmp,ans=LLONG_MAX,L=LLONG_MAX,R=LLONG_MIN;
LL w[MAXN],v[MAXN],l[MAXN],r[MAXN],s[MAXN],cnt[MAXN];
LL read()
{
    LL re=0;
    char ch=getchar();
    while(!isdigit(ch)) ch=getchar();
    while(isdigit(ch)) re=(re<<3)+(re<<1)+ch-'0',ch=getchar();
    return re;
}
LL check(LL lzq)
{
    tmp=0;
    memset(s,0,sizeof s);
    memset(cnt,0,sizeof cnt);
    for(LL i=1;i<=n;i++)
    {
        if(w[i]>=lzq) s[i]=s[i-1]+v[i],cnt[i]=cnt[i-1]+1;
        else s[i]=s[i-1],cnt[i]=cnt[i-1];
    }
    for(LL i=1;i<=m;i++) tmp+=(s[r[i]]-s[l[i]-1])*(cnt[r[i]]-cnt[l[i]-1]);
    sum=tmp-S;
    if(sum<0) sum=-sum;
    return tmp>S;
}
int main()
{
    n=read(),m=read(),S=read();
    for(LL i=1;i<=n;i++) w[i]=read(),v[i]=read(),L=min(L,w[i]),R=max(R,w[i]);
    for(LL i=1;i<=m;i++) l[i]=read(),r[i]=read();
    while(L<=R)
    {
        LL mid=(L+R)>>1;
        if(!check(mid)) R=mid-1;
        else L=mid+1;
        if(ans>sum) ans=sum;
    }
    printf("%lld",ans);
    return 0;
}