当前位置: 首页 > news >正文

[BZOJ 3531][Sdoi2014]旅行(树链剖分+线段树)

Description

S国有N个城市,编号从1到N。城市间用N-1条双向道路连接,满足从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教,如飞天面条神教、隐形独角兽教、绝地教都是常见的信仰。为了方便,我们用不同的正整数代表各种宗教,  S国的居民常常旅行。旅行时他们总会走最短路,并且为了避免麻烦,只在信仰和他们相同的城市留宿。当然旅程的终点也是信仰与他相同的城市。S国政府为每个城市标定了不同的旅行评级,旅行者们常会记下途中(包括起点和终点)留宿过的城市的评级总和或最大值。

在S国的历史上常会发生以下几种事件:
”CC x c”:城市x的居民全体改信了c教;
”CW x w”:城市x的评级调整为w;
”QS x y”:一位旅行者从城市x出发,到城市y,并记下了途中留宿过的城市的评级总和;
”QM x y”:一位旅行者从城市x出发,到城市y,并记下了途中留宿过
的城市的评级最大值。
由于年代久远,旅行者记下的数字已经遗失了,但记录开始之前每座城市的信仰与评级,还有事件记录本身是完好的。请根据这些信息,还原旅行者记下的数字。    为了方便,我们认为事件之间的间隔足够长,以致在任意一次旅行中,所有城市的评级和信仰保持不变。

Solution

只需要对每个宗教建一棵线段树就好啦

