//Accepted    18440 KB    5556 ms
/*
    source:HYSBZ 2243
    time  :2015.5.29
    by    :songt
  */
/*题解:
树链剖分

  */

#include <cstdio>
#include <cstring>

;
struct Edge
{
    int u,v;
    Edge(){}
    Edge(int u,int v):u(u),v(v){}
}edge[*imax_n];
int head[imax_n];
*imax_n];
int tot;
int top[imax_n];
int son[imax_n];
int num[imax_n];
int p[imax_n];
int fp[imax_n];
int deep[imax_n];
int fa[imax_n];
int pos;

void init()
{
    memset(head,-,sizeof(head));
    memset(next,-,sizeof(next));
    tot=;
    memset(son,-,sizeof(son));
    pos=;
}

void addEdge(int u,int v)
{
    edge[tot]=Edge(u,v);
    next[tot]=head[u];
    head[u]=tot++;
}

void dfs1(int u)
{
    num[u]=;
    ;i=next[i])
    {
        int v=edge[i].v;
        if (v!=fa[u])
        {
            fa[v]=u;
            deep[v]=deep[u]+;
            dfs1(v);
            num[u]+=num[v];
             || num[son[u]]<num[v])
                son[u]=v;
        }
    }
}

void dfs2(int u)
{
    p[u]=pos++;
    fp[p[u]]=u;
    ) return ;
    top[son[u]]=top[u];
    dfs2(son[u]);
    ;i=next[i])
    {
        int v=edge[i].v;
        if (v!=son[u] && v!=fa[u])
        {
            dfs2(top[v]=v);
        }
    }
}

struct Tree
{
    int l,r;
    int lc,rc;
    int num;
    int same;
    int color;
}f[imax_n*];
int color[imax_n];

void push_down(int t)
{
    f[*t].same=f[*t+].same=;
    f[*t].color=f[*t+].color=f[t].color;
    f[*t].num=f[*t+].num=;
    f[*t].lc=f[*t].rc=f[t].color;
    f[*t+].lc=f[*t+].rc=f[t].color;
    f[t].same=;
}

void push_up(int t)
{
    f[t].num=f[*t].num+f[*t+].num-(f[*t].rc==f[*t+].lc);
    f[t].lc=f[*t].lc;
    f[t].rc=f[*t+].rc;
}

void build(int t,int l,int r)
{
    f[t].l=l;
    f[t].r=r;
    f[t].same=;
    if (l==r)
    {
        f[t].lc=f[t].rc=color[fp[l]];
        f[t].color=color[fp[l]];
        f[t].num=;
        f[t].same=;
        return ;
    }
    ;
    build(*t,l,mid);
    build(*t+,mid+,r);
    push_up(t);
}

void update(int t,int l,int r,int color)
{
    if (f[t].l==l && f[t].r==r)
    {
        f[t].color=color;
        f[t].num=;
        f[t].lc=f[t].rc=color;
        f[t].same=;
        return ;
    }
    if (f[t].same) push_down(t);
    ;
    *t,l,r,color);
    else
    {
        *t+,l,r,color);
        else
        {
            update(*t,l,mid,color);
            update(*t+,mid+,r,color);
        }
    }
    push_up(t);
}

int query(int t,int l,int r,int &cl,int &rl)
{
    if (f[t].l==l && f[t].r==r)
    {
        cl=f[t].lc;
        rl=f[t].rc;
        return f[t].num;
    }
    if (f[t].same) push_down(t);
    ;
    *t,l,r,cl,rl);
    else
    {
        *t+,l,r,cl,rl);
        else
        {
            int lcl,lcr,rcl,rcr;
            int numl,numr;
            numl=query(*t,l,mid,lcl,lcr);
            numr=query(*t+,mid+,r,rcl,rcr);
            cl=lcl;
            rl=rcr;
            return numl+numr-(lcr==rcl);
        }
    }
}

