/** * */ // header {{{ #include using namespace std; using i8 = int8_t; using u8 = uint8_t; using i16 = int16_t; using u16 = uint16_t; using i32 = int32_t; using u32 = uint32_t; using i64 = int64_t; using u64 = uint64_t; #ifdef __SIZEOF_INT128__ using i128 = __int128; using u128 = unsigned __int128; #endif using f32 = float; using f64 = double; using f80 = __float80; using f128 = __float128; // }}} template constexpr T PROCON_INF(); template<> constexpr i64 PROCON_INF() { return 1'010'000'000'000'000'000LL; } template<> constexpr f64 PROCON_INF() { return 1e100; } constexpr i64 INF = PROCON_INF(); constexpr f64 FINF = PROCON_INF(); constexpr i64 MOD = 1'000'000'007LL; constexpr f64 EPS = 1e-12; constexpr f64 PI = 3.14159265358979323846; // util {{{ #define FOR(i, start, end) for(i64 i = (start), i##_end=(end); i < i##_end; ++i) #define REP(i, n) FOR(i, 0, n) #define ALL(f,c,...) (([&](decltype((c)) cccc) { return (f)(std::begin(cccc), std::end(cccc), ## __VA_ARGS__); })(c)) #define SLICE(f,c,l,r,...) (([&](decltype((c)) cccc, decltype((l)) llll, decltype((r)) rrrr) {\ auto iiii = llll <= rrrr ? std::begin(cccc)+llll : std::end(cccc);\ auto jjjj = llll <= rrrr ? std::begin(cccc)+rrrr : std::end(cccc);\ return (f)(iiii, jjjj, ## __VA_ARGS__);\ })(c,l,r)) #define GENERIC(f) ([](auto&&... args) -> decltype(auto) { return (f)(std::forward(args)...); }) // BoolArray {{{ class BoolArray { public: using value_type = bool; using reference = value_type&; using const_reference = const value_type&; using iterator = value_type*; using const_iterator = const value_type*; using difference_type = ptrdiff_t; using size_type = size_t; using reverse_iterator = std::reverse_iterator; using const_reverse_iterator = std::reverse_iterator; BoolArray() : BoolArray(0) {} explicit BoolArray(size_t n) : BoolArray(n,false) {} BoolArray(size_t n, bool value) : size_(n), data_(new bool[n]) { ALL(fill, *this, value); } BoolArray(initializer_list init) : size_(init.size()), data_(new bool[size_]) { ALL(copy, init, begin()); } template BoolArray(InputIt first, InputIt last) { deque tmp(first, last); size_ = tmp.size(); data_ = new bool[size_]; ALL(copy, tmp, begin()); } BoolArray(const BoolArray& other) : size_(other.size_), data_(new bool[size_]) { ALL(copy, other, begin()); } BoolArray(BoolArray&& other) noexcept : size_(other.size_), data_(other.data_) { other.data_ = nullptr; } BoolArray& operator=(const BoolArray& other) { if(this == &other) return *this; if(!data_ || size_ < other.size_) { delete[] data_; data_ = new bool[other.size_]; } size_ = other.size_; ALL(copy, other, begin()); return *this; } BoolArray& operator=(BoolArray&& other) noexcept { if(this == &other) return *this; size_ = other.size_; data_ = other.data_; other.data_ = nullptr; } BoolArray& operator=(initializer_list init) { if(!data_ || size_ < init.size()) { delete[] data_; data_ = new bool[init.size()]; } size_ = init.size(); ALL(copy, init, begin()); return *this; } void swap(BoolArray& other) noexcept { std::swap(size_, other.size_); std::swap(data_, other.data_); } ~BoolArray() { delete[] data_; data_ = nullptr; } bool empty() const noexcept { return size_ == 0; } size_type size() const noexcept { return size_; } size_type max_size() const noexcept { return 1'010'000'000; } iterator begin() noexcept { return data_; } const_iterator begin() const noexcept { return data_; } const_iterator cbegin() const noexcept { return data_; } iterator end() noexcept { return data_+size_; } const_iterator end() const noexcept { return data_+size_; } const_iterator cend() const noexcept { return data_+size_; } reverse_iterator rbegin() noexcept { return reverse_iterator(end()); } const_reverse_iterator rbegin() const noexcept { return const_reverse_iterator(end()); } const_reverse_iterator crbegin() const noexcept { return const_reverse_iterator(end()); } reverse_iterator rend() noexcept { return reverse_iterator(begin()); } const_reverse_iterator rend() const noexcept { return const_reverse_iterator(begin()); } const_reverse_iterator crend() const noexcept { return const_reverse_iterator(begin()); } reference operator[](size_type pos) { return data_[pos]; } const_reference operator[](size_type pos) const { return data_[pos]; } bool* data() noexcept { return data_; } const bool* data() const noexcept { return data_; } private: size_t size_; bool* data_; }; void swap(BoolArray& lhs, BoolArray& rhs) noexcept { lhs.swap(rhs); } bool operator==(const BoolArray& lhs, const BoolArray& rhs) { return equal(begin(lhs), end(lhs), begin(rhs), end(rhs)); } bool operator!=(const BoolArray& lhs, const BoolArray& rhs) { return !(lhs == rhs); } bool operator<(const BoolArray& lhs, const BoolArray& rhs) { return lexicographical_compare(begin(lhs), end(lhs), begin(rhs), end(rhs)); } bool operator> (const BoolArray& lhs, const BoolArray& rhs) { return rhs < lhs; } bool operator<=(const BoolArray& lhs, const BoolArray& rhs) { return !(rhs < lhs); } bool operator>=(const BoolArray& lhs, const BoolArray& rhs) { return !(lhs < rhs); } // }}} // 多次元 vector {{{ // 最内周が vector になるのを避けるための措置 template struct Array1Container { using type = vector; }; template<> struct Array1Container { using type = BoolArray; }; // イテレート用 template struct is_arrayn_container { static constexpr bool value = false; }; template struct is_arrayn_container> { static constexpr bool value = true; }; template<> struct is_arrayn_container { static constexpr bool value = true; }; template auto arrayn_make(i64 n, T x) { using Cont = typename Array1Container::type; return Cont(n, x); } template = nullptr> auto arrayn_make(i64 n, Args... args) { auto inner = arrayn_make(args...); return vector(n, inner); } template enable_if_t::value> arrayn_foreach(T& e, F f) { f(e); } template enable_if_t::value> arrayn_foreach(T& ary, F f) { for(auto& e : ary) arrayn_foreach(e, f); } template enable_if_t::value> arrayn_fill(T& ary, const U& x) { arrayn_foreach(ary, [&x](auto& e) { e = x; }); } // }}} // 多次元生配列 {{{ template enable_if_t::value==0> CARRAY_FOREACH(T& e, F f) { f(e); } template enable_if_t::value!=0> CARRAY_FOREACH(Array& ary, F f) { for(auto& e : ary) CARRAY_FOREACH(e, f); } template enable_if_t::value!=0> CARRAY_FILL(Array& ary, const U& v) { CARRAY_FOREACH(ary, [&v](auto& e) { e = v; }); } // }}} // メモ化ラッパー (8引数まで) {{{ template class Memoized1 { static_assert(N1 >= 1, ""); public: explicit Memoized1(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1) const { using R = decltype(f_(*this,x1)); static bool done[N1] {}; static R memo[N1]; if(!done[x1]) { memo[x1] = f_(*this,x1); done[x1] = true; } return memo[x1]; } private: const F f_; }; template class Memoized2 { static_assert(N1 >= 1 && N2 >= 1, ""); public: explicit Memoized2(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2) const { using R = decltype(f_(*this,x1,x2)); static bool done[N1][N2] {}; static R memo[N1][N2]; if(!done[x1][x2]) { memo[x1][x2] = f_(*this,x1,x2); done[x1][x2] = true; } return memo[x1][x2]; } private: const F f_; }; template class Memoized3 { static_assert(N1 >= 1 && N2 >= 1 && N3 >= 1, ""); public: explicit Memoized3(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2, i64 x3) const { using R = decltype(f_(*this,x1,x2,x3)); static bool done[N1][N2][N3] {}; static R memo[N1][N2][N3]; if(!done[x1][x2][x3]) { memo[x1][x2][x3] = f_(*this,x1,x2,x3); done[x1][x2][x3] = true; } return memo[x1][x2][x3]; } private: const F f_; }; template class Memoized4 { static_assert(N1 >= 1 && N2 >= 1 && N3 >= 1 && N4 >= 1, ""); public: explicit Memoized4(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2, i64 x3, i64 x4) const { using R = decltype(f_(*this,x1,x2,x3,x4)); static bool done[N1][N2][N3][N4] {}; static R memo[N1][N2][N3][N4]; if(!done[x1][x2][x3][x4]) { memo[x1][x2][x3][x4] = f_(*this,x1,x2,x3,x4); done[x1][x2][x3][x4] = true; } return memo[x1][x2][x3][x4]; } private: const F f_; }; template class Memoized5 { static_assert(N1 >= 1 && N2 >= 1 && N3 >= 1 && N4 >= 1 && N5 >= 1, ""); public: explicit Memoized5(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2, i64 x3, i64 x4, i64 x5) const { using R = decltype(f_(*this,x1,x2,x3,x4,x5)); static bool done[N1][N2][N3][N4][N5] {}; static R memo[N1][N2][N3][N4][N5]; if(!done[x1][x2][x3][x4][x5]) { memo[x1][x2][x3][x4][x5] = f_(*this,x1,x2,x3,x4,x5); done[x1][x2][x3][x4][x5] = true; } return memo[x1][x2][x3][x4][x5]; } private: const F f_; }; template class Memoized6 { static_assert(N1 >= 1 && N2 >= 1 && N3 >= 1 && N4 >= 1 && N5 >= 1 && N6 >= 1, ""); public: explicit Memoized6(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2, i64 x3, i64 x4, i64 x5, i64 x6) const { using R = decltype(f_(*this,x1,x2,x3,x4,x5,x6)); static bool done[N1][N2][N3][N4][N5][N6] {}; static R memo[N1][N2][N3][N4][N5][N6]; if(!done[x1][x2][x3][x4][x5][x6]) { memo[x1][x2][x3][x4][x5][x6] = f_(*this,x1,x2,x3,x4,x5,x6); done[x1][x2][x3][x4][x5][x6] = true; } return memo[x1][x2][x3][x4][x5][x6]; } private: const F f_; }; template class Memoized7 { static_assert(N1 >= 1 && N2 >= 1 && N3 >= 1 && N4 >= 1 && N5 >= 1 && N6 >= 1 && N7 >= 1, ""); public: explicit Memoized7(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2, i64 x3, i64 x4, i64 x5, i64 x6, i64 x7) const { using R = decltype(f_(*this,x1,x2,x3,x4,x5,x6,x7)); static bool done[N1][N2][N3][N4][N5][N6][N7] {}; static R memo[N1][N2][N3][N4][N5][N6][N7]; if(!done[x1][x2][x3][x4][x5][x6][x7]) { memo[x1][x2][x3][x4][x5][x6][x7] = f_(*this,x1,x2,x3,x4,x5,x6,x7); done[x1][x2][x3][x4][x5][x6][x7] = true; } return memo[x1][x2][x3][x4][x5][x6][x7]; } private: const F f_; }; template class Memoized8 { static_assert(N1 >= 1 && N2 >= 1 && N3 >= 1 && N4 >= 1 && N5 >= 1 && N6 >= 1 && N7 >= 1 && N8 >= 1, ""); public: explicit Memoized8(F&& f) : f_(forward(f)) {} decltype(auto) operator()(i64 x1, i64 x2, i64 x3, i64 x4, i64 x5, i64 x6, i64 x7, i64 x8) const { using R = decltype(f_(*this,x1,x2,x3,x4,x5,x6,x7,x8)); static bool done[N1][N2][N3][N4][N5][N6][N7][N8] {}; static R memo[N1][N2][N3][N4][N5][N6][N7][N8]; if(!done[x1][x2][x3][x4][x5][x6][x7][x8]) { memo[x1][x2][x3][x4][x5][x6][x7][x8] = f_(*this,x1,x2,x3,x4,x5,x6,x7,x8); done[x1][x2][x3][x4][x5][x6][x7][x8] = true; } return memo[x1][x2][x3][x4][x5][x6][x7][x8]; } private: const F f_; }; template decltype(auto) MEMOIZE(F&& f) { return Memoized1(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized2(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized3(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized4(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized5(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized6(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized7(forward(f)); } template decltype(auto) MEMOIZE(F&& f) { return Memoized8(forward(f)); } // }}} // lambda で再帰 {{{ template class FixPoint { public: explicit constexpr FixPoint(F&& f) : f_(forward(f)) {} template constexpr decltype(auto) operator()(Args&&... args) const { return f_(*this, forward(args)...); } private: const F f_; }; template decltype(auto) FIX(F&& f) { return FixPoint(forward(f)); } // }}} // tuple {{{ template = nullptr> constexpr auto tuple_head(const tuple& t) { return get<0>(t); } template constexpr auto tuple_tail_helper(const tuple& t, index_sequence) { return make_tuple(get(t)...); } template = nullptr> constexpr auto tuple_tail(const tuple&) { return make_tuple(); } template = nullptr> constexpr auto tuple_tail(const tuple& t) { return tuple_tail_helper(t, make_index_sequence()); } // }}} // FST/SND {{{ template T1& FST(pair& p) { return p.first; } template const T1& FST(const pair& p) { return p.first; } template T2& SND(pair& p) { return p.second; } template const T2& SND(const pair& p) { return p.second; } template = nullptr> auto& FST(tuple& t) { return get<0>(t); } template = nullptr> const auto& FST(const tuple& t) { return get<0>(t); } template = nullptr> auto& SND(tuple& t) { return get<1>(t); } template = nullptr> const auto& SND(const tuple& t) { return get<1>(t); } // }}} template i64 SIZE(const C& c) { return static_cast(c.size()); } template i64 SIZE(const T (&)[N]) { return static_cast(N); } bool is_odd (i64 x) { return x % 2 != 0; } bool is_even(i64 x) { return x % 2 == 0; } template i64 cmp(T x, T y) { return (y i64 sgn(T x) { return cmp(x, T(0)); } // lo:OK, hi:NG template i64 bisect_integer(i64 lo, i64 hi, Pred pred) { assert(lo < hi); while(lo+1 < hi) { i64 mid = (lo+hi) / 2; if(pred(mid)) lo = mid; else hi = mid; } return lo; } template f64 bisect_real(f64 lo, f64 hi, Pred pred, i64 iter=100) { assert(lo < hi); REP(_, iter) { f64 mid = (lo+hi) / 2; if(pred(mid)) lo = mid; else hi = mid; } return lo; } i64 ipow(i64 x, i64 e) { assert(e >= 0); i64 res = 1; REP(_, e) { res *= x; } return res; } i64 isqrt(i64 x) { assert(x >= 0); i64 lo = 0; i64 hi = min(x/2+2, 3037000500LL); return bisect_integer(lo, hi, [x](i64 r) { return r*r <= x; }); } // 0 <= ilog2(x) <= 62 i64 ilog2(i64 x) { assert(x > 0); return 63 - __builtin_clzll(x); } // 0 <= ilog10(x) <= 18 i64 ilog10(i64 x) { assert(x > 0); static constexpr i64 TABLE[18] { 9LL, 99LL, 999LL, 9999LL, 99999LL, 999999LL, 9999999LL, 99999999LL, 999999999LL, 9999999999LL, 99999999999LL, 999999999999LL, 9999999999999LL, 99999999999999LL, 999999999999999LL, 9999999999999999LL, 99999999999999999LL, 999999999999999999LL, }; REP(i, SIZE(TABLE)) { if(x <= TABLE[i]) return i; } return SIZE(TABLE); } // Haskell の divMod と同じ pair divmod(i64 a, i64 b) { i64 q = a / b; i64 r = a % b; if((b>0 && r<0) || (b<0 && r>0)) { --q; r += b; } return {q,r}; } i64 div_ceil(i64 a, i64 b) { i64 q = a / b; i64 r = a % b; if((b>0 && r>0) || (b<0 && r<0)) ++q; return q; } i64 div_floor(i64 a, i64 b) { return divmod(a,b).first; } i64 modulo(i64 a, i64 b) { return divmod(a,b).second; } bool feq(f64 x, f64 y, f64 eps=EPS) { return fabs(x-y) < eps; } template> bool chmax(T& xmax, const U& x, Comp comp={}) { if(comp(xmax, x)) { xmax = x; return true; } return false; } template> bool chmin(T& xmin, const U& x, Comp comp={}) { if(comp(x, xmin)) { xmin = x; return true; } return false; } template> ForwardIt bsearch_find(ForwardIt first, ForwardIt last, const T& x, Comp comp={}) { auto it = lower_bound(first, last, x, comp); if(it == last || comp(x,*it)) return last; return it; } // x 未満の最後の要素 template> BidiIt bsearch_lt(BidiIt first, BidiIt last, const T& x, Comp comp={}) { auto it = lower_bound(first, last, x, comp); if(it == first) return last; return prev(it); } // x 以下の最後の要素 template> BidiIt bsearch_le(BidiIt first, BidiIt last, const T& x, Comp comp={}) { auto it = upper_bound(first, last, x, comp); if(it == first) return last; return prev(it); } // x より大きい最初の要素 template> BidiIt bsearch_gt(BidiIt first, BidiIt last, const T& x, Comp comp={}) { return upper_bound(first, last, x, comp); } // x 以上の最初の要素 template> BidiIt bsearch_ge(BidiIt first, BidiIt last, const T& x, Comp comp={}) { return lower_bound(first, last, x, comp); } template auto SUM(InputIt first, InputIt last) { using T = typename iterator_traits::value_type; return accumulate(first, last, T()); } template ForwardIt transform_self(ForwardIt first, ForwardIt last, UnaryOperation op) { return transform(first, last, first, op); } template void UNIQ(C& c) { c.erase(ALL(unique,c), end(c)); } template auto FLIP(BinaryFunc f) { return [f](const auto& x, const auto& y) { return f(y,x); }; } template auto ON(BinaryFunc bf, UnaryFunc uf) { return [bf,uf](const auto& x, const auto& y) { return bf(uf(x), uf(y)); }; } template auto LT_ON(F f) { return ON(less<>(), f); } template auto GT_ON(F f) { return ON(greater<>(), f); } template auto EQ_ON(F f) { return ON(equal_to<>(), f); } template auto NE_ON(F f) { return ON(not_equal_to<>(), f); } template> auto EQUIV(Comp comp={}) { return [comp](const auto& lhs, const auto& rhs) { return !comp(lhs,rhs) && !comp(rhs,lhs); }; } struct IDENTITY { template constexpr T&& operator()(T&& x) const noexcept { return forward(x); } }; char digit_chr(i64 n) { return static_cast('0' + n); } i64 digit_ord(char c) { return c - '0'; } char lower_chr(i64 n) { return static_cast('a' + n); } i64 lower_ord(char c) { return c - 'a'; } char upper_chr(i64 n) { return static_cast('A' + n); } i64 upper_ord(char c) { return c - 'A'; } // 出力は operator<< を直接使わず、このテンプレート経由で行う // 提出用出力とデバッグ用出力を分けるため template struct Formatter { static ostream& write_str(ostream& out, const T& x) { return out << x; } static ostream& write_repr(ostream& out, const T& x) { return out << x; } }; template ostream& WRITE_STR(ostream& out, const T& x) { return Formatter::write_str(out, x); } template ostream& WRITE_REPR(ostream& out, const T& x) { return Formatter::write_repr(out, x); } template ostream& WRITE_JOIN_STR(ostream& out, InputIt first, InputIt last, const string& sep) { while(first != last) { WRITE_STR(out, *first++); if(first != last) out << sep; } return out; } template ostream& WRITE_JOIN_REPR(ostream& out, InputIt first, InputIt last, const string& sep) { while(first != last) { WRITE_REPR(out, *first++); if(first != last) out << sep; } return out; } template ostream& WRITE_RANGE_STR(ostream& out, InputIt first, InputIt last) { return WRITE_JOIN_STR(out, first, last, " "); } template ostream& WRITE_RANGE_REPR(ostream& out, InputIt first, InputIt last) { out << "["; WRITE_JOIN_REPR(out, first, last, ", "); out << "]"; return out; } template void FROM_STR(const string& s, T& x) { istringstream in(s); in >> x; } template string TO_STR(const T& x) { ostringstream out; WRITE_STR(out, x); return out.str(); } template string TO_REPR(const T& x) { ostringstream out; WRITE_REPR(out, x); return out.str(); } template string RANGE_TO_STR(InputIt first, InputIt last) { ostringstream out; WRITE_RANGE_STR(out, first, last); return out.str(); } template string RANGE_TO_REPR(InputIt first, InputIt last) { ostringstream out; WRITE_RANGE_REPR(out, first, last); return out.str(); } template string JOIN(InputIt first, InputIt last, const string& sep) { ostringstream out; WRITE_JOIN_STR(out, first, last, sep); return out.str(); } template<> struct Formatter { static ostream& write_str(ostream& out, i64 x) { return out << x; } static ostream& write_repr(ostream& out, i64 x) { if(x == INF) return out << "INF"; if(x == -INF) return out << "-INF"; return out << x; } }; template<> struct Formatter { static ostream& write_str(ostream& out, f64 x) { return out << x; } static ostream& write_repr(ostream& out, f64 x) { #pragma GCC diagnostic push #pragma GCC diagnostic ignored "-Wfloat-equal" if(x == FINF) return out << "FINF"; if(x == -FINF) return out << "-FINF"; #pragma GCC diagnostic pop return out << x; } }; template struct Formatter> { static ostream& write_str(ostream& out, const vector& v) { return WRITE_RANGE_STR(out, begin(v), end(v)); } static ostream& write_repr(ostream& out, const vector& v) { out << "vector"; return WRITE_RANGE_REPR(out, begin(v), end(v)); } }; template<> struct Formatter { static ostream& write_str(ostream& out, const BoolArray& a) { return WRITE_RANGE_STR(out, begin(a), end(a)); } static ostream& write_repr(ostream& out, const BoolArray& a) { out << "BoolArray"; return WRITE_RANGE_REPR(out, begin(a), end(a)); } }; template struct Formatter> { static ostream& write_str(ostream& out, const pair& p) { WRITE_STR(out, p.first); out << ' '; WRITE_STR(out, p.second); return out; } static ostream& write_repr(ostream& out, const pair& p) { out << "("; WRITE_REPR(out, p.first); out << ","; WRITE_REPR(out, p.second); out << ")"; return out; } }; template struct Formatter> { template = nullptr> static ostream& write_str_impl(ostream& out, const tuple&) { return out; } template = nullptr> static ostream& write_str_impl(ostream& out, const tuple& t) { if(I != 0) out << ' '; WRITE_STR(out, get(t)); return write_str_impl(out, t); } template = nullptr> static ostream& write_repr_impl(ostream& out, const tuple&) { if(sizeof...(TS) == 0) out << "("; return out << ")"; } template = nullptr> static ostream& write_repr_impl(ostream& out, const tuple& t) { if(I == 0) out << "("; else out << ","; WRITE_REPR(out, get(t)); return write_repr_impl(out, t); } static ostream& write_str(ostream& out, const tuple& t) { return write_str_impl(out, t); } static ostream& write_repr(ostream& out, const tuple& t) { return write_repr_impl(out, t); } }; template void RD(T& x) { cin >> x; #ifdef PROCON_LOCAL assert(cin); #endif } template auto RD_ARRAY(i64 n) { auto res = arrayn_make(n, T()); arrayn_foreach(res, [](T& e) { RD(e); }); return res; } template auto RD_ARRAY2(i64 h, i64 w) { auto res = arrayn_make(h,w, T()); arrayn_foreach(res, [](T& e) { RD(e); }); return res; } template pair RD_PAIR() { T1 x; RD(x); T2 y; RD(y); return { x, y }; } template = nullptr> auto RD_TUPLE() { return make_tuple(); } template auto RD_TUPLE() { T x; RD(x); return tuple_cat(make_tuple(x), RD_TUPLE()); } void PRINT() {} template void PRINT(const T& x, const TS& ...args) { WRITE_STR(cout, x); if(sizeof...(args)) { cout << ' '; PRINT(args...); } } template void PRINTLN(const TS& ...args) { PRINT(args...); cout << '\n'; } [[noreturn]] void EXIT() { #ifdef __OPTIMIZE__ # ifdef PROCON_LOCAL cerr.flush(); # endif cout.flush(); _Exit(0); #else exit(0); #endif } template void DBG_IMPL(i64 line, const char* expr, const T& value) { #ifdef PROCON_LOCAL cerr << "[L " << line << "]: "; cerr << expr << " = "; WRITE_REPR(cerr, value); cerr << "\n"; #endif } template void DBG_CARRAY_IMPL(i64 line, const char* expr, const T (&ary)[N]) { #ifdef PROCON_LOCAL cerr << "[L " << line << "]: "; cerr << expr << " = "; WRITE_RANGE_REPR(cerr, begin(ary), end(ary)); cerr << "\n"; #endif } template void DBG_RANGE_IMPL(i64 line, const char* expr1, const char* expr2, InputIt first, InputIt last) { #ifdef PROCON_LOCAL cerr << "[L " << line << "]: "; cerr << expr1 << "," << expr2 << " = "; WRITE_RANGE_REPR(cerr, first, last); cerr << "\n"; #endif } #define DBG(expr) DBG_IMPL(__LINE__, #expr, (expr)) #define DBG_CARRAY(expr) DBG_CARRAY_IMPL(__LINE__, #expr, (expr)) #define DBG_RANGE(first,last) DBG_RANGE_IMPL(__LINE__, #first, #last, (first), (last)) #define PAIR make_pair #define TUPLE make_tuple // }}} // init {{{ struct ProconInit { static constexpr int IOS_PREC = 15; static constexpr bool AUTOFLUSH = false; ProconInit() { cin.tie(nullptr); ios::sync_with_stdio(false); cout << fixed << setprecision(IOS_PREC); #ifdef PROCON_LOCAL cerr << fixed << setprecision(IOS_PREC); #endif if(AUTOFLUSH) cout << unitbuf; } } PROCON_INIT; // }}} //-------------------------------------------------------------------- template void advance_bounded(ForwardIt first, ForwardIt last, ForwardIt& it, i64 n) { if(n > 0) { auto bound = distance(it, last); advance(it, min(n, bound)); } else if(n < 0) { auto bound = distance(it, first); advance(it, max(n, bound)); } } template struct SuffixArray { using value_type = i64; using const_reference = const value_type&; using const_iterator = typename vector::const_iterator; using difference_type = typename vector::difference_type; using size_type = typename vector::size_type; ForwardIt first_; ForwardIt last_; vector sa_; SuffixArray(ForwardIt first, ForwardIt last) : first_(first), last_(last), sa_(distance(first,last)+1) { i64 n = distance(first_,last_); vector rnk(n+1); ALL(iota, sa_, 0); copy(first_, last_, std::begin(rnk)); rnk[n] = -1; vector rnk_nex(n+1); for(i64 k = 1; k <= n; k *= 2) { // TODO: 条件部の等号いらない気がするんだが… auto comp = LT_ON([n,&rnk,k](i64 i) { i64 snd = i+k <= n ? rnk[i+k] : -1; return make_pair(rnk[i], snd); }); ALL(sort, sa_, comp); rnk_nex[sa_[0]] = 0; FOR(i, 1, n+1) { rnk_nex[sa_[i]] = rnk_nex[sa_[i-1]] + (comp(sa_[i-1],sa_[i]) ? 1 : 0); } swap(rnk, rnk_nex); } } const_iterator begin() const { return std::begin(sa_); } const_iterator end() const { return std::end(sa_); } const_iterator cbegin() const { return std::cbegin(sa_); } const_iterator cend() const { return std::cend(sa_); } template const_iterator find(ForwardIt2 first2, ForwardIt2 last2) const { auto x = make_pair(first2,last2); auto comp = Compare(first_, last_); auto it = ALL(std::lower_bound, sa_, x, comp); if(it == end() || comp(x,*it)) return end(); return it; } template pair equal_range(ForwardIt2 first2, ForwardIt2 last2) const { auto comp = Compare(first_, last_); return ALL(std::equal_range, sa_, make_pair(first2,last2), comp); } template struct Compare { ForwardIt first_; ForwardIt last_; Compare(ForwardIt first, ForwardIt last) : first_(first), last_(last) {} bool operator()(i64 i, const pair& p) const { auto q = get_range(i, distance(FST(p),SND(p))); return lexicographical_compare(FST(q), SND(q), FST(p), SND(p)); } bool operator()(const pair& p, i64 i) const { auto q = get_range(i, distance(FST(p),SND(p))); return lexicographical_compare(FST(p), SND(p), FST(q), SND(q)); } pair get_range(i64 i, i64 n) const { auto first1 = first_ + i; auto last1 = first_ + i; advance_bounded(first_, last_, last1, n); return make_pair(first1, last1); } }; }; template SuffixArray make_suffix_array(ForwardIt first, ForwardIt last) { return SuffixArray(first, last); } void solve() { string S; RD(S); i64 M; RD(M); auto sa = ALL(make_suffix_array, S); i64 ans = 0; REP(_, M) { string t; RD(t); auto p = ALL(sa.equal_range, t); if(distance(FST(p),SND(p)) > 0) assert(*ALL(sa.find,t) == *FST(p)); ans += distance(FST(p), SND(p)); } PRINTLN(ans); } signed main() { solve(); EXIT(); }