結果

問題 No.2575 Almost Increasing Sequence
ユーザー chineristACchineristAC
提出日時 2023-12-03 03:49:05
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 723 ms / 10,000 ms
コード長 1,251 bytes
コンパイル時間 295 ms
コンパイル使用メモリ 81,700 KB
実行使用メモリ 81,392 KB
スコア 10,000
最終ジャッジ日時 2023-12-03 03:49:22
合計ジャッジ時間 16,356 ms
ジャッジサーバーID
(参考情報)
judge11 / judge13
純コード判定しない問題か言語
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 694 ms
81,392 KB
testcase_01 AC 720 ms
81,392 KB
testcase_02 AC 688 ms
81,392 KB
testcase_03 AC 713 ms
81,392 KB
testcase_04 AC 687 ms
81,392 KB
testcase_05 AC 716 ms
81,392 KB
testcase_06 AC 691 ms
81,392 KB
testcase_07 AC 700 ms
81,392 KB
testcase_08 AC 684 ms
81,392 KB
testcase_09 AC 698 ms
81,392 KB
testcase_10 AC 715 ms
81,392 KB
testcase_11 AC 692 ms
81,392 KB
testcase_12 AC 717 ms
81,392 KB
testcase_13 AC 675 ms
81,392 KB
testcase_14 AC 723 ms
81,392 KB
testcase_15 AC 683 ms
81,392 KB
testcase_16 AC 699 ms
81,392 KB
testcase_17 AC 691 ms
81,392 KB
testcase_18 AC 690 ms
81,392 KB
testcase_19 AC 704 ms
81,392 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

import sys
from itertools import permutations
from heapq import heappop,heappush
from collections import deque
import random
import bisect

input = lambda :sys.stdin.readline().rstrip()
mi = lambda :map(int,input().split())
li = lambda :list(mi())

def cmb(n, r, mod):
    if ( r<0 or r>n ):
        return 0
    return (g1[n] * g2[r] % mod) * g2[n-r] % mod


mod = 998244353
N = 2*10**5
g1 = [1]*(N+1)
g2 = [1]*(N+1)
inverse = [1]*(N+1)

for i in range( 2, N + 1 ):
    g1[i]=( ( g1[i-1] * i ) % mod )
    inverse[i]=( ( -inverse[mod % i] * (mod//i) ) % mod )
    g2[i]=( (g2[i-1] * inverse[i]) % mod )
inverse[0]=0

def calc_small(N,K):
    n_pow = [pow(i,K,mod) for i in range(N+1)]

    res = [0] * (N+1)
    for a1 in range(1,N+1):
        for a2 in range(1,a1+1):
            for a3 in range(1,a2+1):
                n = a1 + a2 + a3
                if n > N:
                    continue

                res[n] += (((a1-a2+1) * (a1-a3+2) % mod) * ((a2-a3+1) % mod) % mod * (g2[a1+2] * g2[a2+1] % mod) * g2[a3] % mod)**2 % mod * n_pow[a1] % mod
    
    for n in range(N+1):
        res[n] *= g1[n] * g1[n]
        res[n] %= mod

    return res[1:]

def calc_mid(N,K):
    assert False

K = int(input())
N = 500
print(N)
print(*calc_small(N,K))
0