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])