[LOJ#2980][THUSCH2017]大魔法师(线段树+矩阵)

每个线段树维护一个行向量[A,B,C,len]分别是这个区间的A,B,C区间和与区间长度，转移显然。

以及此题卡常，稍微哪里写丑了就能100->45。

 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<iostream>
 #define ls (x<<1)
 #define rs (ls|1)
 #define lson ls,L,mid
 #define rson rs,mid+1,R
 #define rep(i,l,r) for (int i=(l); i<=(r); i++)
 using namespace std;

 const int N=,mod=;
 int n,T,op,l,r,x,a[N],b[N],c[N];

 int rd(){
     int x=; char ch=getchar(); bool f=;
     while (ch<'' || ch>'') f|=(ch=='-'),ch=getchar();
     while (ch>='' && ch<='') x=x*+ch-'',ch=getchar();
     return f ? -x : x;
 }

 inline void inc(int &x,int y){ x+=y; x-=x>=mod ? mod : ; }

 struct Mat{
     int a[][];
     Mat (){ memset(a,,sizeof(a)); }
 }I,A[N],v[N<<],tag[N<<];

 Mat operator +(const Mat &a,const Mat &b){
     Mat c=a;
     rep(i,,) inc(c.a[][i],b.a[][i]);
     return c;
 }

 void put(int x,const Mat &a){
     Mat nw;
     rep(i,,) rep(k,,) if (a.a[i][k])
         inc(nw.a[][k],1ll*v[x].a[][i]*a.a[i][k]%mod);
     v[x]=nw; nw=Mat();
     rep(i,,) rep(j,,) if (tag[x].a[i][j]) rep(k,,) if (a.a[j][k])
         inc(nw.a[i][k],1ll*tag[x].a[i][j]*a.a[j][k]%mod);
     tag[x]=nw;
 }

 void push(int x){ put(ls,tag[x]); put(rs,tag[x]); tag[x]=I; }

 void build(int x,int L,int R){
     tag[x]=I;
     if (L==R){ v[x]=A[L]; return; }
     int mid=(L+R)>>;
     build(lson); build(rson); v[x]=v[ls]+v[rs];
 }

 void mdf(int x,int L,int R,int l,int r,Mat &k){
     if (L==l && r==R){ put(x,k); return; }
     int mid=(L+R)>>; push(x);
     if (r<=mid) mdf(lson,l,r,k);
     else if (l>mid) mdf(rson,l,r,k);
         else mdf(lson,l,mid,k),mdf(rson,mid+,r,k);
     v[x]=v[ls]+v[rs];
 }

 Mat que(int x,int L,int R,int l,int r){
     if (L==l && r==R){ return v[x]; }
     int mid=(L+R)>>; push(x);
     if (r<=mid) return que(lson,l,r);
     else if (l>mid) return que(rson,l,r);
         else return que(lson,l,mid)+que(rson,mid+,r);
 }

 int main(){
     freopen("magic.in","r",stdin);
     freopen("magic.out","w",stdout);
     n=rd();
     rep (i,,) I.a[i][i]=;
     rep (i,,n){ A[i].a[][]=; rep (j,,) A[i].a[][j]=rd(); }
     build(,,n);
     for (T=rd(); T--; ){
         op=rd(); l=rd(); r=rd(); Mat s=I;
         if (op>= && op<=) x=rd();
         if (op==){
             Mat res=que(,,n,l,r);
             printf("%d %d %d\n",res.a[][],res.a[][],res.a[][]);
             continue;
         }
         if (op==) s.a[][]=;
         if (op==) s.a[][]=;
         if (op==) s.a[][]=;
         if (op==) s.a[][]=x;
         if (op==) s.a[][]=x;
         if (op==) s.a[][]=,s.a[][]=x;
         mdf(,,n,l,r,s);
     }
     return ;
 }