Codeforces 718C 线段树+矩乘

题意：

维护一个序列，支持两种操作：
1.区间[l,r]的权值+x
2.询问区间[l,r]的函数和，即∑fib(x)这里的函数即斐波那契函数
数据范围：1≤n,q≤10⁵

思路：
一般求斐波那契函数的方法可以考虑矩阵乘法，这里也是这样的。
我们不用线段树维护权值，我们用线段树维护线段树维护区间矩阵和。
有一个矩阵乘法的性质：A*B+A*C=A*(B+C)
在求斐波那契数列中，是A*F，A是变换矩阵，F是列矩阵
那么我们用线段树的lazy标记维护A矩阵，然后用sum维护F矩阵
之后在线段树上，就变成了区间更新乘以x。

//By SiriusRen
#include <bits/stdc++.h>
using namespace std;
const int mod=,N=;
int n,m,op,xx,yy,zz;
struct Matrix{
    int a[][];
    void init(){memset(a,,sizeof(a));}
    void dia(){a[][]=a[][]=,a[][]=a[][]=;}
}sum[N<<],lazy[N<<],base;
Matrix operator*(Matrix a,Matrix b){
    Matrix c;c.init();
    for(int i=;i<;i++)
        for(int j=;j<;j++)
            for(int k=;k<;k++)
                c.a[i][j]=(1ll*a.a[i][k]*b.a[k][j]+c.a[i][j])%mod;
    return c;
}
Matrix operator+(Matrix a,Matrix b){
    Matrix c;c.init();
    for(int i=;i<;i++)
        for(int j=;j<;j++)
            c.a[i][j]=(a.a[i][j]+b.a[i][j])%mod;
    return c;
}
Matrix operator^(Matrix a,int b){
    Matrix c;c.dia();
    while(b){
        if(b&)c=c*a;
        a=a*a,b>>=;
    }return c;
}
void push_up(int pos){sum[pos]=sum[pos<<]+sum[pos<<|];}
void build(int l,int r,int pos){
    if(l==r){scanf("%d",&xx);sum[pos]=base^(xx-);return;}
    int mid=(l+r)>>,lson=pos<<,rson=pos<<|;
    build(l,mid,lson),build(mid+,r,rson);
    lazy[pos].dia(),push_up(pos);
}
void push_down(int pos,int lson,int rson){
    lazy[lson]=lazy[lson]*lazy[pos],sum[lson]=sum[lson]*lazy[pos];
    lazy[rson]=lazy[rson]*lazy[pos],sum[rson]=sum[rson]*lazy[pos];
    lazy[pos].dia();
}
Matrix query(int l,int r,int pos,int L,int R){
    if(l>=L&&r<=R)return sum[pos];
    int mid=(l+r)>>,lson=pos<<,rson=pos<<|;
    push_down(pos,lson,rson);
    if(mid<L)return query(mid+,r,rson,L,R);
    else if(mid>=R)return query(l,mid,lson,L,R);
    else return query(l,mid,lson,L,R)+query(mid+,r,rson,L,R);
}
void insert(int l,int r,int pos,int L,int R,Matrix wei){
    if(l>=L&&r<=R){sum[pos]=sum[pos]*wei;lazy[pos]=lazy[pos]*wei;return;}
    int mid=(l+r)>>,lson=pos<<,rson=pos<<|;
    push_down(pos,lson,rson);
    if(mid<L)insert(mid+,r,rson,L,R,wei);
    else if(mid>=R)insert(l,mid,lson,L,R,wei);
    else insert(l,mid,lson,L,R,wei),insert(mid+,r,rson,L,R,wei);
    push_up(pos);
}
int main(){
    base.a[][]=base.a[][]=base.a[][]=;
    scanf("%d%d",&n,&m),build(,n,);
    while(m--){
        scanf("%d%d%d",&op,&xx,&yy);
        if(op==)scanf("%d",&zz),insert(,n,,xx,yy,base^zz);
        else printf("%d\n",query(,n,,xx,yy).a[][]);
    }
}