結果
問題 | No.2503 Typical Path Counting Problem on a Grid |
ユーザー |
![]() |
提出日時 | 2023-10-13 23:56:31 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,213 ms / 2,000 ms |
コード長 | 971 bytes |
コンパイル時間 | 161 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 392,524 KB |
最終ジャッジ日時 | 2024-09-15 19:49:23 |
合計ジャッジ時間 | 10,592 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 10 |
ソースコード
def matrix_exponentiation(a,r,b):n = len(a); m = r.bit_length()res = [a]+[[[0]*n for _ in range(n)] for _ in range(m-1)]for h in range(m-1):for i in range(n):for j in range(n):for k in range(n): res[h+1][i][j] += res[h][i][k]*res[h][k][j]%modres[h+1][i][j] %= modfor h in range(m):if not r>>h&1: continuetmp = [0]*nfor i in range(n):for j in range(n): tmp[i] += res[h][i][j]*b[j]%modtmp[i] %= modb = tmpreturn bmod = 998244353; k = 10**7a = list(range(2,2*k+2,2))b = list(range(1,k+1))dp = [1]+[0]*kfor i in range(k):if i<=k-2: dp[i+2] = dp[i]*b[i]%moddp[i+1] = (dp[i+1]+dp[i]*a[i])%modfor _ in range(int(input())):n,m = map(int,input().split())if n>m: n,m = m,nif n==0: print(1); continuex,y = matrix_exponentiation([[2*n+1,n],[1,0]],m-n,[dp[n],dp[n-1]])print((x*dp[n]+y*n*dp[n-1])%mod)