void swap(int &a,int &b)
{
    int t=a;
    a=b;
    b=t;
}

void OpC(int u,int v,int c)
{
    int f1=top[u],f2=top[v];
    while (f1!=f2)
    {
        //printf("u=%d v=%d top[u]=%d top[v]=%d\n",u,v,f1,f2);
        if (deep[f1]<deep[f2])
        {
            swap(f1,f2);
            swap(u,v);
        }
        update(,p[f1],p[u],c);
        //printf("update %d %d\n",f1,u);
        u=fa[f1];
        f1=top[u];
    }
    if (deep[u]>deep[v]) swap(u,v);
    update(,p[u],p[v],c);
    //printf("update %d %d\n",u,v);
}

int OpQ(int u,int v)
{
    int f1=top[u],f2=top[v];
    ,cv=-;
    ;
    int lc,rc;
    while (f1!=f2)
    {
        if (deep[f1]<deep[f2])
        {
            swap(f1,f2);
            swap(u,v);
            swap(cu,cv);
        }
        ,p[f1],p[u],lc,rc);
        //printf("%d %d num=%d lc=%d rc=%d\n",f1,u,tmp,lc,rc);
        ans+=tmp;
        ans-=(cu==rc);
        cu=lc;
        u=fa[f1];
        f1=top[u];
    }
    if (deep[u]>deep[v])
    {
        swap(u,v);
        swap(cu,cv);
    }
    ,p[u],p[v],lc,rc);
    //printf("%d %d num=%d lc=%d rc=%d\n",u,v,tmp,lc,rc);
    ans+=query(,p[u],p[v],lc,rc);
    ans-=(cu==lc);
    ans-=(cv==rc);
    return ans;
}

int n,m;
];
int u,v,c;

int main()
{
    //while (scanf("%d%d",&n,&m)==2)
    scanf("%d%d",&n,&m);
    {
        init();
        ;i<=n;i++)
        {
            scanf("%d",&color[i]);
        }
        ;i<n-;i++)
        {
            scanf("%d%d",&u,&v);
            addEdge(u,v);
            addEdge(v,u);
        }
        fa[]=;
        deep[]=;
        dfs1();
        dfs2(top[]=);
        build(,,pos-);
        ;i<m;i++)
        {
            scanf("%s",op);
            ]=='Q')
            {
                scanf("%d%d",&u,&v);
                printf("%d\n",OpQ(u,v));
            }
            else
            {
                scanf("%d%d%d",&u,&v,&c);
                OpC(u,v,c);
            }
        }
    }
    ;
}

