結果
問題 | No.895 MESE |
ユーザー |
![]() |
提出日時 | 2019-09-27 21:46:27 |
言語 | C++17(clang) (17.0.6 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 100 ms / 2,000 ms |
コード長 | 2,076 bytes |
コンパイル時間 | 1,054 ms |
コンパイル使用メモリ | 148,124 KB |
実行使用メモリ | 11,264 KB |
最終ジャッジ日時 | 2024-11-30 14:58:09 |
合計ジャッジ時間 | 4,483 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 26 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS #include <stdio.h> #include <algorithm> #include <utility> #include <functional> #include <cstring> #include <queue> #include <stack> #include <math.h> #include <iterator> #include <vector> #include <string> #include <set> #include <math.h> #include <iostream> #include <random> #include<map> #include <iomanip> #include <time.h> #include <stdlib.h> #include <list> #include <typeinfo> #include <list> #include <set> #include <cassert> #include<fstream> #include <unordered_map> #include <cstdlib> #include <complex> #include <cctype> #include <bitset> using namespace std; typedef string::const_iterator State; #define Ma_PI 3.141592653589793 #define eps 1e-5 #define LONG_INF 1000000000000000000LL #define GOLD 1.61803398874989484820458 #define MAX_MOD 1000000007LL #define GYAKU 500000004LL #define MOD 998244353LL #define seg_size 262144*4 #define REP(a,b) for(long long a = 0;a < b;++a) long long powing(long long now, long long b) { long long ans = 1; while (b != 0) { if (b % 2) { ans *= now; ans %= MAX_MOD; } b /= 2; now *= now; now %= MAX_MOD; } return ans; } long long inv(long long now) { return powing(now, MAX_MOD - 2LL); } long long mae[500000]; long long gyaku[500000]; long long comb(long long a, long long b) { long long ans = mae[b]; ans *= gyaku[a]; ans %= MAX_MOD; ans *= gyaku[b - a]; ans %= MAX_MOD; return ans; } int main() { mae[0] = 1; gyaku[0] = 1; for (long long i = 1; i < 500000; ++i) { mae[i] = mae[i - 1] * i; mae[i] %= MAX_MOD; gyaku[i] = inv(mae[i]); } long long a, b, c; cin >> a >> b >> c; a--; b--; long long ans = 0; while (a != -1) { //Choose B and starting! //期待値ゲー long long percent = c * inv(a + b + c); // cを選ぶ確率 long long toori = (comb(a, a + b + c) * comb(b, b + c)) % MAX_MOD;//通り数 long long tmp = (powing(2LL, a + b + c) + MAX_MOD - 1LL) % MAX_MOD; percent %= MAX_MOD; ans += (((percent * toori) % MAX_MOD) * tmp) % MAX_MOD; ans %= MAX_MOD; //choose A and continue a--; } cout << ans << endl; return 0; }