結果
問題 | No.1879 How many matchings? |
ユーザー | Eki1009 |
提出日時 | 2022-03-18 22:44:29 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 43 ms / 2,000 ms |
コード長 | 527 bytes |
コンパイル時間 | 375 ms |
コンパイル使用メモリ | 82,476 KB |
実行使用メモリ | 60,816 KB |
最終ジャッジ日時 | 2024-04-14 10:47:13 |
合計ジャッジ時間 | 1,844 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 36 ms
53,220 KB |
testcase_01 | AC | 37 ms
52,408 KB |
testcase_02 | AC | 36 ms
52,480 KB |
testcase_03 | AC | 37 ms
52,464 KB |
testcase_04 | AC | 37 ms
52,316 KB |
testcase_05 | AC | 36 ms
53,220 KB |
testcase_06 | AC | 37 ms
52,388 KB |
testcase_07 | AC | 37 ms
53,488 KB |
testcase_08 | AC | 36 ms
52,756 KB |
testcase_09 | AC | 37 ms
53,268 KB |
testcase_10 | AC | 43 ms
60,816 KB |
testcase_11 | AC | 43 ms
59,820 KB |
testcase_12 | AC | 41 ms
60,200 KB |
testcase_13 | AC | 37 ms
53,700 KB |
testcase_14 | AC | 43 ms
60,232 KB |
ソースコード
mod = 10**9+7 def mul(A,B): C = [0]*(len(A)+len(B)-1) for i,a in enumerate(A): for j,b in enumerate(B): C[i+j] = (C[i+j]+a*b)%mod return C def bostan_mori(P,Q,n): if n < 0: return 0 while n: X = Q.copy() for i in range(1,len(X),2): X[i] *= -1 P = mul(P,X)[n&1::2] Q = mul(Q,X)[::2] n >>= 1 return P[0] n = int(input()) P = [1] Q = [1,-1,-1] if n%2 == 0: ans = bostan_mori(P,Q,n//2) else: Q = mul(Q,Q) ans = (bostan_mori(P,Q,n//2)+bostan_mori(P,Q,n//2-1))%mod print(ans)