#include using namespace std; #define FOR(i,n) for(int i = 0; i < (n); i++) #define sz(c) ((int)c.size()) #define ten(n) ((int)1e##n) using ll = long long; template static inline int getchar_unlocked(void) { return getchar(); } template static inline void putchar_unlocked(int c) { putchar(c); } #define mygc(c) (c)=getchar_unlocked() #define mypc(c) putchar_unlocked(c) void reader(int& x) { int k, m = 0; x = 0; for (;;) { mygc(k); if (k == '-') { m = 1; break; }if ('0' <= k&&k <= '9') { x = k - '0'; break; } }for (;;) { mygc(k); if (k<'0' || k>'9')break; x = x * 10 + k - '0'; }if (m) x = -x; } void reader(ll& x) { int k, m = 0; x = 0; for (;;) { mygc(k); if (k == '-') { m = 1; break; }if ('0' <= k&&k <= '9') { x = k - '0'; break; } }for (;;) { mygc(k); if (k<'0' || k>'9')break; x = x * 10 + k - '0'; }if (m) x = -x; } int reader(char c[]) { int i, s = 0; for (;;) { mygc(i); if (i != ' '&&i != '\n'&&i != '\r'&&i != '\t'&&i != EOF) break; }c[s++] = i; for (;;) { mygc(i); if (i == ' ' || i == '\n' || i == '\r' || i == '\t' || i == EOF) break; c[s++] = i; }c[s] = '\0'; return s; } int reader(string& c) { int i; for (;;) { mygc(i); if (i != ' '&&i != '\n'&&i != '\r'&&i != '\t'&&i != EOF) break; }c.push_back(i); for (;;) { mygc(i); if (i == ' ' || i == '\n' || i == '\r' || i == '\t' || i == EOF) break; c.push_back(i); }; return sz(c); } template void reader(T& x, S& y) { reader(x); reader(y); } template void reader(T& x, S& y, U& z) { reader(x); reader(y); reader(z); } template void reader(T& x, S& y, U& z, V & w) { reader(x); reader(y); reader(z); reader(w); } void writer(int x, char c) { int s = 0, m = 0; char f[10]; if (x<0)m = 1, x = -x; while (x)f[s++] = x % 10, x /= 10; if (!s)f[s++] = 0; if (m)mypc('-'); while (s--)mypc(f[s] + '0'); mypc(c); } void writer(ll x, char c) { int s = 0, m = 0; char f[20]; if (x<0)m = 1, x = -x; while (x)f[s++] = x % 10, x /= 10; if (!s)f[s++] = 0; if (m)mypc('-'); while (s--)mypc(f[s] + '0'); mypc(c); } void writer(const char c[]) { int i; for (i = 0; c[i] != '\0'; i++)mypc(c[i]); } void writer(const char x[], char c) { int i; for (i = 0; x[i] != '\0'; i++)mypc(x[i]); mypc(c); } template void writerLn(T x) { writer(x, '\n'); } template void writerLn(T x, S y) { writer(x, ' '); writer(y, '\n'); } template void writerLn(T x, S y, U z) { writer(x, ' '); writer(y, ' '); writer(z, '\n'); } template void writerArr(T x[], int n) { if (!n) { mypc('\n'); return; }FOR(i, n - 1)writer(x[i], ' '); writer(x[n - 1], '\n'); } template void writerArr(vector& x) { writerArr(x.data(), (int)x.size()); } template T gcd(T a, T b) { return b ? gcd(b, a % b) : a; } template T lcm(T a, T b) { return a / gcd(a, b) * b; } ll mod_pow(ll a, ll n, ll mod) { ll ret = 1; ll p = a % mod; while (n) { if (n & 1) ret = ret * p % mod; p = p * p % mod; n >>= 1; } return ret; } template T extgcd(T a, T b, T& x, T& y) { for (T u = y = 1, v = x = 0; a;) { T q = b / a; swap(x -= q * u, u); swap(y -= q * v, v); swap(b -= q * a, a); } return b; } template T mod_inv(T a, T m) { T x, y; extgcd(a, m, x, y); return (m + x % m) % m; } using Pii = pair; using Pll = pair; struct PalindromicTree { // // private: struct node { map link; int suffix_link; int len; int count; }; vector c; string s; int active_idx; node* create_node() { c.emplace_back(); node* ret = &c.back(); ret->count = 0; return ret; } // this->s の状態に依存する int find_prev_palindrome_idx(int node_id) { const int pos = int(s.size()) - 1; while (true) { const int opposite_side_idx = pos - 1 - c[node_id].len; if (opposite_side_idx >= 0 && s[opposite_side_idx] == s.back()) break; node_id = c[node_id].suffix_link; // 次の回文に移動 } return node_id; } bool debug_id2string_dfs(int v, int id, vector& charas) { if (v == id) return true; for (auto kv : c[v].link) { if (debug_id2string_dfs(kv.second, id, charas)) { charas.push_back(kv.first); return true; } } return false; } public: PalindromicTree() { node* size_m1 = create_node(); // 長さ-1のノードを作成 size_m1->suffix_link = 0; // -1 の親suffixは自分自身 size_m1->len = -1; node* size_0 = create_node(); // 長さ0のノードを作成 size_0->suffix_link = 0; // 親は長さ-1のノード size_0->len = 0; active_idx = 0; } int get_active_idx() const { return active_idx; } node* get_node(int id) { return &c[id]; } void add(char ch) { s.push_back(ch); // ch + [A] + ch が回文となるものを探す const int a = find_prev_palindrome_idx(active_idx); //新しいノードへのリンクが発生するか試す const auto inserted_result = c[a].link.insert(make_pair(ch, int(c.size()))); active_idx = inserted_result.first->second; // insertの成否に関わらず、iteratorが指す先は新しい回文のindex if (!inserted_result.second) { c[active_idx].count++; // その回文が現れた回数が増加 return; // 既にリンクが存在したので、新しいノードを作る必要がない } // 実際に新しいノードを作成 node* nnode = create_node(); nnode->count = 1; nnode->len = c[a].len + 2; // ch + [A] + ch だから、長さは len(A) + 2 // suffix_linkの設定 if (nnode->len == 1) { // この時だけsuffix_linkはsize 0に伸ばす nnode->suffix_link = 1; } else { // ch + [B] + ch が回文になるものを探す。ただし長さはaより小さいもの const int b = find_prev_palindrome_idx(c[a].suffix_link); nnode->suffix_link = c[b].link[ch]; } } //各文字列が何回現れるか計算する // O(n) vector build_frequency() { vector frequency(c.size()); //常に親ノードのid < 子ノードのidが成り立つので、idを大きい順から回せばよい for (int i = int(c.size()) - 1; i > 0; i--) { frequency[i] += c[i].count; frequency[c[i].suffix_link] += frequency[i]; } return frequency; } // debug用 // idがどのような回文を表しているのかを返す // O(n) string debug_id2string(int id) { if (id == 0) { return "(-1)"; } else if (id == 1) { return "(0)"; } vector charas; debug_id2string_dfs(0, id, charas); debug_id2string_dfs(1, id, charas); string ret(charas.begin(),charas.end()); int start = int(charas.size()) - 1; if (c[id].len % 2 == 1) start--; for (int i = start; i >= 0; i--) ret.push_back(charas[i]); return ret; } void display_frequencies() { auto freq = build_frequency(); printf("frequencies of each palindrome...\n"); for (int i = 0; i < int(c.size()); i++) { const string palindrome = debug_id2string(i); printf(" %s : %d\n", palindrome.c_str(), freq[i]); } } }; int main() { string s, t; reader(s, t); PalindromicTree pt; for (auto c : s) pt.add(c); auto dp = pt.build_frequency(); FOR(i, sz(pt.c)) pt.c[i].count = 0; pt.add('!'); pt.add('?'); for (auto c : t) pt.add(c); auto dp2 = pt.build_frequency(); ll ans = 0; for(int i = 2; i < sz(dp); i++) ans += dp[i] * dp2[i]; writerLn(ans); }