結果
問題 |
No.2445 奇行列式
|
ユーザー |
![]() |
提出日時 | 2023-08-25 23:29:47 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,065 bytes |
コンパイル時間 | 1,532 ms |
コンパイル使用メモリ | 81,920 KB |
実行使用メモリ | 325,376 KB |
最終ジャッジ日時 | 2024-12-24 11:29:57 |
合計ジャッジ時間 | 18,657 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 17 WA * 3 |
ソースコード
N,mod=map(int,input().split()) dp=[0]*(2**(N+1)) z=2**N z*=N u=[0]*z G=[[] for i in range(N+1)] for bit in range(2**N): v=[] for k in range(N): if (bit>>k)&1: v.append(k) count=0 G[len(v)].append(bit) for x in range(N): for j in range(len(v)-1,-1,-1): if v[j]>x: count+=1 else: break count%=2 u[N*bit+x]=count A=[list(map(int,input().split())) for i in range(N)] for i in range(N): for j in range(N): A[i][j]%=mod for j in range(N): bit=2**j if j<100: dp[2*bit]=A[0][j] else: dp[2*bit+1]=A[0][j] for i in range(N-1): for bit in G[i+1]: for k in range(2): if dp[2*bit+k]==0: continue for j in range(N): if (bit>>j)&1: continue w=u[N*bit+j] if w==0: k2=k dp[2*(bit+2**j)+k2]+=dp[2*bit+k]*A[i+1][j] dp[2*(bit+2**j)+k2]%=mod else: k2=(k+1)%2 dp[2*(bit+2**j)+k2]+=dp[2*bit+k]*A[i+1][j] dp[2*(bit+2**j)+k2]%=mod result=dp[2**(N+1)-1] result%=mod print(result)