#include <bits/stdc++.h>

using namespace std;

#ifdef LOCAL
template <class... T> void dbg(T... x) { char e{}; ((cerr << e << x, e = ' '), ...); }
template <class T> void org(T l, T r) { while (l != r) cerr << ' ' << *l++; cerr << '\n'; }
#define debug(x...) dbg(#x, '=', x, '\n')
#define olist(x...) dbg(#x, '='), org(x)
#else
// #pragma GCC optimize("O3,unroll-loops")
// #pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt")
#define debug(...) ((void)0)
#define olist(...) ((void)0)
#endif
#define all(v) (v).begin(), (v).end()
#define rall(v) (v).rbegin(), (v).rend()
#define ff first
#define ss second

using u32 = unsigned int;
using i64 = long long;
using u64 = unsigned long long;
using i128 = __int128;
using u128 = unsigned __int128;
using Pt = pair<int, int>;

template <class T> inline constexpr T inf = numeric_limits<T>::max() / 2;
constexpr int mod = 998244353;

template<class T> bool chmin(T &a, T b) { return (b < a and (a = b, true)); }
template<class T> bool chmax(T &a, T b) { return (a < b and (a = b, true)); }
template<class... T> int add(T... x) { int t{}; return (((t += x) %= mod), ...), t; }
template<class... T> int mul(T... x) { i64 t{1}; return (((t *= x) %= mod), ...), t; }

#include <bits/extc++.h>
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>
#include <ext/pb_ds/hash_policy.hpp>
#include <ext/pb_ds/priority_queue.hpp>
using namespace __gnu_pbds;
template<class T>
using BST = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>;

constexpr u64 P = 131;
int fail[805];

void solve() {
    string S;
    cin >> S;
    
    const int n = S.size();
    vector<u64> pw(n + 1, 1);
    for (int i = 1; i <= n; i++) {
        pw[i] = pw[i - 1] * P;
    }

    auto build = [&](string_view) -> vector<u64> {
        vector<u64> h(n);
        for (int i = 0; i < n; i++) {
            h[i] = S[i] * pw[i];
            if (i) h[i] += h[i - 1];
        }
        return h;
    };

    auto cal = [&](string S, bool same) -> i64 {
        auto h = build(S);
        i64 ans = 0;

        auto get = [&](int l, int r) -> u64 {
            return (h[r] - (l == 0 ? 0 : h[l - 1])) * pw[n - l];
        };
        
        unordered_map<u64, i64> cnt;
        for (int l = 1; l < n; l++) {
            for (int i = 0; i < l; i++) {
                cnt[get(i, l - 1)]++;
            }
            
            int p = -1;
            for (int r = l; r < n; r++) {
                while (p != -1 and S[l + p + 1] != S[r]) {
                    p = fail[p];
                }
                if (S[l + p + 1] == S[r]) p++;
                if (r == l) p = -1;
                fail[r - l] = p;

                for (int g = p; g != -1; g = fail[g]) {
                    if (l + g + 1 <= r - g - 1) {
                        ans += cnt[get(l + g + 1, r - g - 1)];
                    }
                }

                if (same) {
                    ans += cnt[get(l, r)];
                }
            }
        }
        
        return ans;
    };

    i64 ans = cal(S, 1);
    reverse(all(S));
    ans += cal(S, 0);

    cout << ans << '\n';
}

signed main() {
    cin.tie(0)->sync_with_stdio(false);
    cin.exceptions(cin.failbit);
    int T = 1;
    // cin >> T;
    while (T--) {
        solve();
    }
    return 0;
}