#include // created [2019/12/24] 23:17:25 #pragma GCC diagnostic ignored "-Wsign-compare" #pragma GCC diagnostic ignored "-Wsign-conversion" using i32 = int32_t; using i64 = int64_t; using u32 = uint32_t; using u64 = uint64_t; using uint = unsigned int; using usize = std::size_t; using ll = long long; using ull = unsigned long long; using ld = long double; template using arr = T (&)[n]; template using c_arr = const T (&)[n]; template constexpr T popcount(const T u) { return u ? static_cast(__builtin_popcountll(static_cast(u))) : static_cast(0); } template constexpr T log2p1(const T u) { return u ? static_cast(64 - __builtin_clzll(static_cast(u))) : static_cast(0); } template constexpr T msbp1(const T u) { return log2p1(u); } template constexpr T lsbp1(const T u) { return __builtin_ffsll(u); } template constexpr T clog(const T u) { return u ? log2p1(u - 1) : static_cast(u); } template constexpr bool ispow2(const T u) { return u and (static_cast(u) & static_cast(u - 1)) == 0; } template constexpr T ceil2(const T u) { return static_cast(1) << clog(u); } template constexpr T floor2(const T u) { return u == 0 ? static_cast(0) : static_cast(1) << (log2p1(u) - 1); } template constexpr bool btest(const T mask, const usize ind) { return static_cast((static_cast(mask) >> ind) & static_cast(1)); } template void bset(T& mask, const usize ind) { mask |= (static_cast(1) << ind); } template void breset(T& mask, const usize ind) { mask &= ~(static_cast(1) << ind); } template void bflip(T& mask, const usize ind) { mask ^= (static_cast(1) << ind); } template void bset(T& mask, const usize ind, const bool b) { (b ? bset(mask, ind) : breset(mask, ind)); } template constexpr T bcut(const T mask, const usize ind) { return ind == 0 ? static_cast(0) : static_cast((static_cast(mask) << (64 - ind)) >> (64 - ind)); } template bool chmin(T& a, const T& b) { return (a > b ? a = b, true : false); } template bool chmax(T& a, const T& b) { return (a < b ? a = b, true : false); } constexpr unsigned int mod = 1000000007; template constexpr T inf_v = std::numeric_limits::max() / 4; template constexpr Real pi_v = Real{3.141592653589793238462643383279502884}; auto mfp = [](auto&& f) { return [=](auto&&... args) { return f(f, std::forward(args)...); }; }; template T in() { T v; return std::cin >> v, v; } template T in_v(typename std::enable_if<(i == n), c_arr>::type) { return in(); } template auto in_v(typename std::enable_if<(i < n), c_arr>::type& szs) { const usize s = (usize)szs[i]; std::vector(szs))> ans(s); for (usize j = 0; j < s; j++) { ans[j] = in_v(szs); } return ans; } template auto in_v(c_arr szs) { return in_v(szs); } template auto in_t() { return std::tuple...>{in()...}; } struct io_init { io_init() { std::cin.tie(nullptr), std::ios::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(20); } void clear() { std::cin.tie(), std::ios::sync_with_stdio(true); } } io_setting; template int out(const T& v) { return std::cout << v, 0; } template int out(const std::vector& v) { for (usize i = 0; i < v.size(); i++) { if (i > 0) { std::cout << ' '; } out(v[i]); } return std::cout << "\n", 0; } template int out(const std::pair& v) { return out(v.first), std::cout << ' ', out(v.second), 0; } template int out(const T& v, const Args... args) { return out(v), std::cout << ' ', out(args...), 0; } template int outln(const Args... args) { return out(args...), std::cout << '\n', 0; } template void outel(const Args... args) { return out(args...), std::cout << std::endl, 0; } # define SHOW(...) static_cast(0) constexpr ull TEN(const usize n) { return n == 0 ? 1ULL : TEN(n - 1) * 10ULL; } template auto make_v(typename std::enable_if<(i == n), c_arr>::type, const T& v = T{}) { return v; } template auto make_v(typename std::enable_if<(i < n), c_arr>::type szs, const T& v = T{}) { const usize s = (usize)szs[i]; return std::vector(szs, v))>(s, make_v(szs, v)); } template auto make_v(c_arr szs, const T& t = T{}) { return make_v(szs, t); } template class sparse_table { public: using semigroup_type = SemiGroup; using value_type = typename semigroup_type::value_type; template sparse_table(const InIt first, const InIt last) : lg{clog(static_cast(std::distance(first, last) + 1))}, table(lg, std::vector(1 << lg)) { const usize cap = 1 << lg; for (usize d = 0, base = 1; d < lg; d++, base <<= 1) { std::copy(first, last, table[d].begin()); for (usize b = 0; b < cap / base; b++) { if (b % 2 == 1) { for (usize i = 1; i < base; i++) { table[d][b * base + i] = semigroup_type::merge(table[d][b * base + i - 1], table[d][b * base + i]); } } else { for (usize i = 1; i < base; i++) { table[d][(b + 1) * base - i - 1] = semigroup_type::merge(table[d][(b + 1) * base - i - 1], table[d][(b + 1) * base - i]); } } } } } value_type fold(const usize l, const usize r) const { assert(l < r); if (r - l == 1) { return table[0][l]; } const usize d = log2p1(r ^ l) - 1; return semigroup_type::merge(table[d][l], table[d][r - 1]); } private: const usize lg; std::vector> table; }; template, typename Bucket = uint64_t> class static_rmq { public: using value_type = Value; using comparator_type = Comp; template static_rmq(const InIt first, const InIt last) : sz{static_cast(std::distance(first, last))}, bn{wind(sz + bs - 1)}, val{first, last}, bucket_val([&]() { std::vector ans(bn); for (usize i = 0; i < sz; i++) { ans[wind(i)] = i % bs == 0 ? val[i] : std::min(ans[wind(i)], val[i], comparator_type{}); } return ans; }()), masks(sz, 0), st(bucket_val.begin(), bucket_val.end()) { for (usize i = 0; i < bn; i++) { std::vector g(bs, sz); std::stack stack; for (usize j = 0; j < bs and ind(i, j) < sz; j++) { for (; not stack.empty() and not comparator_type{}(val[stack.top()], val[ind(i, j)]); stack.pop()) {} g[j] = stack.empty() ? sz : stack.top(), stack.push(ind(i, j)); } for (usize j = 0; j < bs and ind(i, j) < sz; j++) { masks[ind(i, j)] = g[j] == sz ? static_cast(0) : (masks[g[j]] | static_cast(1) << (g[j] - i * bs)); } } } value_type fold(const usize l, const usize r) const { assert(l < r), assert(r <= sz); const usize lb = (l + bs - 1) / bs, rb = r / bs; if (lb > rb) { return brmq(l, r); } else { return lb < rb ? (l < bs * lb ? (bs * rb < r ? std::min({st.fold(lb, rb), brmq(l, bs * lb), brmq(bs * rb, r)}, comparator_type{}) : std::min(st.fold(lb, rb), brmq(l, bs * lb), comparator_type{})) : (bs * rb < r ? std::min(st.fold(lb, rb), brmq(bs * rb, r), comparator_type{}) : st.fold(lb, rb))) : (l < bs * lb ? (bs * rb < r ? std::min(brmq(l, bs * lb), brmq(bs * rb, r), comparator_type{}) : brmq(l, bs * lb)) : (bs * rb < r ? brmq(bs * rb, r) : value_type{})); } } private: using bucket_type = Bucket; static constexpr usize bs = sizeof(bucket_type) * 8; static constexpr usize bslog = log2p1(bs) - 1; static constexpr usize wind(const usize n) { return n >> (bslog); } static constexpr usize bind(const usize n) { return bcut(n, bslog); } static constexpr usize ind(const usize w, const usize b) { return (w << bslog) | b; } value_type brmq(const usize l, usize r) const { r--; const Bucket w = masks[r] >> (l % bs); return w == 0 ? val[r] : val[l + lsbp1(w) - 1]; } struct min { using value_type = Value; static value_type merge(const value_type& a, const value_type& b) { return std::min(a, b, comparator_type{}); } }; const usize sz, bn; std::vector val, bucket_val; std::vector masks; sparse_table st; }; int main() { const auto N = in(); const auto ss = in_v({N}); SHOW(ss); auto lcp = [](const std::string& s, const std::string& t) { int ans = 0; for (; ans < s.size() and ans < t.size() and s[ans] == t[ans]; ans++) {} return ans; }; std::vector ls(N - 1); for (int i = 0; i < N - 1; i++) { ls[i] = lcp(ss[i], ss[i + 1]); } SHOW(ls); static_rmq rmq(ls.begin(), ls.end()); ll ans = 0; for (const auto& s : ss) { ans += s.size(); } constexpr int L = 500; constexpr int W = 500; for (int i = 0; i <= N - 1; i++) { for (int l = 1; l <= L; l++) { // LNlogN int inf = i, sup = N; while (sup - inf > 1) { const int mid = (inf + sup) / 2; const int len = rmq.fold(i, mid); (len >= l ? inf : sup) = mid; } if (inf == i) { break; } const int width = inf - i + 1; ans += (ll)width * (width + 1) / 2 - 1; } } for (int i = 0; i < N - 1; i++) { for (int w = 1; i + w <= N - 1; w++) { // N^2/L const int len = rmq.fold(i, i + w); if (len > L) { ans += (ll)len * (w + 1); } else { break; } } } outln(ans); return 0; }