#include <bits/stdc++.h>
using namespace std;
const int MOD = 1000000007;
int g[1000010], f[1000010], trans[1000010][26];
int fail[1000010], len[1000010];
char st[1000010];
int c[1000010];
int d[1000010], pos[1000010];
int n, lst, cnt;
inline int mo(int x) {
  if (x >= MOD) return x - MOD;
  return x;
}
int getfail(int i, int x) {
  while (c[i - len[x] - 1] != c[i]) x = fail[x];
  return x;
}
int main() {
  scanf("%s", st + 1);
  n = strlen(st + 1);
  for (int i = 1; i <= n; i++)
    c[i] = (i & 1 ? st[(i + 1) >> 1] : st[n - (i >> 1) + 1]) - 'a';
  len[1] = -1;
  fail[0] = cnt = 1;
  f[0] = 1;
  c[0] = -1;
  for (int i = 1; i <= n; i++) {
    int x = getfail(i, lst);
    if (!trans[x][c[i]]) {
      int v = getfail(i, fail[x]);
      fail[++cnt] = trans[v][c[i]];
      trans[x][c[i]] = cnt;
      len[cnt] = len[x] + 2;
      d[cnt] = len[cnt] - len[fail[cnt]];
      pos[cnt] = d[cnt] == d[fail[cnt]] ? pos[fail[cnt]] : fail[cnt];
    }
    lst = trans[x][c[i]];
    for (x = lst; x; x = pos[x]) {
      g[x] = f[i - len[pos[x]] - d[x]];
      if (pos[x] != fail[x]) g[x] = mo(g[x] + g[fail[x]]);
      if (!(i & 1)) f[i] = mo(f[i] + g[x]);
    }
  }
  printf("%d\n", f[n]);
}
