// compile: make data // run: ./data < data.in #include <bits/stdc++.h> using namespace std; #pragma GCC optimize("O3,unroll-loops") #pragma GCC target("avx2,bmi,bmi2,lzcnt,popcnt") #ifdef LOCAL #include <debug/codeforces.h> #define debug(x...) _debug_print(#x, x); #define Debug(x...) _debug_print_format(#x, x); std::ifstream terminal("/dev/tty"); #define PP cerr<<"\033[1;30mpause...\e[0m",terminal.ignore(); #else #define debug(x...) #define Debug(x...) #define PP #endif template<typename...Args> void print_(Args...args){((cout<<args<<" "),...)<<endl;} #define VI vector<int> #define VII vector<vector<int>> #define VIII vector<vector<vector<int>>> #define rep(i,a,b) for(int i=(a);i<(int)(b);++i) #define sz(v) ((int)(v).size()) #define print(...) print_(__VA_ARGS__); #define FIND(a, x) ((find(a.begin(),a.end(),(x))!=a.end())?1:0) #define cmin(x,...) x=min({(x), __VA_ARGS__}) #define cmax(x,...) x=max({(x), __VA_ARGS__}) #define INTMAX (int)(9223372036854775807) #define INF (int)(1152921504606846976) #define NaN (int)(0x8b88e1d0595d51d1) #define double long double #define int long long #define uint unsigned int #define MAXN 200010 #define P 1000000007 #define MOD0 11635419037240549 #define MOD1 10487639585196761 #define MOD2 10671969493879289 namespace math { constexpr int safe_mod(int x, int m) { x %= m; if (x < 0) x += m; return x; } struct barrett { uint _m, im; explicit barrett(uint m) : _m(m), im((uint)(-1) / m + 1) {} uint umod() const { return _m; } uint mul(unsigned int a, unsigned int b) const { uint z = a; z *= b; uint x = (uint)(((unsigned __int128)(z)*im) >> 64); uint y = x * _m; return (uint)(z - y + (z < y ? _m : 0)); } }; int pow_mod(int x, int n, int m = INF) { if (m == 1) return 0; barrett bt((uint)(m)); uint r = 1, y = (uint)(safe_mod(x, m)); while (n) { if (n & 1) r = bt.mul(r, y); y = bt.mul(y, y); n >>= 1; } return r; } } template <char F, int B, int N> struct strhash { const int M[10] = {999999937, 999999929, 999999893, 999999883, 999999797, 999999761, 999999757, 999999739, 999999733, 999999677}; string s; int len; vector<vector<int>> pw, hs; strhash(string S): s(S), len(S.size()) { assert(N <= 10); pw.assign(len, vector<int>(N)); fill(pw[0].begin(), pw[0].end(), 1); rep(i, 1, len) rep(j, 0, N) pw[i][j] = pw[i-1][j] * B % M[j]; hs.assign(len+1, vector<int>(N, 0)); rep(i, 0, len) rep(j, 0, N) hs[i+1][j] = (hs[i][j] * B + s[i] - F) % M[j]; } vector<int> hash(int l, int r) { if (!l) return hs[r+1]; vector<int> res(N); rep(i, 0, N) res[i] = (hs[r+1][i] - hs[l][i] * pw[r-l+1][i] % M[i] + M[i]) % M[i]; return res; } }; // 11813125 int32_t main() { ios::sync_with_stdio(false); cin.tie(nullptr); cout.tie(nullptr); int n; cin >> n; string s; cin >> s; int m = s.length(); strhash<'A', 26, 2> hs(s); auto getlen = [&](int len, char c) -> int { string cur = s.substr(0, len) + c; int maxlen = cur.length(); strhash<'A', 26, 2> hcur(cur); while (cur.substr(cur.length() - maxlen, maxlen) != s.substr(0, maxlen)) --maxlen; return maxlen; }; vector<vector<int>> nxt(m+1, vector<int>(26, 0)); rep(i, 0, m+1) rep(c, 0, 26) nxt[i][c] = getlen(i, (char)(c + 'A')); vector<vector<int>> dp(n+1, vector<int>(m+1, 0)); dp[0][0] = 1; rep(i, 0, n) { rep(j, 0, min(i+1, m)) { rep(c, 0, 26) { int len = nxt[j][c]; dp[i+1][len] += dp[i][j]; dp[i+1][len] %= P; } } } int sum = math::pow_mod(26, n, P), ans = 0; rep(i, 0, m) ans = (ans + dp[n][i]) % P; ans = (sum - ans + P) % P; cout << ans << endl; return 0; }
No Comments