結果

問題 No.3030 Kruskal-Katona
ユーザー Yotugi
提出日時 2025-02-21 22:20:59
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 2,189 bytes
コンパイル時間 444 ms
コンパイル使用メモリ 82,220 KB
実行使用メモリ 63,404 KB
最終ジャッジ日時 2025-02-21 22:21:10
合計ジャッジ時間 3,147 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15 WA * 12
権限があれば一括ダウンロードができます

ソースコード

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 = 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



def main(n,k):
    if k == 1:
        print(n)
        exit() 

    def solve(n,i,pre):
        global ans 
        if n == 0:
            return 
        if i == 1:
            ans.append(n) 
            return 
        ok = i
        ng = min(pre,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,ok)

    solve(n,k,n+1)
    print(*ans)



def debug(ans,n,k):
    tmp = 0 
    for x in ans:
        tmp += comb(x,k)
        k -= 1 
    
    if tmp != n:
        print(n,k) 
        print(ans)
        exit() 


n,k = mi() 
ans = [] 
main(n,k) 
# from random import randint

# for _ in range(10):
#     k = randint(2,100) 
#     n = randint(k+1,100000000)
#     ans = [] 
#     print(n,k)
#     main(n,k)
#     debug(ans,n,k)
0