結果
問題 | 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 200000long 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;}