#pragma GCC optimize ("O3") #pragma GCC target ("avx") #include "bits/stdc++.h" // define macro "/D__MAI" using namespace std; typedef long long int ll; #define xprintf(fmt,...) fprintf(stderr,fmt,__VA_ARGS__) #define debugv(v) {printf("L%d %s > ",__LINE__,#v);for(auto e:v){cout< ",__LINE__,#m);for(int x=0;x<(w);x++){cout<<(m)[x]<<" ";}cout<\n",__LINE__,#m);for(int y=0;y<(h);y++){for(int x=0;x<(w);x++){cout<<(m)[y][x]<<" ";}cout< ostream& operator <<(ostream &o, const pair p) { o << "(" << p.first << ":" << p.second << ")"; return o; } template inline size_t argmin(iterator begin, iterator end) { return distance(begin, min_element(begin, end)); } template inline size_t argmax(iterator begin, iterator end) { return distance(begin, max_element(begin, end)); } template T& maxset(T& to, const T& val) { return to = max(to, val); } template T& minset(T& to, const T& val) { return to = min(to, val); } void bye(string s, int code = 0) { cout << s << endl; exit(0); } mt19937_64 randdev(8901016); inline ll rand_range(ll l, ll h) { return uniform_int_distribution(l, h)(randdev); } #ifdef __MAI #define getchar_unlocked getchar #define putchar_unlocked putchar #endif #ifdef __VSCC #define getchar_unlocked _getchar_nolock #define putchar_unlocked _putchar_nolock #endif namespace { #define isvisiblechar(c) (0x21<=(c)&&(c)<=0x7E) class MaiScanner { public: template void input_integer(T& var) { var = 0; T sign = 1; int cc = getchar_unlocked(); for (; cc<'0' || '9'>(int& var) { input_integer(var); return *this; } inline MaiScanner& operator>>(long long& var) { input_integer(var); return *this; } inline MaiScanner& operator>>(string& var) { int cc = getchar_unlocked(); for (; !isvisiblechar(cc); cc = getchar_unlocked()); for (; isvisiblechar(cc); cc = getchar_unlocked()) var.push_back(cc); return *this; } template void in(IT begin, IT end) { for (auto it = begin; it != end; ++it) *this >> *it; } }; } MaiScanner scanner; namespace std { template class hash> { public: size_t operator()(const pair& x) const { return hash()(x.first) ^ hash()(x.second); } }; } unordered_map, int> memo; namespace rare { int b_search_low(int low, int high, auto compare) { register int l, h, m; l = low; h = high; while (l < h) { m = (l + h) / 2; if (0 sa; SuffixArray() {} SuffixArray(const string &d) { data = d.c_str(); size = d.size(); } // Suffix Array 生成 O(nlogn) void build() { int i; sa.resize(size); for (i = 0; istrcmp(l, r); }); // (´・ω・`) vector g(size + 1), b(size + 1); vector v(size + 1); repeat(i, size + 1) { v[i] = i, g[i] = data[i]; } b[0] = 0; b[size] = 0; sort(v.begin(), v.end(), [&](int l, int r) {return l == r ? false : g[l] != g[r] ? g[l] < g[r] : g[l] < g[r]; }); for (int h = 1; b[size] != size; h *= 2) { sort(v.begin(), v.end(), [&](int l, int r) {return l == r ? false : g[l] != g[r] ? g[l] < g[r] : g[l + h] < g[r + h]; }); repeat(i, size) b[i + 1] = b[i] + [&](int l, int r) {return l == r ? false : g[l] != g[r] ? g[l] < g[r] : g[l + h] < g[r + h]; }(v[i], v[i + 1]); repeat(i, size + 1) g[v[i]] = b[i]; } repeat(i, size) { sa[i] = data + v[i]; } return; } void build(const string &d) { data = d.c_str(); size = d.size(); build(); } // find keyword // out_low <= out_highの範囲にkeyword文字列が存在する.out_high-out_low+1個見つかった,とも解釈できる. // foundToIdx(out_low..out_high)とすることで,見つかった文字列の先頭インデックスを得る. // O(mlogn) void find(const string &keyword, int &out_low, int &out_high) { int low, high, l, h, i; auto length = keyword.size(); const char* c_keyword = keyword.c_str(); low = 0; high = sa.size(); for (i = 0; i < length; i++) { l = rare::b_search_low(low, high, [&](auto p) {return *(c_keyword + i) - *(sa[p] + i); }); h = rare::b_search_high(low, high, [&](auto p) {return *(c_keyword + i) - *(sa[p] + i); }); if (l == -1 || h == -1 || h revsa; inline bool my_same(int left, int right, int k) { for (int l = 0; l < k; ++l) { if (str[left + l] != str[right - k + l]) return false; } return true; } ll dfs(int left, int right) { //cout << make_pair(left,right) << endl; if (right - left <= 1) return 1; if (memo.count(make_pair(right, left))) return memo[make_pair(right, left)]; // + [left,right]を1つのブロック // + [left,left+k],[right-k,right] で1つのpair ll ans = 1; // (´・ω・`) for (int k = 1; left + k <= right - k; ++k) { //cout << " " << make_pair(left + k, right - k) << endl; if (my_same(left, right, k)) ans = (ans + dfs(left + k, right - k)) % MD; } //cout << make_pair(left,right)<> str; n = str.size(); sa.build(str); { revsa.resize(n); repeat(i, n) { revsa[(sa.sa[i] - sa.data)] = i; } } ll ans = dfs(0, n); cout << ans << endl; return 0; }