【BZOJ 1093】【ZJOI 2007】【最大半连通子图】【tarjan强连通分量】【拓扑排序dp】

1093: [ZJOI2007]最大半连通子图

Time Limit: 30 Sec   Memory Limit: 162 MB
Submit: 3596   Solved: 1419
[ Submit][ Status][ Discuss]

Description

　　一个有向图G=(V,E)称为半连通的(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的余数。

Input

　　第一行包含两个整数N，M，X。N，M分别表示图G的点数与边数，X的意义如上文所述接下来M行，每行两个正整
数a, b，表示一条有向边(a, b)。图中的每个点将编号为1,2,3…N，保证输入中同一个(a,b)不会出现两次。N ≤1
00000, M ≤1000000；对于100%的数据， X ≤10^8

Output

　　应包含两行，第一行包含一个整数K。第二行包含整数C Mod X.

Sample Input

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

Sample Output

3
3

HINT

Source

[ Submit][ Status][ Discuss] 
HOME  Back

此题目是一道比较裸的题目（算是吧），开始先将一定能满足条件的强连通分量查询出来并缩成一点，然后用dp去求出最大半联通分量 ，因为一个半联通分量中点之间只需要任意两点之间一方能到达另外一方即可，而一个半联通分量若可以达到另一个‘点’（在被tarjan缩点意义后的点），且连接该点的是一个‘被走到’的点，那么该点和该半联通分量能形成新的半联通分量；那么我们确定方向，设当某一点A不会走向其他边时的其所在最大强连通分量为dp[A]，则dp[A]会被所有到他的点更新，然而更新完后，A则会去更新其他的点，那么就是用拓扑排序来确定更新顺序即可（后面会发现很多有向无环图上的dp都需要用到拓扑排序）

代码如下：

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define M 120001
using namespace std;
int n , m , mod;
int ind_[M] , nex_[M * 20] , e_[M * 20] , in[M] , cnt_;
int ind[M] , nex[M * 20] , e[M * 20] , cnt;
int dfn[M] , ouf[M] , f[M] , tim , tot;
int dp[M] , s[M] , t[M] , vis[M];
int sta[M] , top , ous[M];
queue<int> qu;
void add(int a , int b){
	nex[++cnt] = ind[a];
	e[ind[a] = cnt] = b;
}
void add_(int a , int b){
	nex_[++cnt_] = ind_[a];
	e_[ind_[a] = cnt_] = b;
	in[b]++;
}
void tarjan(int v){
	dfn[v] = ouf[v] = ++tim;
	sta[++top] = v;	
	for(int i = ind[v] ; i ; i = nex[i]){
		if(ous[e[i]])continue;
		if(!dfn[e[i]])tarjan(e[i]);
		ouf[v] = min(ouf[e[i]] , ouf[v]);
	}
	int now = 0;
	if(dfn[v] == ouf[v]){
		t[++tot] = 1;
		while(now != v){
			now = sta[top--];
			f[now] = tot;
			dp[tot]++;
			ous[now] = 1;
		}
		s[tot] = dp[tot];
	}
}

void dp_(){
	for(int i = 1 ; i <= tot ; ++i)
		if(!in[i])qu.push(i);
	while(!qu.empty()){
		int v = qu.front();
		qu.pop();
		for(int i = ind_[v] ; i ; i = nex_[i]){
			int u = e_[i];
			in[u]--;
			if(!in[u])qu.push(u);
			if(vis[u] == v)continue;
			if(dp[v] + s[u] == dp[u])t[u] = (t[u] + t[v]) % mod;
			if(dp[v] + s[u]  > dp[u])dp[u] = dp[v] + s[u] , t[u] = t[v] % mod;
			vis[u] = v;
		}
	}
}

int main(){
	int anss = 0 , anss_len = 0;
	scanf("%d%d%d",&n,&m,&mod);
	for(int i = 1 ; i <= m ; ++i){
		int u , v;
		scanf("%d%d",&u,&v);
		add(u , v);
	}
	for(int i = 1 ; i <= n ; ++i)
		if(!dfn[i])tarjan(i);
	for(int i = 1 ; i <= n ; ++i)
		for(int j = ind[i] ; j ; j = nex[j])
			if(f[i] != f[e[j]])add_(f[i] , f[e[j]]);
	dp_();
	for(int i = 1 ; i <= tot ; i++){
		if(dp[i] == anss)anss_len = (anss_len + t[i]) % mod;
		if(dp[i] >  anss)anss = dp[i] , anss_len = t[i];
	}
	printf("%d\n%d",anss,anss_len);
}