結果

問題 No.3030 Kruskal-Katona
ユーザー Yotugi
提出日時 2025-02-21 22:05:32
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,706 bytes
コンパイル時間 383 ms
コンパイル使用メモリ 82,552 KB
実行使用メモリ 64,592 KB
最終ジャッジ日時 2025-02-21 22:05:35
合計ジャッジ時間 3,160 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2 WA * 1
other AC * 8 WA * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

from sys import stdin,setrecursionlimit#,set_int_max_str_digits
#import pypyjit
#pypyjit.set_param('max_unroll_recursion=-1')
setrecursionlimit(20000000) # これこどふぉだと無理
#set_int_max_str_digits(200010)

mod = 998244353  
ii = lambda :int(stdin.readline())
mi = lambda :map(int,stdin.readline().split())
li = lambda :list(mi())
gmi = lambda :map(lambda x: int(x) - 1, stdin.readline().split())
gi = lambda :list(map(lambda x: 0 if x == "." else 1,input())) # グリッド入力受け取り
py = lambda :print("Yes")
pn = lambda :print("No")
pf = lambda :print("First") 
ps = lambda :print("Second")
vec = [(1,0),(-1,0),(0,-1),(0,1)]
vec1 = [(1,0),(1,1),(0,1),(-1,1),(-1,0),(-1,-1),(0,-1),(1,-1)] #8方向
inf = 10**18


n,k = mi() 

N = 10**5 
fact = [0 for _ in range(N)]
invfact = [0 for _ in range(N)]
fact[0] = 1
for i in range(1, N):
    fact[i] = fact[i - 1] * i % mod

invfact[N - 1] = pow(fact[N - 1], mod - 2, mod)

for i in range(N - 2, -1, -1):
    invfact[i] = invfact[i + 1] * (i + 1) % mod

def comb(n, k):
    if k < 0 or n < k:
        return 0
    else:
        return (fact[n] * invfact[k] % mod) * invfact[n - k] % mod

def nHk(n, k):
    return comb(n + k - 1, k)

def p(n,k):
    if n < k:
        return 0
    return fact[n]*invfact[n-k]%mod


if k == 1:
    print(n)
    exit() 

ans = [] 

def solve(n,i):
    global ans 
    if n == 0:
        return 
    if i == 1:
        ans.append(n) 
        return 
    
    ok = i
    ng = min(N,n)

    while abs(ok-ng) > 1:
        mid = (ok+ng)//2 
        if comb(mid,i) <= n:
            ok = mid 
        else:
            ng = mid 
    
    ans.append(ok) 
    solve(n-comb(ok,i),i-1)

solve(n,k)
print(*ans)
0