#include #include #include using namespace std; using namespace atcoder; using mint = modint998244353; const mint inv2 = mint(2).inv(), inv3 = mint(3).inv(), inv5 = mint(5).inv(); // 多項式 using F = vector; template using matrix = vector>; // 多項式の和 F operator+=(F &f, const F &g) { if (f.size() < g.size()) f.resize(g.size()); for (int i = 0; i < g.size(); i++) f[i] += g[i]; return f; } // 遷移行列 const matrix P = { {{inv2, 0}, {0, inv2}, {0, 0}}, {{inv3, 0}, {0, inv3}, {0, inv3}}, {{0 , 0}, {0, 2*inv3}, {0, inv3}} }; // 行列積 matrix operator*=(matrix &A, const matrix &B) { int N = A.size(); matrix C(N, vector(N)); for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { for (int k = 0; k < N; k++) { C[i][j] += convolution(A[i][k], B[k][j]); } } } return A = C; } // ドット積 vector operator*=(vector &f, const matrix &A) { int N = f.size(); vector g(N); for (int i = 0; i < N; i++) { for (int j = 0; j < N; j++) { g[i] += convolution(f[j], A[j][i]); } } return f = g; } // 繰り返し二乗法 matrix pow(matrix A, int n) { int N = A.size(); matrix B(N, vector(N)); for (int i = 0; i < N; i++) B[i][i] = {1}; while (n) { if (n & 1) B *= A; A *= A; n >>= 1; } return B; } int main() { int N, H; cin >> N >> H; vector v = {{2*inv5, 0}, {0, 2*inv5}, {0, inv5}}; v *= pow(P, N-1); F f; for (int i = 0; i < v.size(); i++) f += v[i]; mint ans = 0; for (int i = H; i <= N; i++) ans += f[i]; cout << ans.val() << endl; }