首页 > 代码库 > P1314 聪明的质监员
P1314 聪明的质监员
题目描述
小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 node20 {21 lli w,v;22 }kuang[MAXN];23 struct xw24 {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 node20 {21 lli w,v;22 }kuang[MAXN];23 struct xw24 {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 }
P1314 聪明的质监员