题解:AT_abc156_d [ABC156D] Bouquet

此题是一道排列组合题。显然可得答案是 2 n − C n a − C n b 2^n-C^a_n-C^b_n 2nCnaCnb

其中 2 n 2^n 2n 可用快速幂得出,两个组合数可以直接套公式暴力乘即可。

由于求组合数的时候最后要除一下,需要使用逆元,由于 m o d mod mod 是一个质数, a × b m o d − 2 m o d    1 0 9 + 7 a \times b^{mod-2} \mod 10^9+7 a×bmod2mod109+7 就可以求出。

#include <bits/stdc++.h>
using namespace std;
#define int long long
typedef long long ll;
typedef pair<int, int> pr;
#define up(i, l, r) for(int i = (l); i <= (r); i++)
#define down(i, r, l) for(int i = (r); i >= (l); i--)
const int mod = 1000000007;
const int base = 2333;
const double eps = 1e-6;

inline int read() {
	int x = 0, f = 1;
	char ch = getchar();
	while (ch < '0' || ch > '9') {
		if (ch == '-') f = -1;
		ch = getchar();
	}
	while (ch >= '0' && ch <= '9') {
		x = (x << 1) + (x << 3) + (ch ^ 48);
		ch = getchar();
	}
	return x * f;
}

int n, m, k, Q, T, _, ans = 0;
int a, b;

inline int power(int x, int y) {
	int res = 1;
	while (y) {
		if (y & 1) res = (res * x) % mod;
		x = (x * x) % mod;
		y >>= 1;
	}
	return res;
}

inline int C(int x, int y) {
	int res1 = 1, res2 = 1;
	for (int i = 1; i <= x; i++)
		res1 = (res1 * i) % mod;
	for (int i = y - x + 1; i <= y; i++)
		res2 = (res2 * i) % mod;
	return res2 * power(res1, mod - 2) % mod;
}

signed main() {
	n = read(), a = read(), b = read();
	cout << (power(2, n) - 1 - C(a, n) - C(b, n) + mod * 2) % mod;
	return 0;
}
  • 25
    点赞
  • 16
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值