题意
题目描述
小(T)是一名质量监督员，最近负责检验一批矿产的质量。这批矿产共有(n)个矿石，从(1)到(n)逐一编号，每个矿石都有自己的重量(w_i)以及价值(v_i)。检验矿产的流程是：
-
给定(m)个区间([L_i,R_i])
-
选出一个参数(W)；
-
对于一个区间([L_i,R_i])，计算矿石在这个区间上的检验值(Y_i)
这批矿产的检验结果(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)。注意：不同区间可能重合或相互重叠。
输出格式：
一个整数，表示所求的最小值。
输入输出样例
输入样例：
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 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;
}