結果

問題 No.1750 ラムドスウイルスの感染拡大-hard
ユーザー NatsubiSogan
提出日時 2021-11-19 23:33:15
言語 Python3
(3.13.1 + numpy 2.2.1 + scipy 1.14.1)
結果
TLE  
実行時間 -
コード長 1,624 bytes
コンパイル時間 418 ms
コンパイル使用メモリ 12,544 KB
実行使用メモリ 26,320 KB
最終ジャッジ日時 2025-01-01 22:53:32
合計ジャッジ時間 71,853 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 4
other AC * 7 TLE * 23
権限があれば一括ダウンロードができます

ソースコード

diff # 

import typing

class Matrix:
	def __init__(self, n: int, m: int, mat: typing.Union[list, None] = None, mod: int = 998244353) -> None:
		self.n = n
		self.m = m
		self.mat = [[0] * self.m for i in range(self.n)]
		self.mod = mod
		if mat:
			for i in range(self.n):
				self.mat[i] = mat[i]
	
	def is_square(self) -> None:
		return self.n == self.m
	
	def __getitem__(self, key: int) -> int:
		if isinstance(key, slice):
			return self.mat[key]
		else:
			assert key >= 0
			return self.mat[key]

	def __mul__(self, other: typing.Union["Matrix", int]) -> "Matrix":
		if other.__class__ == Matrix:
			res = [[0] * other.m for i in range(self.n)]
			for i in range(self.n):
				for k in range(self.m):
					for j in range(other.m):
						res[i][j] += self[i][k] * other[k][j]
			for i in range(self.n):
				for j in range(other.m): res[i][j] %= self.mod
			return Matrix(self.n, other.m, res)
		else:
			return self.times(other)
	
	def __rmul__(self, other: typing.Union["Matrix", int]) -> "Matrix":
		return self.times(other)

	def __pow__(self, k: int) -> "Matrix":
		tmp = Matrix(self.n, self.n, self.mat)
		res = Matrix(self.n, self.n)
		for i in range(self.n):
			res[i][i] = 1
		while k:
			if k & 1:
				res *= tmp
			tmp *= tmp
			k >>= 1
		return res

def main():
	import sys
	input = sys.stdin.buffer.readline
	n, m, t = map(int, input().split())
	G = [[0] * n for i in range(n)]
	for _ in range(m):
		a, b = map(int, input().split())
		G[a][b] = 1
		G[b][a] = 1
	G = Matrix(n, n, G)
	G **= t
	b = [[0] * n for i in range(n)]
	b[0][0] = 1
	G *= Matrix(n, n, b)
	print(G.mat[0][0])
if __name__ == '__main__':
	main()
0