[CC-SUBWAY]Subway Ride
[CC-SUBWAY]Subway Ride
题目大意:
一棵\(n(n\le5\times10^5)\)个点的含重边的树,总边数为\(m(m\le10^6)\),每条边有一个颜色。\(q(q\le5\times10^5)\)次询问,每次询问两个结点之间的路径(不经过重复结点)最多包含多少连续的同色段。
思路:
一个显然的结论是,对于连接相同两个点的重边,如果不同的颜色数量\(\le 3\),那么至少有一条边和路径两侧的颜色不相同,这时我们只需要发明一种新的颜色来代替这些颜色即可。因此每条边最多只需要保留\(2\)种颜色。
用\(f[i][0/1]\)表示到了第\(i\)个结点,上一条边是第\(0/1\)种颜色的答案。每次询问跑一遍这样的DP,单次的时间复杂度是\(\mathcal O(n)\)。
使用树上倍增的思想,\(f[i][j][k][l]\)表示\(i\sim anc[i][j]\)的路径上,第一条边是第\(k\)种颜色,最后一条边是第\(l\)种颜色,询问时将一整段转移即可。
时间复杂度\(\mathcal O(q\log n)\)。
源代码:
#include<map>
#include<cstdio>
#include<cctype>
#include<vector>
inline int getint() {
register char ch;
while(!isdigit(ch=getchar()));
register int x=ch^'0';
while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
return x;
}
const int N=5e5+1,logN=19;
int tot,anc[N][logN],dep[N],f[N][2],g[N][2];
int seg[N][logN][2][2],w[N][logN][2];
std::vector<int> e[N];
inline void add_edge(const int &u,const int &v) {
e[u].push_back(v);
e[v].push_back(u);
}
std::map<std::pair<int,int>,std::pair<int,int>> map;
inline int lg2(const float &x) {
return ((unsigned&)x>>23&255)-127;
}
inline void up(int &x,const int &y) {
x=std::max(x,y);
}
void dfs(const int &x,const int &par) {
if(par) {
int u=x,v=par;
if(u>v) std::swap(u,v);
const auto &p=map[std::make_pair(u,v)];
if(p.first) {
w[x][0][0]=p.first;
seg[x][0][0][0]=1;
}
if(p.second) {
w[x][0][1]=p.second;
seg[x][0][1][1]=1;
}
}
anc[x][0]=par;
dep[x]=dep[par]+1;
for(register int i=1;i<=lg2(dep[x]);i++) {
anc[x][i]=anc[anc[x][i-1]][i-1];
w[x][i][0]=w[anc[x][i-1]][i-1][0];
w[x][i][1]=w[anc[x][i-1]][i-1][1];
for(register int j=0;j<2;j++) {
if(!w[x][0][j]) continue;
for(register int k=0;k<2;k++) {
if(!w[x][i][k]) continue;
for(register int l=0;l<2;l++) {
if(!w[x][i-1][l]) continue;
for(register int m=0;m<2;m++) {
if(!w[anc[x][i-1]][0][m]) continue;
up(seg[x][i][j][k],seg[x][i-1][j][l]+seg[anc[x][i-1]][i-1][m][k]-(w[x][i-1][l]==w[anc[x][i-1]][0][m]));
}
}
}
}
}
for(auto &y:e[x]) {
if(y!=par) dfs(y,x);
}
}
inline int query(int u,int v) {
if(u==v) return 0;
if(dep[u]<dep[v]) std::swap(u,v);
f[u][0]=f[u][1]=g[u][0]=g[u][1]=0;
f[v][0]=f[v][1]=g[v][0]=g[v][1]=0;
for(register int i=lg2(dep[u]-dep[v]);i>=0;i--) {
const int &p=anc[u][i];
if(dep[p]<dep[v]) continue;
f[p][0]=f[p][1]=g[p][0]=g[p][1]=0;
for(register int j=0;j<2;j++) {
for(register int k=0;k<2;k++) {
if(!w[u][0][k]) continue;
for(register int l=0;l<2;l++) {
if(f[u][j]+seg[u][i][k][l]-(g[u][j]==w[u][0][k])>f[p][l]) {
f[p][l]=f[u][j]+seg[u][i][k][l]-(g[u][j]==w[u][0][k]);
g[p][l]=w[u][i][l];
}
}
}
}
u=p;
}
if(u==v) return std::max(f[u][0],f[u][1])-1;
for(register int i=lg2(dep[u]);i>=0;i--) {
const int &p=anc[u][i],&q=anc[v][i];
if(p==q) continue;
f[p][0]=f[p][1]=g[p][0]=g[p][1]=0;
for(register int j=0;j<2;j++) {
for(register int k=0;k<2;k++) {
if(!w[u][0][k]) continue;
for(register int l=0;l<2;l++) {
if(f[u][j]+seg[u][i][k][l]-(g[u][j]==w[u][0][k])>f[p][l]) {
f[p][l]=f[u][j]+seg[u][i][k][l]-(g[u][j]==w[u][0][k]);
g[p][l]=w[u][i][l];
}
}
}
}
f[q][0]=f[q][1]=g[q][0]=g[q][1]=0;
for(register int j=0;j<2;j++) {
for(register int k=0;k<2;k++) {
if(!w[v][0][k]) continue;
for(register int l=0;l<2;l++) {
if(f[v][j]+seg[v][i][k][l]-(g[v][j]==w[v][0][k])>f[q][l]) {
f[q][l]=f[v][j]+seg[v][i][k][l]-(g[v][j]==w[v][0][k]);
g[q][l]=w[v][i][l];
}
}
}
}
u=p;
v=q;
}
int p=anc[u][0];
f[p][0]=f[p][1]=g[p][0]=g[p][1]=0;
for(register int j=0;j<2;j++) {
for(register int k=0;k<2;k++) {
if(!w[u][0][k]) continue;
if(f[u][j]+seg[u][0][k][k]-(g[u][j]==w[u][0][k])>f[p][k]) {
f[p][k]=f[u][j]+seg[u][0][k][k]-(g[u][j]==w[u][0][k]);
g[p][k]=w[u][0][k];
}
}
}
p=0;
f[p][0]=f[p][1]=g[p][0]=g[p][1]=0;
for(register int j=0;j<2;j++) {
for(register int k=0;k<2;k++) {
if(!w[v][0][k]) continue;
if(f[v][j]+seg[v][0][k][k]-(g[v][j]==w[v][0][k])>f[p][k]) {
f[p][k]=f[v][j]+seg[v][0][k][k]-(g[v][j]==w[v][0][k]);
g[p][k]=w[v][0][k];
}
}
}
int ans=0;
for(register int i=0;i<2;i++) {
for(register int j=0;j<2;j++) {
up(ans,f[anc[u][0]][i]+f[0][j]-(g[anc[u][0]][i]==g[0][j]));
}
}
return ans-1;
}
int main() {
const int n=getint(),m=getint();
for(register int i=0;i<m;i++) {
int u=getint(),v=getint(),c=getint();
if(u>v) std::swap(u,v);
auto &p=map[std::make_pair(u,v)];
if(p.first<0) continue;
if(p.first==c||p.second==c) continue;
if(p.first==0) {
add_edge(u,v);
p.first=c;
continue;
}
if(p.second==0) {
p.second=c;
continue;
}
p.first=--tot;
p.second=0;
}
dfs(1,0);
const int q=getint();
for(register int i=0;i<q;i++) {
printf("%d\n",query(getint(),getint()));
}
return 0;
}
[CC-SUBWAY]Subway Ride的更多相关文章
- java设计模式学习笔记
简介 设计模式可以分为五类 接口型 模式:适配器模式,外观模式,合成模式,桥接模式 职责型 模式:单例模式,观察者模式,调停者模式,代理模式,职责链模式,享元模式 构造型 模式:构建者模式,工厂方法模 ...
- C#中异常:“The type initializer to throw an exception(类型初始值设定项引发异常)”的简单分析与解决方法
对于C#中异常:“The type initializer to throw an exception(类型初始值设定项引发异常)”的简单分析,目前本人分析两种情况,如下: 情况一: 借鉴麒麟.NET ...
- Android开发中常见的设计模式
对于开发人员来说,设计模式有时候就是一道坎,但是设计模式又非常有用,过了这道坎,它可以让你水平提高一个档次.而在android开发中,必要的了解一些设计模式又是非常有必要的.对于想系统的学习设计模式的 ...
- 基于Spring Boot的RESTful API实践(一)
1. RESTful简述 REST是一种设计风格,是一组约束条件及原则,而遵循REST风格的架构就称为RESTful架构,资源是RESTful的核心,一个好的RESTful架构,通过URL就能很 ...
- android 开发设计模式---Strategy模式
假设我们要出去旅游,而去旅游出行的方式有很多,有步行,有坐火车,有坐飞机等等.而如果不使用任何模式,我们的代码可能就是这样子的. 12345678910111213141516171819202122 ...
- Android开发中常见的设计模式(四)——策略模式
策略模式定义了一些列的算法,并将每一个算法封装起来,而且使它们还可以相互替换.策略模式让算法独立于使用它的客户而独立变换. 假设我们要出去旅游,而去旅游出行的方式有很多,有步行,有坐火车,有坐飞机等等 ...
- Angular material mat-icon 资源参考_Maps
ul,li>ol { margin-bottom: 0 } dt { font-weight: 700 } dd { margin: 0 1.5em 1.5em } img { height: ...
- Java_地铁购票系统
定义了两个类,在Subway类中定义三个私有数据变量,线路号,经过站点,换乘站.以及4个方法分别实现从txt文件中导入线路信息:输出线路信息:查询两个站点经过站点数,并输出经过站点以及在某站换乘几号线 ...
- 2019.9.30极限测试 04.JAVA语言课堂测试试卷-极限测试
题目存储在上传的文件当中. 代码实现 Subway 类: package ClassroomTest; public class Subway { private String railway; pr ...
- Thomson Plaza里面的三家店以及水果大会
旅行应该是一个发现的过程,至少我是这么认为的.很多时候并不一定要到什么特别的地方,也可以感受到旅游的乐趣.我觉得只要能看到值得回味的东西就好了.而能回味的东西,往往是需要仔细地来品.像旅行社安排的那样 ...
随机推荐
- μC/OS-Ⅱ在C8051F060上的移植及其应用
嵌入式操作系统是嵌入式应用的基础和核心.随着应用系统的不断复杂化和系统实时性需求的不断提高,对相应软件的逻辑结构.稳定性.实时性也提出了更高的要求,以传统的前后台编程模式编制软件将更加困难,而且容易出 ...
- ASM配置OGG
两种方法:http://blog.sina.com.cn/s/blog_aa84cfe40101lsks.html 使用ACFS配置OGG:http://ylw6006.blog.51cto.com/ ...
- mysql高可用架构 -> MHA配置binlog-server-06
前期准备 1.准备一台新的mysql实例(db03),GTID必须开启. 2.将来binlog接收目录,不能和主库binlog目录一样 停止mha masterha_stop --conf=/etc/ ...
- K/V式枚举
public enum OType { LOGIN { public String getDesc() { return "登录"; } }, ADD { public Strin ...
- mac上安装完成node,就升级好了npm,之后的设置
1.打开终端输入: npm config list 找到npmrc globalconfig /usr/local/etc/npmrc 2.打开npmrc sudo vim /usr/local/et ...
- python基础--hashlib模块
hashlib模块用于加密操作,代替了md5和sha模块, 主要提供SHA1, SHA224, SHA256, SHA384, SHA512 ,MD5 算法. # -*- coding:utf-8 - ...
- TF-图像的深度和通道的概念(转)
图像的深度和通道概念 图像的深度: 图片是由一个个像素点构成的,所有不同颜色的像素点构成了一副完整的图像,计算机存储图片是以二进制来进行的. 1 bit : 用一位来存储,那么这个像素点的取值范围就是 ...
- 20165203 2017-2018-2 《Java程序设计》第一周学习总结
20165203 2017-2018-2<Java程序设计>第一周学习总结 教材学习内容总结 (一)Java的地位 Java是面向对象编程,并涉及网络.多线程等重要的基础知识,是一门很好的 ...
- Phoenix的安装使用与SQL查询HBase
一. Phoenix的简介 1. 什么是phoenix 现有hbase的查询工具有很多如:Hive,Tez,Impala,Shark/Spark,Phoenix等.今天主要说Phoenix.phoen ...
- CF293B. Distinct Paths
B. Distinct Paths time limit per test 2 seconds memory limit per test 256 megabytes input standard i ...