[ZJOI2007]最大半连通子图(拓扑排序 + tarjin +DP)

C-[ZJOI2007]最大半连通子图_2022图论班第二章连通性例题与习题 (nowcoder.com)

题目描述t
一个有向图G=(VE)称为半连通的(Semi-Connected)，如果满足: u,v∈V，满足u→v或v一→u，即对于图中任意两点u，v,存在一条u到v的有向路径或者从v到u的有向路径。
若G'=(V,E')满足V∈V，E'是E中所有跟V有关的边，则称G'是G的一个导出子图。若G'是G的导出子图，且G半连通，则称G'为G的半连通子图。
若G'是G所有半连通子图中包含节点数最多的，则称G'是G的最大半连通子图。
给定一个有向图G，请求出G的最大半连通子图拥有的节点数K，以及不同的最大半连通子图的数目C。由于C可能比较大，仅要求输出C对X的余数。
输入描述:
第一行包含两个整数N，M，X。N，M分别表示图c的点数与边数，x的意义如上文所述接下来M行，每行两个正整数a,b，表示一条有向边(a,b)。
图中的每个点将编号为1,2，3...N，保证输入中同一个(a,b)不会出现两次。N≤100000,M≤ 1000000;对于100%的数据，x≤10^8
输出描述:
应包含两行，第一行包含一个整数K。第二行包含整数c Mod x.
示例1
输入
复制
6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4

输出

3
3

题解:

根据题意,我们发现在一个强联通分量内,点可以互相到达,肯定是满足题意的,所以我们把他考虑为一个点,用tarjin算法找到所有的强联通分量

其次我们可以发现他要我们找的图是一条链(强联通分量看作点),所以我们应该找一条最长的链,并找到有多少这样的链

进行缩点后,利用拓扑排序,我们可以找到最长的链的长度,同时进行DP

DP的思想如下:

f[i]i之前链的最大长度

ff[i]有多少种方式可以达到f[i]

如果f[u] + 此时点(不是真的点,强连通分量的大小)的数目 > f[v]

f[v] = f[u] + siz[v]

ff[v] = ff[u]

如果f[u] + 此时点(不是真的点,强连通分量的大小)的数目 =  f[v]

ff[v] = ff[u] + ff[v]

最终把强连通分量中,最长链数目 = ma 的ff[i]相加即可

#include<iostream>
#include<algorithm>
#include<string>
#include<queue>
#include<vector>
#include<map>
#include<cstring>
#include<cmath>
#include<stack>
#include<set>
using namespace std;
#define int long long
typedef pair<int,int> PII;
int n,m,x;
const int N = 2e5 + 10;
vector<int> p[N];
int dfn[N];
int low[N];
int vis[N];
int idx;
stack<int> s;
int cnt;
int siz[N];
int color[N];
vector<int> ne[N]; 
int rd[N];
int f[N];
int ff[N];
int pre[N];
int ma;
int cd[N];
queue<int> q;
void tarjin(int u)
{
	dfn[u] = low[u] = ++idx;
	vis[u] = 1;
	s.push(u);
	for(auto v:p[u])
	{
		if(!dfn[v])
		{
			tarjin(v);
			low[u] = min(low[u],low[v]);
		}
		else if(vis[v])
		{
			low[u] = min(low[u],dfn[v]);
		}
	}
	if(low[u] == dfn[u])
	{
		cnt++;
		while(1)
		{
			int t = s.top();
			s.pop();	
			color[t] = cnt;
			siz[cnt] ++;
			vis[t] = 0;
			if(t == u)
			break;
		}
	}
}
void toupu()
{
	while(q.size())
	{
		int u = q.front();
		q.pop();
		ma = max(ma,f[u]);
		for(auto v:ne[u])
		{
			rd[v]--;
			if(rd[v] == 0)
			{
				q.push(v);
			}
			if(pre[v] == u)
			{
				continue;
			}
			if(f[u] + siz[v] > f[v])
			{
				f[v] = f[u] +siz[v];
				ff[v] = ff[u];
			}
			else if(f[u] + siz[v] == f[v])
			{
				ff[v] = (ff[u] + ff[v])%x;
			}
			pre[v] = u;
		}
	}
}
void solve()
{
	cin >> n >> m >> x;
	for(int i = 1;i <= m;i++)
	{
		int l,r;
		cin >> l >> r;
		p[l].push_back(r);
	}
	for(int i = 1;i <= n;i++)
	{
		if(!dfn[i])
		tarjin(i);
	}
	for(int i = 1;i <= n;i++)
	{
		for(auto j : p[i])
		{
			int u = color[i];
			int v = color[j];
			if(u!=v)
			{
				ne[u].push_back(v); 
				cd[u]++;
				rd[v]++;
			}
		}
	}
	for(int i = 1;i <= cnt;i++)
	{
		if(!rd[i])
		{
			ff[i] = 1;
			f[i] = siz[i];
			q.push(i);
		}
	}
	toupu();
	int ans = 0;
	for(int i = 1;i <= cnt;i++)
	{
		if(f[i] == ma)
		{
			ans = (ans + ff[i])%x;
		}
	}
	cout <<ma<<"\n"<<ans;
	
	
}
//1 -7 1 -7
signed main(){
//	ios::sync_with_stdio(false);
//	cin.tie(0);
//	cout.tie(0);
	int t = 1;
//	cin >> t;
    while(t--)
	{

		solve();
	} 
}
//3 3 2 2 1 0
//1 2 1 1 1 0