BZOJ 2600: [Ioi2011]ricehub

2600: [Ioi2011]ricehub

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 628  Solved: 325
[Submit][Status][Discuss]

Description

乡间有一条笔直而长的路称为“米道”。沿着这条米道上 R 块稻田,每块稻田的坐标均
为一个 1 到 L 之间(含 1 和 L)的整数。这些稻田按照坐标以不减的顺序给出,即对于 0 ≤ i <
R,稻田 i 的坐标 X[i]满足 1 ≤ X[0] ≤ ... ≤ X[R-1] ≤ L。 
注意:可能有多块稻田位于同一个坐标上。 
我们计划建造一个米仓用于储存尽可能多的稻米。和稻田一样,米仓将建在米道上,其
坐标也是一个 1 到 L 之间的整数(含 1 和 L)。这个米仓可以建在满足上述条件的任一个位
置上,包括那些原来已有一个或多个稻田存在的位置。 
在收获季节,每一块稻田刚好出产一滿货车的稻米。为了将这些稻米运到米仓,需要雇
用一位货车司机来运米。司机的收费是每一满货车运送一个单位的距离收取 1 元。換言之,
将稻米从特定的稻田运到米仓的费用在数值上等于稻田坐标与米仓坐标之差的绝对值。 
不幸的是,今年预算有限,我们至多只能花费 B 元运费。你的任务是要帮我们找出一个
建造米仓的位置,可以收集到尽可能多的稻米。 

Input

第一行 三个整数 R L B
接下来R行 每行一个整数 表示X[i]

Output


一个整数 最多稻米数

Sample Input


5 20 6
1
2
10
12
14

Sample Output


3
HINT
1 ≤ R ≤ 100,000
1 ≤ L ≤ 1,000,000,000
0 ≤ B ≤ 2,000,000,000,000,000

这个题目居然在我们学校讲课时被归类为了二分(没认真听课QAQ不知道是怎么兹磁的啊)%%
显然易得出我们选取的米仓一定收集的一段连续的区间,而米仓一定是在这段区间的中位数上。
我们将右区间r从1到R枚举,每一次维护的左区间l。
记res为剩余的钱数。如果res小于0则左区间右移(即序列最早入队的稻田出队),直到res大于0。中位数可以O(1)的转移。
 
关于复杂度:因为每一个稻田至多只有入队和出队两次操作,所以复杂度为O(n)

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstdlib>
 4 #include<cmath>
 5 #include<algorithm>
 6 #include<cmath>
 7 #include<queue>
 8 #define llg long long
 9 #define yyj(s) freopen(s".in","r",stdin),freopen(s".out","w",stdout);
10 #define maxn 100010
11 using namespace std;
12 
13 llg sc()
14 {
15     llg i=0;char c=getchar();
16     while(c>'9'||c<'0')c=getchar();
17     while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();
18     return i;
19 }
20 
21 llg x,i,j,k,n,l,r,mid,m,a[maxn],ans;
22 double res,b,wz,lastwz;
23 int main()
24 {
25     yyj("a");
26     cin>>n>>m>>b;
27     for (i=1;i<=n;i++) a[i]=sc();
28     res=b;
29     x=1;  
30     ans=1; 
31     l=0; lastwz=a[1];
32     for (r=2;r<=n;r++)
33         {
34             if (a[r]>m) break;
35             x++;
36             if (x % 2) wz=a[l+x/2+1]; else wz=(double)(a[l+x/2]+a[l+x/2+1])/2;
37             if (x % 2==0) res-=(wz-lastwz);
38             res-=a[r]-wz; lastwz=wz;
39             while (res<0)
40                 {
41                     x--; l++; res+=lastwz-a[l];
42                     if (x % 2) wz=a[l+x/2+1]; else wz=(a[l+x/2]+a[l+x/2+1])/2;
43                     if (x % 2) res+=wz-lastwz;
44                     lastwz=wz;
45                 }
46             ans=max(ans,x);
47         }
48     cout<<ans;
49     return 0;
50 }
 
 
 
 
本文作者:xrdog 作者博客:http://www.cnblogs.com/Dragon-Light/ 转载请注明出处,侵权必究,保留最终解释权!
原文地址:https://www.cnblogs.com/Dragon-Light/p/5671105.html