首页 > 其他分享 >OI loves Math(三)——组合数

OI loves Math(三)——组合数

时间:2022-11-10 22:22:59浏览次数:62  
标签:exp OI inv long base loves ans Math mod

又来了(喜)
这次我们聊组合数,也就是 $ C_n^m $ 。。。。。。

何为组合数

没人不知道吧。。。。。。
组合数就是问从 $ n $ 项里选 $ m $ 项有多少种选法,记作 $ C_n^m $ 或 $ \left( \begin{matrix} n \\ m \end{matrix} \right) $ 。
但是,竞赛中一般求 $ C_n^m \bmod 998244353 $ 。

怎么求

注:下文 $ mod $ 为模数。
有很多,但这儿介绍一种 $ O(n \log mod) $ 预处理,$ O(1) $ 单词查询的方法!
其实就是组合数公式。。。。。。

\[C_n^m = \frac{n!}{m!(n - m)!} \]

再加上以前的逆元:

\[\frac{1}{x} \equiv x^{mod - 2} \pmod{mod} \]

预处理 $ x! $ 和 $ x!^{mod - 2} \bmod mod $ 就可以了。
两段小代码:

prod[0] = 1, inv[0] = 1;
for (long long i = 1; i <= maxn; i++) {
	prod[i] = prod[i - 1] * i % mod;
	inv[i] = qpow(prod[i], mod - 2, mod);
}
long long C(long long n, long long m) {
	return prod[n] * inv[n - m] % mod * inv[m] % mod;
}

例题

CF1717D

把选手按编号顺序排列:(注:转载于OIer某罗

我们可以发现,选手 $ x $ 进行 $ x - 1 $ ,然后 $ \operatorname{popcount}(x - 1) $ 就是胜的场数。
那么问题就变成了求满足 $ \operatorname{popcount}(x) \le k \ \ \ (0 \le x \le 2^n - 1) $ 的 $ x $ 的个数,也就是 $ \sum _ {i = 0} ^ {\min \{ n, k \}} C _ {n} ^ {i} $
代码:

#include <bits/stdc++.h>
using namespace std;

long long prod[100005], inv[100005];
const long long mod = 1000000007;

long long qmul(long long base, long long exp, long long mod) {
	long long ans = 0;
	while (exp) {
		if (exp & 1) {
			ans = (ans + base) % mod;
		}
		exp >>= 1;
		base = (base + base) % mod;
	}
	return ans;
}

long long qpow(long long base, long long exp, long long mod) {
	long long ans = 1;
	while (exp) {
		if (exp & 1) {
//			ans = qmul(ans, base, mod);
			ans = ans * base % mod;
		}
		exp >>= 1;
//		base = qmul(base, base, mod);
		base = base * base % mod;
	}
	return ans;
}

long long C(long long n, long long m) {
	return prod[n] * inv[n - m] % mod * inv[m] % mod;
}

int main() {
	long long n, k;
	scanf("%lld %lld", &n, &k);
	long long ans = 0;
	prod[0] = 1, inv[0] = 1;
	k = min(n, k);
	for (long long i = 1; i <= n; i++) {
		prod[i] = prod[i - 1] * i % mod;
		inv[i] = qpow(prod[i], mod - 2, mod);
	}
	for (long long i = 0; i <= k; i++) {
		ans = (ans + C(n, i)) % mod;
	}
	printf("%lld", ans);
	return 0;
}

再见!

标签:exp,OI,inv,long,base,loves,ans,Math,mod
From: https://www.cnblogs.com/AProblemSolver/p/16878977.html

相关文章

  • Android开发Compose版本、Kotlin 版本、KSP版本版本对应关系
    Android开发Compose版本、Kotlin版本、KSP版本版本对应关系是要遵循官方给出的,不然容易出锅甚至编译都不过,即使编译通过也可能导致潜在崩溃ComposeCompiler版本和兼......
  • 探究Android中的注解
    本文系GDGAndroidMeetup分享内容总结文章注解是我们经常接触的技术,Java有注解,Android也有注解,本文将试图介绍Android中的注解,以及ButterKnife和Otto这些基于注解的库......
  • 一些快速提高Android开发的脚本与技巧(终端篇)
    正所谓“工欲善其事必先利其器”,一个好的工具或者技巧能让提升工作效率,起到事半功倍的效果。在这里斗胆列出一些窃以为一些可能快速提高Android日常开发的脚本,希望可以为大......
  • P7737 [NOI2021] 庆典
    题意给定一张有向图,每次询问给出\(s,t\),求从\(s\tot\)的路径上(可以有重复点)可能会经过多少个点,每次询问会临时加入\(k\)条边。其中,题目给出的图有如下特点:若\(x\)......
  • Android基于坐标对View进行模拟点击事件
    在Android中,我们对于View进行模拟点击事件,很容易,比如调用​​View.performClick​​即可。但是有些时候,我们想要更加精细的点击,比如View的某一区域或者某一点进行点击。比如......
  • 一个简单实用的Android调试应用技巧
    在应用开发中,我们常常会进行日志打印或者debug调试,以此来分析运行时的一些信息,便于发现bug和问题。AndroidStudio的Debug功能很好用,但是有时候有些情况下,就显得不是那么快......
  • 关于 Android 应用多进程的整理
    在计算机操作系统中,进程是进行资源分配和调度的基本单位。这对于基于Linux内核的Android系统也不例外。在Android的设计中,一个应用默认有一个(主)进程。但是我们通过配置可......
  • Android中使用ViewStub提高布局性能
    在Android开发中,View是我们必须要接触的用来展示的技术.通常情况下随着View视图的越来越复杂,整体布局的性能也会随之下降.这里介绍一个在某些场景下提升布局性能的View,它......
  • 小姿势 之 Android Studio 3.5 Retry 问题解决
    总会有那么一个人,让你觉得这个世界一切都是值得的。纵使结果不尽人意,曾经拥有即是最好。前言家里的MBP静静地躺了一段时间,某天心血来潮想嗨起来,其实就是瞎折腾一把,结果......
  • Android代码规范利器: Checkstyle
    程序代码向来都不仅仅是用来运行的,写的一手好代码,易读,可维护应该是每个程序员所追求的。每个团队都(应该)有一套优良统一的代码规范,而规范的检查依赖于人工检测就不太现实,好在......