結果

問題 No.2763 Macaron Gift Box
ユーザー sasa8uyauyasasa8uyauya
提出日時 2024-09-19 22:13:00
言語 PyPy3
(7.3.15)
結果
TLE  
実行時間 -
コード長 2,072 bytes
コンパイル時間 245 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 261,692 KB
最終ジャッジ日時 2024-09-19 22:13:49
合計ジャッジ時間 9,279 ms
ジャッジサーバーID
(参考情報)
judge5 / judge2
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 TLE -
testcase_01 TLE -
testcase_02 TLE -
testcase_03 TLE -
testcase_04 TLE -
testcase_05 TLE -
testcase_06 TLE -
testcase_07 TLE -
testcase_08 TLE -
testcase_09 TLE -
testcase_10 TLE -
testcase_11 -- -
testcase_12 -- -
testcase_13 -- -
testcase_14 -- -
testcase_15 -- -
testcase_16 -- -
権限があれば一括ダウンロードができます

ソースコード

diff #

mod = 998244353
R = 3
Rinv = 332748118
W = [pow(R, (mod-1)>>i, mod) for i in range(24)]
Winv = [pow(Rinv, (mod-1)>>i, mod) for i in range(24)]


def fft(k, f):
	for l in range(k, 0, -1):
		d = 1<<l-1
		U = [1]
		for i in range(d):
			U.append(U[-1]*W[l]%mod)
		for i in range(1<<k-l):
			for j in range(d):
				s = i*2*d+j
				f[s], f[s+d] = (f[s]+f[s+d])%mod, U[j]*(f[s]-f[s+d])%mod


def fftinv(k, f):
	for l in range(1, k+1):
		d = 1<<l-1
		for i in range(1<<k-l):
			u = 1
			for j in range(i*2*d, (i*2+1)*d):
				f[j+d] *= u
				f[j], f[j+d] = (f[j]+f[j+d])%mod, (f[j]-f[j+d])%mod
				u *= Winv[l]
				u %= mod


def convolution(a, b):
	le = len(a)+len(b)-1
	k = le.bit_length()
	n = 1<<k
	a = a+[0]*(n-len(a))
	b = b+[0]*(n-len(b))
	fft(k, a)
	fft(k, b)
	for i in range(n):
		a[i] *= b[i]
		a[i] %= mod
	fftinv(k, a)
	ninv = pow(n, mod-2, mod)
	for i in range(le):
		a[i] *= ninv
		a[i] %= mod
	return a[:le]


def FPSinv(H):
  I=[pow(H[0],M-2,M)]
  l=1
  while l<len(H):
    I+=[0]*l
    nI=convolution(H[:l*2],convolution(I,I)[:l*2])[:l*2]
    for i in range(l*2):
      nI[i]=(2*I[i]-nI[i])%M
    I=nI
    l*=2
  return I[:len(H)]

def FPSlog(H):
  H1=[H[i]*i%M for i in range(1,len(H))]+[0]
  H2=FPSinv(H)
  I=convolution(H1,H2)
  I=[0]+[I[i]*fb[i+1]*fa[i]%M for i in range(len(H)-1)]
  return I

def FPSexp(H):
  I=[1]
  l=1
  while l<len(H):
    I+=[0]*l
    I2=FPSlog(I)[:l*2]
    I3=H[:l*2]
    I3[0]+=1
    for i in range(l*2):
      I3[i]-=I2[i]
    nI=convolution(I,I3)[:l*2]
    I=nI
    l*=2
  return I[:len(H)]

n,k=map(int,input().split())

L=1<<18
M=998244353

fa=[1,1]
fb=[1,1]
for i in range(2,L+1):
  fa+=[fa[-1]*i%M]
  fb+=[fb[-1]*(M//i)*fb[M%i]*fa[M%i-1]*(-1)%M]

q1=[1,-1]+[0]*(L-2)
q2=[0]+[-fb[i]*fa[i-1]%M for i in range(1,L)]
q3=[0]*L
for i in range(1,n+1):
  for j in range(L//(i)+1):
    if i*j<L:
      q3[i*j]-=q2[j]
      q3[i*j]%=M
    else:
      break
for i in range(1,n+1):
  for j in range(L//(i*(k+1))+1):
    if i*j*(k+1)<L:
      q3[i*j*(k+1)]+=q2[j]
      q3[i*j*(k+1)]%=M
    else:
      break
q4=FPSexp(q3)

print(*q4[1:n+1])
0