Luogu P3313 [SDOI2014]旅行

最新推荐文章于 2019-08-19 10:53:00 发布
最新推荐文章于 2019-08-19 10:53:00 发布
阅读量333 收藏 1
点赞数 1
分类专栏： Luogu 树链剖分 线段树
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42142540/article/details/81027466
版权
Luogu 同时被 3 个专栏收录
22 篇文章 0 订阅
订阅专栏
线段树
3 篇文章 0 订阅
订阅专栏
树链剖分
1 篇文章 0 订阅
订阅专栏
Luogu P3313 [SDOI2014]旅行

血的教训q-q
错误示范
因为这个我WA了MLE了不知道多少次q-q
我终于自我察觉了，然后删了
然后就AC了
。
我和LHY说了我非常沙雕的错误
LHY非常好人教了我，如果需要赋值时应该这样
正确的打开方式

本题解法：树链剖分+动态开点线段树
对于这么多个宗教，每个宗教开一棵线段树
树链剖分寻找路径
再去线段树上查找，也就不用考虑途中城市是不是同一宗教
*数据保证询问的st和ed同一宗教

代码：

#include<cstdio>
#include<cstring>

int n,m,len=0;
struct nod1{int x,y,next;}b[200010];
struct nod2{int lc,rc,mx,tot;}tr[20000010];
int size[100010],son[100010],newid[100010],first[100010],f[100010];
int deep[100010],top[100010],w[100010],c[100010],root[100010];

int maxx(int x,int y)
{
    return x>y?x:y;
}

int ins(int x,int y)
{
    len++;
    b[len].x=x;
    b[len].y=y;
    b[len].next=first[x];
    first[x]=len;
}

int dfs1(int x)
{
    size[x]=1;
    for(int i=first[x];i;i=b[i].next)
    {
        int y=b[i].y;
        if(y!=f[x])
        {
            deep[y]=deep[x]+1;
            f[y]=x;
            dfs1(y);
            if(size[son[x]]<size[y])son[x]=y;
            size[x]+=size[y];
        }
    }
}

int dfs2(int x,int ttop)
{
    top[x]=ttop;
    len++;newid[x]=len;
    ///id[len]=x;
    if(son[x]!=0)dfs2(son[x],ttop);
    for(int i=first[x];i;i=b[i].next)
    {
        int y=b[i].y;
        if(y!=son[x]&&y!=f[x])
        {
            dfs2(y,y);
        }
    }
}

void update(int &now,int l,int r,int x,int k)
{
    //printf("--%d %d %d %d %d--\n",now,l,r,x,k);
    if(now==0)
    {
        len++;now=len;
    }
    tr[now].tot+=k;
    if(l==r)
    {
        if(k>0)tr[now].mx=k;
        else tr[now].mx=0;
        return ;
    }
    int mid=(l+r)/2;
    if(x<=mid)update(tr[now].lc,l,mid,x,k);
    else update(tr[now].rc,mid+1,r,x,k);

    tr[now].mx=maxx(tr[tr[now].lc].mx,tr[tr[now].rc].mx);
    //printf("|%d %d|\n",tr[now].tot,tr[now].mx);
}

int change1(int x,int y)
{
    update(root[c[x]],1,n,newid[x],-w[x]);
    c[x]=y;
    update(root[c[x]],1,n,newid[x],w[x]);
}

int change2(int x,int y)
{
    update(root[c[x]],1,n,newid[x],-w[x]);
    w[x]=y;
    update(root[c[x]],1,n,newid[x],w[x]);
}

int swap(int &x,int &y)
{
    int t=x;x=y;y=t;
}

int findsum(int now,int l,int r,int x,int y)
{
    if(l==x&&y==r)return tr[now].tot;
    int mid=(l+r)/2,lc=tr[now].lc,rc=tr[now].rc;
    if(y<=mid)return findsum(lc,l,mid,x,y);
    else if(mid+1<=x)return findsum(rc,mid+1,r,x,y);
    else return findsum(lc,l,mid,x,mid)+findsum(rc,mid+1,r,mid+1,y);
}

