結果
問題 | No.891 隣接3項間の漸化式 |
ユーザー | vwxyz |
提出日時 | 2024-04-10 08:07:06 |
言語 | Python3 (3.12.2 + numpy 1.26.4 + scipy 1.12.0) |
結果 |
AC
|
実行時間 | 33 ms / 2,000 ms |
コード長 | 2,292 bytes |
コンパイル時間 | 261 ms |
コンパイル使用メモリ | 13,056 KB |
実行使用メモリ | 11,136 KB |
最終ジャッジ日時 | 2024-10-02 09:50:57 |
合計ジャッジ時間 | 2,925 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 31 ms
10,880 KB |
testcase_01 | AC | 30 ms
11,008 KB |
testcase_02 | AC | 30 ms
11,008 KB |
testcase_03 | AC | 31 ms
11,136 KB |
testcase_04 | AC | 31 ms
10,880 KB |
testcase_05 | AC | 30 ms
10,880 KB |
testcase_06 | AC | 31 ms
10,880 KB |
testcase_07 | AC | 31 ms
10,880 KB |
testcase_08 | AC | 32 ms
10,880 KB |
testcase_09 | AC | 31 ms
11,008 KB |
testcase_10 | AC | 33 ms
11,008 KB |
testcase_11 | AC | 31 ms
10,880 KB |
testcase_12 | AC | 30 ms
10,880 KB |
testcase_13 | AC | 30 ms
11,136 KB |
testcase_14 | AC | 31 ms
10,880 KB |
testcase_15 | AC | 30 ms
10,880 KB |
testcase_16 | AC | 31 ms
11,008 KB |
testcase_17 | AC | 31 ms
11,136 KB |
testcase_18 | AC | 30 ms
10,880 KB |
testcase_19 | AC | 31 ms
10,880 KB |
testcase_20 | AC | 30 ms
10,880 KB |
testcase_21 | AC | 30 ms
10,880 KB |
testcase_22 | AC | 31 ms
11,008 KB |
testcase_23 | AC | 31 ms
11,008 KB |
testcase_24 | AC | 31 ms
11,136 KB |
testcase_25 | AC | 31 ms
10,880 KB |
testcase_26 | AC | 31 ms
11,008 KB |
testcase_27 | AC | 33 ms
10,880 KB |
testcase_28 | AC | 31 ms
10,880 KB |
testcase_29 | AC | 31 ms
11,008 KB |
testcase_30 | AC | 31 ms
10,880 KB |
testcase_31 | AC | 30 ms
11,136 KB |
testcase_32 | AC | 31 ms
10,880 KB |
testcase_33 | AC | 30 ms
10,880 KB |
testcase_34 | AC | 30 ms
10,880 KB |
testcase_35 | AC | 30 ms
11,008 KB |
testcase_36 | AC | 30 ms
11,008 KB |
testcase_37 | AC | 30 ms
10,880 KB |
testcase_38 | AC | 30 ms
10,880 KB |
testcase_39 | AC | 30 ms
11,008 KB |
testcase_40 | AC | 31 ms
11,008 KB |
testcase_41 | AC | 31 ms
11,136 KB |
ソースコード
def Berlekamp_Massey(A): n = len(A) B, C = [1], [1] l, m, p = 0, 1, 1 for i in range(n): d = A[i] for j in range(1, l + 1): d += C[j] * A[i - j] d %= mod if d == 0: m += 1 continue T = C.copy() q = pow(p, mod - 2, mod) * d % mod if len(C) < len(B) + m: C += [0] * (len(B) + m - len(C)) for j, b in enumerate(B): C[j + m] -= q * b C[j + m] %= mod if 2 * l <= i: B = T l, m, p = i + 1 - l, 1, d else: m += 1 res = [-c % mod for c in C[1:]] return res def Bostan_Mori(poly_nume,poly_deno,N,mod=0,convolve=None): #if type(poly_nume)==Polynomial: # poly_nume=poly_nume.polynomial #if type(poly_deno)==Polynomial: # poly_deno=poly_deno.polynomial if convolve==None: def convolve(poly_nume,poly_deno): conv=[0]*(len(poly_nume)+len(poly_deno)-1) for i in range(len(poly_nume)): for j in range(len(poly_deno)): x=poly_nume[i]*poly_deno[j] if mod: x%=mod conv[i+j]+=x if mod: for i in range(len(conv)): conv[i]%=mod return conv while N: poly_deno_=[-x if i%2 else x for i,x in enumerate(poly_deno)] if N%2: poly_nume=convolve(poly_nume,poly_deno_)[1::2] else: poly_nume=convolve(poly_nume,poly_deno_)[::2] poly_deno=convolve(poly_deno,poly_deno_)[::2] if mod: for i in range(len(poly_nume)): poly_nume[i]%=mod for i in range(len(poly_deno)): poly_deno[i]%=mod N//=2 return poly_nume[0] def BMBM(A,N,mod=0): deno=[1]+[-c for c in Berlekamp_Massey(A)] nume=[0]*(len(deno)-1) for i in range(len(A)): for j in range(len(deno)): if i+j<len(nume): nume[i+j]+=A[i]*deno[j] nume[i+j]%=mod return Bostan_Mori(nume,deno,N,mod=mod) A,B,N=map(int,input().split()) mod=10**9+7 X=[0,1] for i in range(2,10): X.append((A*X[i-1]+B*X[i-2])%mod) ans=BMBM(X,N,mod) print(ans)