做砍树时问大佬说,“这是一个数论分块”模板题
我:???
原来只有我没学过数论分块吗?
https://www.cnblogs.com/0xfffe/p/9648943.html
略微理解了理解,写的非常清楚
你说这是向下取整,不是向上取整,砍树要向上取整,那篇博客不适用于砍树?
确实不适用
我们要分块的是等式右面的$\sum_{i}^{n} a[i] +k$除以d
因为C是固定的,所以这是一个分段函数,我们要处理的是不同的右面的值最后再跟左面对应
我们然后f存下这一段具体的值,
r存下具体右端点
然后就完了
代码
![](https://i-blog.csdnimg.cn/blog_migrate/8f900a89c6347c561fdf2122f13be562.gif)
![](https://i-blog.csdnimg.cn/blog_migrate/961ddebeb323a10fe0623af514929fc1.gif)
#include<bits/stdc++.h> using namespace std; #define ll long long #define A 11000000 ll l=1,r,n,m,a[A],dl[A],R[A],f[A],zz=0,num=0,ans=0,sum=0; void precl() { while(1){ if(!(sum/l)) break; r=sum/(sum/l); f[++num]=sum/r; R[num]=r; l=r+1; } } int main(){ scanf("%lld%lld",&n,&m); sum=m; for(ll i=1;i<=n;i++){ scanf("%lld",&a[i]); sum+=a[i]; } precl(); for(ll j=1;j<=num;j++) { ll t=0; for(ll i=1;i<=n;i++){ t+=ceil((double)a[i]/(double)R[j]); } // printf("f=%lld R=%lld\n",f[j],R[j]); if(t<=f[j]) ans=R[j]; } cout<<ans<<endl; }
以下是我完全错误的解释
设$k\times i-p=N$ 向上取整设
$\large{\lceil \frac N{i+d} \rceil}=k$
于是$k\times (i+d)-p2=N$
同样得出p2=p+kd
就是照猫画虎的一个过程
底下我不具体推了,
$\large \left \lceil \frac N{\left \lfloor \frac Ni \right \rfloor } \right \rceil$
所以对砍树这道题来说,这确实是个模板题,分析发现这是一个分段函数,维护每一段大小相同,维护l,r下一个l=r+1
具体来说
$\large \left \lceil \frac {a[i]}{d} \right \rceil$不是为我们具体分块的值
$\large \lfloor \frac Ni \rfloor$才是
然后等式右面是$\sum_{i}^{n} a[i] +k$再除以d
这个N就是$\sum_{i}^{n} a[i] +k$
那么这个题就迎刃而解了。