【转载】【转自AekdyCoin的组合数取模】

原创 2015年07月07日 18:47:38
本篇文章主要介绍了"【组合数求模】 转自AekdyCoin",主要涉及到【组合数求模】 转自AekdyCoin方面的内容,对于【组合数求模】 转自AekdyCoin感兴趣的同学可以参考一下。

这个表示的是从n个元素中选取m个元素的方案数。

(PS.组合数求模似乎只用在信息学竞赛和 ACM竞赛等计算机编程设计大赛中……,求在现实中的运用)

可以知道当n,m 取得比较大的时候,组合数可能很大很大 (天文数字?无法度量?)

例如 C(100, 50) = 100891344545564193334812497256, 于是计算机的 64位整数型已经没法阻止它了!C(1000000000, 500000000) ? C( 2^50000, 2^49999 ) ? (Note:这里^表示次方,你能计算得到2的50000次级别的组合数么?它有多少位?)

看起来似乎高精度神马的都无法阻止这个邪恶的函数的急速扩张了……

庆幸的是,在竞赛中我们能够遇到的规模也就只有10^9级别(显然是mod上某个数字K,否则输出的文件那叫一个大啊……),这是多么的小呀呀呀呀!(Note: 相比较2^50000 -_-)

一.           入门篇:我会暴力!

(1)  K = 1: 今天你学数论了么? 难度系数: 0

(2)  (K> 1) n, m <= 1000 (n * n 是可以接受的) 难度系数: 1

递推!

   c(n,m) =c(n - 1,m) + c(n – 1, m – 1)

某人: 555555 这个公式太复杂, 记忆不能!

c(5,2) = 10 = c(4,2) + c(4,1) = 6 + 4 ……

我们知道mod操作满足加法性质,即

(a + b) mod c = ( (a mod c) + (b mod c) ) mod c

c(n,m) = ( c(n - 1,m) + c(n – 1, m – 1) ) mod K

证明利用模的定义即可……很简单的

于是如此,我们只需要简单的开上一个 f[ N ][ N ],2个循环搞定!

其实我们遇到的大部分情况需要的 组合数 都可以用这个来搞定~

这里唯一可能被邪恶的其实是 K + K 溢出! 所以如果某个邪恶的题目出到 K = 2*10^9,在某些倒霉的场合会出现2个接近K的int相加,那么就溢出了!不要忘记用unsigned int… (我从来没出过这种题的!真的!)

(3)   n 巨大(10^9 级别), m巨小(10^4级别), k 很小,大约10^9

a)     m<= 1: 今天你学数论了么? 难度系数: 0

b)     m<= 10000     难度系数: 2

 

可以发现分子分母的项数都少到可以接受!于是我们可以采取各种方式来通过:

i)                   对于每个数字,分解素因子,合并,二分求幂! (你会数论!)

ii)                 对于每个数字,只分解包含于K的素因子,例如K里面有一个素因子3,那么分解的时候我只考虑3呀,因为其他部分显然与3互质……最后统计3的次数即可……

例子:

计算C(10, 3) mod 36

C(10, 3) = (10 * 9 * 8) / (1 * 2 * 3)

对于分母:

1 : ok 逆元(有区别么?)

2: 没法逆元, (2, 36) = 2

3: 没法逆元, (3,36) = 3

为了神马啊!! 还不让人逆元啊!显然是因为邪恶的2和3,如果他们不存在,那么多么美好呀!

于是我开2个变量,记录2,3的次数

对于分子:

10: 里面只有1个2,去掉了2,剩下的部分是 10 / 2 = 5. 

9: 里面只有2个3,去掉去掉, 剩下的是 9 / (3^2) = 1.

8: 里面只有3个2,去掉去掉,剩下的是 8 / (2^3) = 1

于是啊,分母我们把剩下的部分乘起来,得到了神马?得到了 和 2,3 因子完全无关的 部分mod 36的值!就是 5 * 1 * 1 = 5了。

接下来,还有分母呢

1: 逆元(其实你可以无视它)

2: 一个2,去掉去掉, 剩下1, 逆元继续是1(继续无视)

3: 一个3,同上

接下来发现,2有几个? 分子有4个,分母1个,所以一共只有4 – 1 = 3个

3有几个? 同上的做法,显然只有1个。

于是呢答案就是:

5 * 1 * 2^3 * 3^1 = 12( mod 36)

解释:

5 -> 分子除了因子2,3的积

1 ->分母除了因子2,3 的逆元的积

2^3 -> 最终统计发现有3个2

3^1 ->最终统计发现有1个3

请好好理解本例子,你会发现这个问题是如此的美妙!

经典例题:

http://acm.fzu.edu.cn/problem.php?pid=2020 

c)     m<= n 别想了!我不会!你会了教我!难度系数: -1

二.           基础篇:我会数论!

1)     n,m<= 10^6, K是10^9级别

对于n! 分解素因子,这里就不说了,可以参考各种帖子。

之后保存个数,二分求幂啊啊啊啊啊

2)     n,m<= 10^10, k是素数,并且K 很小(比如几百?)

其实遇到这种情况我都用一个叫Lucas定理的东西。

 

ni,mi 就是把 n,m分解p进制的第i位的值。

例如:

                     计算 C(12, 4) mod 7

                     n = 12  (15)(base_7)

                     m = 4     (4) (base_7)

                     为了对齐,我们前面的部分补0

                     m = 4   (04) (base_7)

                     于是

                     Ans = C(5,4) * C(1,0) mod 7= 5 (mod 7)

                     有人又要问了, 如果mi > ni 怎么办呀?

                     直接为0!!!!!!!!!

              这里不给出证明,证明可以搜索到。同时由于这个应用的区域比较狭窄,显然有更简单,更好理解的算法,于是这里被无视了。

