Codeforces 57C Array

本文介绍了一个使用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;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值