(题面保密,内部人员可览)
首先观察题面,可得出如下公式
∑(ceil(a[i] /d)*d−a[i])≤k
其中,ceil(a[i] /d)表示在需要被砍伐之前所经过的轮数,ceil函数是为了保证一定可以被砍伐,即轮数=(a[i]%d)?(a[i]/d):(a[i]/d+1)
轮数*d表示一共生长的长度,减去a[i]即为需要砍伐的长度,取sigma即为一共需砍伐的长度,与长度上限k比较,若小于等于k,即为符合条件。
推得公式后,大部分OIER的思路都是二分求解,但这个题不满足二分单调性,可以手动模拟一下,有如下几种方法:
1.暴力对拍
2.将check函数返回值输出,理想状态:“1 1 1 1 1 1 1 1 0 0 0 0 0 0”,实际状态:“1 0 1 1 0 1 0 0 1 1 0”(疯狂跳跃)。。
3.将二分左端点值赋作以求得的ans值,右端点不变,再次二分,发现ans值发生变化
于是对于此题而言,我们应该用什么方法求解呢?
首先我们化简公式
∑(ceil(a[i] /d))≤(k+∑(a[i]))/d
我们不难发现 k+∑(a[i]) 是一个定值,我们设之为limit
考虑函数单调性思想,我们发现右侧(k+∑(a[i]))/d的值应向下取整,因为若一个函数的最大值小于另一个函数的最小值,那么这个函数小于另一个函数恒成立,由此可得解
那么总结规律,随着d的不断增加(k+∑(a[i]))/d的值呈分段下降,而∑(ceil(a[i] /d))也同样分段下降,且区间长度远小于前者
于是不难看出,可以采用数论分块向下取整的思想来处理。
例如:floor(107/36)=2,floor(107/53)=2,36~53即为一段区间,那么如何求呢?公式如下
r=k/(k/l);
证明如下:
而这也正是此题的核心内容
由于题目要求求最大的值,那么我们只需要满足每一个区间的右端点满足条件即可,因为∑(ceil(a[i] /d))满足区间单调性,若一个区间,左端点满足条件,那么整个区间一定满足条件,若右端点不满足条件,那整体一定不满足,而一个区间的最优解,一定是右端点,因此不需考虑中间交点。
代码如下:
#include<bits/stdc++.h> #define re register #define ll long long using namespace std; ll n,k,limit,a[103]; inline ll read(){ re ll a=0,b=1;re char ch=getchar(); while(ch<'0'||ch>'9') b=(ch=='-')?-1:1,ch=getchar(); while(ch>='0'&&ch<='9') a=(a<<3)+(a<<1)+(ch^48),ch=getchar(); return a*b; } signed main() { n=read(),k=read();limit=k; for(re ll i=1;i<=n;++i) a[i]=read(),limit+=a[i]; for(re ll d=limit;d>=0;d--) { re ll flag=limit/d,tot=0; for(re ll i=1;i<=n;++i) tot+=(a[i]-1)/d+1; if(tot<=flag) {printf("%lld\n",d);return 0;} d=(limit/(limit/d+1))+1; } return 0; }