#line 1 "test/aoj_0355.test.cpp"
#define PROBLEM "http://judge.u-aizu.ac.jp/onlinejudge/description.jsp?id=0355"
#include <cstdint>
#include <cstdio>
#include <algorithm>
#include <string>
#include <utility>
#include <vector>
#line 1 "ModularArithmetic/modint.cpp"
/**
* @brief 合同算術用クラス
* @author えびちゃん
*/
#line 10 "ModularArithmetic/modint.cpp"
#include <limits>
#include <type_traits>
#line 13 "ModularArithmetic/modint.cpp"
template <intmax_t Modulo>
class modint {
public:
using value_type = typename std::conditional<
(0 < Modulo && Modulo < std::numeric_limits<int>::max() / 2), int, intmax_t
>::type;
private:
static constexpr value_type S_cmod = Modulo; // compile-time
static value_type S_rmod; // runtime
value_type M_value = 0;
static constexpr value_type S_inv(value_type n, value_type m) {
value_type x = 0;
value_type y = 1;
value_type a = n;
value_type b = m;
for (value_type u = y, v = x; a;) {
value_type q = b / a;
std::swap(x -= q*u, u);
std::swap(y -= q*v, v);
std::swap(b -= q*a, a);
}
if ((x %= m) < 0) x += m;
return x;
}
static value_type S_normalize(intmax_t n, value_type m) {
if (n >= m) {
n %= m;
} else if (n < 0) {
if ((n %= m) < 0) n += m;
}
return n;
}
public:
modint() = default;
modint(intmax_t n): M_value(S_normalize(n, get_modulo())) {}
modint& operator =(intmax_t n) {
M_value = S_normalize(n, get_modulo());
return *this;
}
modint& operator +=(modint const& that) {
if ((M_value += that.M_value) >= get_modulo()) M_value -= get_modulo();
return *this;
}
modint& operator -=(modint const& that) {
if ((M_value -= that.M_value) < 0) M_value += get_modulo();
return *this;
}
modint& operator *=(modint const& that) {
intmax_t tmp = M_value;
tmp *= that.M_value;
M_value = tmp % get_modulo();
return *this;
}
modint& operator /=(modint const& that) {
intmax_t tmp = M_value;
tmp *= S_inv(that.M_value, get_modulo());
M_value = tmp % get_modulo();
return *this;
}
modint& operator ++() {
if (++M_value == get_modulo()) M_value = 0;
return *this;
}
modint& operator --() {
if (M_value-- == 0) M_value = get_modulo()-1;
return *this;
}
modint operator ++(int) { modint tmp(*this); ++*this; return tmp; }
modint operator --(int) { modint tmp(*this); --*this; return tmp; }
friend modint operator +(modint lhs, modint const& rhs) { return lhs += rhs; }
friend modint operator -(modint lhs, modint const& rhs) { return lhs -= rhs; }
friend modint operator *(modint lhs, modint const& rhs) { return lhs *= rhs; }
friend modint operator /(modint lhs, modint const& rhs) { return lhs /= rhs; }
modint operator +() const { return *this; }
modint operator -() const {
if (M_value == 0) return *this;
return modint(get_modulo() - M_value);
}
friend bool operator ==(modint const& lhs, modint const& rhs) {
return lhs.M_value == rhs.M_value;
}
friend bool operator !=(modint const& lhs, modint const& rhs) {
return !(lhs == rhs);
}
value_type get() const { return M_value; }
static value_type get_modulo() { return ((S_cmod > 0)? S_cmod: S_rmod); }
template <int M = Modulo, typename Tp = typename std::enable_if<(M <= 0)>::type>
static Tp set_modulo(value_type m) { S_rmod = m; }
};
template <intmax_t N>
constexpr typename modint<N>::value_type modint<N>::S_cmod;
template <intmax_t N>
typename modint<N>::value_type modint<N>::S_rmod;
#line 1 "utility/monoid/rolling_hash.cpp"
/**
* @brief ロリハの演算のモノイド
* @author えびちゃん
*/
#line 10 "utility/monoid/rolling_hash.cpp"
template <typename ModInt>
class rolling_hash_monoid {
public:
using value_type = ModInt;
private:
value_type M_x = 0, M_p = 1;
public:
rolling_hash_monoid() = default; // identity
rolling_hash_monoid(value_type const& x, value_type const& p): M_x(x), M_p(p) {};
rolling_hash_monoid& operator +=(rolling_hash_monoid const& that) {
M_x = M_x * that.M_p + that.M_x;
M_p *= that.M_p;
return *this;
}
friend bool operator ==(rolling_hash_monoid const& lhs, rolling_hash_monoid const& rhs) {
return lhs.M_x == rhs.M_x && lhs.M_p == rhs.M_p;
}
friend rolling_hash_monoid operator +(rolling_hash_monoid lhs, rolling_hash_monoid const& rhs) {
return lhs += rhs;
}
friend bool operator !=(rolling_hash_monoid const& lhs, rolling_hash_monoid const& rhs) {
return !(lhs == rhs);
}
value_type const& get() const { return M_x; }
value_type const& coefficient() const { return M_p; }
};
#line 1 "utility/monoid/set.cpp"
/**
* @brief モノイドクラス
* @author えびちゃん
*/
template <typename Tp>
class set_monoid {
public:
using value_type = Tp;
private:
bool M_empty = true;
value_type M_x;
public:
set_monoid() = default; // identity
set_monoid(value_type const& x): M_empty(false), M_x(x) {}
set_monoid& operator +=(set_monoid const& that) {
M_empty = that.M_empty;
if (!that.M_empty) M_x = that.M_x;
return *this;
}
friend bool operator ==(set_monoid const& lhs, set_monoid const& rhs) {
if (lhs.M_empty && rhs.M_empty) return true;
if (lhs.M_empty != rhs.M_empty) return false;
return lhs.M_x == rhs.M_x;
}
friend set_monoid operator +(set_monoid lhs, set_monoid const& rhs) { return lhs += rhs; }
friend bool operator !=(set_monoid const& lhs, set_monoid const& rhs) {
return !(lhs == rhs);
}
bool empty() const noexcept { return M_empty; }
value_type const& get() const { return M_x; }
};
#line 1 "DataStructure/segment_tree.cpp"
/**
* @brief 区間作用・区間和セグメント木
* @author えびちゃん
*/
#include <cstddef>
#line 11 "DataStructure/segment_tree.cpp"
#line 1 "integer/bit.cpp"
/**
* @brief ビット演算
* @author えびちゃん
*/
// XXX integral promotion 関連の注意をあまりしていません
#include <climits>
#include <type_traits>
template <typename Tp>
constexpr auto countl_zero(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, int>::type
{
using value_type = typename std::make_unsigned<Tp>::type;
int bits = (sizeof(value_type) * CHAR_BIT);
if (n == 0) return bits;
int res = 0;
for (int i = bits / 2; i > 0; i /= 2) {
value_type mask = ((static_cast<value_type>(1) << i) - 1) << i;
if (n & mask) n >>= i;
else res += i;
}
return res;
}
template <typename Tp>
constexpr auto countl_one(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, int>::type
{
using value_type = typename std::make_unsigned<Tp>::type;
return countl_zero(static_cast<value_type>(~n));
}
template <typename Tp>
constexpr auto countr_zero(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, int>::type
{
using value_type = typename std::make_unsigned<Tp>::type;
int bits = (sizeof(value_type) * CHAR_BIT);
if (n == 0) return bits;
int res = 0;
for (int i = bits / 2; i > 0; i /= 2) {
value_type mask = ((static_cast<value_type>(1) << i) - 1);
if (!(n & mask)) res += i, n >>= i;
}
return res;
}
template <typename Tp>
constexpr auto countr_one(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, int>::type
{
using value_type = typename std::make_unsigned<Tp>::type;
return countr_zero(static_cast<value_type>(~n));
}
constexpr unsigned long long half_mask[] = {
0x5555555555555555uLL, 0x3333333333333333uLL, 0x0F0F0F0F0F0F0F0FuLL,
0x00FF00FF00FF00FFuLL, 0x0000FFFF0000FFFFuLL, 0x00000000FFFFFFFFuLL
};
template <typename Tp>
constexpr auto popcount(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, int>::type
{
int bits = static_cast<int>((sizeof n) * CHAR_BIT);
for (int i = 0, j = 1; j < bits; ++i, j *= 2) {
if (j <= 8) n = (n & half_mask[i]) + ((n >> j) & half_mask[i]);
else n += n >> j;
}
return n & 0xFF;
}
template <typename Tp>
constexpr auto parity(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, int>::type
{ return popcount(n) & 1; }
template <typename Tp>
int clz(Tp n) { return countl_zero(static_cast<typename std::make_unsigned<Tp>::type>(n)); }
template <typename Tp>
int ctz(Tp n) { return countr_zero(static_cast<typename std::make_unsigned<Tp>::type>(n)); }
template <typename Tp>
int ilog2(Tp n) {
return (CHAR_BIT * sizeof(Tp) - 1) - clz(static_cast<typename std::make_unsigned<Tp>::type>(n));
}
template <typename Tp>
bool is_pow2(Tp n) { return (n > 0) && ((n & (n-1)) == 0); }
template <typename Tp>
Tp floor2(Tp n) { return is_pow2(n)? n: static_cast<Tp>(1) << ilog2(n); }
template <typename Tp>
Tp ceil2(Tp n) { return is_pow2(n)? n: static_cast<Tp>(2) << ilog2(n); }
template <typename Tp>
constexpr auto reverse(Tp n)
-> typename std::enable_if<std::is_unsigned<Tp>::value, Tp>::type
{
int bits = static_cast<int>((sizeof n) * CHAR_BIT);
for (int i = 0, j = 1; j < bits; ++i, j *= 2) {
n = ((n & half_mask[i]) << j) | ((n >> j) & half_mask[i]);
}
return n;
}
#line 13 "DataStructure/segment_tree.cpp"
template <typename Operation>
class segment_tree {
public:
using size_type = size_t;
using operation = Operation;
using operand_type = typename operation::operand_type;
using action_type = typename operation::action_type;
using value_type = operand_type;
private:
size_type M_n;
std::vector<operand_type> M_c;
std::vector<action_type> M_d; // deferred
void M_build(size_type i) {
while (i > 1) {
i >>= 1;
M_c[i] = (M_c[i<<1|0] + M_c[i<<1|1]);
operation::act(M_c[i], M_d[i]);
}
}
void M_resolve(size_type i) {
size_type h = ilog2(M_n) + 2; // ilog2p1(M_n*2)
for (size_type s = h; s > 0; --s) {
size_type p = i >> s;
action_type id{};
if (M_d[p] != id) {
M_apply(p<<1|0, M_d[p]);
M_apply(p<<1|1, M_d[p]);
M_d[p] = id;
}
}
}
void M_apply(size_type i, action_type const& x) {
operation::act(M_c[i], x);
if (i < M_n) M_d[i] += x;
}
public:
segment_tree() = default;
explicit segment_tree(size_type n):
M_n(n), M_c(n+n, operand_type{}), M_d(n, action_type{}) {}
segment_tree(size_type n, operand_type const& x):
M_n(n), M_c(n+n, x), M_d(n, action_type{})
{
for (size_type i = n; i--;) M_c[i] = M_c[i<<1|0] + M_c[i<<1|1];
}
template <typename InputIt>
segment_tree(InputIt first, InputIt last): M_c(first, last) {
M_n = M_c.size();
M_d.assign(M_n, action_type{});
M_c.insert(M_c.begin(), M_n, operand_type{});
for (size_type i = M_n; i--;) M_c[i] = M_c[i<<1|0] + M_c[i<<1|1];
}
void assign(size_type n) {
M_n = n;
M_c(n+n, operand_type{});
M_d(n, action_type{});
}
void assign(size_type n, operand_type const& x) {
M_n = n;
M_c(n+n, x);
M_d(n, action_type{});
for (size_type i = n; i--;) M_c[i] = M_c[i<<1|0] + M_c[i<<1|1];
}
template <typename InputIt>
void assign(InputIt first, InputIt last) {
M_c.assign(first, last);
M_n = M_c.size();
M_d.assign(M_n, action_type{});
M_c.insert(M_c.begin(), M_n, operand_type{});
for (size_type i = M_n; i--;) M_c[i] = M_c[i<<1|0] + M_c[i<<1|1];
}
void act(size_type l, size_type r, action_type const& x) {
if (l == r) return;
l += M_n;
r += M_n;
size_type l0 = l;
size_type r0 = r;
M_resolve(l0);
M_resolve(r0-1);
while (l < r) {
if (l & 1) M_apply(l++, x);
if (r & 1) M_apply(--r, x);
l >>= 1;
r >>= 1;
}
M_build(l0);
M_build(r0-1);
}
operand_type fold(size_type l, size_type r) {
operand_type resl{}, resr{};
if (l == r) return resl;
l += M_n;
r += M_n;
M_resolve(l);
M_resolve(r-1);
while (l < r) {
if (l & 1) resl += M_c[l++];
if (r & 1) resr = M_c[--r] + std::move(resr);
l >>= 1;
r >>= 1;
}
return resl += resr;
}
operand_type operator [](size_type i) {
i += M_n;
M_resolve(i);
return M_c[i];
}
};
#line 14 "test/aoj_0355.test.cpp"
constexpr intmax_t mod = 998244353;
using mi = modint<mod>;
template <typename ModInt>
struct action_set_to_rolling_hash {
using operand_type = rolling_hash_monoid<ModInt>;
using action_type = set_monoid<std::pair<ModInt, ModInt>>;
static void act(operand_type& op, action_type const& a) {
if (a.empty()) return;
operand_type tmp(a.get().first, a.get().second);
if (op.coefficient() == 1) {
op = tmp;
return;
}
while (tmp.coefficient() != op.coefficient()) tmp += tmp;
op = tmp;
}
};
int main() {
size_t n;
char buf[262144];
int q;
scanf("%zu %s %d", &n, buf, &q);
std::string s = buf;
intmax_t b1 = 2352983, b2 = 4986917;
segment_tree<action_set_to_rolling_hash<mi>> s1(n), s2(n);
auto f = [](auto x, auto y) {
return set_monoid<std::pair<mi, mi>>({x, y});
};
for (size_t i = 0; i < n; ++i) {
s1.act(i, i+1, f(s[i], b1));
s2.act(i, i+1, f(s[i], b2));
}
for (int i = 0; i < q; ++i) {
scanf("%s", buf);
std::string type = buf;
if (type == "set") {
size_t x, y;
char z;
scanf("%zu %zu %c", &x, &y, &z);
--x, --y;
s1.act(x, y+1, f(z, b1));
s2.act(x, y+1, f(z, b2));
} else if (type == "comp") {
size_t a, b, c, d;
scanf("%zu %zu %zu %zu", &a, &b, &c, &d);
--a, --b, --c, --d;
size_t sl = b-a+1, tl = d-c+1;
auto eq = [&](size_t l) {
return s1.fold(a, a+l) == s1.fold(c, c+l) && s2.fold(a, a+l) == s2.fold(c, c+l);
};
if (eq(std::min(sl, tl))) {
printf("%c\n", ((sl < tl)? 's': (sl > tl)? 't': 'e'));
continue;
}
size_t lb = 0, ub = std::min(sl, tl);
while (ub-lb > 1) {
size_t mid = (lb+ub) >> 1;
(eq(mid)? lb: ub) = mid;
}
int sc = s1[a+lb].get().get(), tc = s1[c+lb].get().get();
printf("%c\n", ((sc < tc)? 's': 't'));
}
}
}
test/aoj_0355.test.cpp