zoukankan      html  css  js  c++  java
  • luogu P2272 [ZJOI2007]最大半连通子图

    题目描述

    一个有向图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的余数。

    输入格式

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

    输出格式

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

    说明/提示

    对于100%的数据,N<=100000, M<=1000000, X<=1e8

    #include<queue>
    #include<cstdio>
    #include<cstring>
    #include<iostream>
    #include<algorithm>
    using namespace std;
    int mod;
    const int N=2e5+10,M=2e6+10;
    int next[M],head[N],go[M],tot;
    inline void add(int u,int v){
    	next[++tot]=head[u];head[u]=tot;go[tot]=v;
    }
    int dfn[N],low[N],co[N],st[N],sz[N],top,num,col;
    inline void Tarjan(int u){
    	dfn[u]=low[u]=++num;
    	st[++top]=u;
    	for(int i=head[u];i;i=next[i]){
    		int v=go[i];
    		if(!dfn[v]){
    			Tarjan(v);
    			low[u]=min(low[u],low[v]);
    		}else if(!co[v]){
    			low[u]=min(low[u],dfn[v]);
    		}
    	}
    	if(dfn[u]==low[u]){
    		co[u]=++col;
    		sz[col]=1;
    		while(st[top]!=u){
    			co[st[top]]=col;
    			sz[col]++;
    			--top;
    		}
    		--top;
    	}
    }
    
    struct E{
    	int u,v;
    }e[M];
    int in[N];
    int n,m;
    struct node{
    	int u,dat;
    };
    int dis[N],dp[N];
    inline void topsort(){
    	queue<int>q;
    	for(int i=1;i<=col;i++){
    		if(!in[i])q.push(i);
    		dis[i]=sz[i],dp[i]=1;
    	}
    	while(q.size()){
    		int u=q.front();q.pop();
    		for(register int i=head[u];i;i=next[i]){
    			int v=go[i];
    			if(dis[v]<dis[u]+sz[v]){
    				dis[v]=dis[u]+sz[v];
    				dp[v]=dp[u];
    			}else if(dis[v]==dis[u]+sz[v]){
    				dp[v]=(dp[v]+dp[u])%mod;
    			}
    			--in[v];
    			if(in[v]==0)q.push(v);
    		}
    	}
    	int ans1=0,ans2=0;
    	for(int i=1;i<=col;i++){
    		if(dis[i]==ans1)
    		ans2=(ans2+dp[i])%mod;
    		if(dis[i]>ans1)
    		ans1=dis[i],ans2=dp[i];
    	}
    	cout<<ans1<<endl<<ans2;
    }
    bool cmp(E t1,E t2){
    	if(co[t1.u]==co[t2.u])return co[t1.v]<co[t2.v];
    	else return co[t1.u]<co[t2.u];
    }
    signed main(){
    	
    	cin>>n>>m>>mod;
    	for(int i=1,u,v;i<=m;i++){
    		scanf("%d%d",&u,&v);
    		e[i]=(E){u,v};
    		add(u,v);
    	}
    	for(int i=1;i<=n;i++)
    	if(!dfn[i])Tarjan(i);
    	memset(next,0,sizeof(next));
    	memset(head,0,sizeof(head));
    	memset(go,0,sizeof(go)),tot=0;
    	sort(e+1,e+1+m,cmp);
    	for(int i=1;i<=m;i++){
    		if(co[e[i].u]==co[e[i-1].u]&&co[e[i].v]==co[e[i-1].v])continue;
    		if(co[e[i].u]!=co[e[i].v])
    		add(co[e[i].u],co[e[i].v]),in[co[e[i].v]]++;
    	}
    	topsort();
    }
    
  • 相关阅读:
    nodejs第二天之Buffer类
    2019-2020-1 20175334 20175322 20175315 实验四 《外设驱动程序设计》实验报告
    2019-2020-1 20175315 20175322 20175334 实验三 实时系统
    2019-2020-2-20175315 20175322 20175334-实验二固件程序设计
    2019-2020-1 实验一 开发环境的熟悉
    2018-2019-20175315 实验五《网络编程与安全》实验报告
    实验四《Android程序设计》
    20175315Mycp
    20175315 实验二《Java面向对象程序设计》实验报告
    《Java程序设计》第七周学习总结
  • 原文地址:https://www.cnblogs.com/naruto-mzx/p/11656872.html
Copyright © 2011-2022 走看看