图中异色点对最短距离(最小生成树+线段树)

problem

给定一个 nnn 个点,mmm 条边的无向连通图,图有边权,每个点有一个颜色。

qqq 次操作,每次操作可更改某一个点颜色。

求每次操作后图中不同颜色点之间的最短距离。

若图中点颜色全相同,输出 000

一行给定 n,m,c,qn,m,c,qn,m,c,q。接下来 mmm 行形如 (u,v,w)(u,v,w)(u,v,w),表示 u,vu,vu,v 间有一条边长 www

接下来一行 nnn 个数表示颜色。最后 qqq 行,每行形如 x,yx,yx,y,把 xxx 点颜色改为 yyy

c≤n≤2e5,m≤3e5,q≤2e5,1≤w≤1e6c\le n\le 2e5,m\le 3e5,q\le 2e5,1\le w\le 1e6cn2e5,m3e5,q2e5,1w1e6

solution

observation:因为 www 是正的,所以答案一定是某条边的权值。

observation:这条边一定在图的随便一棵 MST\text{MST}MST(最小生成树)上。

在这里插入图片描述

只要发现了这两个性质,这道题就转化成了求树上任意两异色相邻点的最短距离。

我们考虑每个点对答案的贡献,即求出这个点儿子中距离最近且与之异色的儿子。

怎么快速求得?

  • 将该点的所有儿子扔到线段树上,线段树合并颜色信息。

    具体而言,如果两段相邻颜色不同,合并后颜色赋零,否则赋成同样颜色。

    并且儿子要按照距离该点的距离有序排列重编号后再以重编号建线段树。

    直接在线段树上走与该点颜色不同的区间,先走左儿子(距离最近)。

    动态开点哦~

然后将所有点对答案的贡献放到一个优先队列里面,最短距离就从这中间产生。

现在考虑修改某个点颜色后的变化。显然只会影响该点的贡献以及该点父亲的贡献。

  • 该点贡献。直接在其对应线段树上重新以新颜色为标准找即可。

  • 该点父亲贡献。直接在父亲对应线段上修改该点的信息,然后父亲再查一遍。

  • 将两者产生的新贡献直接覆盖原来的贡献。

    具体而言:可以开一个数组记录每个点实时贡献,最后输出最短距离时候判断一下存的贡献和对应点现在提供的贡献是否相同。不同扔掉即可。

时间复杂度就只带个 log⁡\loglog

code

