結果
問題 | 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 sysinput = sys.stdin.readlineN,M,T=map(int,input().split())mod=1+(1<<30)# 行列の計算(numpyを使えないとき,modを使用)def prod(A,B,k,l,m):# A:k*l,B:l*mC=[[None for i in range(m)] for j in range(k)]for i in range(k):for j in range(m):ANS=0for pl in range(l):ANS=(ANS+A[i][pl]*B[pl][j])%modC[i][j]=ANSreturn Cdef plus(A,B,k,l):# a,B:k*lC=[[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])%modreturn CA=[[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=0for i in X[0]:if i>0:ANS+=1print(ANS)