結果
問題 | No.2369 Some Products |
ユーザー |
|
提出日時 | 2023-06-30 22:34:15 |
言語 | PyPy3 (7.3.15) |
結果 |
AC
|
実行時間 | 1,584 ms / 2,500 ms |
コード長 | 1,553 bytes |
コンパイル時間 | 205 ms |
コンパイル使用メモリ | 82,412 KB |
実行使用メモリ | 470,368 KB |
最終ジャッジ日時 | 2024-07-07 10:25:20 |
合計ジャッジ時間 | 10,677 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 14 |
ソースコード
import sysinput = lambda :sys.stdin.readline()[:-1]ni = lambda :int(input())na = lambda :list(map(int,input().split()))yes = lambda :print("yes");Yes = lambda :print("Yes");YES = lambda : print("YES")no = lambda :print("no");No = lambda :print("No");NO = lambda : print("NO")#######################################################################mod = 998244353nn = 10**5fact = [1] * nnfor i in range(nn - 1):fact[i + 1] = fact[i] * (i + 1) % modinvfact = [1] * nninvfact[nn - 1] = pow(fact[nn - 1], mod - 2, mod)for i in range(nn - 1)[::-1]:invfact[i] = invfact[i + 1] * (i + 1) % moddef binom(x, y):if x < 0 or y < 0 or x - y < 0:return 0return fact[x] * invfact[y] % mod * invfact[x - y] % modn = ni()a = na()mod = 998244353dp = [[0 for j in range(n+1)]for i in range(n+1)]dp[0][0] = 1for i in range(n):for j in range(n):if dp[i][j] == 0:continuedp[i+1][j] += dp[i][j]dp[i+1][j] %= modif j + 1 <= n:dp[i+1][j+1] += dp[i][j] * a[i] % moddp[i+1][j+1] %= moddp2 = [[0 for j in range(n+1)]for i in range(n+1)]dp2[0][0] = 1for i in range(1, n+1):for j in range(n+1):if j:dp2[i][j] = (dp2[i-1][j] + dp2[i][j-1] * (-a[i-1])%mod)%modelse:dp2[i][j] = dp2[i-1][j]q = ni()for _ in range(q):a,b,k = na()a -= 1ans = 0#print(dp[b])#print(dp2[a])for i in range(k+1):ans += dp[b][i] * dp2[a][k-i] % modans %= modprint(ans)