#include using namespace std; #define rep(i,a,b) for (int i=(int)(a);i<(int)(b);i++) #define all(p) p.begin(),p.end() namespace po167{ template struct Binomial{ std::vector fact_vec, fact_inv_vec; void extend(int m = -1){ int n = fact_vec.size(); if (m == -1) m = n * 2; if (n >= m) return; fact_vec.resize(m); fact_inv_vec.resize(m); for (int i = n; i < m; i++){ fact_vec[i] = fact_vec[i - 1] * T(i); } fact_inv_vec[m - 1] = T(1) / fact_vec[m - 1]; for (int i = m - 1; i > n; i--){ fact_inv_vec[i - 1] = fact_inv_vec[i] * T(i); } } Binomial(int MAX = 0){ fact_vec.resize(1, T(1)); fact_inv_vec.resize(1, T(1)); extend(MAX + 1); } T fact(int i){ if (i < 0) return 0; while (int(fact_vec.size()) <= i) extend(); return fact_vec[i]; } T invfact(int i){ if (i < 0) return 0; while (int(fact_inv_vec.size()) <= i) extend(); return fact_inv_vec[i]; } }; } #include using mint = atcoder::modint998244353; // DEAR MYSTERIES / TOMOO int main() { ios::sync_with_stdio(false); cin.tie(nullptr); int N, K; cin >> N >> K; po167::Binomial table; vector ans(N), res(N); vector g(N); rep(i, 0, N) g[i] = (mint(N - 1 - i)).pow(K); /* * v : 左側の (1, 2) 成分の多項式 * 返り値 : (a, b) * a : (1, 2) 成分 * b : (2, 2) 成分 */ auto f = [&](auto self, int l, int r, vector v) -> pair, vector> { /*cout << "# " << l << " " << r << endl; for (auto x : v) cout << x.val() << " "; cout << endl;*/ if (l + 1 == r){ ans[l] = v.back(); return {{l + 1}, {l + 1, 2}}; } int m = (l + r) / 2; vector nv; // v のサイズ調整をする // m - l しかいらない nv.resize(m - l); for (int i = 0; i < m - l; i++) nv[i] = v[i + r - m]; auto [LA, LB] = self(self, l, m, nv); // v の更新をする // nv = v * LB + g * fact(l) * LA { nv.resize(r - m); auto nv1 = atcoder::convolution(v, LB); vector ng(r - l); rep(i, 0, r - l) ng[i] = g[i + N - r + l]; auto nv2 = atcoder::convolution(ng, LA); rep(i, 0, r - m) nv[i] = nv1[i + m - l] + nv2[i + m - l] * table.fact(l); } auto [RA, RB] = self(self, m, r, nv); // A, B を返す // A = LA * RB + fact(l, m) * RA // B = LB * RB vector nA, nB; nA = atcoder::convolution(LA, RB); rep(i, 0, RA.size()){ nA[i] += RA[i] * table.fact(m) * table.invfact(l); } nB = atcoder::convolution(LB, RB); return {nA, nB}; }; f(f, 0, N, g); rep(i, 1, N){ res[i] = (ans[i] - ans[i - 1] * i) * table.fact(N - 1) * table.invfact(i); } mint inv2 = (mint(2)).inv(); rep(i, 0, N){ res[i] *= inv2; cout << res[i].val() << "\n"; } }