codechef FIBTREE 码农题 线段树 树剖 标记永久化

好烦啊，调了半天

线段树部分标记比较多，手抖打错了一个

剩下的都是取模的问题

我自己瞎jb推的公式里保留了abs，但是在模意义下是gg的，所以必须把正负区分开

调试的时候一定要注意构造各种形状的树，不要只做随机树

随机树深度只有log，很难体现一些链上的性质

我用随机树拍了一下午没出错，一掏出直链就秒秒钟出错

最后找到了那个该死的abs

还是逻辑不够严谨啊

#include <bits/stdc++.h>
#define DEBUG 0
#define mid (l+r>>1)
#define MOD 1000000009
#define Que(rt,x,y) que(rt,1,n,x,y)
#define Len (dep[x]-dep[top[x]]+1)
using namespace std;
struct node
{
    int f,s,F,S,sum;
    node(long long a,long long b,long long c,long long d,long long e)
    {
        f=a%MOD;s=b%MOD;F=c%MOD;S=d%MOD;sum=e%MOD;
    }
    node()
    {
    }
} tr[40000001];
int ls[40000001],rs[40000001];
int TIME,E,NODE,n,m,p,q,x,y;char ch;bool rev;
int pos[200001],start[200001],en[200001];
int dep[200001],top[200001],fa[200001],root[200001];
int to[400001],nex[400001],fir[200001],size[200001];
long long f[200001];
int son(int x,int y)
{
    while(dep[fa[top[y]]]>dep[x]) y=fa[top[y]];
    if(fa[top[y]]==x) return top[y];
    else return pos[start[x]+1];
} 
void add(int p,int q)
{
    to[++E]=q;nex[E]=fir[p];fir[p]=E;
}
int build(int now,int fat)
{
    fa[now]=fat;
    dep[now]=dep[fat]+1;
    size[now]=1;
    for(int i=fir[now];i;i=nex[i])
    if(to[i]!=fat) 
        size[now]+=build(to[i],now);
    return size[now];
}
void pou(int now,int Top)
{
    top[now]=Top;start[now]=++TIME;pos[TIME]=now;
    int id=0;
    for(int i=fir[now];i;i=nex[i])
    if(to[i]!=fa[now])
        if(!id || size[to[i]]>size[id]) id=to[i];
    if(id) pou(id,Top);
    for(int i=fir[now];i;i=nex[i])
    if(to[i]!=fa[now] && to[i]!=id)
        pou(to[i],to[i]);
    en[now]=TIME;
}
int lca(int x,int y)
{
    for(;top[x]!=top[y];x=fa[top[x]])
        if(fa[top[x]]<fa[top[y]]) swap(x,y);
    return dep[x]>dep[y]?y:x;
}
long long fib(int F,int S,int x)
{
    if(x==1) return F;
    if(x==2) return S;
    return (f[x-2]*F+f[x-1]*S)%MOD;
}
long long getsum(int now,int x,int y)
{
    if(y==1) return (tr[now].f+tr[now].F)%MOD;
    if(x==1 && y==2) return ((tr[now].f+tr[now].F)%MOD+(tr[now].s+tr[now].S)%MOD)%MOD;
    if(x==2 && y==2) return (tr[now].s+tr[now].S)%MOD;
    int F=(((y&1?1:-1)*f[y-2]*tr[now].F-(y&1?1:-1)*f[y-1]*tr[now].S)%MOD+MOD)%MOD,S=(((y&1?-1:1)*f[y-3]*tr[now].F+(y&1?1:-1)*f[y-2]*tr[now].S)%MOD+MOD)%MOD;
    return (fib(tr[now].f,tr[now].s,y+2)-fib(tr[now].f,tr[now].s,x+1)+MOD+fib(F,S,y-x+3)-S+MOD)%MOD;
}
int que(int now,int l,int r,int x,int y)
{
    if(!now) return 0;
    if(l==x && r==y)
        return tr[now].sum;
    int ret=getsum(now,x-l+1,y-l+1);
    if(x<=mid) ret=(ret+que(ls[now],l,mid,x,min(y,mid)))%MOD;
    if(y>mid) ret=(ret+que(rs[now],mid+1,r,max(mid+1,x),y))%MOD;
    return ret; 
}
int add(int acc,int l,int r,int x,int y,int st,bool rev)
{
    int now=++NODE;
    if(NODE>40000000) while(1); 
    tr[now]=tr[acc];
    if(l==x && r==y)
    {
        if(rev) tr[now]=node(tr[acc].f,tr[acc].s,tr[acc].F+f[st],tr[acc].S+f[st-1],tr[acc].sum+f[st+2]-f[st+l-r+1]+MOD);
        else tr[now]=node(tr[acc].f+f[st],tr[acc].s+f[st+1],tr[acc].F,tr[acc].S,tr[acc].sum+f[st+r-l+2]-f[st+1]+MOD);
        ls[now]=ls[acc];rs[now]=rs[acc];
        return now;
    }
    if(x<=mid && y<=mid)
        ls[now]=add(ls[acc],l,mid,x,y,st,rev),rs[now]=rs[acc];
    else
    if(y>mid && x>mid)
        ls[now]=ls[acc],rs[now]=add(rs[acc],mid+1,r,x,y,st,rev);
    else
    {
        ls[now]=add(ls[acc],l,mid,x,mid,st,rev);
        rs[now]=add(rs[acc],mid+1,r,mid+1,y,st+(rev?-1:1)*(mid-x+1),rev);
    }
    tr[now].sum=(tr[ls[now]].sum+tr[rs[now]].sum+getsum(now,1,r-l+1))%MOD;
    return now;
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<n;i++)
        scanf("%d%d",&p,&q),
        add(p,q),add(q,p);
    build(1,0);pou(1,1);
    f[1]=1;f[2]=1;
    for(int i=3;i<=n+10;i++)
    {
        f[i]=f[i-2]+f[i-1];
        if(f[i]>MOD) f[i]-=MOD;
    }
    root[0]=1;NODE=1;
    int lastans=0;
    for(int i=1;i<=m;i++)
    {
        root[i]=root[i-1];
        char S[10];
        scanf("%s",&S);
        if(S[0]=='A')
        {
            scanf("%d%d",&x,&y);
            if(!DEBUG)
            x^=lastans;
            if(x>n) while(1); 
            bool rev=0;int len=1,leng=dep[x]+dep[y]-2*dep[lca(x,y)]+1;
            for(;top[x]!=top[y];x=fa[top[x]])
            {
                if(dep[top[x]]<dep[top[y]]) swap(x,y),rev^=1;
                if(rev) leng-=Len;else len+=Len;
                root[i]=add(root[i],1,n,start[top[x]],start[x],rev?leng+1:len-1,!rev);
            }
            if(dep[x]>dep[y]) swap(x,y),rev^=1;
            root[i]=add(root[i],1,n,start[x],start[y],rev?leng:len,rev);
        } 
        if(S[0]=='R')
        {
            scanf("%d",&x);
            if(!DEBUG)
            x^=lastans;
            if(x>=i) while(1);
            root[i]=root[x];
        }
        if(S[0]=='Q')
        {
            scanf("%d%d",&x,&y);
            if(!DEBUG)
            x^=lastans;
            if(x>n) while(1);
            int t=lca(x,y);
            if(S[1]=='S')
            if(x!=y) 
            if(t==y)
            {
                int tem=son(y,x);
                lastans=(Que(root[i],1,n)-Que(root[i],start[tem],en[tem])+MOD)%MOD;
            }
            else
                lastans=Que(root[i],start[y],en[y]);
            else
                lastans=Que(root[i],1,n);
            else
            {
                lastans=0;
                for(;top[x]!=top[y];x=fa[top[x]])
                {
                    if(dep[top[x]]<dep[top[y]]) swap(x,y);
                    lastans+=Que(root[i],start[top[x]],start[x]);
                    lastans%=MOD;
                }
                if(dep[x]>dep[y]) swap(x,y);
                lastans+=Que(root[i],start[x],start[y]);
                lastans%=MOD;
            }
            printf("%d
",lastans);
        }
    }
    return 0;
}