【CF1252K】Addition Robot（线段树，矩阵乘法）

题意：

 

思路：因为线段树上每一段的矩阵之积只有两种，预处理一下，翻转的时候下传tag然后把另一种可能性换上来就好

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned int uint;
typedef unsigned long long ull;
typedef long double ld;
typedef pair<int,int> PII;
typedef pair<ll,ll> Pll;
typedef vector<int> VI;
typedef vector<PII> VII;
//typedef pair<ll,ll>P;
#define N  200010
//#define M  200010
#define INF 1e9
#define fi first
#define se second
#define MP make_pair
#define pb push_back
#define pi acos(-1)
#define mem(a,b) memset(a,b,sizeof(a))
#define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
#define per(i,a,b) for(int i=(int)a;i>=(int)b;i--)
#define lowbit(x) x&(-x)
#define Rand (rand()*(1<<16)+rand())
#define id(x) ((x)<=B?(x):m-n/(x)+1)
#define ls p<<1
#define rs p<<1|1

const ll MOD=1e9+7,inv2=(MOD+1)/2;
      double eps=1e-6;
      int dx[4]={-1,1,0,0};
      int dy[4]={0,0,-1,1};

struct Ma
{
    int n=2,m=2;
    ll a[2][2];
    void init()
    {
        n=m=0;
        mem(a,0);
    }
    Ma operator +(Ma b) const
    {
        Ma c;
        c.n=n; c.m=m;
        rep(i,0,n-1)
         rep(j,0,m-1) c.a[i][j]=a[i][j]+b.a[i][j];
        return c;
    }
    Ma operator *(Ma b) const
    {
        Ma c;
        c.init();
        c.n=n; c.m=b.m;
        rep(i,0,n-1)
        rep(j,0,b.m-1)
        rep(k,0,m-1) c.a[i][j]=(c.a[i][j]+a[i][k]*b.a[k][j]%MOD)%MOD;
        return c;
    }
};

Ma t[N<<2][2],M[2],one,ans;
int b[N],tag[N];
char s[N];

int read()
{
   int v=0,f=1;
   char c=getchar();
   while(c<48||57<c) {if(c=='-') f=-1; c=getchar();}
   while(48<=c&&c<=57) v=(v<<3)+v+v+c-48,   c=getchar();
   return v*f;
}

void pushdown(int p)
{
    if(tag[p])
    {
        swap(t[ls][0],t[ls][1]);
        swap(t[rs][0],t[rs][1]);
        tag[ls]^=1;
        tag[rs]^=1;
        tag[p]=0;
    }
}

void pushup(int p)
{
    t[p][0]=t[ls][0]*t[rs][0];
    t[p][1]=t[ls][1]*t[rs][1];
}

void build(int l,int r,int p)
{
    if(l==r)
    {
        t[p][0]=M[b[l]];
        t[p][1]=M[b[l]^1];
        tag[p]=0;
        return;
    }
    int mid=(l+r)>>1;
    build(l,mid,ls);
    build(mid+1,r,rs);
    pushup(p);

}

void update(int l,int r,int x,int y,int p)
{
    if(x<=l&&r<=y)
    {
        tag[p]^=1;
        swap(t[p][0],t[p][1]);
        return;
    }
    pushdown(p);
    int mid=(l+r)>>1;
    if(x<=mid) update(l,mid,x,y,ls);
    if(y>mid) update(mid+1,r,x,y,rs);
    pushup(p);
}

Ma query(int l,int r,int x,int y,int p)
{
    if(x<=l&&r<=y) return t[p][0];
    pushdown(p);
    int mid=(l+r)>>1;
    Ma res=one;
    if(x<=mid) res=res*query(l,mid,x,y,ls);
    if(y>mid) res=res*query(mid+1,r,x,y,rs);
    return res;
}

int main()
{
    //freopen("1.in","r",stdin);
    int n=read(),q=read();
    scanf("%s",s+1);
    rep(i,1,n) b[i]=s[i]-'A';
    M[0].n=M[0].m=M[1].n=M[1].m=one.n=one.m=ans.n=ans.m=2;
    M[0].a[0][0]=M[0].a[1][0]=M[0].a[1][1]=1;
    M[1].a[0][0]=M[1].a[0][1]=M[1].a[1][1]=1;
    one.a[0][0]=one.a[1][1]=1;
    build(1,n,1);

    while(q--)
    {
        int op=read();
        if(op==1)
        {
            int x=read(),y=read();
            update(1,n,x,y,1);
        }
         else
         {
             int x=read(),y=read(),A=read(),B=read();
             Ma t=query(1,n,x,y,1);
            // printf("t=\n");
             //rep(i,0,1)
             //{
              //  rep(j,0,1) printf("%I64d ",t.a[i][j]);
              //  printf("\n");
             //}
             ans.a[0][0]=A;
             ans.a[0][1]=B;
             ans.a[1][0]=ans.a[1][1]=0;
             ans=ans*t;
             printf("%I64d %I64d\n",ans.a[0][0],ans.a[0][1]);
         }
    }
    return 0;
}