HYSBZ 2243的更多相关文章

  1. HDU 3966 & POJ 3237 & HYSBZ 2243 树链剖分

    树链剖分是一个很固定的套路 一般用来解决树上两点之间的路径更改与查询 思想是将一棵树分成不想交的几条链 并且由于dfs的顺序性 给每条链上的点或边标的号必定是连着的 那么每两个点之间的路径都可以拆成几 ...

  2. HDU 3966 & POJ 3237 & HYSBZ 2243 & HRBUST 2064 树链剖分

    树链剖分是一个很固定的套路 一般用来解决树上两点之间的路径更改与查询 思想是将一棵树分成不想交的几条链 并且由于dfs的顺序性 给每条链上的点或边标的号必定是连着的 那么每两个点之间的路径都可以拆成几 ...

  3. hysbz 2243 染色(树链剖分)

    题目链接:hysbz 2243 染色 题目大意:略. 解题思路:树链剖分+线段树的区间合并,可是区间合并比較简单,节点仅仅要记录左右端点的颜色就可以. #include <cstdio> ...

  4. HYSBZ 2243(树链剖分)

    题目连接:http://acm.hust.edu.cn/vjudge/contest/view.action?cid=28982#problem/D 题意:给定一棵有n个节点的无根树及点权和m个操作, ...

  5. HYSBZ 2243 染色 (树链拆分)

    主题链接~~> 做题情绪:这题思路好想.调试代码调试了好久.第一次写线段树区间合并. 解题思路: 树链剖分 + 线段树区间合并 线段树的端点记录左右区间的颜色.颜色数目.合并的时候就用区间合并的 ...

  6. HYSBZ 2243(染色)

    题目链接:传送门 题目大意:中文题,略 题目思路:树链剖分,区间更新,区间查询. 闲谈:      只想说这道题做的好苦逼..去长春现场赛之前就没A,回来后又做了2天才A掉,蒟蒻太菜了 这道题也没有想 ...

  7. HYSBZ - 2243 染色 (树链剖分+线段树)

    题意:树上每个结点有自己的颜色,支持两种操作:1.将u到v路径上的点颜色修改为c; 2.求u到v路径上有多少段不同的颜色. 分析:树剖之后用线段树维护区间颜色段数.区间查询区间修改.线段树结点中维护的 ...

  8. HYSBZ - 2243 树链剖分 + 线段树 处理树上颜色段数

    用线段树处理颜色段数 记录区间内的颜色段数,区间右端点的颜色,区间右端点的颜色. int tr[maxn<<2], lc[maxn<<2], rc[maxn<<2] ...

  9. HYSBZ 2038 莫队算法

    小Z的袜子(hose) Time Limit:20000MS     Memory Limit:265216KB     64bit IO Format:%lld & %llu Submit  ...

随机推荐

  1. iOS 动画组

    其实早在一个多月以前就已经实现了动作组播放的功能,不过当时感觉好像没有什么难度并没有放在心上,今天突然要用到动画组,发现已经忘记了,所以又将原来的代码翻出来看了下.最后决定还是写下来,以备不时之需.动 ...

  2. javascript实现简单的轮播图片

    方法一: <script> var curIndex=0;//时间间隔(单位毫秒),每秒钟显示一张,数组共有5张图片放在Photos文件夹下. var timeInterval=1000; ...

  3. Leetcode--Generate Parentheses

    主要考察栈的理解 static vector<string> generateParenthesis(int n) { vector<string> res; addingpa ...

  4. python——赋值与深浅拷贝

    初学编程的小伙伴都会对于深浅拷贝的用法有些疑问,今天我们就结合python变量存储的特性从内存的角度来谈一谈赋值和深浅拷贝~~~ 预备知识一——python的变量及其存储 在详细的了解python中赋 ...

  5. linux 挂载光盘:mount: you must specify the filesystem type

    尝试挂载光盘镜像时出现mount: you must specify the filesystem type 使用-t auto -t iso9660 或不加参数都搞不定,最后在以下链接找到解决办法: ...

  6. FileFilter 遍历某个目录下文件名含有某个字符的文件

    由于IIS版本的升级,造成了文件名中含有“+”的特殊字符的文件(多数是图片)在网页中不能被访问,于是必须查找当前目录下含有多少这样的文件,从而制定最佳的解决方案. 废话少说,直接上核心代码: publ ...

  7. WeCenter二次开发教程(一):熟悉模板结构

    <1>程序文件目录介绍: app – 应用目录 models – 模型目录 plugins – 插件目录 static – 静态文件 system – 系统目录 views – 模板目录 ...

  8. html slelect 标签默认值

    <select name="channelCode" id="channelCode" class="all_input" style ...

  9. 关于小组所要做的APP的想法

    关于小组所要做的app,我们敲定下来是做关于在线做题的app,但是,纯粹的做题目的app我认为并没有什么大的吸引力,尤其是拿手机做题.所以,我们考虑准备在以下几个方面做功夫以增加吸引力.第一,我们的题 ...

  10. Transport Block Size, Throughput and Code rate-----http://www.simpletechpost.com/2012/12/transport-block-size-code-rate-protocol.html

    Transport Block Size, Throughput and Code rate   Since the size of transport block is not fixed, oft ...