【GDOI2014模拟】雨天的尾巴 题解+代码

2024-05-29 03:32

本文主要是介绍【GDOI2014模拟】雨天的尾巴 题解+代码,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

Description

深绘里一直很讨厌雨天。
灼热的天气穿透了前半个夏天,后来一场大雨和随之而来的洪水,浇灭了一切。
虽然深绘里家乡的小村落对洪水有着顽固的抵抗力,但也倒了几座老房子,几棵老树被连
根拔起,以及田地里的粮食被弄得一片狼藉。
无奈的深绘里和村民们只好等待救济粮来维生。
不过救济粮的发放方式很特别。
首先村落里的一共有n 座房屋,并形成一个树状结构。然后救济粮分m 次发放,每次选
两个房屋(x,y),然后对于x 到y 的路径上(含x 和y) 每座房子里发放一袋z 类型的救济粮。
然后深绘里想知道,当所有的救济粮发放完毕后,每座房子里存放的最多的是哪种救济粮。

Input

第一行两个正整数n,m,含义如题目所示。
接下来n-1 行,每行两个数(a,b),表示(a,b) 间有一条边。
再接下来m 行,每行三个数(x,y,z),含义如题目所示。

Output

n行,第i 行一个整数,表示第i 座房屋里存放的最多的是哪种救济粮,如果有多种救济粮
存放次数一样,输出编号最小的。
如果某座房屋里没有救济粮,则对应一行输出0。

Sample Input

5 3
1 2
3 1
3 4
5 4
2 3 3
1 5 2
3 3 3

Sample Output

2
3
3
2
2

Data Constraint

对于20% 的数据,1<= n,m <= 100
对于50% 的数据,1 <= n,m <= 2000
对于100% 的数据,1 <= n;m <= 100000; 1 <= a, b, x, y <= n; 1 <= z <= 10^9

Solutions

一颗树,区间改变,可以考虑树链剖分。
但是,这题不用树链剖分。
对于每个节点,开一颗动态开节点的权值线段树,存的是每种救济粮的数目。
但是显然不能每次操作都把x到y之间的所有点维护一次线段树,那么可以用到一个很神奇的东西叫做:线段树合并。
如果一个点x,设它的父亲为fa[x],那么把x与fa[x]合并,就可以把x的状态加到fa[x]里面即x的每个救济粮的数目可以复制到fa[x]中,那么每次操作的时候就只用把点x和y对应位置+1,lca(x,y)和fa[lca(x,y)]对应位置-1,就行了。
那么怎么线段树合并呢?
简单来说:对于两个线段树,如果一个为空,就返回另一个,否则就两个线段树分别左儿子合并,右儿子合并就行了,在叶子节点操作。在这道题中,在叶子节点的操作就是加起来。
线段树合并代码

void merge(int v,int v1,int i,int j)
{if (i==j) {tree[v].data+=tree[v1].data;return;}int mid=(i+j)/2;if (tree[v1].l!=0) if(tree[v].l==0) tree[v].l=tree[v1].l;else merge(tree[v].l,tree[v1].l,i,mid);if (tree[v1].r!=0)if(tree[v].r==0) tree[v].r=tree[v1].r;else merge(tree[v].r,tree[v1].r,mid+1,j);tree[v].data=max(tree[tree[v].l].data,tree[tree[v].r].data);
}

code

#include<cstdio>
#include<cstring>
#include<cmath>
#include<algorithm>
#define MAX 202000
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define N 101000
using namespace std;
int deep[N],fa[N],last[N*10],next[N*10],to[N*10],n,hash[N*10],tot=0,a[N],ans[N],g[N],f[N][17];
bool bz[N];
struct abcde{int l,r,data;
};
abcde tree[N*100];
void putin(int x,int y)
{next[++tot]=last[x];last[x]=tot;to[tot]=y;
}
int lsh(int x)
{int y=x%MAX;while (hash[y]!=0 && hash[y]!=x) y=(y+1)%MAX;hash[y]=x;return y;}
void dg1(int x)
{for(int i=last[x];i;i=next[i]){int k=to[i];if (deep[k]!=0) continue;fa[k]=x;deep[k]=deep[x]+1;dg1(k);}
}
void insert(int v,int i,int j,int x,int y)
{if (i==j){tree[v].data+=y;return;}int mid=(i+j)/2;if (x<=mid){if (tree[v].l==0) tree[v].l=++tot;insert(tree[v].l,i,mid,x,y);}else{if (tree[v].r==0) tree[v].r=++tot;insert(tree[v].r,mid+1,j,x,y);}tree[v].data=max(tree[tree[v].l].data,tree[tree[v].r].data);
}
void merge(int v,int v1,int i,int j)
{if (i==j) {tree[v].data+=tree[v1].data;return;}int mid=(i+j)/2;if (tree[v1].l!=0) if(tree[v].l==0) tree[v].l=tree[v1].l;else merge(tree[v].l,tree[v1].l,i,mid);if (tree[v1].r!=0)if(tree[v].r==0) tree[v].r=tree[v1].r;else merge(tree[v].r,tree[v1].r,mid+1,j);tree[v].data=max(tree[tree[v].l].data,tree[tree[v].r].data);
}
int lca(int x,int y)
{fd(i,16,0) if(deep[f[x][i]]>=deep[y]) x=f[x][i];fd(i,16,0) if(deep[f[y][i]]>=deep[x]) y=f[y][i];fd(i,16,0) if(f[x][i]!=f[y][i]) x=f[x][i],y=f[y][i];if (x!=y) x=f[x][0],y=f[y][0];return x;
}
int find(int v,int i,int j)
{if (i==j) return i;int mid=(i+j)/2;if (tree[tree[v].l].data>=tree[tree[v].r].data) return find(tree[v].l,i,mid);else return find(tree[v].r,mid+1,j);
}
void dfs(int x,int fa)
{for(int i=last[x];i;i=next[i]){int k=to[i];if (k==fa) continue;dfs(k,x);merge(x,k,1,MAX);}ans[x]=find(x,1,MAX);
}
int main()
{int m;scanf("%d%d",&n,&m); fo(i,1,n-1){int x,y;scanf("%d%d",&x,&y);putin(x,y);putin(y,x);}deep[1]=1;dg1(1);tot=n;fo(i,1,n) f[i][0]=fa[i];fo(j,1,16)fo(i,1,n) f[i][j]=f[f[i][j-1]][j-1];fo(i,1,m) {int x,y,z;scanf("%d%d%d",&x,&y,&z);z=lsh(z);int lc=lca(x,y);insert(x,1,MAX,z,1);insert(y,1,MAX,z,1);insert(lc,1,MAX,z,-1);insert(fa[lc],1,MAX,z,-1);}dfs(1,0);fo(i,1,n) printf("%d\n",hash[ans[i]]);
}

