4830: [Hnoi2017]抛硬币

题意:A投a次硬币,B投b次硬币,a比b正面朝上次数多的方案数,模\(10^k\)。

\(b \le a \le b+10000 \le 10^{15}, k \le 9\)


几乎一下午和一晚上杠这道题...中间各种翻《具体数学》各种卡常

有两种做法,这里只说我认为简单的一种。


题目就是要求

\[\sum_{i=0}^a \sum_{j=0}^b [i>j] \binom{a}{i} \binom{b}{j}
\]

化一化得到

\[\sum_{i=1}^a \sum_{j=0}^{a-i} \binom{a}{i+j} \binom{b}{b-j}
\]

根据范德蒙德卷积,得到

\[\sum_{i=b+1}^{a+b} \binom{a+b}{i}
\]

然而并不会计算...根据题目范围显然应该与a,b的差相关...


参考了这里

把这个组合数拆成两部分计算:

\[\sum_{i=\lceil \frac{a+b}{2} \rceil}^{a+b} \binom{a+b}{i}+ \sum_{i=b+1}^{\lfloor \frac{a+b}{2} \rfloor} \binom{a+b}{i}
\]

\(a+b\)为奇数时前一部分就是\(2^{a+b-1}\)。

注意判断\(a+b\)为偶数的时候,要减去\(\frac{1}{2} \binom{a+b}{(a+b)/2}\)


因为模数特殊需要使用扩展lucas定理

但是\(10^9\)直接用会T...

10只有2和5两个质因子,预处理不能整除2或5的数的阶乘就没问题了。


关于模意义再提两点:

  1. \(a \equiv b \pmod {md} \rightarrow a\equiv b \pmod m\),所以可以放在\(10^9\)下做
  2. 组合数除以2,对于\(2^a\)这个质因子2的逆元不存在,我们只能把指数减1.然后2的指数有可能本来就为0,所以不是所有组合数都能除以2(有的本来就是奇数啊),可以证明\(\binom{2n}{n}\)为偶数


最后还有一个卡常技巧

如果\(p\)指数大于9不要继续计算了...节省好多常数

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <ctime>
using namespace std;
typedef long long ll;
const int N=2e6+5;
int mo, p2, p5; ll a, b; int k;
inline ll Pow(ll a, ll b, ll mo) {
ll ans=1;
for(; b; b>>=1, a=a*a%mo)
if(b&1) ans=ans*a%mo;
return ans;
}
inline ll Inv(ll a, ll n) {return Pow(a, ((n&1) ? n/5*4 : n/2) - 1, n);} ll fac2[N], fac5[N];
inline ll Fac(ll n, ll p, ll pr) {
if(n == 0) return 1;
ll ans = 1;
if(p == 2) ans = fac2[pr] %pr; else ans = fac5[pr] %pr;
ans = Pow(ans, n/pr, pr);
if(p == 2) ans = ans * fac2[n%pr] %pr; else ans = ans * fac5[n%pr] %pr;
return ans * Fac(n/p, p, pr) %pr;
}
ll cou(ll n, ll p) { ll ans = 0; for(ll i=n; i; i/=p) ans += i/p; return ans; } ll C(ll n, ll m, bool div2=0) {
if(n < m) return 0; ll ans = 0, a = 0, t;
{
ll p = 2, pr = p2;
ll c = cou(n, p) - cou(m, p) - cou(n-m, p);
if(div2) t = Pow(p, --c, pr);
else t = Pow(p, c, pr); if(c<9) {
a = Fac(n, p, pr) * Inv(Fac(m, p, pr), pr) %pr * Inv(Fac(n-m, p, pr), pr) %pr * t %pr;
a = a * (mo/pr) %mo * Inv(mo/pr, pr) %mo;
ans = (ans + a) %mo;
}
}
{
ll p = 5, pr = p5;
ll c = cou(n, p) - cou(m, p) - cou(n-m, p);
t = Pow(p, c, pr);
if(div2) t = t * Inv(2, pr) %pr; if(c<9) {
a = Fac(n, p, pr) * Inv(Fac(m, p, pr), pr) %pr * Inv(Fac(n-m, p, pr), pr) %pr * t %pr;
a = a * (mo/pr) %mo * Inv(mo/pr, pr) %mo;
ans = (ans + a) %mo;
}
} return ans;
} char c[50];
void print(ll ans) {
for(int i=0; i<20; i++) c[i] = '0';
int p=0;
while(ans) c[++p] = ans%10 + '0', ans /= 10;
for(int i=k; i>=1; i--) putchar(c[i]); puts("");
} void init(int k) {
fac2[0] = 1; p2 = Pow(2, k, 1e9);
for(int i=1; i<=p2; i++) fac2[i] = fac2[i-1] * (i%2 ? i : 1) %p2; fac5[0] = 1; p5 = Pow(5, k, 1e9);
for(int i=1; i<=p5; i++) fac5[i] = fac5[i-1] * (i%5 ? i : 1) %p5;
} int main() {
freopen("in", "r", stdin);
init(9); mo = 1e9;
while(scanf("%lld %lld %d", &a, &b, &k) != EOF) {
ll ans = Pow(2, a+b-1, mo);
for(ll i=b+1; i<=(a+b)/2; i++) ans = (ans + C(a+b, i)) %mo;
if(~(a+b)&1) ans = (ans - C(a+b, (a+b)>>1, 1) + mo) %mo;
int t = 1; for(int i=1; i<=k; i++) t *= 10; ans %= t;
print(ans);
}
//printf("time %lf\n", (double) clock()/CLOCKS_PER_SEC);
}

