結果
問題 | No.766 金魚すくい |
ユーザー | nebukuro09 |
提出日時 | 2018-12-14 10:09:43 |
言語 | D (dmd 2.106.1) |
結果 |
AC
|
実行時間 | 698 ms / 1,500 ms |
コード長 | 1,828 bytes |
コンパイル時間 | 1,035 ms |
コンパイル使用メモリ | 121,856 KB |
実行使用メモリ | 10,284 KB |
最終ジャッジ日時 | 2024-06-13 02:20:50 |
合計ジャッジ時間 | 19,133 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 164 ms
6,784 KB |
testcase_01 | AC | 164 ms
6,912 KB |
testcase_02 | AC | 164 ms
6,912 KB |
testcase_03 | AC | 165 ms
6,784 KB |
testcase_04 | AC | 169 ms
6,784 KB |
testcase_05 | AC | 174 ms
6,784 KB |
testcase_06 | AC | 168 ms
6,784 KB |
testcase_07 | AC | 168 ms
6,784 KB |
testcase_08 | AC | 166 ms
6,784 KB |
testcase_09 | AC | 166 ms
6,912 KB |
testcase_10 | AC | 168 ms
6,784 KB |
testcase_11 | AC | 165 ms
6,784 KB |
testcase_12 | AC | 165 ms
6,784 KB |
testcase_13 | AC | 175 ms
6,912 KB |
testcase_14 | AC | 179 ms
6,912 KB |
testcase_15 | AC | 176 ms
6,912 KB |
testcase_16 | AC | 172 ms
6,912 KB |
testcase_17 | AC | 170 ms
6,912 KB |
testcase_18 | AC | 168 ms
6,912 KB |
testcase_19 | AC | 180 ms
6,912 KB |
testcase_20 | AC | 176 ms
7,040 KB |
testcase_21 | AC | 178 ms
6,912 KB |
testcase_22 | AC | 177 ms
7,040 KB |
testcase_23 | AC | 640 ms
10,108 KB |
testcase_24 | AC | 657 ms
10,096 KB |
testcase_25 | AC | 662 ms
10,144 KB |
testcase_26 | AC | 632 ms
10,056 KB |
testcase_27 | AC | 653 ms
10,148 KB |
testcase_28 | AC | 633 ms
9,948 KB |
testcase_29 | AC | 628 ms
9,952 KB |
testcase_30 | AC | 679 ms
9,944 KB |
testcase_31 | AC | 698 ms
10,284 KB |
testcase_32 | AC | 639 ms
10,084 KB |
testcase_33 | AC | 364 ms
10,104 KB |
testcase_34 | AC | 364 ms
10,240 KB |
testcase_35 | AC | 347 ms
6,912 KB |
testcase_36 | AC | 333 ms
6,784 KB |
testcase_37 | AC | 647 ms
10,012 KB |
testcase_38 | AC | 655 ms
10,100 KB |
testcase_39 | AC | 638 ms
10,136 KB |
testcase_40 | AC | 638 ms
10,068 KB |
testcase_41 | AC | 639 ms
9,476 KB |
testcase_42 | AC | 623 ms
9,336 KB |
testcase_43 | AC | 166 ms
6,912 KB |
testcase_44 | AC | 165 ms
6,784 KB |
ソースコード
import std.stdio, std.array, std.string, std.conv, std.algorithm; import std.typecons, std.range, std.random, std.math, std.container; import std.numeric, std.bigint, core.bitop, core.stdc.string; immutable long MOD = 10^^9 + 7; immutable long MAX = 3 * 10^^5; long[] F, G; void main() { auto s = readln.split.map!(to!int); auto N = s[0]; auto M = s[1]; auto P = s[2]; auto V = readln.split.map!(to!long).array; V.sort!"a > b"; F = new long[](MAX+1); G = new long[](MAX+1); F[0] = 1; foreach (i; 1..MAX+1) F[i] = F[i-1] * i % MOD; foreach (i; 0..MAX+1) G[i] = powmod(F[i], MOD-2, MOD); long ans = 0; long v = 0; foreach (i; 0..N) { // i回成功, M回失敗する long ok = i; long ng = M; long tmp = F[ok+ng-1] * G[ok] % MOD * G[ng-1] % MOD; long p1 = powmod(100-P, ok, MOD) * powmod(powmod(100, MOD-2, MOD), ok, MOD) % MOD; long p2 = powmod(P, ng, MOD) * powmod(powmod(100, MOD-2, MOD), ng, MOD) % MOD; ans += v * p1 % MOD * p2 % MOD * tmp % MOD; ans %= MOD; v = (v + V[i]) % MOD; } foreach (i; 0..M) { // N回成功, i回失敗する long ok = N; long ng = i; long tmp = F[ok+ng-1] * G[ok-1] % MOD * G[ng] % MOD; long p1 = powmod(100-P, ok, MOD) * powmod(powmod(100, MOD-2, MOD), ok, MOD) % MOD; long p2 = powmod(P, ng, MOD) * powmod(powmod(100, MOD-2, MOD), ng, MOD) % MOD; ans += v * p1 % MOD * p2 % MOD * tmp % MOD; ans %= MOD; } ans.writeln; } long comb(long n, long r) { if (n < r) return 0; return F[n] * G[n-r] % MOD * G[r] % MOD; } long powmod(long a, long x, long m) { long ret = 1; while (x) { if (x % 2) ret = ret * a % m; a = a * a % m; x /= 2; } return ret; }