【BZOJ4999】This Problem Is Too Simple!(线段树)

39 篇文章 0 订阅

题面

BZOJ

题解

对于每个值,维护一棵线段树就好啦

动态开点,否则空间开不下

剩下的就是很简单的问题啦

当然了,对于数值要离散化

没必要离线吧,在线用 map 维护就行了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define RG register
#define MAX 120000
#define lson (t[x].ls)
#define rson (t[x].rs)
inline int read()
{
    RG int x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
struct Line{int v,next;}e[MAX<<1];
int h[MAX],cnt=1;
inline void Add(int u,int v){e[cnt]=(Line){v,h[u]};h[u]=cnt++;}
map<int,int> M;
int tim,dfn[MAX],hson[MAX],size[MAX],fa[MAX],top[MAX],C[MAX],dep[MAX];
void dfs1(int u,int ff)
{
    fa[u]=ff;size[u]=1;dep[u]=dep[ff]+1;
    for(int i=h[u];i;i=e[i].next)
    {
        int v=e[i].v;
        if(v==ff)continue;
        dfs1(v,u);
        if(size[v]>size[hson[u]])hson[u]=v;
        size[u]+=size[v];
    }
}
void dfs2(int u,int tp)
{
    top[u]=tp;dfn[u]=++tim;
    if(hson[u])dfs2(hson[u],tp);
    for(int i=h[u];i;i=e[i].next)
        if(e[i].v!=fa[u]&&e[i].v!=hson[u])
            dfs2(e[i].v,e[i].v);
}
struct Node{int ls,rs,v;}t[MAX<<7];
int tot,sum,rt[MAX<<2];
int n,Q;
void Modify(int &x,int l,int r,int p,int w)
{
    if(!x)x=++tot;t[x].v+=w;
    if(l==r)return;
    int mid=(l+r)>>1;
    if(p<=mid)Modify(lson,l,mid,p,w);
    else Modify(rson,mid+1,r,p,w);
}
int Query(int x,int l,int r,int L,int R)
{
    if((!t[x].v)||(!x))return 0;
    if(L<=l&&r<=R)return t[x].v;
    int mid=(l+r)>>1,ret=0;
    if(L<=mid)ret+=Query(lson,l,mid,L,R);
    if(R>mid)ret+=Query(rson,mid+1,r,L,R);
    return ret;
}
int Route(int u,int v,int x)
{
    int ret=0;
    while(top[u]!=top[v])
    {
        if(dep[top[u]]<dep[top[v]])swap(u,v);
        ret+=Query(rt[x],1,n,dfn[top[u]],dfn[u]);
        u=fa[top[u]];
    }
    if(dep[u]<dep[v])swap(u,v);
    ret+=Query(rt[x],1,n,dfn[v],dfn[u]);
    return ret;
}
int main()
{
    n=read();Q=read();
    for(int i=1;i<=n;++i)C[i]=read();
    for(int i=1;i<n;++i)
    {
        int u=read(),v=read();
        Add(u,v);Add(v,u);
    }
    dfs1(1,0);dfs2(1,1);
    for(int i=1;i<=n;++i)
    {
        if(!M[C[i]])M[C[i]]=++sum;
        Modify(rt[M[C[i]]],1,n,dfn[i],1);
    }
    char ch[5];
    while(Q--)
    {
        scanf("%s",ch);
        if(ch[0]=='C')
        {
            int u=read(),x=read();
            Modify(rt[M[C[u]]],1,n,dfn[u],-1);
            if(!M[x])M[x]=++sum;
            Modify(rt[M[x]],1,n,dfn[u],1);
            C[u]=x;
        }
        else
        {
            int u=read(),v=read(),x=read();
            if(!M[x])puts("0");
            else printf("%d\n",Route(u,v,M[x]));
        }
    }
    return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值