码迷,mamicode.com
首页 > 其他好文 > 详细

聪明的质监员[NOIP2011]

时间:2017-08-04 19:39:35      阅读:178      评论:0      收藏:0      [点我收藏+]

标签:return   ems   函数   今天   frame   noi   个数   空格   开始   

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

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

这批矿产的检验结果Y为各个区间的检验值之和。即:Y=i=1mYi

若这批矿产的

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

 

【输入输出样例】

qc.in

5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3

qc.out

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。

 

【题解】

      对于我来说今天的题只有这一道相当可做,但是当时并没有意识到这一点……考试的时候想得太复杂,刚开始画绝对值函数想三分+线段树,打完了才发现线段树行不通,坚信这是一道线段树套Splay或无旋Treap区间求和的题,发现一来不太会打二来数据范围根本承受不了,然后万念俱灰直接弃疗打了个O(n^2)的暴力查询+三分。

       虽然这是个单峰函数,但是并不是必须三分。|y-s|与w的绝对值函数是单峰的,但是绝对值函数是根据y与w反相关的单调函数得来的,所以二分也能解决这个问题,只要寻找与最小值最接近的值即可。至于check,甚至不用什么数据结构,只要前缀和就能解决。把[l,r]转化成[1,r]-[1,l-1],再分别统计到每个位置大于w的矿石个数、价值和,可以在O(n)的时间内计算出取此w的答案。

       一道NOIP题,十分清楚十分简单,却被我生生做成了“强行数据结构,自己还不会打”。前缀和、区间 转化这种方法在ad爷的数数题里面已经出现过,当时没想到,现在还是没想到。要充分利用已有的知识,以之为工具解决问题,这里的工具不仅有算法和数据结构,更有许多转化变通的思想。

 

技术分享
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
const int sj=200010;
int n,m,jd,w[sj],v[sj],l[sj],r[sj],zj,yj,cnt[sj];
long long s,jg,vh[sj];
long long check(int x)
{
     memset(cnt,0,sizeof(cnt));
     memset(vh,0,sizeof(vh));
     for(int i=1;i<=n;i++)
     {
        cnt[i]=cnt[i-1]+((w[i]>=x)?1:0);
        vh[i]=vh[i-1]+((w[i]>=x)?v[i]:0);
     }
     long long ans=0;
     for(int i=1;i<=m;i++)
        ans+=(vh[r[i]]-vh[l[i]-1])*(cnt[r[i]]-cnt[l[i]-1]);
     return ans;
}
int main()
{
    scanf("%d%d%lld",&n,&m,&s);
    jg=s;
    for(int i=1;i<=n;i++)
    {
       scanf("%d%d",&w[i],&v[i]);
       if(w[i]>jd) jd=w[i];
    }
    for(int i=1;i<=m;i++)
       scanf("%d%d",&l[i],&r[i]);
    yj=jd;
    zj=1;
    while(zj<yj-1)
    {
        int mid=(zj+yj)>>1;
        long long ans;
        ans=check(mid);
        if(ans>=s) zj=mid;
        else yj=mid;
        if(abs(ans-s)<jg) jg=abs(ans-s);
        if(jg==0) break;
    }
    long long ans=check(zj);
    if(abs(ans-s)<jg) jg=abs(ans-s);
    ans=check(yj);
    if(abs(ans-s)<jg) jg=abs(ans-s);
    printf("%lld",jg);
    return 0;
}
qc

 

聪明的质监员[NOIP2011]

标签:return   ems   函数   今天   frame   noi   个数   空格   开始   

原文地址:http://www.cnblogs.com/moyiii-/p/7286406.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!