結果
| 問題 |
No.1693 Invasion
|
| コンテスト | |
| ユーザー |
terasa
|
| 提出日時 | 2022-06-14 21:29:24 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 178 ms / 2,000 ms |
| コード長 | 2,135 bytes |
| コンパイル時間 | 196 ms |
| コンパイル使用メモリ | 82,096 KB |
| 実行使用メモリ | 84,480 KB |
| 最終ジャッジ日時 | 2024-10-02 22:32:42 |
| 合計ジャッジ時間 | 3,349 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 21 |
ソースコード
import sys
import pypyjit
import itertools
import heapq
import math
from collections import deque, defaultdict
import bisect
input = sys.stdin.readline
sys.setrecursionlimit(10 ** 6)
pypyjit.set_param('max_unroll_recursion=-1')
def index_lt(a, x):
'return largest index s.t. A[i] < x or -1 if it does not exist'
return bisect.bisect_left(a, x) - 1
def index_le(a, x):
'return largest index s.t. A[i] <= x or -1 if it does not exist'
return bisect.bisect_right(a, x) - 1
def index_gt(a, x):
'return smallest index s.t. A[i] > x or len(a) if it does not exist'
return bisect.bisect_right(a, x)
def index_ge(a, x):
'return smallest index s.t. A[i] >= x or len(a) if it does not exist'
return bisect.bisect_left(a, x)
class Combination:
def __init__(self, N, p):
self.N = N
self.p = p
self.f = [None] * (N + 1)
self.finv = [None] * (N + 1)
self.inv = [None] * (N + 1)
self.f[0] = 1
self.f[1] = 1
self.finv[0] = 1
self.finv[1] = 1
self.inv[1] = 1
for i in range(2, N + 1):
self.f[i] = self.f[i - 1] * i % p
self.inv[i] = p - self.inv[p % i] * (p // i) % p
self.finv[i] = self.finv[i - 1] * self.inv[i] % p
def P(self, n, k):
if n < k:
return 0
if n < 0 or k < 0:
return 0
return self.f[n] * self.finv[n - k] % self.p
def C(self, n, k):
if n < k:
return 0
if n < 0 or k < 0:
return 0
return self.f[n] * (self.finv[k] * self.finv[n - k] % self.p) % self.p
# 重複組合せ
# n種類のものからk個選ぶ
def H(self, n, k):
return self.C(k + n - 1, k)
N, M = map(int, input().split())
A = list(map(int, input().split()))
mod = 998244353
comb = Combination(M, mod)
INF = 1 << 30
dpmin = [INF] * (M + 1)
dpmin[0] = 0
for i in range(M):
for a in A:
if i + a <= M:
dpmin[i + a] = min(dpmin[i + a], dpmin[i] + 1)
ans = 0
for i in range(M + 1):
ans += comb.C(M - dpmin[i], i - dpmin[i])
ans %= mod
print(ans)
terasa