int sumsum(int x,int y)
{
    int sum=0;
    int st=x;
    int topx=top[x],topy=top[y];
    while(topx!=topy)
    {
        if(deep[topx]>deep[topy])
        {
            swap(x,y);swap(topx,topy);
        }
        sum+=findsum(root[c[st]],1,n,newid[topy],newid[y]);
        y=f[topy];topy=top[y];
    }
    if(deep[x]>deep[y])swap(x,y);
    sum+=findsum(root[c[st]],1,n,newid[x],newid[y]);
    printf("%d\n",sum);
}

int findmax(int now,int l,int r,int x,int y)
{
    if(l==x&&y==r)return tr[now].mx;
    int mid=(l+r)/2,lc=tr[now].lc,rc=tr[now].rc;
    if(y<=mid)return findmax(lc,l,mid,x,y);
    if(mid+1<=x)return findmax(rc,mid+1,r,x,y);
    else return maxx(findmax(lc,l,mid,x,mid),findmax(rc,mid+1,r,mid+1,y));
}

int mxmx(int x,int y)
{
    int ma=0;
    int st=x;
    int topx=top[x],topy=top[y];
    while(topx!=topy)
    {
        if(deep[topx]>deep[topy])
        {
            swap(x,y);swap(topx,topy);
        }
        int g=findmax(root[c[st]],1,n,newid[topy],newid[y]);
        if(g>ma)ma=g;
        y=f[topy];topy=top[y];
    }
    if(deep[x]>deep[y])swap(x,y);
    int g=findmax(root[c[st]],1,n,newid[x],newid[y]);
    if(g>ma)ma=g;
    printf("%d\n",ma);
}

int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d %d",&w[i],&c[i]);
    for(int i=1;i<=n-1;i++)
    {
        int x,y;
        scanf("%d %d",&x,&y);
        ins(x,y);ins(y,x);
    }
    len=0;
    deep[1]=1;son[1]=0;
    size[1]=0;f[1]=0;
    dfs1(1);len=0;
    dfs2(1,1);len=0;

    /*for(int i=1;i<=n;i++)
    {
        printf("|%d |\n",root[c[i]]);
    }*/

    for(int i=1;i<=n;i++)
        update(root[c[i]],1,n,newid[i],w[i]);
        //ok

    for(int i=1;i<=m;i++)
    {
        int x,y;
        char ss[10];
        scanf("%s",ss+1);
        scanf("%d %d",&x,&y);
        if(ss[2]=='C')change1(x,y);
        else if(ss[2]=='W')change2(x,y);
        else if(ss[2]=='S')sumsum(x,y);
        else if(ss[2]=='M')mxmx(x,y);
    }
}

代码+注释：

#include<cstdio>
#include<cstring>

int n,m,len=0;
struct nod1{int x,y,next;}b[200010];
struct nod2{int lc,rc,mx,tot;}tr[20000010];
int size[100010],son[100010],newid[100010],first[100010],f[100010];
int deep[100010],top[100010],w[100010],c[100010],root[100010];

int maxx(int x,int y)
{
    return x>y?x:y;
}

int ins(int x,int y)
{
    len++;
    b[len].x=x;
    b[len].y=y;
    b[len].next=first[x];
    first[x]=len;
}

int dfs1(int x)
{
    size[x]=1;//以x为根的子树至少一个结点
    for(int i=first[x];i;i=b[i].next)
    {
        int y=b[i].y;
        if(y!=f[x])
        {//不是父亲才能是自己儿子 
            deep[y]=deep[x]+1;
            f[y]=x;
            dfs1(y);
            if(size[son[x]]<size[y])son[x]=y;
            //如果你比我当前重儿子还厉害，那我让你当重儿子 
            size[x]+=size[y];
            //你的后代也是我的后代 
        }
    }
}

