#include #include using namespace std; using namespace atcoder; typedef modint998244353 mint; typedef long long ll; // isprime bool isprime(ll n){ if (n == 1) return false; for (ll i = 2; i * i <= n; i++){ if (n % i == 0) return false; } return true; } // importrandom // don't forget randomjumon // and also isprime ll RandomMod(ll l, ll r){ ll ret = l + rand() % (r - l); while (!isprime(ret)) ret = l + rand() % (r - l); return ret; } ll randrange(ll l, ll r){ return l + rand() % (r - l); } // ----- int main(){ ios_base::sync_with_stdio(false); cin.tie(NULL); srand((unsigned)time(NULL)); // Rolling Hash // don't forget randomjumon and importrandom typedef dynamic_modint<0> mint1; typedef dynamic_modint<1> mint2; int mod1 = RandomMod(7e8, 1e9); int mod2 = RandomMod(7e8, 1e9); mint1::set_mod(mod1); mint2::set_mod(mod2); mint1 b1 = randrange(100, 200); mint2 b2 = randrange(100, 200); int MAX_m = 10000; vector b1l(MAX_m+1); b1l[0] = 1; for (int i=0; i b2l(MAX_m+1); b2l[0] = 1; for (int i=0; i> s; int n = s.size(); vector> dp(n+1, vector(n+1, -INF)); dp[0][0] = INF; mint1 r1 = 0; mint2 r2 = 0; vector s1(n+1); vector s2(n+1); vector t1(n+1); vector t2(n+1); for (int i=0; i=0; i--){ r1 = r1 * b1 + s[i] - 'a' + 1; r2 = r2 * b2 + s[i] - 'a' + 1; t1[i] = r1; t2[i] = r2; } mint1 x1, y1; mint2 x2, y2; for (int i=0; i