首页 > 其他 > 详细

求组合数小结

时间:2019-02-10 20:15:35      阅读:176      评论:0      收藏:0      [点我收藏+]

今天学了一天数学,觉得自己都要转竞了23333

题目链接https://vjudge.net/contest/282927#problem/E

这里说一说求组合数的方法吧

其实就是求阶乘及其逆元的方法:

规定mod为模数,n为数据规模

1.mod为素数

   费马小定理:nlogn

   线性求逆元(n较小)

2.mod不为素数

  欧拉筛出两个阶乘的素数,再计算出每个素数的次幂,最后快速幂乘起来即可

  详见

	n = read(),p = read();
	ans = 1;
	for(int i = 2;i <= 2 * n;i++)
	{
		if(isprime[i] == 0) prime[++tot] = i,minm[i] = i;
		for(int j = 1;j <= tot && prime[j] * i <= 2 * n;j++)
			{
				isprime[i * prime[j]] = 1;
				minm[i * prime[j]] = prime[j];
				if(i % prime[j] == 0) break;
			}
	}
	for(int i = n + 2;i <= 2 * n;i++) cnt[i] = 1;
	for(int i = 1;i <= n;i++) cnt[i] = -1;
	for(int i = 2 * n;i > 1;i--)
	{
		if(isprime[i])
		{
			cnt[minm[i]] += cnt[i];
			cnt[i / minm[i]] += cnt[i];
		}
	}
	for(int i = 2;i <= 2 * n;i++)
		if(isprime[i] == 0)
			ans = ans * power(i,cnt[i]) % p;
	printf("%lld",ans);

  数学使人憔悴啊qaq

求组合数小结

原文:https://www.cnblogs.com/LM-LBG/p/10360041.html

(0)
(0)
   
举报
评论 一句话评论(0
关于我们 - 联系我们 - 留言反馈 - 联系我们:wmxa8@hotmail.com
© 2014 bubuko.com 版权所有
打开技术之扣,分享程序人生!