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
1 2
2 1
1 3
2 4
5 6
6 4
Sample Output
3
3
3
Solution
这一题可以看出强联通分量一定是半联通分量,所以先用Tarjan缩点,然后重新构图,这样就变成了一个DAG(有向无环图),然后拓扑排序一边顺便dp记录结果就好了。注意联通块缩成点之后可能会有重边,拓扑排序的时候要注意不能算多次。
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
inline int read(){
int xx=0,f=1;char ch=getchar();
for(;ch<'0'||ch>'9';ch=getchar())if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=getchar())xx=xx*10+ch-'0';
return xx*f;
}
const int maxn=100010,maxm=1000010;
struct edge{
int to,next;
}e[maxm],e2[maxm];
int n,m,p,num,num2,cnt,head[maxn],head2[maxn],dfn[maxn],low[maxn],len[maxn],sum[maxn];
int sta[maxn],top,ltk[maxn],size[maxn],ans1,ans2,du[maxn],last[maxn];
void add(int u,int v){
e[++num].to=v;e[num].next=head[u];head[u]=num;
}
void add2(int u,int v){
du[v]++;
e2[++num2].to=v;e2[num2].next=head2[u];head2[u]=num2;
}
void dfs(int x){
dfn[x]=low[x]=++cnt;sta[++top]=x;
for(int i=head[x];i;i=e[i].next){
if(!dfn[e[i].to]){
dfs(e[i].to);
low[x]=min(low[x],low[e[i].to]);
}
else if(!ltk[e[i].to]){
low[x]=min(low[x],dfn[e[i].to]);
}
}
if(low[x]==dfn[x]){
ltk[0]++;
while(top){
ltk[sta[top]]=ltk[0];
size[ltk[0]]++;
if(sta[top--]==x)break;
}
}
}
void Rebuild(){
for(int i=1;i<=n;i++){
for(int j=head[i];j;j=e[j].next){
if(ltk[i]!=ltk[e[j].to]){
add2(ltk[i],ltk[e[j].to]);
}
}
}
}
void Topsort(){
memset(sta,0,sizeof sta);
int h=0;top=0;
for(int i=1;i<=ltk[0];i++){
if(!du[i]){
len[i]=size[i];
sum[i]=1;
sta[++top]=i;
}
}
while(h<top){
int x=sta[++h];
for(int i=head2[x];i;i=e2[i].next){
if(!--du[e2[i].to]){
sta[++top]=e2[i].to;
}
if(last[e2[i].to]==x)continue; //处理重边
last[e2[i].to]=x;
if(len[e2[i].to]==len[x]+size[e2[i].to]){
sum[e2[i].to]=(sum[e2[i].to]+sum[x])%p;
}
else if(len[e2[i].to]<len[x]+size[e2[i].to]){
len[e2[i].to]=len[x]+size[e2[i].to];
sum[e2[i].to]=sum[x];
}
}
}
}
int main(){
n=read();m=read();p=read();
for(int i=1,u,v;i<=m;i++){
u=read();v=read();
add(u,v);
}
for(int i=1;i<=n;i++){
if(!dfn[i])dfs(i);
}
Rebuild();Topsort();
for(int i=1;i<=ltk[0];i++){
ans1=max(ans1,len[i]);
}
for(int i=1;i<=ltk[0];i++){
if(len[i]==ans1){
ans2=(ans2+sum[i])%p;
}
}
printf("%d\n%d\n",ans1,ans2);
return 0;
}