#include <bits/stdc++.h>
using namespace std;
#define maxn 200005
#define maxm 300005
#define inf 0x3f3f3f3f
#define Pair pair < int, int >
struct edge { int u, v, w; }E[maxm];
vector < Pair > G[maxn];
priority_queue < Pair, vector < Pair >, greater < Pair > > ans;
int n, m, C, Q, cnt;
int c[maxn];
int head[maxn], to[maxn], nxt[maxn], len[maxn];
int f[maxn], dis[maxn], siz[maxn], id[maxn];int find( int x ) {return x == f[x] ? x : f[x] = find( f[x] );
}void addedge( int u, int v, int w ) {len[cnt] = w, to[cnt] = v, nxt[cnt] = head[u], head[u] = cnt ++;
}void dfs( int u, int fa ) {f[u] = fa;for( int i = head[u];~ i;i = nxt[i] )if( to[i] ^ fa ) {dfs( to[i], u );G[u].push_back( make_pair( len[i], to[i] ) );siz[u] ++;}sort( G[u].begin(), G[u].end() );for( int i = 1;i <= siz[u];i ++ ) id[G[u][i - 1].second] = i;//根据排序 对于u而言 id越小的儿子离u越近//所以在线段树上尽可能往左走
}int root[maxn];
namespace sgt {struct node { int lson, rson, color; }t[maxn * 30];#define lson t[now].lson#define rson t[now].rson#define mid (l + r >> 1)void pushup( int now ) {if( ! lson ) t[now].color = t[rson].color;else if( ! rson ) t[now].color = t[lson].color;else {if( t[lson].color ^ t[rson].color ) t[now].color = 0;else t[now].color = t[rson].color;}}void modify( int &now, int l, int r, int pos, int col ) {if( ! now ) now = ++ cnt;if( l == r ) { t[now].color = col; return; }if( pos <= mid ) modify( lson, l, mid, pos, col );else modify( rson, mid + 1, r, pos, col );pushup( now );}int query( int now, int l, int r, int col ) {if( t[now].color == col ) return inf;if( l == r ) return l;if( t[lson].color ^ col ) return query( lson, l, mid, col );//先找左else return query( rson, mid + 1, r, col );  }
}int main() {scanf( "%d %d %d %d", &n, &m, &C, &Q );for( int i = 1, u, v, w;i <= m;i ++ ) {scanf( "%d %d %d", &u, &v, &w );if( u == v ) continue;E[i] = (edge){ u, v, w };}for( int i = 1;i <= n;i ++ ) scanf( "%d", &c[i] );sort( E + 1, E + m + 1, []( edge x, edge y ) { return x.w < y.w; } );iota( f + 1, f + n + 1, 1 );memset( head, -1, sizeof( head ) );for( int i = 1;i <= m;i ++ ) {int u = E[i].u, v = E[i].v, w = E[i].w;int fu = find( u ), fv = find( v );if( fu ^ fv ) {f[fv] = fu;addedge( u, v, w );addedge( v, u, w );}}dfs( 1, 0 ); cnt = 0;//对每个点 将其所有儿子建成一棵线段树 上面储存颜色信息for( int i = 2;i <= n;i ++ ) sgt :: modify( root[f[i]], 1, siz[f[i]], id[i], c[i] );for( int i = 1;i <= n;i ++ )if( siz[i] ) {int w = sgt :: query( root[i], 1, siz[i], c[i] );if( w == inf ) dis[i] = inf;else dis[i] = G[i][w - 1].first, ans.push( make_pair( dis[i], i ) );//儿子重编号是1~siz[i]对应vector的下标都要-1}for( int i = 1, x, y;i <= Q;i ++ ) {scanf( "%d %d", &x, &y );c[x] = y;if( siz[x] ) {//改和儿子间的答案int w = sgt :: query( root[x], 1, siz[x], c[x] );if( w == inf ) dis[x] = inf;else if( dis[x] != G[x][w - 1].first ) {dis[x] = G[x][w - 1].first;ans.push( make_pair( dis[x], x ) );}}if( f[x] ) {//改和父亲的答案sgt :: modify( root[f[x]], 1, siz[f[x]], id[x], c[x] ); //在父亲的线段树上修改x的信息int w = sgt :: query( root[f[x]], 1, siz[f[x]], c[f[x]] );if( w == inf ) dis[f[x]] = inf;else if( dis[f[x]] != G[f[x]][w - 1].first ) {dis[f[x]] = G[f[x]][w - 1].first;ans.push( make_pair( dis[f[x]], f[x] ) );}}while( ! ans.empty() ) {//有些点颜色改变信息不再正确 通过一直是正确的dis来弹出错误备选答案if( dis[ans.top().second] ^ ans.top().first ) ans.pop();else break;}if( ans.empty() ) puts("0");else printf( "%d\n", ans.top().first );}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/316486.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Defend Your Country

Defend Your Country 题意&#xff1a; n个点&#xff0c;m条边的简单无向连通图&#xff0c;每个点一个权值ai,一个连通块的贡献&#xff1a;(−1)块内点数∗∑ai[点i在该连通块内](-1)^{块内点数}*\sum a_{i}[点i在该连通块内](−1)块内点数∗∑ai​[点i在该连通块内] 可以…

SkyWalking Liunx 环境搭建NetCore接入

背景前两天看见有小哥介绍windows下安装skywalking的介绍地址。正好最近也在搭建linux环境的SkyWalking&#xff0c;顺便把linux环境搭建的经验分享下&#xff0c;帮助下使用linux部署DotNetCore项目的同学。介绍SkyWalking是开源的apm工具&#xff0c;服务器端使用java编写&am…

.NET 基金会完成第一次全面改选

.NET基金会是一个独立的组织&#xff0c;支持.NET社区和开源&#xff0c;旨在拓宽和加强.NET生态系统和社区。这可以通过多种方式完成&#xff0c;包括项目指导&#xff0c;指导&#xff0c;法律和营销帮助&#xff0c;技术和财务支持设置等&#xff0c;2014年微软组织成立.NET…

xay loves trees

xay loves trees 题意&#xff1a; 有两棵树&#xff0c;现在让你找到一个最大的点集合S&#xff0c;要求S中的点在第一棵树中任意两点存在祖先儿子关系且所有点是连接的&#xff0c;在第二棵树中任意两点都不存在祖先儿子关系&#xff0c;问S集合的最大是多少&#xff1f; …

特来电混沌工程实践

一、导语随着大型分布式系统架构的演进和广泛应用&#xff0c;软件工程的最佳实践也随之改变。我们通过分布式、服务化、DevOps、敏捷开发&#xff0c;快速响应业务的需求变化&#xff0c;支持大规模分布式应用。但这些做法带来效益的同时&#xff0c;也带来了另一个紧迫问题&a…

[ZJOI2007] 棋盘制作(单调栈 / DP悬线法)

problem 洛谷链接 solution1-单调栈 很容易想到&#xff0c;预处理出每个点向上最大能延伸的长度&#xff0c;然后对每个点求一个矩阵面积。 然后思考优化&#xff0c;不难想到每次对一行进行求解。 每一行的所有列一起构成了一个直方图。 直方图直接经典笛卡尔树。笛卡尔…

构建现代Web应用时究竟是选择传统web应用还是SPA

在大前端盛行的今天&#xff0c;似乎前后端分离的开发模式才是大势所趋&#xff0c;而SPA的概念更是应运而生。现在随便构建一个web应用程序如果你不是使用SPA的话&#xff0c;就会感觉有点low&#xff0c;但是真的是这样吗&#xff1f;今天这篇文章我们就来一起探讨下&#xf…

你所不知道的ASP.NET Core MVC/WebApi基础系列(二)

冒个泡&#xff0c;算起来估计有很长时间没更新公众号了&#xff0c;估计是我第一次停更如此之久&#xff0c;人总有懒惰的时候&#xff0c;时间越长越懒惰&#xff0c;但是呢&#xff0c;不学又不行&#xff0c;持续的惰性是不行dei&#xff0c;要不然会被时光所抛弃&#xff…

POJ-2069 Super Star(最小球覆盖)

POJ-2069 Super Star 题意&#xff1a; 给你n个点&#xff0c;求覆盖所有点的最小球的半径 4<n<30 题解&#xff1a; 求最小球覆盖的步骤&#xff1a; &#xff08;1&#xff09;对于一个点&#xff1a;球心就是这个点&#xff0c;且半径无穷小。 &#xff08;2&…

ASP.NET Core 3.0 上的gRPC服务模板初体验(多图)

早就听说ASP.NET Core 3.0中引入了gRPC的服务模板&#xff0c;正好趁着家里电脑刚做了新系统&#xff0c;然后装了VS2019的功夫来体验一把。同时记录体验的过程。如果你也想按照本文的步骤体验的话&#xff0c;那你得先安装.NET Core3.0预览版的SDK。至于开发工具我用的时VS201…

结合使用 Draft 与 Tencent Kubernetes Engine (TKE)

Draft 是一种开源工具&#xff0c;有助于在 Kubernetes 群集中打包和部署应用程序容器&#xff0c;让你专注于开发周期 - 专注开发的“内部循环”。 在开发代码期间&#xff0c;但尚未将代码提交到版本控制之前&#xff0c;Draft 将会运行。 借助 Draft&#xff0c;可在代码发生…

[ZJOI2010] 贪吃的老鼠(二分+差分+神仙建图网络流)

problem luogu-P2570 solution 卧槽网络流尼玛神题 首先这个最小延长时间 TTT &#xff0c;套路地考虑二分&#xff0c;将问题转化为判定性问题。 其次 n,mn,mn,m 和奶酪存在时间 [l,r][l,r][l,r] 的量级差很大&#xff0c;我们肯定会猜想一段时间内选择吃奶酪的老鼠是一样…

基于IdentityServer的系统对接微信公众号

业务需求公司有两个业务系统&#xff0c;A和B&#xff0c;AB用户之间属于多对一的关系&#xff0c;数据库里面也就是两张表&#xff0c;A表有个外键指向B。现在需要实现以下几个功能。A用户扫描B的二维码&#xff0c;填写相关的注册信息&#xff0c;注册完成之后自动属于B。也就…

Ocelot Api网关教程(9)- QoS

本文介绍Ocelot中的QoS(Quality of Service)&#xff0c;其使用了Polly对超时等请求下游失败等情况进行熔断。1、添加Nuget包添加 Ocelot.Provider.Polly 到OcelotGetway项目中2、修改 Startup.ConfigureServices 如下来添加Polly&#xff1a;services .AddOcelot(new Confi…

cf1553F. Pairwise Modulo

cf1553F. Pairwise Modulo 题意&#xff1a; 给你一个数组a&#xff0c;a由n个不同的数组成,让你求出对应的数组p 数组p的定义为&#xff1a; pk∑1≤i,j≤kaimodajp_{k}\sum_{1\leq i,j\leq k}a_{i} \mod a_{j}pk​∑1≤i,j≤k​ai​modaj​ 题解&#xff1a; 官方题解 首…

Abp框架准备加入.NET Foundation

Abp团队正准备尝试将Abp项目加入.NET Foundation!既然申请了,必然有很大的自信能够成功,请大家等待好消息!Abp中文网会第一时间跟进.点击阅读原文查看最新进展.

[ZJOI2014] 璀璨光华(bfs建图 + dfs搜索)

problem luogu-P3342 solution 你感觉这道题没考什么&#xff0c;又感觉考了什么 通过样例以及题面&#xff0c;我们并未获取到『立方体每个小方块的编号是按一定规则命名』的信息。 也就是说&#xff0c;我们需要通过输入的每个小方块相邻的编号的信息来建出这个立方体的…

SignalR 中丰富多彩的消息推送方式

在上一篇 SignalR 文章中&#xff0c;演示了如何通过 SignalR 实现了简单的聊天室功能&#xff1b;本着简洁就是美的原则&#xff0c;这一篇我们也来聊聊在 SignalR 中的用户和组的概念&#xff0c;理解这些基础知识有助于更好的开发基于 SignalR 的应用&#xff0c;通过对用户…

Visual Studio 2019 使用 Live Share

一.前言Visual Studio 2019 在今天发布&#xff08;北京时间&#xff09;了&#xff0c;这次带来了一个比较有趣的 Live Share 功能&#xff0c;使用它可以进行更好的协作开发。主要功能&#xff1a;更多资料可看官方介绍&#xff1a;Visual Studio 实时共享什么是Visual Studi…

cf1523B. Lord of the Values

cf1523B. Lord of the Values 题意&#xff1a; 给你一个数组&#xff0c;有n个数&#xff0c;n为偶数&#xff0c;a1&#xff0c;a2…an 现在有两个操作&#xff1a; 对于i<j 操作1&#xff1a;aiaiaj 操作2&#xff1a;ajaj-ai 把原数组转换为-a1,-a2,-a3… 题解&#…