結果
問題 |
No.1126 SUM
|
ユーザー |
|
提出日時 | 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 |
ソースコード
#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; }