#include<bits/stdc++.h>

using namespace std;

const int mod = 1e9 + 7;

short LCP[10001][10001];
int dp[10000];

int main()
{
  string S;
  cin >> S;
  for(int i = S.size() - 1; i >= 0; i--) {
    for(int j = S.size() - 1; j >= 0; j--) {
      LCP[i][j] = S[i] == S[j] ? LCP[i + 1][j + 1] + 1 : 0;
    }
  }
  dp[0] = 1;
  for(int i = 0; i < S.size() / 2; i++) {
    for(int k = i; k >= 0; k--) {
      if(LCP[k][S.size() - i - 1] >= i - k + 1) {
        (dp[i + 1] += dp[k]) %= mod;
      }
    }
  }
  cout << accumulate(dp, dp + S.size(), 0LL) % mod << endl;
}