#include <iostream>
#include <stdio.h>
#include <cstring>
#include <math.h>
#include <algorithm>
#include <vector>
#include <string>
#include <stdlib.h>
#include <queue>
#include <stack>
#include <utility>
#include <fstream>
#include <random>
#include <map>
#include <unordered_map>
#define rep(i,n) for(int i=0;i<n;i++)
#define rrep(i,n) for(int i=n-1;i>=0;i--)
#define FOR(i,a,b) for(int i=a;i<b;i++)
#define ll long long
#define INF 1000000001
#define mod 1000000007
#define p pair<int,int>
int dx[4]={0,1,0,-1};
int dy[4]={1,0,-1,0};
using namespace std;
int POW(int x,int y){return int(pow(double(x),double(y)));}

vector<vector<int>> a;
string s;
int f(int x,int y){
	if(a[x][y]>0){
		return a[x][y];
	}
	if(x==y){
		return 1;
	}
	a[x][y]=1;
	rep(i,(y-x+1)/2){
		if(s.substr(x,i+1) == s.substr(y-i,i+1)){
			if((x+i+1)-(y-i-1)==1){
				a[x][y]++;
				a[x][y]=a[x][y]%mod;
			}else{
				a[x][y]+=f(x+i+1,y-i-1);
				a[x][y]=a[x][y]%mod;
			}
		}
	}
	return a[x][y];
}




int main(){
	
	
	cin>>s;
	int ss=s.size();
	a.resize(ss);
	rep(i,ss){
		a[i].resize(ss);
		rep(j,ss){
			a[i][j]=0;
		}
	}
	cout<<f(0,ss-1)<<endl;

}