// Too many mind, no mind.
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <vector>
#include <set>
#include <map>
#include <unordered_set>
#include <unordered_map>
#include <queue>
#include <ctime>
#include <cassert>
#include <complex>
#include <string>
#include <cstring>
#include <chrono>
#include <random>
#include <bitset>
#include <array>
#include <iomanip>

using namespace std;

#define mp make_pair
#define fi first
#define se second
#define pb push_back
#define ALL(x) (x).begin(), (x).end()
#define RALL(x) (x).rbegin(), (x).rend()
#define COMP(x) sort(ALL(x)); x.resize(unique(ALL(x)) - (x).begin())
#define forn(i, n) for (int i = 0; i < (int)(n); ++i)
#define fore(i, a, b) for (int i = (int)(a); i <= (int)(b); ++i)
#define ford(i, n) for (int i = (int)(n) - 1; i >= 0; --i)

using pii = pair <int, int>;
using vi = vector <int>;
using vpi = vector <pii>;
using ll = long long;
using pll = pair<ll, ll>;
using vl = vector<ll>;
using ld = long double;
using vld = vector<ld>;

const int maxn = 810, P=31; 
const ll mod[2] = {(ll)1e9+7, (ll)1e9+9};
ll add(ll a, ll b, int m) { return a + b >= mod[m] ? a + b - mod[m] : a + b; }
ll sub(ll a, ll b, int m) { return a >= b ? a - b : a - b + mod[m]; }
ll mul(ll a, ll b, int m) { return a*b%mod[m]; }
ll val[maxn];
ll pw[2][maxn];
ll lcp[maxn][maxn];
ll h[2][maxn];
int n;
char s[maxn];
void roll() {
	pw[0][0] = pw[1][0] = 1;
	fore(i, 1, maxn-1) {
		forn(j, 2) pw[j][i] = mul(pw[j][i-1], P, j);
	}
	forn(i, n) {
		forn(j, 2) {
			h[j][i] = mul(s[i] - 'a' + 1, pw[j][i], j);
			if (i) h[j][i] = add(h[j][i], h[j][i-1], j);
		}
	}
}
ll get(int l, int r) {
	vl hash(2);
	forn(j, 2) {
		hash[j] = h[j][r];
		if (l) hash[j] = sub(hash[j], h[j][l-1], j);
		hash[j] = mul(hash[j], pw[j][maxn-1-r], j);
	}
	return hash[0] << 32 | hash[1];
}
ll hsh[maxn][maxn];
int main() {
	ios_base::sync_with_stdio(0), cin.tie(0);
	//freopen("input.txt", "r", stdin);
	//freopen("output.txt", "w", stdout);
	cin >> s;
	n = strlen(s);
	roll();
	int ans = 0;
	forn(i, n) fore(j, i, n-1) hsh[i][j] = get(i, j);
	forn(i, n) fore(j, i, n-1) {
		int len=j-i+1;
		for(int k=0; k<n; k++) {
			if(k+len-1<n) val[k] = hsh[k][k+len-1] == hsh[i][j];
			else val[k]=0;
		}

		fore(i, 1, n-1) val[i] += val[i-1];
		if(j+1<n) ans += val[n-1] - val[j]; 

		for (int w=1; w <= min(i, n-j-1); w++) {
			if (get(i-w, i-1) != get(j+1, j+w)) continue;

			if (j+1+w<n) ans += val[n-1] - val[j+w];
			if (i-len-w>=0) ans += val[i-len-w];
		}
	}
	cout << ans << endl;
}