結果

問題 No.1126 SUM
ユーザー ooaiu
提出日時 2025-02-12 19:31:07
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 5 ms / 1,000 ms
コード長 889 bytes
コンパイル時間 3,517 ms
コンパイル使用メモリ 275,860 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2025-02-12 19:31:13
合計ジャッジ時間 5,043 ms
ジャッジサーバーID
(参考情報)
judge3 / judge2
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
other AC * 26
権限があれば一括ダウンロードができます

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
#ifdef LOCAL
#include "algo/debug.h"
#else
#define debug(...) (void(0))
#endif
void run_case();
int main() {
	std::ios_base::sync_with_stdio(false);
	std::cin.tie(nullptr);
	int T = 1;
	// cin >> T;
	while (T--) run_case();
	return 0;
}
ll mpow(ll a, ll b, ll m) {
	ll r = 1;
	while(b) {
		if(b & 1) {
			r *= a;
			r %= m;
		}
		a *= a;
		a %= m;
		b >>= 1;
	}
	return r;
}
void run_case() {
	int N, M;
	cin >> N >> M;
	const int MOD = 1e9 + 7;
	const int K = 1 << 17;
	vector<ll> fact(K, 1), inv(K, 1);
	for(int i = 2; i < K; i++) {
		fact[i] = fact[i - 1] * i;
		fact[i] %= MOD;
	}
	inv[K - 1] = mpow(fact[K - 1], MOD - 2, MOD);
	for(int i = K - 2; i >= 2; i--) {
		inv[i] = inv[i + 1] * (i + 1);
		inv[i] %= MOD;
	}
	ll a = fact[M + 1], b = inv[N + 1], c = inv[M - N];
	cout << (((a * b) % MOD) * c) % MOD << endl;
}
0