結果
問題 | No.2511 Mountain Sequence |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
import sysinput = sys.stdin.readlinemod=998244353FACT=[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]%modelse:return 0N,M=map(int,input().split())ANS=0if 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 1if x-y//2>=0:A1=FACT[x]*pow(FACT[x-y//2],mod-2,mod)else:A1=0if 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=0C1=FACT[(y-1)//2]C2=F2[y//2*2-1]#print(A1,A2,C1,C2)if y%2==0:A1*=2return A1*A2*pow(C1*C2,mod-2,mod)%modfor 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)