Codeforces 57C Array

原创于 2019-04-04 12:53:56 发布 · 205 阅读

0 ·

CC 4.0 BY-SA版权

ACM 专栏收录该内容

245 篇文章

订阅专栏

本文介绍了一个使用Lucas定理计算特定形式组合数模特定值的应用案例。通过定义阶乘、扩展欧几里得算法求逆元等辅助函数，实现了高效计算C(2n-1,n) mod p的功能。

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

#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <cstdio>
#include <math.h>
#include <string>
#include <algorithm>
#include <vector>
#define ll long long

using namespace std;
const int mod = 1000000007;

ll fact(ll n, ll p) {//n的阶乘求余p 
	ll ret = 1;
	for (ll i = 1; i <= n; i++) ret = ret * i % p;
	return ret;
}
void ex_gcd(ll a, ll b, ll &x, ll &y, ll &d) {
	if (!b) { d = a, x = 1, y = 0; }
	else {
		ex_gcd(b, a % b, y, x, d);
		y -= x * (a / b);
	}
}
ll inv(ll t, ll p) {//如果不存在，返回-1 //求逆元 
	ll d, x, y;
	ex_gcd(t, p, x, y, d);
	return d == 1 ? (x % p + p) % p : -1;
}
ll comb(ll n, ll m, ll p) {//C(n, m) % p
	if (m < 0 || m > n) return 0;
	return fact(n, p) * inv(fact(m, p), p) % p * inv(fact(n - m, p), p) % p;
}

ll Lucas(ll n, ll m, ll p) {//Lucas定理求C(n, m) % p 
	return m ? Lucas(n / p, m / p, p) * comb(n%p, m%p, p) % p : 1;
}

int n;

int main() {
	scanf("%d", &n);
	printf("%lld\n", (2 * Lucas(2 * n - 1, n, mod) % mod - n));
	//system("pause");
	return 0;
}