bzoj 4999: This Problem Is Too Simple！ 树链剖分+线段树

题意

给您一颗树，每个节点有个初始值。
现在支持以下两种操作：
1. C i x(0<=x<2^31) 表示将i节点的值改为x。
2. Q i j x(0<=x<2^31) 表示询问i节点到j节点的路径上有多少个值为x的节点。

n<=100000,q<=200000

分析

像这种题，一般的套路就是对每种权值开一棵线段树然后直接树链剖分即可。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

const int N=100005;

int n,m,cnt,tim,sz,root[N*3],last[N],val[N],w[N*3],size[N],dep[N],fa[N],top[N],pos[N];
struct edge{int to,next;}e[N*2];
struct data{int op,i,j,x;}q[N*2];
struct tree{int s,l,r;}t[N*60];

int read()
{
    int x=0,f=1;char ch=getchar();
    while (ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while (ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}

void addedge(int u,int v)
{
    e[++cnt].to=v;e[cnt].next=last[u];last[u]=cnt;
    e[++cnt].to=u;e[cnt].next=last[v];last[v]=cnt;
}

void dfs1(int x)
{
    dep[x]=dep[fa[x]]+1;size[x]=1;
    for (int i=last[x];i;i=e[i].next)
    {
        if (e[i].to==fa[x]) continue;
        fa[e[i].to]=x;
        dfs1(e[i].to);
        size[x]+=size[e[i].to];
    }
}

void dfs2(int x,int chain)
{
    pos[x]=++tim;top[x]=chain;int k=0;
    for (int i=last[x];i;i=e[i].next)
        if (e[i].to!=fa[x]&&size[e[i].to]>size[k]) k=e[i].to;
    if (!k) return;
    dfs2(k,chain);
    for (int i=last[x];i;i=e[i].next)
        if (e[i].to!=fa[x]&&e[i].to!=k) dfs2(e[i].to,e[i].to);
}

void ins(int &d,int l,int r,int x,int y)
{
    if (!d) d=++sz;
    t[d].s+=y;
    if (l==r) return;
    int mid=(l+r)/2;
    if (x<=mid) ins(t[d].l,l,mid,x,y);
    else ins(t[d].r,mid+1,r,x,y);
}

int Sig_query(int d,int l,int r,int x,int y)
{
    if (x>y||!d) return 0;
    if (l==x&&r==y) return t[d].s;
    int mid=(l+r)/2;
    return Sig_query(t[d].l,l,mid,x,min(y,mid))+Sig_query(t[d].r,mid+1,r,max(x,mid+1),y);
}

int query(int x,int y,int id)
{
    int ans=0;
    while (top[x]!=top[y])
    {
        if (dep[top[x]]<dep[top[y]]) swap(x,y);
        ans+=Sig_query(root[id],1,n,pos[top[x]],pos[x]);
        x=fa[top[x]];
    }
    if (dep[x]<dep[y]) swap(x,y);
    ans+=Sig_query(root[id],1,n,pos[y],pos[x]);
    return ans;
}

int main()
{
    n=read();m=read();
    for (int i=1;i<=n;i++) val[i]=read(),w[i]=val[i];
    for (int i=1;i<n;i++)
    {
        int x=read(),y=read();
        addedge(x,y);
    }
    dfs1(1);
    dfs2(1,1);
    for (int i=1;i<=m;i++)
    {
        char ch[2];
        scanf("%s",ch);
        if (ch[0]=='C') q[i].op=0,q[i].i=read(),q[i].x=read(),w[i+n]=q[i].x;
        else q[i].op=1,q[i].i=read(),q[i].j=read(),q[i].x=read(),w[i+n]=q[i].x;
    }
    sort(w+1,w+n+m+1);
    int w1=unique(w+1,w+n+m+1)-w-1;
    for (int i=1;i<=n;i++) val[i]=lower_bound(w+1,w+w1+1,val[i])-w;
    for (int i=1;i<=m;i++) q[i].x=lower_bound(w+1,w+w1+1,q[i].x)-w;
    for (int i=1;i<=n;i++) ins(root[val[i]],1,n,pos[i],1);
    for (int i=1;i<=m;i++)
        if (!q[i].op)
        {
            ins(root[val[q[i].i]],1,n,pos[q[i].i],-1);
            val[q[i].i]=q[i].x;
            ins(root[val[q[i].i]],1,n,pos[q[i].i],1);
        }
        else
        {
            printf("%d\n",query(q[i].i,q[i].j,q[i].x));
        }
    return 0;
}