#include <bits/stdc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/rope>

#define ll long long
#define pb push_back
#define f first
#define s second
#define sz(x) (int)(x).size()
#define all(x) x.begin(), x.end()
#define rall(x) x.rbegin(), x.rend()
#define ios ios_base::sync_with_stdio(false);cin.tie(NULL)
#define ld long double
#define li __int128

using namespace std;
using namespace __gnu_pbds;
using namespace __gnu_cxx;

mt19937 rng(chrono::steady_clock::now().time_since_epoch().count());
template<class T> using ordered_set = tree<T, null_type, less<T>, rb_tree_tag,tree_order_statistics_node_update>; ///find_by_order(),order_of_key()
template<int D, typename T>struct vec : public vector<vec<D - 1, T>> {template<typename... Args>vec(int n = 0, Args... args) : vector<vec<D - 1, T>>(n, vec<D - 1, T>(args...)) {}};
template<typename T>struct vec<1, T> : public vector<T> {vec(int n = 0, T val = T()) : vector<T>(n, val) {}};
template<class T1,class T2> ostream& operator<<(ostream& os, const pair<T1,T2>& a) { os << '{' << a.f << ", " << a.s << '}'; return os; }
template<class T> ostream& operator<<(ostream& os, const vector<T>& a){os << '{';for(int i=0;i<sz(a);i++){if(i>0&&i<sz(a))os << ", ";os << a[i];}os<<'}';return os;}
template<class T> ostream& operator<<(ostream& os, const deque<T>& a){os << '{';for(int i=0;i<sz(a);i++){if(i>0&&i<sz(a))os << ", ";os << a[i];}os<<'}';return os;}
template<class T> ostream& operator<<(ostream& os, const set<T>& a) {os << '{';int i=0;for(auto p:a){if(i>0&&i<sz(a))os << ", ";os << p;i++;}os << '}';return os;}
template<class T> ostream& operator<<(ostream& os, const set<T,greater<T> >& a) {os << '{';int i=0;for(auto p:a){if(i>0&&i<sz(a))os << ", ";os << p;i++;}os << '}';return os;}
template<class T> ostream& operator<<(ostream& os, const multiset<T>& a) {os << '{';int i=0;for(auto p:a){if(i>0&&i<sz(a))os << ", ";os << p;i++;}os << '}';return os;}
template<class T> ostream& operator<<(ostream& os, const multiset<T,greater<T> >& a) {os << '{';int i=0;for(auto p:a){if(i>0&&i<sz(a))os << ", ";os << p;i++;}os << '}';return os;}
template<class T1,class T2> ostream& operator<<(ostream& os, const map<T1,T2>& a) {os << '{';int i=0;for(auto p:a){if(i>0&&i<sz(a))os << ", ";os << p;i++;}os << '}';return os;}
int ri(){int x;scanf("%i",&x);return x;}
void rd(int&x){scanf("%i",&x);}
void rd(long long&x){scanf("%lld",&x);}
void rd(double&x){scanf("%lf",&x);}
void rd(long double&x){scanf("%Lf",&x);}
void rd(string&x){cin>>x;}
void rd(char*x){scanf("%s",x);}
template<typename T1,typename T2>void rd(pair<T1,T2>&x){rd(x.first);rd(x.second);}
template<typename T>void rd(vector<T>&x){for(T&p:x)rd(p);}
template<typename C,typename...T>void rd(C&a,T&...args){rd(a);rd(args...);}
//istream& operator>>(istream& is,__int128& a){string s;is>>s;a=0;int i=0;bool neg=false;if(s[0]=='-')neg=true,i++;for(;i<s.size();i++)a=a*10+s[i]-'0';if(neg)a*=-1;return is;}
//ostream& operator<<(ostream& os,__int128 a){bool neg=false;if(a<0)neg=true,a*=-1;ll high=(a/(__int128)1e18);ll low=(a-(__int128)1e18*high);string res;if(neg)res+='-';if(high>0){res+=to_string(high);string temp=to_string(low);res+=string(18-temp.size(),'0');res+=temp;}else res+=to_string(low);os<<res;return os;}
const int mod=1e9+7; //998244353
int add(int a,int b){
    a+=b;
    if(a>=mod)
        a-=mod;
    return a;
}
int sub(int a,int b){
    a-=b;
    if(a<0)
        a+=mod;
    return a;
}
int mul(int a,int b){
    return (long long)a*b%mod;
}
int pwrmod(int x,int k){
    int ans=1;
    for(;k;k>>=1,x=mul(x,x))
        if(k&1)
            ans=mul(ans,x);
    return ans;
}
const int D=29;
const int N=805;
vector<int> st;
vector<vector<int>> st2;
int main()
{
	st.pb(1);
	for(int i=0;i<N;i++){
		st.pb(mul(st.back(),D));
	}
	for(int i=0;i<N;i++){
		st2.pb({});
		for(int j=0;j<D;j++){
			st2.back().pb(mul(st[i],j));
		}
	}
	string s;
	cin >> s;
	int n=sz(s);
	vector<gp_hash_table<int,int>> cnt(n+1);
	ll total=0;
	for(int i=0;i<n;i++){
		for(int j=i;j<n;j++){
			int mxLen=(j-i+1);
			int st=(mxLen&1)?1:2;
			vector<bool> ok={1};
			for(int k=0;2*k+2<mxLen;k++){
				if(!ok.back()){
					ok.pb(0);
				}
				else
					ok.pb(s[i+k]==s[j-k]);
			}
			int md=i+(mxLen-1)/2;
			int tr=0;
			for(int k=0;k<st;k++){
				tr=add(mul(tr,D),s[md+k]-'a');
			}
			int ind=sz(ok)-1;
			int oo=1;
			for(int x=st;x<=mxLen;x+=2){
				if(ok[ind]){
					//printf("%i - %i  %i: %i\n",i,j,x,cnt[x][tr]);
					total+=cnt[x][tr];
				}
				ind--;
				tr=add(tr,st2[x][md-oo]-'a');
				tr=add(mul(tr,D),s[md+oo+st-1]-'a');
				oo++;
			}
		}
		int tr=0;
		for(int j=i;j>=0;j--){
			tr=add(tr,st2[i-j][s[j]-'a']);
			cnt[i-j+1][tr]++;
		}
	}
	for(int i=0;i<=n;i++){
		cnt[i].clear();
	}
	reverse(all(s));
	for(int i=0;i<n;i++){
		for(int j=i;j<n;j++){
			int mxLen=(j-i+1);
			int st=(mxLen&1)?1:2;
			vector<bool> ok={1};
			for(int k=0;2*k+2<mxLen;k++){
				if(!ok.back()){
					ok.pb(0);
				}
				else
					ok.pb(s[i+k]==s[j-k]);
			}
			int md=i+(mxLen-1)/2;
			int tr=0;
			for(int k=0;k<st;k++){
				tr=add(mul(tr,D),s[md+k]-'a');
			}
			int ind=sz(ok)-1;
			int oo=1;
			for(int x=st;x<mxLen;x+=2){
				if(ok[ind]){
					//printf("2 %i - %i  %i: %i\n",i,j,x,cnt[x][tr]);
					total+=cnt[x][tr];
				}
				ind--;
				tr=add(tr,st2[x][md-oo]-'a');
				tr=add(mul(tr,D),s[md+oo+st-1]-'a');
				oo++;
			}
		}
		int tr=0;
		for(int j=i;j>=0;j--){
			tr=add(tr,st2[i-j][s[j]-'a']);
			cnt[i-j+1][tr]++;
		}
	}
	printf("%lld\n",total);
	return 0;
}