結果
問題 | No.1035 Color Box |
ユーザー | maspy |
提出日時 | 2020-04-25 16:45:00 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 631 ms / 2,000 ms |
コード長 | 1,181 bytes |
コンパイル時間 | 100 ms |
コンパイル使用メモリ | 12,800 KB |
実行使用メモリ | 48,900 KB |
最終ジャッジ日時 | 2024-04-24 22:34:04 |
合計ジャッジ時間 | 24,109 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 619 ms
48,636 KB |
testcase_01 | AC | 548 ms
45,872 KB |
testcase_02 | AC | 545 ms
45,588 KB |
testcase_03 | AC | 547 ms
45,592 KB |
testcase_04 | AC | 550 ms
45,972 KB |
testcase_05 | AC | 546 ms
45,848 KB |
testcase_06 | AC | 546 ms
45,844 KB |
testcase_07 | AC | 594 ms
45,848 KB |
testcase_08 | AC | 631 ms
47,992 KB |
testcase_09 | AC | 554 ms
45,852 KB |
testcase_10 | AC | 541 ms
45,464 KB |
testcase_11 | AC | 567 ms
45,460 KB |
testcase_12 | AC | 573 ms
45,464 KB |
testcase_13 | AC | 581 ms
47,236 KB |
testcase_14 | AC | 589 ms
48,132 KB |
testcase_15 | AC | 631 ms
48,628 KB |
testcase_16 | AC | 577 ms
45,588 KB |
testcase_17 | AC | 580 ms
45,840 KB |
testcase_18 | AC | 552 ms
45,468 KB |
testcase_19 | AC | 591 ms
48,900 KB |
testcase_20 | AC | 542 ms
45,592 KB |
testcase_21 | AC | 556 ms
45,460 KB |
testcase_22 | AC | 551 ms
45,456 KB |
testcase_23 | AC | 549 ms
45,848 KB |
testcase_24 | AC | 547 ms
45,844 KB |
testcase_25 | AC | 552 ms
45,844 KB |
testcase_26 | AC | 547 ms
46,224 KB |
testcase_27 | AC | 553 ms
45,588 KB |
testcase_28 | AC | 558 ms
45,840 KB |
testcase_29 | AC | 544 ms
45,848 KB |
testcase_30 | AC | 556 ms
45,456 KB |
testcase_31 | AC | 547 ms
45,944 KB |
testcase_32 | AC | 540 ms
46,112 KB |
testcase_33 | AC | 546 ms
45,840 KB |
testcase_34 | AC | 545 ms
46,116 KB |
testcase_35 | AC | 569 ms
45,360 KB |
testcase_36 | AC | 567 ms
45,844 KB |
testcase_37 | AC | 571 ms
45,588 KB |
ソースコード
import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines import numpy as np MOD = 10 ** 9 + 7 def cumprod(A, MOD=MOD): L = len(A) Lsq = int(L**.5 + 1) A = np.resize(A, Lsq**2).reshape(Lsq, Lsq) for n in range(1, Lsq): A[:, n] *= A[:, n - 1] A[:, n] %= MOD for n in range(1, Lsq): A[n] *= A[n - 1, -1] A[n] %= MOD return A.ravel()[:L] def make_fact(U, MOD=MOD): x = np.arange(U, dtype=np.int64) x[0] = 1 fact = cumprod(x, MOD) x = np.arange(U, 0, -1, dtype=np.int64) x[0] = pow(int(fact[-1]), MOD - 2, MOD) fact_inv = cumprod(x, MOD)[::-1] fact.flags.writeable = False fact_inv.flags.writeable = False return fact, fact_inv fact, fact_inv = make_fact(10 ** 5 + 10) N, M = map(int, read().split()) def power(A, n): if n == 1: return A B = power(A, n // 2) B = B * B % MOD return A * B % MOD if n & 1 else B comb = fact[M] * fact_inv[: M + 1] % MOD * fact_inv[: M + 1][::-1] % MOD power = power(np.arange(M + 1, dtype=np.int64)[::-1], N) x = comb * power % MOD x[1::2] *= (-1) print(x.sum() % MOD)