結果

問題 No.2511 Mountain Sequence
ユーザー titia
提出日時 2023-11-19 03:20:47
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 674 ms / 2,000 ms
コード長 1,242 bytes
コンパイル時間 510 ms
コンパイル使用メモリ 82,304 KB
実行使用メモリ 166,400 KB
最終ジャッジ日時 2024-09-26 06:14:19
合計ジャッジ時間 16,190 ms
ジャッジサーバーID
(参考情報)
judge2 / judge3
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 32
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
input = sys.stdin.readline

mod=998244353

FACT=[1]
for i in range(1,3*10**5+1):
    FACT.append(FACT[-1]*i%mod)

FACT_INV=[pow(FACT[-1],mod-2,mod)]
for i in range(3*10**5,0,-1):
    FACT_INV.append(FACT_INV[-1]*i%mod)

FACT_INV.reverse()

def Combi(a,b):
    if 0<=b<=a:
        return FACT[a]*FACT_INV[b]%mod*FACT_INV[a-b]%mod
    else:
        return 0

N,M=map(int,input().split())
ANS=0

if N==M==1:
    print(1)
    exit()

"""
for x in range(1,M):
    for j in range(N+1):
        ANS+=Combi(x,j)*Combi(x,N-1-j)

    #print(ANS)

print(ANS%mod)
"""

F2=[1,1]

for i in range(2,6*10**5+1):
    F2.append(F2[-2]*i%mod)

def calc(x,y):
    if y==x*2+1:
        return 1
    if x-y//2>=0:
        A1=FACT[x]*pow(FACT[x-y//2],mod-2,mod)
    else:
        A1=0

    if 2*x-1-(y-1)//2*2>=0:
        A2=F2[2*x-1]*pow(F2[2*x-1-(y-1)//2*2],mod-2,mod)
    else:
        A2=0

    C1=FACT[(y-1)//2]
    C2=F2[y//2*2-1]

    #print(A1,A2,C1,C2)

    if y%2==0:
        A1*=2

    return A1*A2*pow(C1*C2,mod-2,mod)%mod

for x in range(1,M):
    ANS+=calc(x,N)

    #if calc(x,N)>0:
    #    print(x,calc(x,N))

    #X=0

    #for j in range(N+1):
    #    X+=Combi(x,j)*Combi(x,N-1-j)

    #print((x,N),calc(x,N),X)


print(ANS%mod)
    
0