// Original Python code: // n,k=map(int,input().split()) // M=998244353 // dp=[[0]*3 for _ in range(n+1)] // dp[0][0]=1 // for i in range(1,2*n+1): // ndp=[[0]*3 for _ in range(n+1)] // for j in range(n+1): // if j+1<=n: // if j+1==k: // ndp[j+1][1]+=sum(dp[j][:2])%M // elif j+1>k: // ndp[j+1][2]+=sum(dp[j])%M // else: // ndp[j+1][0]+=dp[j][0]%M // ndp[j+1][1]+=dp[j][1]%M // ndp[j+1][2]+=dp[j][2]%M // if j-1>=0: // ndp[j-1][0]+=dp[j][0]%M // ndp[j-1][1]+=dp[j][1]%M // ndp[j-1][2]+=dp[j][2]%M // dp=ndp // print(dp[0][1]%M) #include #include using namespace std; int main() { int n, k; cin >> n >> k; const int M = 998244353; vector> dp(n + 1, vector(3, 0)); dp[0][0] = 1; for (int i = 1; i <= 2 * n; ++i) { vector> ndp(n + 1, vector(3, 0)); for (int j = 0; j <= n; ++j) { if (j + 1 <= n) { if (j + 1 == k) { int s = (dp[j][0] + dp[j][1]) % M; ndp[j + 1][1] = (ndp[j + 1][1] + s) % M; } else if (j + 1 > k) { int s = ((dp[j][0] + dp[j][1]) % M + dp[j][2]) % M; ndp[j + 1][2] = (ndp[j + 1][2] + s) % M; } else { ndp[j + 1][0] = (ndp[j + 1][0] + dp[j][0]) % M; ndp[j + 1][1] = (ndp[j + 1][1] + dp[j][1]) % M; ndp[j + 1][2] = (ndp[j + 1][2] + dp[j][2]) % M; } } if (j - 1 >= 0) { ndp[j - 1][0] = (ndp[j - 1][0] + dp[j][0]) % M; ndp[j - 1][1] = (ndp[j - 1][1] + dp[j][1]) % M; ndp[j - 1][2] = (ndp[j - 1][2] + dp[j][2]) % M; } } dp = move(ndp); } cout << dp[0][1] % M << endl; return 0; }