結果

問題 No.2105 Avoid MeX
ユーザー taiga0629kyoprotaiga0629kyopro
提出日時 2022-10-21 22:12:17
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 612 ms / 2,000 ms
コード長 874 bytes
コンパイル時間 163 ms
コンパイル使用メモリ 82,140 KB
実行使用メモリ 88,504 KB
最終ジャッジ日時 2024-07-01 06:47:09
合計ジャッジ時間 5,664 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 6
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 



#############################
#############
cnb_max=10**5
mod=998244353
#############

kai=[1]*(cnb_max+1)
rkai=[1]*(cnb_max+1)
for i in range(cnb_max):
    kai[i+1]=kai[i]*(i+1)%mod

rkai[cnb_max]=pow(kai[cnb_max],mod-2,mod)
for i in range(cnb_max):
    rkai[cnb_max-1-i]=rkai[cnb_max-i]*(cnb_max-i)%mod

def cnb(x,y):
    if y>x:
        return 0
    if x<0:return 0
    if y<0:return 0
    return (kai[x]*rkai[y]%mod)*rkai[x-y]%mod


def inv(n):
    return kai[n-1]*rkai[n]%mod

##################################

c,x=map(int,input().split())
if x==0:
    print(inv(c+1))
    exit()
p=[0]*(x+3)
p[x]=1
while 1:
    if c>=x:break
    q=[0]*(x+3)
    for i in range(x+1):
        p[i]%=mod
        use=x-i
        q[i]+=p[i]*use*inv(c+1)
        q[i-1]+=p[i]*(c+1-use)*inv(c+1)
    p=q[:]
    c+=1
ans=0
for i in range(x+1):
    ans+=p[i]*(1-inv(i+1))
print(ans%mod)
0