#include using namespace std; const long long MOD = 998244353; unordered_map mp; long long modpow(long long a, long long b){ if (mp.count(b)){ return mp[b]; } else { long long ans = 1; while (b > 0){ if (b % 2 == 1){ ans *= a; ans %= MOD; } a *= a; a %= MOD; b /= 2; } mp[b] = ans; return mp[b]; } } unordered_map mp2; long long modinv(long long a){ if (mp2.count(a)){ return mp2[a]; } else { mp2[a] = modpow(a, MOD - 2); return mp2[a]; } } vector ntt(vector A, bool inv){ int N = A.size(); long long r = 3; if (inv){ r = modinv(r); } vector B(N); for (int i = N / 2; i > 0; i /= 2){ long long z = modpow(r, (MOD - 1) / (N / i)); long long z2 = 1; for (int j = 0; j < N; j += i * 2){ for (int k = 0; k < i; k++){ A[i + j + k] *= z2; A[i + j + k] %= MOD; B[j / 2 + k] = (A[j + k] + A[i + j + k]) % MOD; B[N / 2 + j / 2 + k] = (A[j + k] - A[i + j + k] + MOD) % MOD; } z2 = z2 * z % MOD; } swap(A, B); } if (inv){ int Ninv = modinv(N); for (int i = 0; i < N; i++){ A[i] = A[i] * Ninv % MOD; } } return A; } vector convolution(vector A, vector B, int d = -1){ int deg = A.size() + B.size() - 1; int N = 1; while (N < deg){ N *= 2; } A.resize(N); B.resize(N); A = ntt(A, false); B = ntt(B, false); vector ans(N); for (int i = 0; i < N; i++){ ans[i] = A[i] * B[i] % MOD; } ans = ntt(ans, true); if (d == -1){ ans.resize(deg); } else { ans.resize(d); } return ans; } vector add(vector A, vector B){ int N = max(A.size(), B.size()); while (A.size() < N){ A.push_back(0); } while (B.size() < N){ B.push_back(0); } vector ans(N, 0); for (int i = 0; i < N; i++){ ans[i] = A[i] + B[i]; if (ans[i] > MOD){ ans[i] -= MOD; } } return ans; } vector inverse(vector &P){ int N = P.size(); vector ans(1); ans[0] = modinv(P[0]); for (int i = 1; i < N; i *= 2){ vector ans2 = ans; ans2.resize(i * 4); vector f(min(N, i * 2)); for (int j = 0; j < min(N, i * 2); j++){ f[j] = P[j]; } f.resize(i * 4); ans2 = convolution(ans2, ans2); ans2 = convolution(ans2, f); ans2.resize(i * 2); for (int j = 0; j < i; j++){ ans2[j] = MOD - ans2[j] + ans[j] * 2; ans2[j] %= MOD; } swap(ans, ans2); } ans.resize(N); return ans; } int main(){ int N, M; cin >> N >> M; vector A(N); for (int i = 0; i < N; i++){ cin >> A[i]; } int N2 = 1; while (N2 < N){ N2 *= 2; } while (A.size() < N2){ A.push_back(0); } N = N2; vector, int>> ST1(N2 * 2 - 1); for (int i = 0; i < N; i++){ vector f = {1, MOD - A[i]}; ST1[N2 - 1 + i] = make_pair(f, 1); } for (int i = N2 - 2; i >= 0; i--){ int deg = min(M + 1, ST1[i * 2 + 1].second + ST1[i * 2 + 2].second + 1); ST1[i] = make_pair(convolution(ST1[i * 2 + 1].first, ST1[i * 2 + 2].first, deg), deg); } while (ST1[0].first.size() <= M * 2){ ST1[0].first.push_back(0); } /* for (int i = 0; i <= M; i++){ cout << ST1[0].first[i] << ' '; } cout << endl; */ vector, int>> ST2(N2 * 2 - 1); for (int i = 0; i < N; i++){ vector f = {1}; ST2[N2 - 1 + i] = make_pair(f, 1); } for (int i = N2 - 2; i >= 0; i--){ if (ST2[i * 2 + 2].second == 0){ ST2[i] = ST2[i * 2 + 1]; } else { int deg = min(M + 1, ST2[i * 2 + 1].second + ST2[i * 2 + 2].second); ST2[i] = make_pair(add(convolution(ST2[i * 2 + 1].first, ST1[i * 2 + 2].first, deg), convolution(ST1[i * 2 + 1].first, ST2[i * 2 + 2].first, deg)), deg); } } while (ST2[0].first.size() <= M * 2){ ST2[0].first.push_back(0); } /* for (int i = 0; i <= M; i++){ cout << ST2[0].first[i] << ' '; } cout << endl; */ vector g = ST1[0].first; vector f = ST2[0].first; //f/g vector g_inv = inverse(g); /* vector h = convolution(g, g_inv); for (int i = 0; i <= M; i++){ cout << h[i] << ' '; } cout << endl; */ vector ans = convolution(f, g_inv); for (int i = 1; i <= M; i++){ cout << ans[i]; if (i < M){ cout << ' '; } } cout << endl; }