#!/usr/bin/env python3.8 import sys read = sys.stdin.buffer.read readline = sys.stdin.buffer.readline readlines = sys.stdin.buffer.readlines import numpy as np S = np.array(list(read().rstrip()), np.int8) N = len(S) is_palind = np.zeros((N + 1, N + 1), np.bool_) is_palind[0] = 1 is_palind[1, :-1] = 1 for n in range(2, N + 1): is_palind[n, : N + 1 - n] = is_palind[n - 2, 1: N + 2 - n] is_palind[n, : N + 1 - n] &= (S[:N + 1 - n] == S[n - 1:N]) P = np.zeros(N, np.int64) P += is_palind[1:, 0] * 1 PP = np.zeros_like(P) for i, x in enumerate(P): PP[i + 1:] += x * is_palind[1: N - i, i + 1] PPA = np.zeros_like(PP) PPA[1:] = np.cumsum(PP[:-1]) PPAP = np.zeros_like(P) for i, x in enumerate(PPA): PPAP[i + 1:] += x * is_palind[1: N - i, i + 1] print(PPAP[-1])