#include using namespace std; using ll=long long; using ull=unsigned long long; using P=pair; templateusing minque=priority_queue,greater>; templatebool chmax(T &a,const T &b){return (abool chmin(T &a,const T &b){return (a>b?(a=b,true):false);} templateistream &operator>>(istream &is,pair&p){is>>p.first>>p.second;return is;} templateistream &operator>>(istream &is,tuple&a){is>>std::get<0>(a)>>std::get<1>(a)>>std::get<2>(a);return is;} templateistream &operator>>(istream &is,array&a){for(auto&i:a)is>>i;return is;} templateistream &operator>>(istream &is,vector &a){for(auto &i:a)is>>i;return is;} templatevoid operator++(pair&a,int n){a.first++,a.second++;} templatevoid operator--(pair&a,int n){a.first--,a.second--;} templatevoid operator++(vector&a,int n){for(auto &i:a)i++;} templatevoid operator--(vector&a,int n){for(auto &i:a)i--;} #define overload3(_1,_2,_3,name,...) name #define rep1(i,n) for(int i=0;i<(int)(n);i++) #define rep2(i,l,r) for(int i=(int)(l);i<(int)(r);i++) #define rep(...) overload3(__VA_ARGS__,rep2,rep1)(__VA_ARGS__) #define reps(i,l,r) rep2(i,l,r) #define all(x) x.begin(),x.end() #define pcnt(x) __builtin_popcountll(x) #define fin(x) return cout<<(x)<<'\n',static_cast(0) #define yn(x) cout<<((x)?"Yes\n":"No\n") #define uniq(x) sort(all(x)),x.erase(unique(all(x)),x.end()) template inline int fkey(vector&z,T key){return lower_bound(z.begin(),z.end(),key)-z.begin();} ll myceil(ll a,ll b){return (a+b-1)/b;} template auto vec(const int (&d)[n],const T &init=T()){ if constexpr (id(d,init)); else return init; } #ifdef LOCAL #include #define SWITCH(a,b) (a) #else #define debug(...) static_cast(0) #define debugg(...) static_cast(0) #define SWITCH(a,b) (b) templateostream &operator<<(ostream &os,const pair&p){os<>testcase; for(int i=0;i #include template constexpr std::enable_if_t::digits<=32,int>msb(T n){return n==0?-1:31-__builtin_clz(n);} template constexpr std::enable_if_t<(std::numeric_limits::digits>32),int>msb(T n){return n==0?-1:63-__builtin_clzll(n);} template constexpr std::enable_if_t::digits<=32,int>lsb(T n){return n==0?-1:__builtin_ctz(n);} template constexpr std::enable_if_t<(std::numeric_limits::digits>32),int>lsb(T n){return n==0?-1:__builtin_ctzll(n);} template constexpr std::enable_if_t,T>floor_pow2(T n){return n==0?0:T(1)< constexpr std::enable_if_t,T>ceil_pow2(T n){return n<=1?1:T(1)<<(msb(n-1)+1);} template constexpr T safe_div(T a,T b){return a/b-(a%b&&(a^b)<0);} template constexpr T safe_ceil(T a,T b){return a/b+(a%b&&(a^b)>0);} template struct SparseTable{ using S=typename M::S; private: std::vectordat,prefix,suffix; std::vector>sp; public: SparseTable(){} SparseTable(std::vectora):dat(a){ int n=a.size(); n=(n+(1<d2(n>>L,M::e()); for(int i=0;i>L);i++){ for(int j=1;j<(1<>L)-1;i>=0;i--){ for(int j=(1<=1;j--)suffix[(i<=j-w;k--)sp[i][k]=M::op(d2[k],sp[i][k+1]); int r=std::min(d2.size(),j+w); for(int k=j+1;k>L,rid=r>>L; if(lid==rid){ S ret=M::e(); for(int i=l;i<=r;i++)ret=M::op(ret,dat[i]); return ret; } else{ lid++; rid--; S mid=M::e(); if(lid==rid)mid=sp[0][lid]; else if(lid struct ConstantMatrix{ private: std::array,N>dat; public: ConstantMatrix(){ dat.fill([](){std::arrayres;res.fill(T());return res;}()); } static ConstantMatrix E(){ static_assert(N==M); ConstantMatrix res; for(int i=0;ifriend ConstantMatrix operator*(const ConstantMatrix&lhs,const ConstantMatrix&rhs){ ConstantMatrixres; for(int i=0;i &operator[](int i){return dat[i];} const std::array &operator[](int i)const{return dat[i];} friend bool operator==(const ConstantMatrix&lhs,const ConstantMatrix&rhs){ for(int i=0;i ConstantMatrix pow(U k)const{ static_assert(N==M); ConstantMatrix res=ConstantMatrix::E(),a(*this); while(k){ if(k&1)res*=a; a*=a; k>>=1; } return res; } std::vector>to_vector()const{ std::vector>res(N,std::vector(M)); for(int i=0;imat,allmat; ll cntr[2]={0,0}; // ll al=0; ll ans=0; int rop=-1; }; static S op(S x,S y){ if(x.rop==-1)return y; if(y.rop==-1)return x; S res; res.rop=y.rop; res.mat=x.mat+y.mat*x.allmat; res.allmat=y.allmat*x.allmat; rep(i,2){ res.cntr[i]+=y.cntr[i]; } // res.cntl[operate(x.al,0,x.rop)]+=y.cntl[0]; // res.cntl[operate(x.al,1,x.rop)]+=y.cntl[1]; // debug(res.cntr[0]); // debug(y.allmat); debug(y.allmat); rep(i,2)rep(j,2)if(y.allmat[j][i])res.cntr[j]+=x.cntr[i]; // res.cntr[operate(0,y.al,x.rop)]+=x.cntr[0]; // res.cntr[operate(1,y.al,x.rop)]+=x.cntr[1]; // res.al=operate(x.al,y.al,x.rop); res.ans=x.ans+y.ans; // debug(res.cntr[0]); // debug(x.ans,y.ans,res.ans); rep(i,2){ res.ans+=x.cntr[i]*y.mat[1][i]; // if(y.allmat[1][i])res.ans+=x.cntr[i]; // debug(i,res.ans,x.cntr[i]); } // rep(i,2)rep(j,2){ // if(operate(i,j,x.rop))res.ans+=(ll)x.cntr[i]*y.cntl[j]; // } return res; } static S e(){return S();} }; ConstantMatrixgetmat(char c,int v){ ConstantMatrixres; if(c=='+'){ if(v==0)res[0][0]=res[1][1]=1; else res[1][0]=res[1][1]=1; } else if(c=='*'){ if(v==1)res[0][0]=res[1][1]=1; else res[0][0]=res[0][1]=1; } else{ if(v==0)res[0][0]=res[1][1]=1; else res[0][1]=res[1][0]=1; } return res; } void SOLVE(){ int n,q; cin>>n>>q; n=(n+1)/2; string x; cin>>x; // x+='+'; vectorinit(n); rep(i,n){ // init[i].cntl[x[i*2]=='T']++; init[i].cntr[x[i*2]=='T']++; // init[i].al=x[i*2]=='T'; if(x[i*2]=='T')init[i].ans++; init[i].rop=0; init[i].mat=init[i].allmat=getmat(i==0?'+':x[i*2-1],x[i*2]=='T'); init[i].allmat=init[i].mat; // debug(i,init[i]); debug(i,init[i].cntr[0]); } SparseTablesp(init); debug('U'); rep(i,q){ int l,r; cin>>l>>r; l--,r--; l/=2,r/=2; cout<