#include <bits/stdc++.h>
#define MOD 1000000007LL
using namespace std;
typedef long long ll;
using namespace std;
int n;
string s;
vector<int> vec[26];
ll dp[2005][2005];

int main(void){
	cin >> s;
	n=s.size();
	for(int i=0;i<n;i++){
		int a=s[i]-'a';
		vec[a].push_back(i);
	}
	memset(dp,0,sizeof(dp));
	ll ans=0;
	dp[0][n]=1;
	for(int i=0;i<n;i++){
		for(int j=n;j>i;j--){
			if(dp[i][j]==0)continue;
			for(int k=0;k<26;k++){
				int v=lower_bound(vec[k].begin(),vec[k].end(),i)-vec[k].begin();
				if(v==vec[k].size())continue;
				int f=vec[k][v];
				//printf("%d %d %d %d\n",i,j,f,k);
				if(f>=j)continue;
				//printf("%d %d %d %d\n",i,j,f,k);
				ans+=dp[i][j];
				ans%=MOD;
				int v2=lower_bound(vec[k].begin(),vec[k].end(),j)-vec[k].begin();
				v2--;
				if(v2<0)continue;
				//printf("%d %d %d %d\n",i,j,k,v,v2);
				if(vec[k][v2]>f){
					dp[f+1][vec[k][v2]]+=dp[i][j];
					dp[f+1][vec[k][v2]]%=MOD;
				}
			}
		}
	}
	for(int i=1;i<n;i++){
		for(int j=0;j<=n;j++){
			ans+=dp[i][j];
			ans%=MOD;
		}
	}
	printf("%lld\n",ans);
	return 0;
}
