結果

問題 No.1771 A DELETEQ
ユーザー LyricalMaestro
提出日時 2024-12-15 00:16:00
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 198 ms / 3,500 ms
コード長 2,011 bytes
コンパイル時間 312 ms
コンパイル使用メモリ 82,016 KB
実行使用メモリ 62,592 KB
最終ジャッジ日時 2024-12-15 00:16:06
合計ジャッジ時間 4,347 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 26 WA * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

# https://yukicoder.me/problems/no/895

MOD = 998244353


class CombinationCalculator:
    """
    modを考慮したPermutation, Combinationを計算するためのクラス
    """    
    def __init__(self, size, mod):
        self.mod = mod
        self.factorial = [0] * (size + 1)
        self.factorial[0] = 1
        for i in range(1, size + 1):
            self.factorial[i] = (i * self.factorial[i - 1]) % self.mod
        
        self.inv_factorial = [0] * (size + 1)
        self.inv_factorial[size] = pow(self.factorial[size], self.mod - 2, self.mod)

        for i in reversed(range(size)):
            self.inv_factorial[i] = ((i + 1) * self.inv_factorial[i + 1]) % self.mod

    def calc_combination(self, n, r):
        if n < 0 or n < r or r < 0:
            return 0

        if r == 0 or n == r:
            return 1
        
        ans = self.inv_factorial[n - r] * self.inv_factorial[r]
        ans %= self.mod
        ans *= self.factorial[n]
        ans %= self.mod
        return ans
    
    def calc_permutation(self, n, r):
        if n < 0 or n < r:
            return 0

        ans = self.inv_factorial[n - r]
        ans *= self.factorial[n]
        ans %= self.mod
        return ans
        

def main():
    x, y = map(int, input().split())

    if not (x <= 4000 and y <= 4000):
        return

    combi = CombinationCalculator(x + y, MOD)
    answer = 0

    pow2 = [0] * (x + y + 1)
    p = 1
    for k in range(x + y + 1):
        pow2[k] = p
        p *= 2
        p %= MOD

    for z1 in range(min(x, y) + 1):
        for z0 in range(z1 + 1):

            ans = combi.factorial[x + y - 2 * z1 + z0]
            ans *= combi.inv_factorial[x - z1]
            ans %= MOD
            ans *= combi.inv_factorial[y - z1]
            ans %= MOD
            ans *= combi.inv_factorial[z0]
            ans %= MOD
            ans *= pow2[z0]
            ans %= MOD
            answer += ans
            answer %= MOD
    print(answer)




if __name__ == "__main__":
    main()
0