hdu5068 线段树

因为楼梯的状态只在两个门之间有效，对楼梯两侧未合并前的方法书无影响，因此很容易想到区间合并，然后就是记录区间的不同状态的方法数，以及向后延伸的接口状态，然后单点更新区间合并就行了。

附代码：

#include <iostream>
#include <cstring>
#include <cstdio>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define delf int m=(l+r)>>1
#define ll long long int

using namespace std;

const int MOD=1000000007;
const int MAX=50050;
ll sum[MAX<<2][2][2];       //区间起点是i，终点是t的方法数
ll st[MAX<<2][2][2];        //区间末端的接口，末端为i连接下一个门为t的楼梯的状态

struct node
{
    ll sum[2][2];
    ll st[2][2];
};

void pushup(int l,int r,int rt)
{
    delf;
    for (int i=0;i<=1;i++)
    {
        for (int t=0;t<=1;t++)
        {
            sum[rt][i][t]=0;
            sum[rt][i][t]+=(((sum[rt<<1][i][0]*st[rt<<1][0][0])%MOD)*sum[rt<<1|1][0][t])%MOD;
            sum[rt][i][t]+=(((sum[rt<<1][i][0]*st[rt<<1][0][1])%MOD)*sum[rt<<1|1][1][t])%MOD;
            sum[rt][i][t]+=(((sum[rt<<1][i][1]*st[rt<<1][1][0])%MOD)*sum[rt<<1|1][0][t])%MOD;
            sum[rt][i][t]+=(((sum[rt<<1][i][1]*st[rt<<1][1][1])%MOD)*sum[rt<<1|1][1][t])%MOD;
            st[rt][i][t]=st[rt<<1|1][i][t];
        }
    }
    return ;
}

void build(int l,int r,int rt)
{
    sum[rt][0][0]=sum[rt][0][1]=sum[rt][1][0]=sum[rt][1][1]=0;
    if (l==r)
    {
        for (int i=0;i<=1;i++)
            for (int t=0;t<=1;t++)
                sum[rt][i][i]=st[rt][i][t]=1;
        return ;
    }
    delf;
    build (lson);
    build (rson);
    pushup(l,r,rt);
    return ;
}

void update(int x,int k,int t,int l,int r,int rt)
{
    if (l==r)
    {
        st[rt][k][t]=1-st[rt][k][t];
        return ;
    }
    delf;
    if (x<=m)
        update(x,k,t,lson);
    else
        update(x,k,t,rson);
    pushup(l,r,rt);
    return ;
}

node query(int L,int R,int l,int r,int rt)
{
    if (L<=l&&r<=R)
    {
        node a;
        for (int i=0;i<=1;i++)
            for (int t=0;t<=1;t++)
            {
                a.st[i][t]=st[rt][i][t];
                a.sum[i][t]=sum[rt][i][t];
            }
        return a;
    }
    delf;
    if (R<=m)
        return query(L,R,lson);
    if (L>m)
        return query(L,R,rson);
    node a1=query(L,R,lson);
    node a2=query(L,R,rson);
    node a;
    for (int i=0;i<=1;i++)
    {
        for (int t=0;t<=1;t++)
        {
            a.sum[i][t]=0;
            a.sum[i][t]+=(((a1.sum[i][0]*a1.st[0][0])%MOD)*a2.sum[0][t])%MOD;
            a.sum[i][t]+=(((a1.sum[i][0]*a1.st[0][1])%MOD)*a2.sum[1][t])%MOD;
            a.sum[i][t]+=(((a1.sum[i][1]*a1.st[1][0])%MOD)*a2.sum[0][t])%MOD;
            a.sum[i][t]+=(((a1.sum[i][1]*a1.st[1][1])%MOD)*a2.sum[1][t])%MOD;
            a.sum[i][t]%=MOD;
            a.st[i][t]=a2.st[i][t];
        }
    }
    return a;
}

int main()
{
    int n,m;
    while (~scanf("%d%d",&n,&m))
    {
        build(1,n,1);
        for (int i=1;i<=m;i++)
        {
            int op,l,r,v;
            scanf("%d",&op);
            if (op==0)
            {
                scanf("%d%d",&l,&r);
                node a=query(l,r,1,n,1);
                ll ans=0;
                for (int i=0;i<=1;i++)
                    for (int t=0;t<=1;t++)
                        ans+=a.sum[i][t];
                ans%=MOD;
                printf("%I64d\n",ans);
            }
            else
            {
                scanf("%d%d%d",&l,&r,&v);
                update(l,r-1,v-1,1,n,1);
            }
        }
    }
}