#include using namespace std; struct Manacher { vector rads; // O(|s|) Manacher(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 = s.size(); assert(n >= 4 && n <= 500000); for (char c : s) assert(c >= 'a' && c <= 'z'); Manacher mana(s); vector pl(n); vector gpl(n + 1); vector> G; for (int j = 1; j <= n; j++) { // G -> G' vector> H; for (const tuple &idk : G) { int i, d, k; tie(i, d, k) = idk; if (i > 1 && s[i - 2] == s[j - 1]) { H.emplace_back(i - 1, d, k); } } // G' -> G'' deque> I; int r = -j; for (const tuple &idk : H) { int i, d, k; tie(i, d, k) = idk; if (i - r != d) { I.emplace_back(i, i - r, 1); if (k > 1) { I.emplace_back(i + d, d, k - 1); } } else { I.emplace_back(idk); } r = i + (k - 1) * d; } if (j > 1 && s[j - 2] == s[j - 1]) { I.emplace_back(j - 1, j - 1 - r, 1); r = j; } I.emplace_back(j, j - r, 1); G.clear(); // G'' -> next G int ip, dp, kp; tie(ip, dp, kp) = I.front(); I.pop_front(); for (const tuple &idk : I) { int i, d, k; tie(i, d, k) = idk; if (dp == d) { kp += k; } else { G.emplace_back(ip, dp, kp); ip = i; dp = d; kp = k; } } G.emplace_back(ip, dp, kp); // calc pl for (const tuple &idk : G) { int i, d, k; tie(i, d, k) = idk; r = i + (k - 1) * d; int m = 0; if (r - 2 >= 0) m += mana.is_palindrome(0, r - 2); if (k > 1 && i - d >= 0) { m += gpl[i - d]; } if (d <= i) { gpl[i - d] = m; } pl[j - 1] += m; } } vector sum_p2(n + 1); for (int i = 1; i < n; i++) sum_p2[i + 1] = sum_p2[i] + pl[i]; long long ans = 0; for (int k = 0; k < n - 1; k++) { if (mana.is_palindrome(k + 1, n - 1)) { ans += sum_p2[k]; } } cout << ans << endl; }