結果
問題 | No.2207 pCr検査 |
ユーザー | dachengz |
提出日時 | 2023-02-06 12:46:15 |
言語 | PyPy3 (7.3.15) |
結果 |
RE
|
実行時間 | - |
コード長 | 1,254 bytes |
コンパイル時間 | 1,026 ms |
コンパイル使用メモリ | 86,724 KB |
実行使用メモリ | 176,708 KB |
最終ジャッジ日時 | 2023-09-18 00:34:45 |
合計ジャッジ時間 | 17,353 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge14 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | RE | - |
testcase_01 | AC | 70 ms
71,476 KB |
testcase_02 | RE | - |
testcase_03 | RE | - |
testcase_04 | RE | - |
testcase_05 | RE | - |
testcase_06 | RE | - |
testcase_07 | RE | - |
testcase_08 | RE | - |
testcase_09 | RE | - |
testcase_10 | RE | - |
testcase_11 | RE | - |
testcase_12 | AC | 432 ms
123,020 KB |
testcase_13 | AC | 310 ms
106,192 KB |
testcase_14 | AC | 155 ms
84,572 KB |
testcase_15 | AC | 337 ms
108,356 KB |
testcase_16 | AC | 299 ms
103,248 KB |
testcase_17 | AC | 327 ms
107,744 KB |
testcase_18 | RE | - |
testcase_19 | AC | 230 ms
92,360 KB |
testcase_20 | AC | 181 ms
86,436 KB |
testcase_21 | AC | 509 ms
135,200 KB |
testcase_22 | RE | - |
testcase_23 | RE | - |
testcase_24 | AC | 137 ms
79,644 KB |
testcase_25 | AC | 309 ms
100,020 KB |
testcase_26 | AC | 236 ms
89,820 KB |
testcase_27 | AC | 302 ms
96,932 KB |
testcase_28 | RE | - |
testcase_29 | RE | - |
testcase_30 | RE | - |
testcase_31 | RE | - |
ソースコード
from math import log, lgamma def deg(n, p): ret = 0 n //= p while n: ret += n n //= p return ret def isqrt(n): if n <= 0: return 0 x = int((n ** 0.5) * (1 + 1e-14)) while True: y = (x + n // x) // 2 if y >= x: return x x = y k = int(input()) pe = [] for _ in range(k): pe.append([int(x) for x in input().split()]) P = pe[-1][0] P1 = P + 1 needlog = 0.0 for p, e in pe: if e > deg(P, p): print(-1, -1) raise ValueError needlog += log(p) * e lgP = lgamma(P1) lo1, hi1 = 1, P // 2 if lgP - lgamma(hi1) - lgamma(P1 - hi1) - needlog < -1e-6: print(-1, -1) raise ValueError while hi1 - lo1 > 1: md = lo1 + (hi1 - lo1) // 2 lgcomb = lgP - lgamma(md + 1) - lgamma(P1 - md) if lgcomb - needlog < -1e-6: lo1 = md else: hi1 = md lo2, hi2 = 1, P // 2 while hi2 - lo2 > 1: md = lo2 + (hi2 - lo2) // 2 lgcomb = lgP - lgamma(md + 1) - lgamma(P1 - md) if lgcomb - needlog > 1e-6: hi2 = md else: lo2 = md for r in range(lo1, hi2 + 1): if all(deg(P, p) - deg(r, p) - deg(P - r, p) == e for p, e in pe): print(P, r) raise ValueError else: print(-1, -1)