結果
問題 | No.766 金魚すくい |
ユーザー | nebukuro09 |
提出日時 | 2018-12-14 10:09:43 |
言語 | D (dmd 2.105.2) |
結果 |
AC
|
実行時間 | 699 ms / 1,500 ms |
コード長 | 1,828 bytes |
コンパイル時間 | 2,644 ms |
コンパイル使用メモリ | 108,352 KB |
実行使用メモリ | 12,744 KB |
最終ジャッジ日時 | 2023-09-03 21:41:15 |
合計ジャッジ時間 | 21,494 ms |
ジャッジサーバーID (参考情報) |
judge15 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 176 ms
8,528 KB |
testcase_01 | AC | 177 ms
9,028 KB |
testcase_02 | AC | 176 ms
9,740 KB |
testcase_03 | AC | 176 ms
8,508 KB |
testcase_04 | AC | 177 ms
9,000 KB |
testcase_05 | AC | 176 ms
7,984 KB |
testcase_06 | AC | 176 ms
8,240 KB |
testcase_07 | AC | 176 ms
8,716 KB |
testcase_08 | AC | 175 ms
8,728 KB |
testcase_09 | AC | 176 ms
8,200 KB |
testcase_10 | AC | 177 ms
9,044 KB |
testcase_11 | AC | 176 ms
8,512 KB |
testcase_12 | AC | 176 ms
8,728 KB |
testcase_13 | AC | 182 ms
8,464 KB |
testcase_14 | AC | 184 ms
8,244 KB |
testcase_15 | AC | 183 ms
8,484 KB |
testcase_16 | AC | 183 ms
9,016 KB |
testcase_17 | AC | 184 ms
9,272 KB |
testcase_18 | AC | 180 ms
8,500 KB |
testcase_19 | AC | 184 ms
8,508 KB |
testcase_20 | AC | 183 ms
8,240 KB |
testcase_21 | AC | 183 ms
8,508 KB |
testcase_22 | AC | 187 ms
8,976 KB |
testcase_23 | AC | 665 ms
11,256 KB |
testcase_24 | AC | 682 ms
11,604 KB |
testcase_25 | AC | 699 ms
11,052 KB |
testcase_26 | AC | 672 ms
11,440 KB |
testcase_27 | AC | 696 ms
11,748 KB |
testcase_28 | AC | 671 ms
11,252 KB |
testcase_29 | AC | 667 ms
12,000 KB |
testcase_30 | AC | 673 ms
11,200 KB |
testcase_31 | AC | 693 ms
11,768 KB |
testcase_32 | AC | 673 ms
11,944 KB |
testcase_33 | AC | 384 ms
11,888 KB |
testcase_34 | AC | 386 ms
11,008 KB |
testcase_35 | AC | 363 ms
8,728 KB |
testcase_36 | AC | 352 ms
9,776 KB |
testcase_37 | AC | 682 ms
12,444 KB |
testcase_38 | AC | 687 ms
12,744 KB |
testcase_39 | AC | 677 ms
11,244 KB |
testcase_40 | AC | 672 ms
11,628 KB |
testcase_41 | AC | 676 ms
10,652 KB |
testcase_42 | AC | 663 ms
10,788 KB |
testcase_43 | AC | 177 ms
7,960 KB |
testcase_44 | AC | 176 ms
8,524 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; }