zoukankan      html  css  js  c++  java
  • 点双连通分量

    它是什么?

    对于一个无向图,如果它没有割点,则称其为“点双联通图”

    无向图的极大点双连通子图成为“点双连通分量”

    它可以用来做什么?

    如果对无向图中的所有点双连通分量缩点,可以得到一颗树,可以比较方便地将一些路径问题转化为树上问题

    怎么求?

    我们可以在(Tarjan)求割点时,顺便把所有(v-DCC)求出来,流程如下:

    1. 当一个节点第一次被访问时,入栈
    2. (dfn[x]leq low[v])成立(即割点判定法则)时
      从弹栈,直至节点(v)被弹出,所有弹出的节点与(x)共同构成一个(v-DCC)

    缩点相对复杂一些这就是我一直不会写点双的原因其实好像也不是很难,因为一个割点可能属于多个(v-DCC)

    设图中一共有(k)个割点,(t)(v-DCC),建立一张包含(k+t)个节点的新图,把每个(v-DCC)和割点作为新图中的节点,并在在每个割点和所有包含它的(v-DCC)之间连边

    最终获得的图即是一颗树(或是森林)


    例题

    BZOJ3331 压力

    TLDR一句话题意:
    给定(n)个点(m)条边的无向连通图,以及(q)个点对
    对于每个点对,它们的路径上有一些点是必经点(包括起点和终点)
    对于图上的每个点,求出它们是多少个点对路径的必经点

    Description

    如今,路由器和交换机构建起了互联网的骨架。处在互联网的骨干位置的
    核心路由器典型的要处理(100Gbit/s)的网络流量。他们每天都生活在巨大的压力
    之下。
    小强建立了一个模型。这世界上有(N)个网络设备,他们之间有M个双向的
    链接。这个世界是连通的。在一段时间里,有(Q)个数据包要从一个网络设备发
    送到另一个网络设备。
    一个网络设备承受的压力有多大呢?很显然,这取决于(Q)个数据包各自走
    的路径。不过,某些数据包无论走什么路径都不可避免的要通过某些网络设
    备。
    你要计算:对每个网络设备,必须通过(包括起点、终点)他的数据包有
    多少个?

    Input

    第一行包含3个由空格隔开的正整数N,M,Q。
    接下来M行,每行两个整数u,v,表示第u个网络设备(从1开始编号)和
    第v个网络设备之间有一个链接。u不会等于v。两个网络设备之间可能有多个
    链接。
    接下来Q行,每行两个整数p,q,表示第p个网络设备向第q个网络设备发
    送了一个数据包。p不会等于q。

    Output

    输出N行,每行1个整数,表示必须通过某个网络设备的数据包的数量。

    Sample Input

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

    Sample Output

    2
    1
    1
    2

    HINT

    【样例解释】
    设备(1,2,3)之间两两有链接,(4)只和(1)有链接。(4)想向(2)(3)各发送一个数据包。显然,这两个数据包必须要经过它的起点、终点和(1)

    【数据规模和约定】
    对于(40\%)的数据,(N,M,Q≤2000)
    对于(60\%)的数据,(N,M,Q≤40000)
    对于(100\%)的数据,(N≤100000,M,Q≤200000)

    Solution

    显然对于原图的一个点,如果它不是割点,那么除了作为起点和和终点,它不会成为任何一条路径的必经点,否则它一定是经过它的路径的必经点

    对原图求点双并进行缩点,获得一棵树,新图中的点只有割点和各个点双
    对于非割点,统计它是多少条路径的起点和终点
    对于割点,计算它被多少条路径覆盖,LCA求路径,树上差分即可
    这样,我们就成功地把原问题转化为树上问题

    Code

    #include <cstdio>
    #include <algorithm>
    #include <cstdlib>
    #include <cstring>
    #include <cmath>
    #include <iostream>
    #include <vector>
    #define maxn 200005//不知道为什么,如果只开100000在BZOJ上会REemmmmm
    #define maxm 200005
    using namespace std;
    typedef long long ll;
    
    int n,m,q;
    int ans[maxn],c[maxn];
    
    struct edge
    {
    	int u,v,nxt;
    }g[maxm*2];
    
    int nn;//强连通分量编号
    vector<int> dcc[maxn];//使用vector记录每个v-DCC所包含的节点
    
    int head[maxn],ecnt;
    void eADD(int u,int v)
    {
    	g[++ecnt].u=u;
    	g[ecnt].v=v;
    	g[ecnt].nxt=head[u];
    	head[u]=ecnt;
    }
    
    int dfn[maxn],low[maxn],sign;
    int stk[maxn],top;
    int bel[maxn],cut[maxn];
    int root;
    void Tarjan(int u,int fa)
    {
    	dfn[u]=low[u]=++sign;
        if(u==root && head[u]==0)
        {
        	dcc[++nn].push_back(u);
            return;
        }//孤立点
    	stk[++top]=u;
    	int flag=0;
    	for(register int i=head[u];i;i=g[i].nxt)
    	{
    		int v=g[i].v;
    		if(v==fa)
    			continue;
    		if(!dfn[v])
    		{
    			Tarjan(v,u);
    			low[u]=min(low[u],low[v]);
    			if(low[v]>=dfn[u])
    			{
    				++flag;
    				if(u!=root || flag>1)
    					cut[u]=1;
    				++nn;
    				while(stk[top+1]!=v)
    					dcc[nn].push_back(stk[top--]);
    				dcc[nn].push_back(u);
    			}
    		}
    		else
    			low[u]=min(low[u],dfn[v]);
    	}
    }
    
    int prt[maxn][21],dep[maxn];
    void dfs0(int u,int fa,int depth)
    {
    	prt[u][0]=fa;
    	dep[u]=depth;
    	for(register int i=1;(1<<i)<=depth;++i)
    		prt[u][i]=prt[prt[u][i-1]][i-1];
    	for(register int i=head[u];i;i=g[i].nxt)
    	{
    		int v=g[i].v;
    		if(v!=fa)
    			dfs0(v,u,depth+1);
    	}
    }
    
    int LCA(int x,int y)
    {
    	if(dep[x]<dep[y])
    		swap(x,y);
    	for(register int i=20;i>=0;--i)
    		if(dep[prt[x][i]]>=dep[y])
    			x=prt[x][i];
    	if(x==y)
    		return x;
    	for(register int i=20;i>=0;--i)
    		if(prt[x][i]!=prt[y][i])
    		{
    			x=prt[x][i];
    			y=prt[y][i];
    		}
    	return prt[x][0];
    }
    
    void dfs1(int u,int fa)
    {
    	for(register int i=head[u];i;i=g[i].nxt)
    	{
    		int v=g[i].v;
    		if(v!=fa)
    		{
    			dfs1(v,u);
    			c[u]+=c[v];
    		}
    	}
    }
    
    int main()
    {
    	scanf("%d%d%d",&n,&m,&q);
    	for(register int i=1;i<=m;++i)
    	{
    		int u,v;
    		scanf("%d%d",&u,&v);
    		eADD(u,v),eADD(v,u);
    	}
    	for(register int i=1;i<=n;++i)
    		if(!dfn[i])
    			Tarjan(root=i,0);
    	int new_id=nn;
    	for(register int i=1;i<=n;++i)
    		if(cut[i])
    			bel[i]=++new_id;//为每个割点分配单独的bel
    	memset(g,0,sizeof(g));
    	memset(head,0,sizeof(head));
    	ecnt=0;//清空原图,建立新图
    	for(register int i=1;i<=nn;++i)
    		for(register int j=0;j<dcc[i].size();++j)
    			if(cut[dcc[i][j]])
    				eADD(bel[dcc[i][j]],i),eADD(i,bel[dcc[i][j]]);
    			else
    				bel[dcc[i][j]]=i;//非割点只会属于一个点双
    	dfs0(1,0,1);
    	while(q--)
    	{
    		int u,v;
    		scanf("%d%d",&u,&v);
    		if(bel[u]==bel[v])//特判
    		{
    			++ans[u],++ans[v];
    			continue;
    		}
    		int lca=LCA(bel[u],bel[v]);
    		++c[bel[u]],++c[bel[v]];
    		if(!cut[u])++ans[u];//若为割点,则不需要单独增加答案
    		if(!cut[v])++ans[v];
    		c[lca]--,c[prt[lca][0]]--;//差分
    	}
    	dfs1(1,0);
    	for(register int i=1;i<=n;++i)
    	{
    		if(cut[i])
    			printf("%d
    ",c[bel[i]]+ans[i]);
    		else
    			printf("%d
    ",ans[i]);
    	}
    	return 0;
    }
    
  • 相关阅读:
    python爬取斗图网中的 “最新套图”和“最新表情”
    SpringBoot (1) idea下的环境搭建及demo
    python爬取视频网站m3u8视频,下载.ts后缀文件,合并成整视频
    微信小程序—day05
    从零起步做到Linux运维经理, 你必须管好的23个细节
    前后端分离原理
    图文并茂|为你揭开微服务架构的“神秘面纱”!
    swarm集群日常部分操作
    OpenStack 部署运维实战
    京东618:Docker扛大旗,弹性伸缩成重点 (2015-06-23)
  • 原文地址:https://www.cnblogs.com/lizbaka/p/v-DCC.html
Copyright © 2011-2022 走看看