这篇关于【GDOI2014模拟】雨天的尾巴 题解+代码的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1012523

相关文章

利用Python调试串口的示例代码

《利用Python调试串口的示例代码》在嵌入式开发、物联网设备调试过程中,串口通信是最基础的调试手段本文将带你用Python+ttkbootstrap打造一款高颜值、多功能的串口调试助手,需要的可以了... 目录概述:为什么需要专业的串口调试工具项目架构设计1.1 技术栈选型1.2 关键类说明1.3 线程模

Python Transformers库(NLP处理库)案例代码讲解

《PythonTransformers库(NLP处理库)案例代码讲解》本文介绍transformers库的全面讲解,包含基础知识、高级用法、案例代码及学习路径,内容经过组织,适合不同阶段的学习者,对... 目录一、基础知识1. Transformers 库简介2. 安装与环境配置3. 快速上手示例二、核心模

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

使用Java将DOCX文档解析为Markdown文档的代码实现

《使用Java将DOCX文档解析为Markdown文档的代码实现》在现代文档处理中,Markdown(MD)因其简洁的语法和良好的可读性,逐渐成为开发者、技术写作者和内容创作者的首选格式,然而,许多文... 目录引言1. 工具和库介绍2. 安装依赖库3. 使用Apache POI解析DOCX文档4. 将解析

C++使用printf语句实现进制转换的示例代码

《C++使用printf语句实现进制转换的示例代码》在C语言中,printf函数可以直接实现部分进制转换功能,通过格式说明符(formatspecifier)快速输出不同进制的数值,下面给大家分享C+... 目录一、printf 原生支持的进制转换1. 十进制、八进制、十六进制转换2. 显示进制前缀3. 指

使用Python实现全能手机虚拟键盘的示例代码

《使用Python实现全能手机虚拟键盘的示例代码》在数字化办公时代,你是否遇到过这样的场景:会议室投影电脑突然键盘失灵、躺在沙发上想远程控制书房电脑、或者需要给长辈远程协助操作?今天我要分享的Pyth... 目录一、项目概述:不止于键盘的远程控制方案1.1 创新价值1.2 技术栈全景二、需求实现步骤一、需求

Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码

《Java中Date、LocalDate、LocalDateTime、LocalTime、时间戳之间的相互转换代码》:本文主要介绍Java中日期时间转换的多种方法,包括将Date转换为LocalD... 目录一、Date转LocalDateTime二、Date转LocalDate三、LocalDateTim

jupyter代码块没有运行图标的解决方案

《jupyter代码块没有运行图标的解决方案》:本文主要介绍jupyter代码块没有运行图标的解决方案,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录jupyter代码块没有运行图标的解决1.找到Jupyter notebook的系统配置文件2.这时候一般会搜索到

Python通过模块化开发优化代码的技巧分享

《Python通过模块化开发优化代码的技巧分享》模块化开发就是把代码拆成一个个“零件”,该封装封装,该拆分拆分,下面小编就来和大家简单聊聊python如何用模块化开发进行代码优化吧... 目录什么是模块化开发如何拆分代码改进版:拆分成模块让模块更强大:使用 __init__.py你一定会遇到的问题模www.

springboot循环依赖问题案例代码及解决办法

《springboot循环依赖问题案例代码及解决办法》在SpringBoot中,如果两个或多个Bean之间存在循环依赖(即BeanA依赖BeanB,而BeanB又依赖BeanA),会导致Spring的... 目录1. 什么是循环依赖?2. 循环依赖的场景案例3. 解决循环依赖的常见方法方法 1:使用 @La