結果
問題 | No.1964 sum = length |
ユーザー |
![]() |
提出日時 | 2022-06-03 21:48:29 |
言語 | PyPy3 (7.3.15) |
結果 |
WA
|
実行時間 | - |
コード長 | 1,596 bytes |
コンパイル時間 | 300 ms |
コンパイル使用メモリ | 82,048 KB |
実行使用メモリ | 79,488 KB |
最終ジャッジ日時 | 2024-09-21 02:37:10 |
合計ジャッジ時間 | 5,943 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 19 WA * 21 |
ソースコード
#!/usr/bin/env PyPy3from collections import Counter, defaultdict, dequeimport itertoolsimport reimport mathfrom functools import reduceimport operatorimport bisectfrom heapq import *import functoolsmod=998244353import sysinput=sys.stdin.readline#nCrをmodで割ったあまりを求めるdef ncr(n, r):if ( r<0 or r>n ):return 0return g1[n] * g2[r] % mod * g2[n-r] % moddef npr(n, r):if ( r<0 or r>n ):return 0return g1[n] * g2[n-r] % modN = 10 ** 4g1 = [1, 1] # 元テーブルg2 = [1, 1] #逆元テーブルinverse = [0, 1] #逆元テーブル計算用テーブルfor i in range( 2, N + 1 ):g1.append( ( g1[-1] * i ) % mod )inverse.append( ( -inverse[mod % i] * (mod//i) ) % mod )g2.append( (g2[-1] * inverse[-1]) % mod )n = int(input())ans = 0for i in range(n):leng = n - 1 + 2 * i + n - itmp = 10 * i + n - iif tmp > leng:breakrem = leng - tmpdp = [[0] * (rem + 1) for _ in range(n+1-i)]dp[0][0] = 1for j in range(1,n+1-i):cumsum = [0]tmp = 0for k in range(rem+1):tmp += dp[j-1][k]tmp %= modcumsum.append(tmp)for k in range(rem+1):dp[j][k] = cumsum[k+1] - cumsum[max(0,k-8)]dp[j][k] %= mod#print(dp)if i != 0:for j in range(rem+1):#print(dp[-1][j] , ncr(i+rem-j-1,i-1) % mod , ncr(n,i) % mod)ans += dp[-1][j] * ncr(i+rem-j-1,i-1) % mod * ncr(n,i) % modans %= modelse:ans += dp[-1][-1]print(ans)