int dfs2(int x,int ttop)
//x和当前链的顶端 
{
    top[x]=ttop;//更新我的top 
    len++;newid[x]=len;//x在树链剖分维护的线段树上的新编号 
    if(son[x]!=0)dfs2(son[x],ttop);
    //如果我有重儿子，他和我一条重链，他的top和我一样，先更新他 
    for(int i=first[x];i;i=b[i].next)
    {
        int y=b[i].y;
        if(y!=son[x]&&y!=f[x])
        {
            dfs2(y,y);//轻儿子们自己是重链 
        }
    }
}

void update(int &now,int l,int r,int x,int k)
{
    if(now==0)
    {
        len++;now=len;//动态开点 
    }
    tr[now].tot+=k;
    if(l==r)
    {
        if(k>0)tr[now].mx=k;
        else tr[now].mx=0;
        return ;
    }
    int mid=(l+r)/2;
    if(x<=mid)update(tr[now].lc,l,mid,x,k);
    else update(tr[now].rc,mid+1,r,x,k);

    tr[now].mx=maxx(tr[tr[now].lc].mx,tr[tr[now].rc].mx);
}

int change1(int x,int y)
{
    update(root[c[x]],1,n,newid[x],-w[x]);
    c[x]=y;
    update(root[c[x]],1,n,newid[x],w[x]);
}

int change2(int x,int y)
{
    update(root[c[x]],1,n,newid[x],-w[x]);
    w[x]=y;
    update(root[c[x]],1,n,newid[x],w[x]);
}

int swap(int &x,int &y)
{
    int t=x;x=y;y=t;
}

int findsum(int now,int l,int r,int x,int y)
{//线段树查找 
    if(l==x&&y==r)return tr[now].tot;
    int mid=(l+r)/2,lc=tr[now].lc,rc=tr[now].rc;
    if(y<=mid)return findsum(lc,l,mid,x,y);
    else if(mid+1<=x)return findsum(rc,mid+1,r,x,y);
    else return findsum(lc,l,mid,x,mid)+findsum(rc,mid+1,r,mid+1,y);
}

int sumsum(int x,int y)
{
    int sum=0;
    int st=x;
    int topx=top[x],topy=top[y];
    while(topx!=topy)//还没跳到一条链上 
    {
        if(deep[topx]>deep[topy])
        {//我们跳y所以要保证y的深度是比x小的 
            swap(x,y);swap(topx,topy);
        }
        sum+=findsum(root[c[st]],1,n,newid[topy],newid[y]);
        //累加跳掉（离开）的这条链的值 
        y=f[topy];topy=top[y];
        //把y跳上另一条链 
    }
    if(deep[x]>deep[y])swap(x,y);
    sum+=findsum(root[c[st]],1,n,newid[x],newid[y]);
    //因为他们在一条链上之后就退出来了，还没有累加这条链 
    printf("%d\n",sum);
}

int findmax(int now,int l,int r,int x,int y)
{//线段树查找 
    if(l==x&&y==r)return tr[now].mx;
    int mid=(l+r)/2,lc=tr[now].lc,rc=tr[now].rc;
    if(y<=mid)return findmax(lc,l,mid,x,y);
    if(mid+1<=x)return findmax(rc,mid+1,r,x,y);
    else return maxx(findmax(lc,l,mid,x,mid),findmax(rc,mid+1,r,mid+1,y));
}

int mxmx(int x,int y)
{//注释见sumsum 差不多la——— 
    int ma=0;
    int st=x;
    int topx=top[x],topy=top[y];
    while(topx!=topy)
    {
        if(deep[topx]>deep[topy])
        {
            swap(x,y);swap(topx,topy);
        }
        int g=findmax(root[c[st]],1,n,newid[topy],newid[y]);
        if(g>ma)ma=g;
        y=f[topy];topy=top[y];
    }
    if(deep[x]>deep[y])swap(x,y);
    int g=findmax(root[c[st]],1,n,newid[x],newid[y]);
    if(g>ma)ma=g;
    printf("%d\n",ma);
}

