結果
問題 | 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 PyPy3 from collections import Counter, defaultdict, deque import itertools import re import math from functools import reduce import operator import bisect from heapq import * import functools mod=998244353 import sys input=sys.stdin.readline #nCrをmodで割ったあまりを求める def ncr(n, r): if ( r<0 or r>n ): return 0 return g1[n] * g2[r] % mod * g2[n-r] % mod def npr(n, r): if ( r<0 or r>n ): return 0 return g1[n] * g2[n-r] % mod N = 10 ** 4 g1 = [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 = 0 for i in range(n): leng = n - 1 + 2 * i + n - i tmp = 10 * i + n - i if tmp > leng: break rem = leng - tmp dp = [[0] * (rem + 1) for _ in range(n+1-i)] dp[0][0] = 1 for j in range(1,n+1-i): cumsum = [0] tmp = 0 for k in range(rem+1): tmp += dp[j-1][k] tmp %= mod cumsum.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) % mod ans %= mod else: ans += dp[-1][-1] print(ans)