1175. 最大半连通子图（tarjan算法,dp,拓扑图，有向图强连通分量）

1175. 最大半连通子图 - AcWing题库

一个有向图 G=(V,E) 称为半连通的 (Semi-Connected)，如果满足：∀u,v∈V，满足 u→v 或 v→u，即对于图中任意两点 u,v，存在一条 u 到 v 的有向路径或者从 v 到 u 的有向路径。

若 G′=(V′,E′) 满足，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 分别表示图 G 的点数与边数，X 的意义如上文所述；

接下来 M 行，每行两个正整数 a,b，表示一条有向边 (a,b)。

图中的每个点将编号为 1 到 N，保证输入中同一个 (a,b) 不会出现两次。

输出格式

应包含两行。

第一行包含一个整数 K，第二行包含整数 C mod X。

数据范围

1≤N≤105
1≤M≤106
1≤X≤108

输入样例：

6 6 20070603
1 2
2 1
1 3
2 4
5 6
6 4

输出样例：

3
3

解析： 

性质：选择强连通分量内的所有节点一定最优（强连通的一定是半连通的，且两点之间可以互达），在经过缩点之后原图变成拓扑图，则原题就是在这个拓扑图上找一个最长链，每个点的权值就是原图中该连通分量内点的个数。

由于缩点后是拓扑图，所以再求最长链的时候就可以使用dp来处理（按照拓扑序处理）。

f[i]表示：拓扑图中以 i 节点为终点的最长链的最大距离。

g[i]表示：f[i]取最大只得方案数

当 f[j]+s[i]>f[i] 时：

f[i]=f[j]+s[i];

g[i]=g[j];

当 f[j]+s[i]==f[i] 时：

g[i]+=g[j];

算法步骤：

1）tarjan算法

2）缩点，建图同时给边判重

3）dp

#include<iostream>
#include<string>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<utility>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<math.h>
#include<map>
#include<sstream>
#include<deque>
#include<unordered_map>
#include<unordered_set>
using namespace std;
typedef long long LL;
typedef pair<int, int> PII;
const int N = 1e5 + 5, M = 2e6 + 5, INF = 0x3f3f3f3f;
int n, m, mod;
int h[N], hs[N], e[M], ne[M], idx;
int dfn[N], low[N], timestamp;
int stk[N], top;
bool in_stk[N];
int Size[N];
int id[N], scc_cnt;
int f[N], g[N];

void add(int H[], int a, int b) {
	e[idx] = b, ne[idx] = H[a], H[a] = idx++;
}

void tarjan(int u) {
	dfn[u] = low[u] = ++timestamp;
	stk[++top] = u, in_stk[u] = 1;
	for (int i = h[u]; i != -1; i = ne[i]) {
		int j = e[i];
		if (!dfn[j]) {
			tarjan(j);
			low[u] = min(low[u], low[j]);
		}
		else if (in_stk[j])low[u] = min(low[u], dfn[j]);
	}
	if (dfn[u] == low[u]) {
		int y;
		scc_cnt++;
		do {
			y = stk[top--];
			in_stk[y] = 0;
			id[y] = scc_cnt;
			Size[scc_cnt]++;
		} while (y != u);
	}
}

int main() {
	cin >> n >> m >> mod;
	memset(h, -1, sizeof h);
	memset(hs, -1, sizeof hs);
	for (int i = 1,a,b; i <= m; i++) {
		scanf("%d%d", &a, &b);
		add(h,a, b);
	}
	for (int i = 1; i <= n; i++) {
		if (!dfn[i])
			tarjan(i);
	}
	unordered_set<LL>s;
	for (int i = 1; i <= n; i++) {
		for (int j = h[i]; j != -1; j = ne[j]) {
			int k = e[j];
			LL a = id[i], b = id[k];
			LL hash = a * 100000 + b;
			if (a != b&&!s.count(hash)) {
				add(hs, a, b);
				s.insert(hash);
			}
		}
	}
	for (int i = scc_cnt; i ; i--) {
		if (!f[i]) {
			f[i] = Size[i];
			g[i] = 1;
		}
		for (int j = hs[i]; j != -1; j = ne[j]) {
			int k = e[j];
			if (f[i] + Size[k] > f[k]) {
				f[k] = f[i] + Size[k];
				g[k] = g[i];
			}
			else if (f[i] + Size[k] == f[k]) {
				g[k] = (g[i]+g[k])%mod;
			}
		}
	}
	int mx = 0, sum = 0;
	for (int i = 1; i <= scc_cnt; i++) {
		if (mx < f[i]) {
			mx = f[i];
			sum = g[i];
		}
		else if (f[i] == mx)sum = (sum + g[i]) % mod;
	}
	cout << mx << endl;
	cout << sum << endl;
	return 0;
}