hdu 4612 Warm up - 双连通 + 树形DP

15 篇文章 0 订阅
6 篇文章 0 订阅
/*
hdu 4612 Warm up - 双连通 + 树形DP

题意很简单

通过双连通缩图,构成一棵树,然后求最长链,连接首尾即可。剩下的连通块就是所求

https://www.byvoid.com/blog/biconnect
该算法是R.Tarjan发明的。对图深度优先搜索,定义DFS(u)为u在搜索树(以下简称为树)中被遍历到的次序号。定义Low(u)为u或u的子树中能通过非父子边追溯到的最早的节点,即DFS序号最小的节点。根据定义,则有:

Low(u)=Min { DFS(u) DFS(v) (u,v)为后向边(返祖边) 等价于 DFS(v)<DFS(u)且v不为u的父亲节点 Low(v) (u,v)为树枝边(父子边) }

一个顶点u是割点,当且仅当满足(1)或(2) (1) u为树根,且u有多于一个子树。 (2) u不为树根,且满足存在(u,v)为树枝边(或称父子边,即u为v在搜索树中的父亲),使得DFS(u)<=Low(v)。

一条无向边(u,v)是桥,当且仅当(u,v)为树枝边,且满足DFS(u)<Low(v)。

*/
#pragma comment(linker,"/STACK:102400000,102400000")
#include<stdio.h>  
#include<string.h>  
const int N=200000+10,M=1000000+10;
int head[N],low[N],fang[N],du[N];  
struct node  
{  
    int v,next;  
}e[M*2];//这地儿老郁闷了,没有注意到无向图存边需要的是边数量的2倍,之前老是WA  
int yong,n,r,m,dp[N][2],result;  
bool visit[M];
int min(int a,int b)  
{  
    if(a<b)  
        return a;  
    return b;  
}  
void tarjan(int u)  //之前的通过传递父节点来规避直接走回去的情况存在漏洞,当1-2,2-3,3-4,4-2时,4-2-1就退回去了
{  
    int i,j;  
    fang[u]=low[u]=yong++;  
	dp[u][0] = dp[u][1] = 0;
    for(i=head[u];i!=-1;i=e[i].next)  
    {  
        j=e[i].v;  
        if(!visit[i >> 1])
		{  
			visit[i >> 1] = true;
            if(fang[j]==0)//跟强连通一样这往下  
            {  
                tarjan(j);  
				result += fang[u] < low[j];//统计连通块,比实际数目少一个,就是回溯的时候的最后一个
				int temp = dp[j][0] + (fang[u] < low[j]);//dp[i][0]第i个节点的子树中最长的链,dp[i][1]第i个节点的子树中第二长的链,
                if(temp > dp[u][0]){
                    dp[u][1] = dp[u][0];
                    dp[u][0] = temp;
                }else if(temp > dp[u][1]){
                    dp[u][1] = temp;
                }

                low[u]=min(low[u],low[j]);  
            }else low[u]=min(low[u],fang[j]);  
        }  
    }  
}  
void adde(int u,int v)
{
	e[yong].v=v;
	e[yong].next=head[u];
	head[u]=yong++;
}
inline int max(int a,int b)
{
	if(a>b) return a;
	return b;
}
int main()
{
	int i,a,b;
	while(scanf("%d%d",&n,&m),n+m)
	{
		yong=0;  
		result=0;
        memset(fang,0,sizeof(fang));  
        memset(low,0,sizeof(low));  
        memset(head,-1,sizeof(head));  
		for(i=0;i<m;++i)
		{
			scanf("%d%d",&a,&b);
			adde(a,b);
			adde(b,a);
		}
		yong=1;  
		memset(visit, 0, sizeof(visit));
        tarjan(1);
		int temp = 0;
		for(i = 1; i <=n; ++ i){
            temp = max(temp, dp[i][0] + dp[i][1]);//+的时候没有算当前点所在的块,同时因为result也少算一个,扯平了
        }
        printf("%d\n", result - temp);
	}
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值