【BZOJ1856】[Scoi2010]字符串【组合数】【神题】

最新推荐文章于 2019-07-27 15:29:00 发布

原创最新推荐文章于 2019-07-27 15:29:00 发布 · 839 阅读

0 ·

CC 4.0 BY-SA版权

组合数学专栏收录该内容

11 篇文章

订阅专栏

本文提供了一种使用组合数学原理解决特定编程竞赛问题的方法。通过快速幂运算和逆元计算，该算法有效地计算了两个整数集合的所有可能组合，避免了直接使用阶乘带来的大数问题。此解决方案适用于那些需要高效计算组合数的场景。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

【题目链接】

太神了【将狼踩尽的题解】

/* Telekinetic Forest Guard */
#include <cstdio>
#include <algorithm>

using namespace std;

typedef long long LL;

const int maxn = 2000001, p = 20100403;

int n, m, fact[maxn];

inline int mul(int a, int b) {
	return (LL)a * b % p;
}

inline int qpow(int x, int n) {
	int s = 1;
	for(int t = x; n; n >>= 1, t = mul(t, t)) if(n & 1) s = mul(s, t);
	return s;
}

int main() {
	fact[0] = 1;
	for(int i = 1; i < maxn; i++) fact[i] = mul(fact[i - 1], i);
	scanf("%d%d", &n, &m);
	printf("%d\n", (mul(mul(fact[n + m], qpow(fact[n], p - 2)), qpow(fact[m], p - 2)) - mul(mul(fact[n + m], qpow(fact[n + 1], p - 2)), qpow(fact[m - 1], p - 2)) + p) % p);
	return 0;
}