Tarjan缩点+拓扑排序求最长链 - 最大半连通子图「BZOJ1093」[ZJOI2007]

最大半连通子图

描述

一个有向图 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,2,3,⋯,N，保证输入中同一个 ((a,b) 不会出现两次。

输出

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

样例输入

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

样例输出

3
3

提示

对于 20% 的数据，N≤18；
对于 60%的数据，N≤10⁴；
对于 100% 的数据，1≤N≤10⁵,1≤M≤10⁶,X≤10⁸

Analysis

妙啊
Tarjan缩点+拓扑排序求最长链
注意拓扑排序的时候处理一下两个联通块之间有重边的情况，因为要统计方案数，所以重边要特殊处理
还有就是不可以一边拓扑乱搞，一边统计答案
因为一个点可能被多次算入答案

Code

#include<bits/stdc++.h>
#define in read()
#define N 200009
#define M 2000009
using namespace std;
inline int read(){
    int data=0;int w=1; char ch=0;
    ch=getchar();
    while(ch!='-' && (ch<'0' || ch>'9')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0' && ch<='9') data=(data<<3)+(data<<1)+ch-'0',ch=getchar();
    return data*w;
}
int n,m,X,vis[N],f[N],du[N],g[N];
int nxt[M],to[M],head[N],ecnt=0;
inline void add(int x,int y){nxt[++ecnt]=head[x];head[x]=ecnt;to[ecnt]=y;}
int Nxt[M],To[M],Head[N],Ecnt=0;
inline void readd(int x,int y){Nxt[++Ecnt]=Head[x];Head[x]=Ecnt;To[Ecnt]=y;}
int res=0,ans=-1;
int dfn[N],low[N],be[N],sze[N],dfs=0,num=0;
bool insta[N];
stack<int> S;
inline void tarjan(int u){
	S.push(u);insta[u]=1;
	dfn[u]=low[u]=++dfs;
	for(int e=head[u];e;e=nxt[e]){
		int v=to[e];
		if(!dfn[v]) {tarjan(v);low[u]=min(low[v],low[u]);}
		else if(insta[v]) low[u]=min(low[u],dfn[v]);
	}
	if(low[u]==dfn[u]){
		++num;int x;
		do{
			x=S.top();S.pop();insta[x]=0;
			be[x]=num;sze[num]++;
		}while(x!=u);
	}
}
int main(){
	n=in;m=in;X=in;
	int i,j,a,b;
	for(i=1;i<=m;++i){
		a=in;b=in;
		add(a,b);
	}
	for(i=1;i<=n;++i)
		if(!dfn[i]) tarjan(i);
	for(i=1;i<=n;++i){
		for(int e=head[i];e;e=nxt[e]){
			j=to[e];
			if(be[i]!=be[j]) readd(be[i],be[j]),du[be[j]]++;
		}
	}
	memset(f,-1,sizeof(f));
	queue<int> q;
	for(i=1;i<=num;++i)	if(!du[i]) f[i]=sze[i],g[i]=1,q.push(i);
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int e=Head[u];e;e=Nxt[e]){
			int v=To[e];--du[v];
			if(!du[v]) q.push(v);
			if(vis[v]==u) continue;
			if(f[u]+sze[v]>f[v]){
				f[v]=f[u]+sze[v];
				g[v]=g[u];
			}
			else if(f[u]+sze[v]==f[v]) g[v]=(g[v]+g[u])%X;
			vis[v]=u;
		}
	}
	for(i=1;i<=num;++i)
	{
		if(f[i]>ans){ans=f[i];res=g[i];}
		else if(f[i]==ans) res=(res+g[i])%X;
	}
	cout<<ans<<'\n'<<res;
	return 0;
}