#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#define MAXN 100005
using namespace std;
int n,q,w[MAXN],c[MAXN],head[MAXN],cnt=0;
int father[MAXN],top[MAXN],maxv[MAXN],siz[MAXN],deep[MAXN],pos[MAXN],sz=0;
int tot=0,rt[MAXN],ls[MAXN*100],rs[MAXN*100],maxn[MAXN*100],sum[MAXN*100];
int read()
{
    int x=0,f=1;char c=getchar();
    while(c<'0'||c>'9'){if(c=='-')f=-1;c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0';c=getchar();}
    return x*f;
}
struct Node1
{
    int next,to;
}Edges[MAXN*2];
void addedge(int u,int v)
{
    Edges[++cnt].next=head[u];
    head[u]=cnt;
    Edges[cnt].to=v;
}
void dfs1(int u)
{
    siz[u]=1;
    for(int i=head[u];~i;i=Edges[i].next)
    {
        int v=Edges[i].to;
        if(v==father[u])continue;
        father[v]=u,deep[v]=deep[u]+1;
        dfs1(v);
        siz[u]+=siz[v];
        if(siz[v]>siz[maxv[u]])maxv[u]=v;
    }
}
void dfs2(int u,int t)
{
    top[u]=t;
    ++sz,pos[u]=sz;
    if(maxv[u])dfs2(maxv[u],t);
    for(int i=head[u];~i;i=Edges[i].next)
    {
        int v=Edges[i].to;
        if(v==maxv[u]||v==father[u])continue;
        dfs2(v,v);
    }
}
void update(int idx)
{
    sum[idx]=sum[ls[idx]]+sum[rs[idx]];
    maxn[idx]=max(maxn[ls[idx]],maxn[rs[idx]]);
}
void change(int &idx,int l,int r,int p,int w)
{
    if(!idx)++tot,idx=tot;
    if(l==r){maxn[idx]=sum[idx]=w;return;}
    int mid=(l+r)>>1;
    if(p<=mid)change(ls[idx],l,mid,p,w);
    else change(rs[idx],mid+1,r,p,w);
    update(idx);
}
int querysum(int idx,int l,int r,int L,int R)
{
    if(!idx)return 0;
    if(L==l&&R==r)return sum[idx];
    int mid=(l+r)>>1;
    if(R<=mid)return querysum(ls[idx],l,mid,L,R);
    else if(L>mid)return querysum(rs[idx],mid+1,r,L,R);
    else return querysum(ls[idx],l,mid,L,mid)+querysum(rs[idx],mid+1,r,mid+1,R);
}
int querymax(int idx,int l,int r,int L,int R)
{
    if(!idx)return 0;
    if(L==l&&R==r)return maxn[idx];
    int mid=(l+r)>>1;
    if(R<=mid)return querymax(ls[idx],l,mid,L,R);
    else if(L>mid)return querymax(rs[idx],mid+1,r,L,R);
    else return max(querymax(ls[idx],l,mid,L,mid),querymax(rs[idx],mid+1,r,mid+1,R));
}
int QS(int x,int y)
{
    int t=c[x],res=0;
    while(top[x]!=top[y])
    {
        if(deep[top[x]]<deep[top[y]])swap(x,y);
        res+=querysum(rt[t],1,sz,pos[top[x]],pos[x]);
        x=father[top[x]];
    }
    if(deep[x]<deep[y])swap(x,y);
    res+=querysum(rt[t],1,sz,pos[y],pos[x]);
    return res;
}

int QM(int x,int y)
{
    int t=c[x],res=0;
    while(top[x]!=top[y])
    {
        if(deep[top[x]]<deep[top[y]])swap(x,y);
        res=max(res,querymax(rt[t],1,sz,pos[top[x]],pos[x]));
        x=father[top[x]];
    }
    if(deep[x]<deep[y])swap(x,y);
    res=max(res,querymax(rt[t],1,sz,pos[y],pos[x]));
    return res;
}
int main()
{
    memset(head,-1,sizeof(head));
    n=read(),q=read();
    for(int i=1;i<=n;i++)w[i]=read(),c[i]=read();
    for(int i=1;i<n;i++)
    {
        int x=read(),y=read();
        addedge(x,y),addedge(y,x);
    }
    deep[1]=1;dfs1(1),dfs2(1,1);
    for(int i=1;i<=n;i++)change(rt[c[i]],1,sz,pos[i],w[i]);
    for(int i=1;i<=q;i++)
    {
        char opt[5];
        scanf("%s",opt);
        int x=read(),y=read();
        if(opt[1]=='C')
        {
            change(rt[c[x]],1,sz,pos[x],0);
            change(rt[y],1,sz,pos[x],w[x]);
            c[x]=y;
        }
        else if(opt[1]=='W')
        {
            change(rt[c[x]],1,sz,pos[x],y);
            w[x]=y;
        }
        else if(opt[1]=='S')
        printf("%d\n",QS(x,y));
        else printf("%d\n",QM(x,y));
    }
    return 0;
}

 

转载于:https://www.cnblogs.com/Zars19/p/6937529.html

相关文章:

  • TeamCity 2017.1.2 发布,持续集成工具
  • 微软:请叫我 Android 预装服务提供商
  • arcgis engine 获取高亮Feature、element
  • Java -- JDBC 学习--PreparedStatement
  • Vijos P1303 导弹拦截【最长上升子序列+DP】
  • Ta-Lib用法介绍 !
  • 负载均衡方案对比表
  • select标签中option内容加链接
  • kibana创建新的index patterns
  • 使用docker搭建android源码编译环境
  • 使用Android Studio运行NDK程序
  • Foreman 1.15.1 发布,数据中心生命周期管理工具
  • C学习-枚举(九)
  • 关于网易云音乐爬虫的api接口?
  • 设计模式(一)----工厂模式
  • [iOS]Core Data浅析一 -- 启用Core Data
  • Android开发 - 掌握ConstraintLayout(四)创建基本约束
  • windows下如何用phpstorm同步测试服务器
  • 程序员该如何有效的找工作?
  • 对超线程几个不同角度的解释
  • 翻译--Thinking in React
  • 聊聊spring cloud的LoadBalancerAutoConfiguration
  • 模仿 Go Sort 排序接口实现的自定义排序
  • 排序(1):冒泡排序
  • 融云开发漫谈:你是否了解Go语言并发编程的第一要义?
  • 试着探索高并发下的系统架构面貌
  • 想写好前端,先练好内功
  • 学习笔记TF060:图像语音结合,看图说话
  • 用Node EJS写一个爬虫脚本每天定时给心爱的她发一封暖心邮件
  • 用Visual Studio开发以太坊智能合约
  • ​香农与信息论三大定律
  • # .NET Framework中使用命名管道进行进程间通信
  • (¥1011)-(一千零一拾一元整)输出
  • (07)Hive——窗口函数详解
  • (11)MATLAB PCA+SVM 人脸识别
  • (6)设计一个TimeMap
  • (70min)字节暑假实习二面(已挂)
  • (Git) gitignore基础使用
  • (附源码)springboot金融新闻信息服务系统 毕业设计651450
  • (四)Tiki-taka算法(TTA)求解无人机三维路径规划研究(MATLAB)
  • (四)七种元启发算法(DBO、LO、SWO、COA、LSO、KOA、GRO)求解无人机路径规划MATLAB
  • (五)Python 垃圾回收机制
  • (循环依赖问题)学习spring的第九天
  • (转)Mysql的优化设置
  • ******之网络***——物理***
  • .NET 4.0中的泛型协变和反变
  • .net core 6 集成和使用 mongodb
  • .net 重复调用webservice_Java RMI 远程调用详解,优劣势说明
  • @Bean, @Component, @Configuration简析
  • @ConditionalOnProperty注解使用说明
  • @WebService和@WebMethod注解的用法
  • @四年级家长,这条香港优才计划+华侨生联考捷径,一定要看!
  • [ C++ ] STL priority_queue(优先级队列)使用及其底层模拟实现,容器适配器,deque(双端队列)原理了解
  • [2019.3.20]BZOJ4573 [Zjoi2016]大森林
  • [android] 手机卫士黑名单功能(ListView优化)