BZOJ 4085:[Sdoi2015]quality(round 2 音质检测)(数据结构)

居然在考场上把这道题打出来了觉得自己也是有点吊啊（虽然后面就没时间做其他题了囧而且还被卡常数了。。。）

题解自己写了一份TEX的就直接放上来吧。。。。

好啦，在谈点什么别的

什么？你在bz上TLE了？注意一下你的矩阵乘法，这个程序的大部分时间几乎都是跑矩阵乘法的，我是从900000次到600000次在到300000次最后预处理那些2的次幂才过的。把OJ卡了好久真是对不起啊QAQ

CODE：

 #include<cstdio>
 #include<iostream>
 #include<cstring>
 #include<algorithm>
 using namespace std;
 #define maxn 300010
 #define mod 1000000007
 typedef long long ll;
 inline ll power(ll x,int y) {
     ll ans=;
     for (;y;y>>=) {
     if (y&) (ans*=x)%=mod;
     (x*=x)%=mod;
     }
     return ans;
 }
 struct node{
     int l,r,ch[][],lz[],s[][];
     inline void print() {
     printf("%d %d\n",l,r);
     for (int i=;i<;i++,puts(""))
         for (int j=;j<;j++) printf("%d ",ch[i][j]);
     for (int i=;i<;i++,puts(""))
         for (int j=;j<;j++) printf("%d ",s[i][j]);
     printf("%d %d\n",lz[],lz[]);
     }
 }t[maxn*];
 #define lc (x<<1)
 #define rc (lc^1)
 #define mid ((l+r)>>1)
 inline void update(int x){
     for (int i=;i<;i++)
     for (int j=;j<;j++) t[x].ch[i][j]=(t[lc].ch[i][j]+t[rc].ch[i][j])%mod;
     for (int i=;i<;i++)
     for(int j=;j<;j++) t[x].s[i][j]=(t[lc].s[i][j]+t[rc].s[i][j])%mod;
 }
 int a,b,inv;
 inline void add(int x,int l){
     for (int i=;i<;i++) t[x].s[l][i]=t[x].s[l][i+];
     t[x].s[l][]=(t[x].s[l][]*1ll+t[x].s[l][]*1ll*a%mod+b*1ll*(t[x].r-t[x].l+)%mod)%mod;
     if (l==) {
     for (int i=;i<;i++)
         for (int j=;j<;j++) t[x].ch[i][j]=t[x].ch[i+][j];
     for (int i=;i<;i++)
         t[x].ch[][i]=(t[x].ch[][i]*1ll+t[x].ch[][i]*1ll*a%mod+b*1ll*t[x].s[][i]%mod)%mod;
     }
     if (l==) {
     for (int i=;i<;i++)
         for (int j=;j<;j++) t[x].ch[j][i]=t[x].ch[j][i+];
     for (int i=;i<;i++)
         t[x].ch[i][]=(t[x].ch[i][]*1ll+t[x].ch[i][]*1ll*a%mod+b*1ll*t[x].s[][i]%mod)%mod;
     }
 }
 inline void dec(int x,int l){
     if (a==) {
     for (int i=;i+;i--) t[x].s[l][i+]=t[x].s[l][i];
     t[x].s[l][]=(t[x].s[l][]-b*1ll*(t[x].r-t[x].l+)%mod+mod)%mod;
         if (l==) {
         for (int i=;i+;i--)
         for (int j=;j<;j++) t[x].ch[i+][j]=t[x].ch[i][j];
         for (int i=;i<;i++)
         t[x].ch[][i]=(t[x].ch[][i]-b*1ll*t[x].s[][i]%mod+mod)%mod;
     }
     if (l==) {
         for (int i=;i+;i--)
         for (int j=;j<;j++) t[x].ch[j][i+]=t[x].ch[j][i];
         for (int i=;i<;i++)
         t[x].ch[i][]=(t[x].ch[i][]-b*1ll*t[x].s[][i]%mod+mod)%mod;
     }
     return;
     }
     for (int i=;i+;i--) t[x].s[l][i+]=t[x].s[l][i];
     t[x].s[l][]=((t[x].s[l][]*1ll-t[x].s[l][]-b*1ll*(t[x].r-t[x].l+)%mod)%mod*inv%mod+mod)%mod;
     if (l==) {
         for (int i=;i+;i--)
         for (int j=;j<;j++) t[x].ch[i+][j]=t[x].ch[i][j];
     for (int i=;i<;i++)
         t[x].ch[][i]=((t[x].ch[][i]*1ll-t[x].ch[][i]-b*1ll*t[x].s[][i]%mod)%mod*inv%mod+mod)%mod;
     }
     if (l==) {
         for (int i=;i+;i--)
         for (int j=;j<;j++) t[x].ch[j][i+]=t[x].ch[j][i];
     for (int i=;i<;i++)
         t[x].ch[i][]=((t[x].ch[i][]*1ll-t[x].ch[i][]-b*1ll*t[x].s[][i]%mod)%mod*inv%mod+mod)%mod;
     }
 }
 inline void pushback(int x,int y,int z) {
     if (z>)
     for (int i=;i<=z;i++) add(x,y);
     if (z<)
     for (int i=-;i>=z;i--) dec(x,y);
 }
 inline void pb(int x) {
     for (int l=;l<=;l++){
     pushback(lc,l,t[x].lz[l]);
     pushback(rc,l,t[x].lz[l]);
     t[lc].lz[l]+=t[x].lz[l];
     t[rc].lz[l]+=t[x].lz[l];
     t[x].lz[l]=;
     }
 }
 int A[maxn][];
 inline void build (int x,int l,int r) {
     t[x].l=l,t[x].r=r;
     if (l==r) {
     for (int i=;i<;i++) {
         t[x].s[][i]=A[l-][i+];
         t[x].s[][i]=A[l+][i+];
     }
     for (int i=;i<;i++)
         for (int j=;j<;j++)
         t[x].ch[i][j]=t[x].s[][i]*1ll*t[x].s[][j]%mod;
     return ;
     }
     build(lc,l,mid);build(rc,mid+,r);
     update(x);
 }
 inline void add(int x,int x1,int y1,int y,int z) {
     int l=t[x].l,r=t[x].r;
     if (l>y1||r<x1) return ;
     if (x1<=l&&r<=y1) {
     t[x].lz[y]+=z;
     pushback(x,y,z);
     return ;
     }
     pb(x);
     add(lc,x1,y1,y,z);add(rc,x1,y1,y,z);
     update(x);
 }
 inline int que(int x,int x1,int y1) {
     int l=t[x].l,r=t[x].r;
     if (l>y1||r<x1) return ;
     if (x1<=l&&r<=y1) return t[x].ch[][];
     pb(x);
     return (que(lc,x1,y1)+que(rc,x1,y1))%mod;
 }
 struct marix{
     int r,c,a[][];
     marix(){r=c=,memset(a,,sizeof(a));}
     void init(){r=c=;for (int i=;i<;i++) a[i][i]=;}
     void print(){
     printf("%d %d\n",r,c);
     for (int i=;i<r;i++,puts(""))
         for (int j=;j<c;j++) printf("%d ",a[i][j]);
     }
 };
 marix operator *(marix x,marix y) {
     marix ans;
     ans.r=x.r;ans.c=y.c;
     for (int i=;i<ans.r;i++)
     for (int j=;j<;j++)
         for (int k=;k<;k++)
         (ans.a[i][j]+=x.a[i][k]*1ll*y.a[k][j]%mod)%=mod;
     ans.a[][]=;
     return ans;
 }
 marix f[];
 marix pow(int y) {
     marix ans;
     int x=;
     ans.init();
     for (;y;y>>=,x++)
     if (y&) ans=ans*f[x];
     return ans;
 }
 marix mx,my;
 int get(int x){
     marix tmp=my*pow(A[x][]-);
     A[x][]=tmp.a[][];
     A[x][]=tmp.a[][];
 }
 int main(){
     int n,q;
     scanf("%d%d%d%d",&n,&q,&a,&b);
     inv=power(a,mod-);
     mx.r=mx.c=;mx.a[][]=,mx.a[][]=a,mx.a[][]=b;
     mx.a[][]=,mx.a[][]=;
     for (int i=;i<=;i++) {
     f[i]=mx;
     mx=mx*mx;
     }
     my.r=;my.c=;
     my.a[][]=,my.a[][]=,my.a[][]=;
     for (int i=;i<=n;i++) {
     scanf("%d",A[i]);
     get(i);
     A[i][]=(A[i][]*1ll+A[i][]*1ll*a+b)%mod;
     }
     build(,,n-);
     char opt[];
     while (q--) {
     int l,r;
     scanf("%s",opt);
     switch(opt[]) {
         case 'p':
         scanf("%d%d",&l,&r);
         add(,l+,r+,,);
         add(,l-,r-,,);
         break;
         case 'm':
         scanf("%d%d",&l,&r);
         add(,l+,r+,,-);
         add(,l-,r-,,-);
         break;
         case 'q':
         scanf("%d%d",&l,&r);
         printf("%d\n",que(,l+,r-));
         break;
     }
     }
     return ;
 }