CodeForces830C-Round#424 Div1 C Solution:分块导学&膜QC霸霸

题意：有n个竹子，每个竹子有一个理想高度，竹子刚开始高度都是0，每过一天就增高1，你每隔Q天会来查看一次竹子的情况，如果竹子的高度超过（或等于）他的理想高度，那么就给他在理想高度那里来一刀，并且这个竹子停止生长。于是你最终会得到T单位的废弃竹子，现在要求在保证T<=K的情况下，Q的最大值。题目给出n,K和a[i]（理想高度）。

题解：我们考虑：假设我现在想q天来查看一次，那么每次查看的时候，那些等待我砍刀的竹子高度肯定是q的倍数q*s，如果q*s>=a[i]那么我就砍这一刀，以后不再考虑他。并且我们可以容忍k单位的下脚料。那么如果q这个答案是合法的，设砍下来的废弃竹子有P单位，则P<=K，我们把这P单位的竹子接回去，所有的竹子的高度都变成了q的倍数，设这些竹子的自然高度总和是SUM，有SUM<=K+∑a[i]，那么我们其实明白了，我们先求出MAX=K+∑a[i]，然后我们这个q把MAX分成MAX/q块（不足q的那个块没用），每个竹子占据了(a[i]+q-1)/q块，我们只需要判断在当前使用q来分块的方案下，能不能把所有的理想高度都放到块里边就行了。

但是我们如果暴力的从1-SUM来枚举q会爆炸的。这里用到了分块的根号优化。（膜QC霸霸）

先来看k/（k/i）是个什么东西（QC霸霸写的好哇）：https://zhuanlan.zhihu.com/p/26020672

我们考虑有块大小a和b，a能够分SUM为Na=SUM/a块，b为Nb=SUM/b块，如果说Na==Nb，也就是说分出来的块数相同（不足的那一块都没用），显然单块更大的要优于单块小的，因为单块更大，有可能使得竹子占的块数变少，更可能使得解合法。

Code：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
int n, a[111];
ll ans, k;
bool chk(ll d, ll k){
    ll sum = 0;
    for(int i = 1; i <= n; ++i) sum += (a[i] + d - 1) / d;
    return sum <= k;
}

int main(){
    scanf("%d%I64d",&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){
        r = k / (k / l);
        if (chk(r, k / r)) ans = max(ans, r);
    }
    printf("%I64d",ans);
    return 0;
}