【数学分析+除法分块】Codeforces830C[Bamboo Partition]题解

题目概述

有 $n$ 个竹子，每个竹子刚开始高度为 $0$ ，每过一天会长高 $1$ ，现在每过 $d$ 天来检查一次，如果 $i$ 竹子的高度 $\ge a_i$ 则将该竹子高度砍为 $a_i$ 且该竹子不再生长，求砍掉高度不超过 $k$ 时 $d$ 最大为多少。

解题报告

除法分块性质： $\lfloor \frac{n}{i}\rfloor$ 排序去重之后是 $O(\sqrt{n})$ 的。

简单说明：一个数的因子是 $O(\sqrt n)$ 的。

---

先推式子喽：

$$\sum_{i=1}^{n}\{d-[(a_i-1)\ mod \ d+1]\}\le k\\ \sum_{i=1}^{n}[d-(a_i-1-d\lfloor \frac{a_i-1}{d}\rfloor+1)]\le k\\ nd-\sum_{i=1}^{n}(a_i-d\lfloor \frac{a_i-1}{d}\rfloor)\le k\\ nd+d\sum_{i=1}^{n}\lfloor \frac{a_i-1}{d}\rfloor\le k+\sum_{i=1}^{n}a_i$$

根据除法分块的性质， $\lfloor \frac{a_i-1}{d}\rfloor$ 只有 $\sqrt{a_i}$ 种情况，而 $n$ 个则有 $n\sqrt{a_{max}}$ 种情况，所以暴力枚举然后 $O(n)$ 验证就可以了，效率 $O(n^2\sqrt{a_{max}})$ ，好像有点慢……而且有较多的细节处理……

---

不妨继续推：

$$d\sum_{i=1}^{n}\lfloor \frac{a_i+d-1}{d}\rfloor\le k+\sum_{i=1}^{n}a_i\\ \sum_{i=1}^{n}\lfloor \frac{a_i+d-1}{d}\rfloor\le \frac{k+\sum_{i=1}^{n}a_i}{d}\\ \sum_{i=1}^{n}\lfloor \frac{a_i+d-1}{d}\rfloor\le \lfloor \frac{k+\sum_{i=1}^{n}a_i}{d}\rfloor$$
由于右边只有 $\sqrt{k+\sum_{i=1}^{n}a_i}​$ 块，而每块中 $d​$ 越大，左边的式子越小（右边不变），所以每块中有价值的只有最大的 $d​$ ，那么还是暴力枚举即可。效率 $O(n\sqrt{k+\sum_{i=1}^{n}a_i})​$ ，比之前的要好一些，也没什么细节处理。

示例程序

我真的没有在水博客……

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long LL;
const int maxn=100;

int n,a[maxn+5];LL K,ans;

bool check(LL d) {LL sum=0;for (int i=1;i<=n;i++) sum+=(a[i]-1)/d+1;return sum<=K/d;}
int main()
{
    freopen("program.in","r",stdin);
    freopen("program.out","w",stdout);
    scanf("%d%lld",&n,&K);for (int i=1;i<=n;i++) scanf("%d",&a[i]),K+=a[i];
    for (LL l=1,r;l<=K;l=r+1) if (check(r=K/(K/l))) ans=r;
    return printf("%lld\n",ans),0;
}