#include <bits/stdc++.h>
using namespace std;
template <class T>
bool umin(T& a, const T& b) {
  return a > b ? a = b, true : false;
}
template <class T>
bool umax(T& a, const T& b) {
  return a < b ? a = b, true : false;
}
template <int sz>
using tut = array<int, sz>;
void usaco(string s) {
  freopen((s + ".in").c_str(), "r", stdin);
  freopen((s + ".out").c_str(), "w", stdout);
  ios::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
}
const int N = 3e3 + 5;
const int mod = (119 << 23) + 1;
const long long inf = 1e18;
const long double Pi = acos(-1);
void add(int& a, int b, int mod = ::mod) { a = (a + b) % mod; }
void sub(int& a, int b, int mod = ::mod) { a = ((a - b) % mod + mod) % mod; }
void mul(int& a, int b, int mod = ::mod) { a = (a * 1ll * b) % mod; }
int n, m, k, t, q, ans, res, a[N];
int dp[N][N];
void solve(int t_case) {
  string s, t;
  cin >> s >> t;
  n = (int)s.size(), s = "#" + s;
  m = (int)t.size(), t = "#" + t;
  for (int j = 1; j <= m; j++)
    if (s[1] == t[j]) dp[1][j] = 2;
  dp[1][m + 1] = 2;
  for (int i = 2; i <= n; i++) {
    for (int j = 1; j <= m + 1; j++) {
      if (j == m + 1) {
        add(dp[i][j], dp[i - 1][j]);
      }
      if (j > 1 && t[j - 1] == s[i]) {
        add(dp[i][j - 1], dp[i - 1][j]);
      }
      if (j + i - 1 <= m && t[j + i - 1] == s[i]) {
        add(dp[i][j], dp[i - 1][j]);
      } else if (j + i - 1 > m) {
        add(dp[i][j], dp[i - 1][j]);
      }
    }
  }
  for (int i = m; i <= n; i++) add(res, dp[i][1]);
  cout << res << "\n";
}
signed main() {
  ios::sync_with_stdio(0);
  cin.tie(0);
  cout.tie(0);
  if (0) {
    int t;
    cin >> t;
    for (int t_case = 1; t_case <= t; t_case++) solve(t_case);
  } else
    solve(1);
  return 0;
}
