noip2011 聪明的质监员（二分+前缀和处理+读入优化）

最新推荐文章于 2021-03-08 08:05:32 发布

yuyanggo

最新推荐文章于 2021-03-08 08:05:32 发布

阅读量1.0k

点赞数

文章标签： noip2011 聪明的质监员二分预处理

本文链接：https://blog.csdn.net/yuyanggo/article/details/48787977

版权

P1740聪明的质检员

Accepted

标签：其他二分查找 NOIP提高组2011

描述

小 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的绝对值最小。请你帮忙求出这个最小值。

格式

输入格式

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

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

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

输出格式

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

样例1

样例输入1[复制]

样例输出1[复制]

限制

提示

样例说明：当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。

来源

NOIp2011提高组Day2第二题

解析：二分标准w，然后计算标准为w时的检验值（这里要用到预处理）。

s的读入弄成%d，wa了好多次。。。。。。

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cctype>
#include<climits>
using namespace std;

typedef long long LL;
const int maxn=2e5;
int n,m;
int w[maxn+20],v[maxn+20];
LL s,sw[maxn+20],sv[maxn+20];
int L[maxn+20],r[maxn+20];

int getin()
{
  int ans=0;char tmp;
  while(!isdigit(tmp=getchar()));
  do ans=(ans<<3)+(ans<<1)+tmp-'0';
  while(isdigit(tmp=getchar()));
  return ans;
}

long long get(int k)
{
  LL i,ans=0;
  for(i=1;i<=n;i++)
    {
      sw[i]=sw[i-1],sv[i]=sv[i-1];
      if(w[i]>=k)sv[i]+=v[i],sw[i]++;
	}
  for(i=1;i<=m;i++)
    ans+=(sw[r[i]]-sw[L[i]-1])*(sv[r[i]]-sv[L[i]-1]);
  return ans;
}

int main()
{
  int i;
  scanf("%d%d%I64d",&n,&m,&s);
  for(i=1;i<=n;i++)w[i]=getin(),v[i]=getin();
  for(i=1;i<=m;i++)L[i]=getin(),r[i]=getin();
  int mid,x=w[1],y=w[1];
  for(i=2;i<=n;i++)x=min(x,w[i]),y=max(y,w[i]);
  LL k,ans=min(abs(get(x)-s),(LL)s);
  while(x<=y)
    {
      if(ans==0)break;
      mid=(x+y)/2,k=get(mid);
      ans=min(ans,abs(k-s));
      if(ans<0)printf("%d \n",mid);
      if(k<s)y=mid-1;
      else x=mid+1;
	}
  printf("%I64d\n",ans);
  return 0;
}