#include #include #include using namespace std; struct Manachar { vector rads; // O(|s|) Manachar(const string &s) : rads(s.size() * 2 - 1){ int size = (int) rads.size(); int i = 0, j = 0; while (i < size) { while (i - j >= 0 && i + j < size && get(s, i - j) == get(s, i + j)) { ++j; } rads[i] = j; int k = 1; while (i - k >= 0 && i + k < size && k + rads[i - k] < j) { rads[i + k] = rads[i - k], ++k; } i += k; j -= k; } } // s[l, r] is palindrome? // O(1) bool is_palindrome(int l, int r) { assert(l >= 0); assert(r >= l); assert(r * 2 <= (int) rads.size()); return rads[l + r] >= r - l + 1; } private: static char get(const string &s, int i) { if (i & 1) return '^'; else return s[i >> 1]; } }; int main() { cin.tie(nullptr); ios::sync_with_stdio(false); string s; cin >> s; int n = (int) s.size(); assert(n >= 4); assert(n <= 5000); for (int i = 0; i < n; i++) { assert('a' <= s[i] && s[i] <= 'z'); } Manachar mana(s); vector p1; for (int i = 0; i < n - 3; i++) { if (mana.is_palindrome(0, i)) { p1.push_back(i); } } vector p2(n - 1); for (int i : p1) { for (int j = i + 1; j < n - 2; j++) { if (mana.is_palindrome(i + 1, j)) { p2[j]++; } } } vector sum_p2(n); for (int i = 0; i < n - 1; i++) sum_p2[i + 1] = sum_p2[i] + p2[i]; long long ans = 0; for (int k = 2; k < n - 1; k++) { if (mana.is_palindrome(k + 1, n - 1)) { ans += sum_p2[k]; } } cout << ans << endl; }