bzoj 4830: [Hnoi2017]抛硬币 [范德蒙德卷积 扩展lucas]的更多相关文章

  1. bzoj 4830: [Hnoi2017]抛硬币

    Description 小A和小B是一对好朋友,他们经常一起愉快的玩耍.最近小B沉迷于**师手游,天天刷本,根本无心搞学习.但是 已经入坑了几个月,却一次都没有抽到SSR,让他非常怀疑人生.勤勉的小A ...

  2. 【刷题】BZOJ 4830 [Hnoi2017]抛硬币

    Description 小A和小B是一对好朋友,他们经常一起愉快的玩耍.最近小B沉迷于**师手游,天天刷本,根本无心搞学习.但是已经入坑了几个月,却一次都没有抽到SSR,让他非常怀疑人生.勤勉的小A为 ...

  3. 【题解】幼儿园篮球题(范德蒙德卷积+斯特林+NTT)

    [题解]幼儿园篮球题(NTT+范德蒙德卷积+斯特林数) 题目就是要我们求一个式子(听说叫做超几何分布?好牛逼的名字啊) \[ \sum_{i=1}^{S}\dfrac 1 {N \choose n_i ...

  4. 【BZOJ4830】[HNOI2017]抛硬币(组合计数,拓展卢卡斯定理)

    [BZOJ4830][HNOI2017]抛硬币(组合计数,拓展卢卡斯定理) 题面 BZOJ 洛谷 题解 暴力是啥? 枚举\(A\)的次数和\(B\)的次数,然后直接组合数算就好了:\(\display ...

  5. luogu P3726 [AH2017/HNOI2017]抛硬币

    传送门 我是真的弱,看题解都写了半天,,, 这题答案应该是\(\sum_{i=1}^{a}\binom{a}{i}\sum_{j=0}^{min(b,i-1)}\binom{b}{j}\) 上面那个式 ...

  6. 题解 洛谷 P3726 【[AH2017/HNOI2017]抛硬币】

    可以分别枚举两人正面朝上的次数来统计答案,所求即为 \[\sum_{i=0}^{a}\sum_{j=0}^{b} \binom{a}{i} \binom{b}{j} [i>j] \] 将\(i\ ...

  7. [AH2017/HNOI2017]抛硬币

    传送门 这个题的暴力比较好想--然后用一些组合的知识就可以变成正解了. 首先我们考虑a=b的情况.我们把扔出来的硬币看成是一个01序列,那么对于一个b获胜的序列,他在每一位都按位异或1之后必然是一个a ...

  8. Codeforces 785D - Anton and School - 2 - [范德蒙德恒等式][快速幂+逆元]

    题目链接:https://codeforces.com/problemset/problem/785/D 题解: 首先很好想的,如果我们预处理出每个 "(" 的左边还有 $x$ 个 ...

  9. 浅谈范德蒙德(Vandermonde)方阵的逆矩阵的求法以及快速傅里叶变换(FFT)中IDFT的原理

    浅谈范德蒙德(Vandermonde)方阵的逆矩阵与拉格朗日(Lagrange)插值的关系以及快速傅里叶变换(FFT)中IDFT的原理 标签: 行列式 矩阵 线性代数 FFT 拉格朗日插值 只要稍微看 ...

随机推荐

  1. C语言dsPIC / PIC24 serial bootloader和C#语言bootloader PC端串口通信程序

    了解更多关于bootloader 的C语言实现,请加我QQ: 1273623966 (验证信息请填 bootloader),欢迎咨询或定制bootloader(在线升级程序). 新dsPIC/PIC2 ...

  2. (转)CVPR 2016 Visual Tracking Paper Review

    CVPR 2016 Visual Tracking Paper Review  本文摘自:http://blog.csdn.net/ben_ben_niao/article/details/52072 ...

  3. 运维工作中sed常规操作命令梳理

    sed是一个流编辑器(stream editor),一个非交互式的行编辑器.它一次处理一行内容.处理时,把当前处理的行存储在临时缓冲区中,称为"模式空间",接着用sed命令处理缓冲 ...

  4. CPU是怎么制造的

    大概的过程就是,先选一堆好沙子(纯净的沙子),初步加工一般在沿海,然而都是初加工,因为技术不行,所以一般用比较污染环境的方法加工大99.9%纯度的硅,然后低价卖给国外企业,用高精尖技术加工到99.99 ...

  5. Java和.NET在开发中的不同盘点

    我是用VS2008和VS2010开发.NET程序,通过MyEclipse8.5开发JAVA程序,下面从IDE.语言.插件的不同点来做下简单的说明.但由于经验知识还有限,本篇文章只能从比较表面的以及自己 ...

  6. keepalived VS zookeeper

    转载请标明出处http://www.cnblogs.com/haozhengfei/p/e3db73cb83afb213a3bff43a850d56c4.html keepalived VS zook ...

  7. Java 读取 Json格式的 内容

    一.Json 报文格式如下: 二.获取 Json 报文中字段的内容 import java.io.IOException; import com.fasterxml.jackson.core.Json ...

  8. learning coap protocol

    reference: http://coap.technology/spec.html

  9. 21.QT-QTreeWidget,QTabWidget

    QTreeWidget树形列表 设置标签相关函数 void QTreeWidget::setHeaderItem (QTreeWidgetItem * item ); void QTreeWidget ...

  10. Python-Django&#160;整合Django和jquery-easyui

    整合Django和jquery-easyui by:授客 QQ:1033553122 测试环境 win7 64 Python 3.4.0 jquery-easyui-1.5.1 下载地址1:http: ...