三.           究极篇

n,m <= 10^9, p <= 10^5

是不是怎么看怎么不可做呢?

第一次见到这种题目是不是觉得作者NC了,出个不可做题 >_<

第一次交发现一坨人全部WA,是不是觉得作者的数据搞疵?!!!!

首先要知道,这题其实等价是求:

 

求完直接合并一个模方程即可。(CRT)

p^c 的规模大约是10^5。

c 不是1,lucas阻止不了它。

n,m太大,因子分解也阻止不了它。

下面介绍我的做法:

假设 p = 3, c = 2,也就是mod 9

假设n = 19

n! = 1 * 2 * 3 * 4 * 5 * 6 * 7 * 8 * …… * 19

要是可以快速得到 n! 中除掉3 以后 mod 9的结果,那么多好呀!

看3多讨厌,直接砍

type cal( int n) :

n! = [ 1 * 2 * 4 * 5 * 7 * 8 * …  * 16 * 17 * 19 ] * (3 * 6 * 9 * 12 * 15 * 18)= [ 1 * 2 * 4 * 5 * 7 * 8 * …  * 16 * 17* 19 ] * 3^6( 1 * 2 * 3 * 4 * 5 * 6)

然后发现后面的一坨实际上是 cal( n / p) !!!!

再看前半部分,尼玛是以 p^c 为周期的啊!!!

[1 * 2 * 4 * 5 * 7 * 8 ] = [10 * 11 * 13 * 14 * 16 *17 ] = (mod 9)

于是说白了,对于前面的部分,由于周期,都是浮云了

下面是 孤立出来的19

可以知道孤立出来的 长度 不超过 p^c ,于是暴力啊,暴力啊!

于是完美解决n! 中和 p无关的项 mod p^c的值!!!

接下来是分母部分,一模一样,无非多了一个求逆元(因为都和p没关系了,逆元必然存在)

我们来分析一下,这样的复杂度是如何的呢

每次递归,规模变为原来的 1/p

logp N的啊!!!

当然是层数= =

于是问题完美解决!


版权声明:本文为博主原创文章,未经博主允许不得转载。

换BLOG了~http://hi.baidu.com/aekdycoin/blog

http://hi.baidu.com/aekdycoin/blog
  • tao_tao_bu_jue
  • tao_tao_bu_jue
  • 2009年02月10日 21:15
  • 1090

【AC大牛陈鸿的ACM总结贴】【ID AekdyCoin】人家当初也一样是菜鸟

acm总结帖_By AekdyCoin       各路大牛都在中国大陆的5个赛区结束以后纷纷发出了退役帖,总结帖,或功德圆满,或死不瞑目,而这或许又会造就明年的各种“炸尸”风波。为了考虑在发退役贴以...
  • lsgqjh
  • lsgqjh
  • 2015年03月24日 18:39
  • 1212

组合数取模

Lucas定理: m,n是非负整数,p是质数。m,n分别写成p进制: 则组合数C(n,m)与C(n[k],m[k])*C(n[k-1],m[k-1])*...*C(n[0],m[0]) mod p同...
  • ACdreamers
  • ACdreamers
  • 2012年10月03日 12:41
  • 23410

组合数取模终极版

以前讲述过很多组合数取模问题,详见:http://blog.csdn.net/acdreamers/article/details/8037918   今天,我们继续学习一些稍有难度的组合数取模问题,...
  • ACdreamers
  • ACdreamers
  • 2015年02月21日 20:08
  • 5961

FZU 2020 组合 -- (大组合数取模 Lucas定理)

Problem 2020 组合 Accept: 892    Submit: 2120 Time Limit: 1000 mSec    Memory Limit : 32768 KB  Pro...
  • qq_27717967
  • qq_27717967
  • 2016年05月24日 21:49
  • 498

利用扩展欧几里得算法快速计算组合数取模

  • 2011年08月17日 12:44
  • 1KB
  • 下载

组合数取模之逆元方法+模板

参自: http://www.cnblogs.com/liziran/p/6804803.html https://baike.baidu.com/item/%E8%B4%B9%E9%A9%A...
  • xiaoli_nu
  • xiaoli_nu
  • 2017年10月08日 22:15
  • 210

组合数取模常用函数与方法总结

组合数取模常用函数与方法总结 这一类题目一般是先手算出答案的组合表达式,然后套用组合数取模的一些模板实现 一般是两种思路,一种是硬算(牺牲时间),一种打出阶乘表(牺牲空间) 有时候需要的话还可以打出...
  • MartaYang
  • MartaYang
  • 2017年01月26日 12:00
  • 230

【快速幂模和组合数取模】

【快速幂模和组合数取模】 记得在我第一天学习编程的时候老师就和我们说:为什么要学计算机呢?因为它速度快啊。一句玩笑话,不过确实快啊。这几天水题的时候就遇到了这样一个问题a^b%c,想了好几个方法都一...
  • mengzhongsharen
  • mengzhongsharen
  • 2017年03月13日 22:56
  • 281

组合数取模

学习总结
  • AaronGZK
  • AaronGZK
  • 2016年02月12日 11:11
  • 1800
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【转载】【转自AekdyCoin的组合数取模】
举报原因:
原因补充:

(最多只允许输入30个字)