小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。
输出只有一行,包含一个整数,表示所求的最小值。
输入样例#1:
5 3 15
1 5
2 5
3 5
4 5
5 5
1 5
2 4
3 3
输出样例#1:
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 }