結果
問題 | No.1340 おーじ君をさがせ |
ユーザー |
![]() |
提出日時 | 2021-01-15 23:15:19 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,530 ms / 2,000 ms |
コード長 | 1,071 bytes |
コンパイル時間 | 354 ms |
コンパイル使用メモリ | 82,008 KB |
実行使用メモリ | 88,832 KB |
最終ジャッジ日時 | 2024-11-27 23:21:45 |
合計ジャッジ時間 | 42,293 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 59 |
ソースコード
import sys input = sys.stdin.readline N,M,T=map(int,input().split()) mod=1+(1<<30) # 行列の計算(numpyを使えないとき,modを使用) def prod(A,B,k,l,m):# A:k*l,B:l*m C=[[None for i in range(m)] for j in range(k)] for i in range(k): for j in range(m): ANS=0 for pl in range(l): ANS=(ANS+A[i][pl]*B[pl][j])%mod C[i][j]=ANS return C def plus(A,B,k,l):# a,B:k*l C=[[None for i in range(l)] for j in range(k)] for i in range(k): for j in range(l): C[i][j]=(A[i][j]+B[i][j])%mod return C A=[[0]*N for i in range(N)] for i in range(M): a,b=map(int,input().split()) A[a][b]=1 # 漸化式を行列累乗で求める(ダブリング) POWA=[A] for i in range(60): POWA.append(prod(POWA[-1],POWA[-1],N,N,N)) # ベキを求めて X=[[0]*N] X[0][0]=1 #print(X) while T: X=prod(X,POWA[T.bit_length()-1],1,N,N) # n乗の場合 T-=1<<(T.bit_length()-1) #print(X) ANS=0 for i in X[0]: if i>0: ANS+=1 print(ANS)