【题解】洛谷P3225 矿场搭建（割点 tarjan）

最新推荐文章于 2022-03-23 14:58:35 发布

原创最新推荐文章于 2022-03-23 14:58:35 发布 · 326 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#题解 #洛谷 #tarjan #割点

题解同时被 3 个专栏收录

140 篇文章

订阅专栏

tarjan

9 篇文章

订阅专栏

割点

3 篇文章

订阅专栏

本文介绍了一种用于寻找图中的割点和割边的算法，并通过具体代码实现了基于Tarjan算法的割点和割边检测。文章还讨论了如何通过深度优先搜索(DFS)来确定图中割点和割边的数量，以及如何基于这些信息进一步计算特定类型的子图数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

https://www.luogu.org/blog/cjyyb/solution-p3225

思路还是比较好理解的，不过实现起来挺复杂。。。

虽然数据水但还是注意初始化别搞错了。。我的Cut数组初始化为0时sizeof写的是cut，看了半天没找出来。。。

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#define ll long long
using namespace std;
const int maxn=610;
int dfn[maxn],low[maxn],color[maxn];
bool Cut[maxn];
ll head[maxn],nnext[maxn*4],to[maxn*4];
int tot;
ll num,cut,id,n,m,rs,root,ans1,ans2,group,Case;
void add(ll x,ll y)
{
	tot++;
	nnext[tot]=head[x];
	head[x]=tot;
	to[tot]=y;
}
void ycl()
{
	memset(Cut,0,sizeof(Cut));
	memset(head,0,sizeof(head));
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
	memset(color,0,sizeof(color));
	ans1=group=n=tot=id=0;
	ans2=1;
}
void tarjan(int x,int fa)
{
	dfn[x]=low[x]=++id;
	for(int i=head[x];i;i=nnext[i])
	{
		int y=to[i];
		if(dfn[y]==0)
		{
			tarjan(y,x);
			if(y==fa) continue;
			low[x]=min(low[x],low[y]);
			if(low[y]>=dfn[x])
			{
				if(x!=root) Cut[x]=true;
				else rs++;
			}
		}
		else
		{
			low[x]=min(low[x],low[y]);
		}
	}
}
void dfs(int x)
{
	color[x]=group;
	num++;
	for(int i=head[x];i;i=nnext[i])
	{
		int y=to[i];
		if(Cut[y]&&color[y]!=group)
		{
			cut++;
			color[y]=group;
		}
		if(!color[y]) dfs(y);
	}
}
int main()
{
	Case=1;
	ll x,y;
	while(cin>>m&&m)
	{
		
		ycl();
		for(int i=1;i<=m;i++)
		{
			
			cin>>x>>y;
			add(x,y);
			add(y,x);
			n=max(n,x);
			n=max(n,y);
		}
		for(int i=1;i<=n;i++)
		{
			if(!dfn[i])
			{
				root=i;
				rs=0;
				tarjan(i,i);
				if(rs>=2) Cut[i]=true;
			}
		}
		for(int i=1;i<=n;i++)
		{
			if(!color[i]&&!Cut[i])
			{
				++group;
				num=cut=0;
				dfs(i);
				
				if(cut==0)
				{
					ans1+=2;
					ans2*=(num-1)*num/2; 
				}
				if(cut==1)
				{
					ans1+=1;
					ans2*=num;
				}
				if(cut>=2)
				{
					;
				}
			}
		}
		cout<<"Case "<<Case++<<": "<<ans1<<" "<<ans2<<endl;
	}
	return 0;
}