int main()
{
    scanf("%d %d",&n,&m);
    for(int i=1;i<=n;i++)
        scanf("%d %d",&w[i],&c[i]);
    for(int i=1;i<=n-1;i++)
    {
        int x,y;
        scanf("%d %d",&x,&y);
        ins(x,y);ins(y,x);
    }
    len=0;
    deep[1]=1;son[1]=0;
    size[1]=0;f[1]=0;
    dfs1(1);len=0;
    dfs2(1,1);len=0;

    for(int i=1;i<=n;i++)
        update(root[c[i]],1,n,newid[i],w[i]);
        //先把初始信息添加 

    for(int i=1;i<=m;i++)
    {
        int x,y;
        char ss[10];
        scanf("%s",ss+1);
        scanf("%d %d",&x,&y);
        if(ss[2]=='C')change1(x,y);
        else if(ss[2]=='W')change2(x,y);
        else if(ss[2]=='S')sumsum(x,y);
        else if(ss[2]=='M')mxmx(x,y);
    }
}

树链剖分好像也不是很坏（还是比较友善吧）qvq

ilcely
关注
专栏目录
[SDOI2014]旅行
C202044zxy的博客
08-04 182
一、题目及数据范围 点此看题 二、解法 这题一看就是树链剖分（可是菜鸡作者不会）。 点此自学树链剖分 现在好了，默认大家都会树链剖分。 我们对于每一种宗教，维护一个线段树，存tottottot和maxmaxmax。 但是这题有1e51e51e5，这样开不是开玩笑？ 所以我们要动态开点。 由于我们要维护的是每一条树链的信息，再经过合并输出答案，所以我们不需要去处理树链中没有出现过的宗教，且对于每一宗...
「SDOI2014」旅行
BFZD的博客
03-23 361
「SDOI2014」旅行 谨以此题纪念第一次动态开点的题 大意 给定一棵树和每一个节点的颜色，和权值。 要求能够修改，并且可以查询一条路径上同种颜色的节点的信息。 思路 看到查询一条路径可以想到用树链剖分和线段树维护，但是要求要查询同一种颜色的节点信息肿么办？ 按颜色来建线段树？ 空间会非常爆炸，那么就动态开点。 修改颜色，将原来颜色线段树上的节点的权值赋成0 并且将当前颜色的线段...
AC日记——旅行 洛谷 P3313
weixin_33743703的博客
02-17 167
题目描述 S国有N个城市，编号从1到N。城市间用N-1条双向道路连接，满足从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教，如飞天面条神教、隐形独角兽教、绝地教都是常见的信仰。 为了方便，我们用不同的正整数代表各种宗教， S国的居民常常旅行。旅行时他们总会走最短路，并且为了避免麻烦，只在信仰和他们相同的城市留宿。当然旅程的终点也是信仰与他相同的城市。S国政府为每个城市标定了不同的旅行...
luogu_P3313 [SDOI2014]旅行
weixin_33979363的博客
12-17 112
传送门 Solution 第二次学习可持久化线段树 打了一道裸题来练习一下…… 对于每个宗教都可以开一个主席树 基础操作 树剖lca Code #include<bits/stdc++.h> #define ll long long #define max(a,b) ((a)>(b)?(a):(b)) #define min(a,b) ((a)<(b)?(...
洛谷 P3313 [SDOI2014]旅行 解题报告
weixin_30832143的博客
06-17 189
P3313 [SDOI2014]旅行 题目描述 S国有N个城市，编号从1到N。城市间用N-1条双向道路连接，满足从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教，如飞天面条神教、隐形独角兽教、绝地教都是常见的信仰。 为了方便，我们用不同的正整数代表各种宗教， S国的居民常常旅行。旅行时他们总会走最短路，并且为了避免麻烦，只在信仰和他们相同的城市留宿。当然旅程的终点也是信仰与他相同的城市。...
【SDOI2014】旅行
Facico的博客
04-02 1286
DescriptionSolution用什么很明显，这是一个树上操作，但是需要类似线段树的维护，一个很裸的树链剖分加线段树。怎么做因为有c个宗教，那么就开c棵线段树。然后就可以了。空间太大了我用了动态开空间，类似主席树那样打。在线大法好！ 也可以打一个离散化。离线大法好！Code#include<iostream> #include<cstring> #include<cstdio> #inclu
锣鼓声声迎新春——新年工作计划ppt模板.zip
09-08
标题中的“锣鼓声声迎新春——新年工作计划ppt模板.zip”表明这是一个为新年工作计划设计的PPT模板，主题风格融入了中国传统的新年元素，如锣鼓和红色的喜庆氛围，暗示了激励和活力的开端。描述进一步强调了模板的...
锣鼓队表演解说词精选.doc
10-11
锣鼓队表演解说词精选 锣鼓队表演讲解词是一种特殊的艺术形式，它通过锣鼓的敲击和音乐的演奏，展现出激情四射的歌、豪迈粗狂的舞、刚毅奔放的魂等多种气概和情感。这种艺术形式可以激发人们的情感，激励人们...
2016湘艺版音乐四年级下册第10课《神奇的民间锣鼓“蛇蜕壳”》课件
11-23
### 2016湘艺版音乐四年级下册第10课《神奇的民间锣鼓“蛇蜕壳”》课件知识点解析 #### 一、课程简介与目标 本节课主要介绍了湖南地区的一种传统民间音乐——“蛇蜕壳”。通过本课的学习，学生能够初步了解这种...
锣鼓队表演解说词.doc
09-18
锣鼓队表演解说词.doc
丰收锣鼓 欣赏PPT学习教案.pptx
10-02
丰收锣鼓 欣赏PPT学习教案.pptx
【bzoj3531】[Sdoi2014]旅行
Joky_2002的博客
03-06 282
3531: [Sdoi2014]旅行Time Limit: 40 Sec  Memory Limit: 512 MBSubmit: 2833  Solved: 1227[Submit][Status][Discuss]Description S国有N个城市，编号从1到N。城市间用N-1条双向道路连接，满足从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教，如飞天面条神教、隐形独角兽教、绝地...
BZOJ3531: [Sdoi2014]旅行
weixin_30340617的博客
06-29 121
3531: [Sdoi2014]旅行 Time Limit: 20 SecMemory Limit: 512 MBSubmit: 323Solved: 192[Submit][Status] Description S国有N个城市，编号从1到N。城市间用N-1条双向道路连接，满足 从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教，如飞天面条神教、隐形独角兽教、绝...
SDOI2014旅行
Ninaye的博客
08-02 238
题意：有一棵树，每个节点有两个值（一个是种类，一个是价值） 有四个操作 1、修改某点的种类。 2、修改某点的价值。 3、给出两个点（种类相同），求这两点之间路径上和他们种类相同的点的最大价值 4、给出两个点（种类相同），求这两点之间路径上和他们种类 这个题敲了整整一天啊，一开始不太明白动态开点什么意思，最后终于明白了，然后无情的被T了一片一片的，终于A了，好开心！！！
[bzoj 3531]sdoi2014 旅行
Zeyu_King的专栏
06-22 669
一眼树剖题 对每种颜色开一颗线段树，动态开点，记录最大值与和，查询直接查找，修改的话只会修改单点，所以最多增加2logn个节点 （原来把题目看错了，以为修改也是修改连续一段的，似乎这样空间会大很多）
「SDOI2014」旅行（信息学奥赛一本通 1564）（洛谷 3313）
anglanjing7414的博客
08-19 623
题目描述 S国有N个城市，编号从1到N。城市间用N-1条双向道路连接，满足从一个城市出发可以到达其它所有城市。每个城市信仰不同的宗教，如飞天面条神教、隐形独角兽教、绝地教都是常见的信仰。 为了方便，我们用不同的正整数代表各种宗教， S国的居民常常旅行。旅行时他们总会走最短路，并且为了避免麻烦，只在信仰和他们相同的城市留宿。当然旅程的终点也是信仰与他相同的城市。S国政府为每个城市标定了...
南锣鼓巷高德poi下载
最新发布
11-11
南锣鼓巷是北京著名的胡同之一，吸引着大量的游客和游览者。如果你想在南锣鼓巷游览时不迷路，可以下载高德...总之，通过下载高德地图，你可以更方便地探索南锣鼓巷，了解周边环境，规划行程，让你的旅行更加顺利愉快。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值