[HDU6155]Subsequence Count(线段树+矩阵)

DP式很容易得到，发现是线性递推形式，于是可以矩阵加速。又由于是区间形式，所以用线段树维护。

https://www.cnblogs.com/Miracevin/p/9124511.html

关键在于证明区间操作中，可以直接在打标记的位置翻转矩阵两行两列。

上面网址用代数形式证了一遍，这里考虑从矩阵本身解释。

由线代内容可知，将一个矩阵作初等行变换，相当于将其左乘一个作了相应初等列变换的单位矩阵。同理将一个矩阵作初等列变换，相当于将其又乘一个作了相应初等行变换的单位矩阵。

这里，左乘的矩阵$T=egin{bmatrix}0 & 1 & 0 \ 1 & 0 & 0 \ 0 & 0 & 1end{bmatrix}$，右乘的矩阵$T'$同样也是这个。

我们发现，$T imes T'$就是单位矩阵，也就是说$T$的逆矩阵就是自己。

于是有，$T imes A imes T' imes T imes B imes T'=T imes A imes B imes T'$。

这就说明中间的所有矩乘操作都可以被省略，只留下首尾的$T$和$T'$。

这也就证明了，对区间矩阵积直接做变换是正确的。

#include<cstdio>
#include<cstring>
#include<algorithm>
#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=100010,mod=1e9+7;
bool tag[N<<2];
int n,Q,T,op,l,r;

struct Mat{ int a[3][3]; }v[N<<2];
const Mat D[2]={(Mat){1,0,0,1,1,0,1,0,1},(Mat){1,1,0,0,1,0,0,1,1}};
int calc(Mat a){ return (a.a[2][0]+a.a[2][1])%mod; }

Mat operator *(const Mat &a,const Mat &b){
    Mat c; memset(c.a,0,sizeof(c.a));
    rep(i,0,2) rep(j,0,2) rep(k,0,2)
        c.a[i][k]=(c.a[i][k]+1ll*a.a[i][j]*b.a[j][k])%mod;
    return c;
}

void put(int x){
    tag[x]^=1;
    swap(v[x].a[0][0],v[x].a[1][0]);
    swap(v[x].a[0][1],v[x].a[1][1]);
    swap(v[x].a[0][2],v[x].a[1][2]);
    swap(v[x].a[0][0],v[x].a[0][1]);
    swap(v[x].a[1][0],v[x].a[1][1]);
    swap(v[x].a[2][0],v[x].a[2][1]);
}

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

void build(int x,int L,int R){
    tag[x]=0;
    if (L==R){ char t; scanf(" %c",&t); v[x]=D[t-'0']; return; }
    int mid=(L+R)>>1; build(lson); build(rson); v[x]=v[ls]*v[rs];
}

void mdf(int x,int L,int R,int l,int r){
    if (L==l && r==R){ put(x); return; }
    int mid=(L+R)>>1; push(x);
    if (r<=mid) mdf(lson,l,r);
    else if (l>mid) mdf(rson,l,r);
        else mdf(lson,l,mid),mdf(rson,mid+1,r);
    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)>>1; 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+1,r);
}

int main(){
    freopen("hdu6155.in","r",stdin);
    freopen("hdu6155.out","w",stdout);
    for (scanf("%d",&T); T--; ){
        scanf("%d%d",&n,&Q); build(1,1,n);
        while (Q--){
            scanf("%d%d%d",&op,&l,&r);
            if (op==1) mdf(1,1,n,l,r); else printf("%d
",calc(que(1,1,n,l,r)));
        }
    }
    return 0;
}