P1314 聪明的质监员
2018-06-17 22:19:58来源:未知 阅读 ()
题目描述
小T 是一名质量监督员,最近负责检验一批矿产的质量。这批矿产共有 n 个矿石,从 1到n 逐一编号,每个矿石都有自己的重量 wi 以及价值vi 。检验矿产的流程是:
1 、给定m 个区间[Li,Ri];
2 、选出一个参数 W;
3 、对于一个区间[Li,Ri],计算矿石在这个区间上的检验值Yi:
这批矿产的检验结果Y 为各个区间的检验值之和。即:Y1+Y2...+Ym
若这批矿产的检验结果与所给标准值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 。
一开始写了个莫队。。发现莫队只能过样例。。。
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<cmath> 5 #include<queue> 6 #include<algorithm> 7 #define lli long long int 8 using namespace std; 9 const lli MAXN=200001; 10 void read(lli &n) 11 { 12 char c='+';lli x=0;bool flag=0; 13 while(c<'0'||c>'9') 14 {c=getchar();if(c=='-')flag=1;} 15 while(c>='0'&&c<='9') 16 {x=x*10+(c-48);c=getchar();} 17 flag==1?n=-x:n=x; 18 } 19 struct node 20 { 21 lli w,v; 22 }kuang[MAXN]; 23 struct xw 24 { 25 lli l,r,id; 26 }q[MAXN]; 27 lli n,m,s; 28 lli ans=0x7ffffff; 29 lli pos[MAXN]; 30 lli base; 31 lli num,va; 32 lli comp(const xw &a,const xw &b) 33 { 34 if(pos[a.l]==pos[b.l]) 35 return a.r<b.r; 36 else 37 return pos[a.l]<pos[b.l]; 38 } 39 void add(lli p,lli need) 40 { 41 if(kuang[p].w>=need) 42 { 43 num++; 44 va+=kuang[p].v; 45 } 46 } 47 void dele(lli p,lli need) 48 { 49 if(kuang[p].w>=need) 50 { 51 num--; 52 va-=kuang[p].v; 53 } 54 } 55 lli pd(lli need) 56 { 57 lli ll=1,rr=0; 58 lli now=0; 59 num=0;va=0; 60 for(lli i=1;i<=m;i++) 61 { 62 for(;ll<q[i].l;ll++) 63 add(ll-1,need); 64 for(;ll>q[i].l;ll--) 65 dele(ll,need); 66 for(;rr>q[i].r;rr--) 67 dele(rr,need); 68 for(;rr<q[i].r;rr++) 69 add(rr+1,need); 70 now+=num*va; 71 } 72 if(now-s<ans&&now-s>0) 73 { 74 ans=now-s; 75 return 1; 76 } 77 else return 0; 78 } 79 int main() 80 { 81 read(n);read(m);read(s); 82 base=sqrt(n); 83 for(lli i=1;i<=n;i++) 84 pos[i]=(i-1)/base+1; 85 for(lli i=1;i<=n;i++) 86 {read(kuang[i].w);read(kuang[i].v);} 87 for(lli i=1;i<=m;i++) 88 { 89 read(q[i].l);read(q[i].r); 90 } 91 sort(q+1,q+m+1,comp); 92 93 lli ll=0,rr=5; 94 for(lli i=-500;i<=500;i++) 95 pd(i); 96 printf("%lld",ans); 97 return 0; 98 }
后来听大佬说直接暴力然后用前缀和优化一下就好
1 #include<iostream> 2 #include<cstdio> 3 #include<cstring> 4 #include<cmath> 5 #include<queue> 6 #include<algorithm> 7 #define lli long long int 8 using namespace std; 9 const lli MAXN=300001; 10 void read(lli &n) 11 { 12 char c='+';lli x=0;bool flag=0; 13 while(c<'0'||c>'9') 14 {c=getchar();if(c=='-')flag=1;} 15 while(c>='0'&&c<='9') 16 {x=x*10+(c-48);c=getchar();} 17 flag==1?n=-x:n=x; 18 } 19 struct node 20 { 21 lli w,v; 22 }kuang[MAXN]; 23 struct xw 24 { 25 lli l,r,id; 26 }q[MAXN]; 27 lli n,m,s; 28 lli ans=0x7fffffff; 29 lli sum[MAXN],num[MAXN]; 30 lli tot=0; 31 lli pd(lli W) 32 { 33 memset(num,0,sizeof(num)); 34 memset(sum,0,sizeof(sum)); 35 for (lli i=1;i<=n;i++) 36 { 37 sum[i]=sum[i-1]; 38 num[i]=num[i-1]; 39 if (kuang[i].w>=W) 40 { 41 sum[i]+=kuang[i].v; 42 num[i]++; 43 } 44 } 45 tot=0; 46 for (lli i=1;i<=m;i++) 47 tot+=(sum[q[i].r]-sum[q[i].l-1])*(num[q[i].r]-num[q[i].l-1]); 48 if (tot<=s) return true; 49 else return false; 50 } 51 int main() 52 { 53 //freopen("qc.in","r",stdin); 54 //freopen("qc.out","w",stdout); 55 read(n);read(m);read(s); 56 for(lli i=1;i<=n;i++) 57 { 58 read(kuang[i].w);read(kuang[i].v); 59 } 60 for(lli i=1;i<=m;i++) 61 { 62 read(q[i].l);read(q[i].r); 63 } 64 65 lli ll=0,rr=1000000; 66 while(ll<=rr) 67 { 68 lli mid=(ll+rr)>>1; 69 if(pd(mid)) 70 { 71 rr=mid-1; 72 ans=mid; 73 } 74 else ll=mid+1; 75 } 76 ll=pd(ans); 77 lli a=tot; 78 rr=pd(ans-1); 79 lli b=tot; 80 printf("%lld",min(s-a,b-s)); 81 return 0; 82 }
标签:
版权申明:本站文章部分自网络,如有侵权,请联系:west999com@outlook.com
特别注意:本站所有转载文章言论不代表本站观点,本站所提供的摄影照片,插画,设计作品,如需使用,请与原作者联系,版权归原作者所有
- P1816 忠诚 2018-06-17
IDC资讯: 主机资讯 注册资讯 托管资讯 vps资讯 网站建设
网站运营: 建站经验 策划盈利 搜索优化 网站推广 免费资源
网络编程: Asp.Net编程 Asp编程 Php编程 Xml编程 Access Mssql Mysql 其它
服务器技术: Web服务器 Ftp服务器 Mail服务器 Dns服务器 安全防护
软件技巧: 其它软件 Word Excel Powerpoint Ghost Vista QQ空间 QQ FlashGet 迅雷
网页制作: FrontPages Dreamweaver Javascript css photoshop fireworks Flash