結果
問題 | No.766 金魚すくい |
ユーザー |
![]() |
提出日時 | 2018-12-14 01:54:13 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 31 ms / 1,500 ms |
コード長 | 1,989 bytes |
コンパイル時間 | 837 ms |
コンパイル使用メモリ | 98,628 KB |
実行使用メモリ | 8,544 KB |
最終ジャッジ日時 | 2024-09-25 04:54:39 |
合計ジャッジ時間 | 2,984 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 42 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS #include<iostream> #include<sstream> #include<cstdio> #include<cstdlib> #include<cstring> #include<climits> #include<cmath> #include<string> #include<vector> #include<set> #include<map> #include<queue> #include<numeric> #include<functional> #include<algorithm> #include<bitset> #include<tuple> #include<unordered_set> #include<unordered_map> #include<random> #include<array> #include<cassert> using namespace std; #define INF (1<<29) #define rep(i,n) for(int i=0;i<(int)(n);i++) #define all(v) v.begin(),v.end() #define uniq(v) v.erase(unique(all(v)),v.end()) #define MOD 1000000007 #define MAXN 200000 long long inv[MAXN + 1];//MODを法とする乗法の逆元 long long fact[MAXN + 1];//階乗 long long ifact[MAXN + 1];//階乗の逆元 void init() { inv[1] = 1; for (int i = 2; i <= MAXN; i++) inv[i] = inv[MOD%i] * (MOD - MOD / i) % MOD; fact[0] = ifact[0] = 1; for (int i = 1; i <= MAXN; i++) { fact[i] = i * fact[i - 1] % MOD; ifact[i] = ifact[i - 1] * inv[i] % MOD; } } long long C(int n, int r) { if (n < 0 || r < 0 || r > n)return 0; if (r > n / 2)r = n - r; return fact[n] * ifact[n - r] % MOD*ifact[r] % MOD; } int n, m, p; int v[100000]; int main() { ios::sync_with_stdio(0); cin.tie(0); init(); cin >> n >> m >> p; rep(i, n)cin >> v[i]; sort(v, v + n); reverse(v, v + n); rep(i, n - 1) (v[i + 1] += v[i]) %= MOD; long long ans = 0; long long suc = 1, fail = 1; rep(i, m) { fail *= p; fail %= MOD; fail *= inv[100]; fail %= MOD; } rep(i, n - 1) { suc *= 100 - p; suc %= MOD; suc *= inv[100]; suc %= MOD; ans += v[i] * C(i + 2 + m - 1 - 1, m - 1) % MOD * suc%MOD * fail%MOD; ans %= MOD; } suc = 1; fail = 1; rep(i, n) { suc *= 100 - p; suc %= MOD; suc *= inv[100]; suc %= MOD; } rep(i, m) { ans += v[n-1] * C(n + i - 1, i) % MOD * suc%MOD * fail%MOD; ans %= MOD; fail *= p; fail %= MOD; fail *= inv[100]; fail %= MOD; } cout << ans << endl; return 0; }