[分块][数学][瞎搞]Codeforces Round #424 .C Bamboo Partition

刚开始的做法是对的…然而越做头越昏，竟然用错误的例子把自己叉掉…
还好最后半个小时清醒了

先列出式子

要求最大的d满足

∑i=1nd−((ai−1)%d+1)≤k

其中ai要先减一再取模是为了防止ai是d的倍数的情况，

再推一推

n×d−∑i=1nai−⌊ai−1d⌋×d≤k

n×d+d∑i=1n⌊ai−1d⌋≤k+∑i=1nai

大概就是这样
然后发现，⌊ai−1d⌋总共有O(a√i)种取值
那么所有⌊ai−1d⌋的取值就是O(nm−−√)，其中m是ai的最大值
然后从大到小枚举每一种取值，显然答案在是第一种满足条件的取值。
然后把这个取值转化成最大的d就可以了，具体看代码
这样复杂度是O(n2m−−√)，看起来很大的样子…跑得过就行

#include <cstdio>#include <iostream>#include <algorithm>#include <vector>#include <map>#include <set>#include <cmath>#define fi first#define se second#define pb push_backusing namespace std;typedef long long ll;typedef pair<int,int> par;inline char nc(){    static char buf[100000],*p1=buf,*p2=buf;    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;}inline void rea(int &x){    char c=nc(); x=0; int f=1;    for(;c>'9'||c<'0';c=nc())if(c=='-') f=-1;    for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc()); x*=f;}inline void rea(ll &x){    char c=nc(); x=0; int f=1;    for(;c>'9'||c<'0';c=nc())if(c=='-') f=-1;    for(;c>='0'&&c<='9';x=x*10+c-'0',c=nc()); x*=f;}const int N=110;int n,a[N];vector<int> v;ll goal,k;int main(){    rea(n); rea(k);    for(int i=1;i<=n;i++)        rea(a[i]),goal+=a[i],a[i]--;    goal+=k;    for(int i=1;i<=n;i++)        for(int j=1,d;j<=a[i];j=d+1)            d=a[i]/(a[i]/j),v.pb(j);    for(int i=1;i<=n;i++) v.pb(a[i]+1);    sort(v.begin(),v.end());     //for(int i=0;i<v.size();i++) printf("%d ",v[i]); putchar('\n');    v.resize(unique(v.begin(),v.end())-v.begin());    //for(int i=0;i<v.size();i++) printf("%d ",v[i]); putchar('\n');    int td,last=1<<30;    for(int s=v.size()-1;~s;s--){        int d=v[s]; ll cur=0;        for(int i=1;i<=n;i++){            cur+=d*(a[i]/d);        }        if(1LL*n*d+cur<=goal){            td=d; break;        }        last=d;    }    ll cur=0;    for(int i=1;i<=n;i++) cur+=(a[i]/td);    printf("%lld\n",goal/(n+cur));    return 0;}