标签:src dash blank 应用 组合 http 简单 png 复杂
额,前两天刚讲了数据结构,今天我来讲讲组合数学中的一种奇妙优化——Lucas
先看这样一个东西
没学过lucas的肯定会说:还不简单?处理逆元,边乘边膜呗
是,可以,但注意一下数据范围
你算这一次,你需要跑25000下
那么你如果求C199999 1~C199999 52222 呢?
你会发现你的复杂度上天了
所以我们会用到一个神奇的定理:Lucas定理
定理内容如下:
Lucas(n,m,p)=c(n%p,m%p)*Lucas(n/p,m/p,p)
不好玩,是吗?
那么我来证明一下
由二项式定理可得,cmn等于(x+1)^m中n次项的系数
那么我们按Lucas展开
原式=(x+1)^(p^k*ak)*(x+1)^(p^(k-1)*a(k-1))*……*(x+1)^(p*a1)+(x+1)^(1*a0)
=[(x+1)^(p^k)]^ak……
由费马小定理可知,其%p后可转化为 (x+1)^ak*(x+1)^a(k-1)*……*(x+1)^a1*(x+1)^a0
原题转化为求上式的n次项的系数
同理,由于每一项已经消去了p^k次方 故即求每一项中的bk次项系数 即为Lucas
以前的组合数,我们一位一位地算(累死了)
现在的组合数,我们只用算%p出来的(这不就是log p m次吗?)
O(n)到O(logn)
大有长进啊
大家应该都会
下面看几道好题
裸的板子,不说什么
带了预处理阶乘的Lucas
lucas合并
礼物,古代猪文
扩展Lucas
标签:src dash blank 应用 组合 http 简单 png 复杂
原文地址:https://www.cnblogs.com/ztz11/p/9174296.html