結果

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

ソースコード

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 = -1
        tmp = 0 
        for ok in range(i,min(pre,N)):
            if comb(ok,i) <= n:
                OK = ok 
            else:
                break
        
        ok = OK
        ans.append(ok) 
        solve(n-comb(ok,i),i-1,ok)

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


n,k = mi() 

ans = [] 
main(n,k) 

print(*ans)
nokori = n 
for i in ans:
    nokori -= comb(i,k)
    #print(i,k,nokori) 
    k -= 1






def debug(ans,n,k):
    assert ans == sorted(ans)[::-1]
    for i in ans:
        assert ans.count(i) == 1
    tmp = 0 
    K = k
    for x in ans:
        if not x > k:
            print(n,K)
            #print(x,k)
            assert 0
        tmp += comb(x,k)
        k -= 1 
    
    if tmp != n:
        print(n,k) 
        print(ans)
        exit() 

deb = 0
if deb:
    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)
        print(ans)
0