#include <iostream>
using namespace std;

uint64_t DP[5002];

uint64_t mod_pow(uint64_t a, uint64_t e, uint64_t p){
    uint64_t res = 1;
    for(;e>0;e>>=1){
        if(e&1)res=(res*a)%p;
        a=(a*a)%p;
    }
    return res;
}

int main(){
    string S;
    cin>>S;
    uint64_t ans = 1;
    DP[0]=1;
    uint64_t base=1009,mod=1000000007,hl,hr;
    for(int i=0;i<S.size()/2;i++){
        hl=0;hr=0;
        for(int j=i;j>=0;j--){
            hl=(hl*base+S[j])%mod;
            hr=(hr+mod_pow(base,i-j,mod)*S[S.size()-j-1])%mod;
            // cout<<S.substr(j,i-j+1)<<' '<<S.substr(S.size()-i-1,i-j+1)<<' '<<hl<<' '<<hr<<endl;
            if(hl==hr){
                DP[i+1]+=DP[j];
                DP[i+1]%=mod;
            }
        }
        ans+=DP[i+1];
        ans%=mod;
    }
    cout<<ans<<endl;
    return 0;
}