week8 C - 班长竞选（强联通分量，缩点）

一、题目描述

大学班级选班长，N 个同学均可以发表意见 若意见为 A B 则表示 A 认为 B 合适，意见具有传递性，即 A 认为 B 合适，B 认为 C 合适，则 A 也认为 C 合适 勤劳的 TT 收集了M条意见，想要知道最高票数，并给出一份候选人名单，即所有得票最多的同学，你能帮帮他吗？

Input

本题有多组数据。第一行 T 表示数据组数。每组数据开始有两个整数 N 和 M (2 <= n <= 5000, 0 <m <= 30000)，接下来有 M 行包含两个整数 A 和 B(A != B) 表示 A 认为 B 合适。

Output

对于每组数据，第一行输出 “Case x: ”，x 表示数据的编号，从1开始，紧跟着是最高的票数。 接下来一行输出得票最多的同学的编号，用空格隔开，不忽略行末空格！

Sample Input

2
4 3
3 2
2 0
2 1
3 3
1 0
2 1
0 2

Sample Output

Case 1: 2
0 1
Case 2: 2
0 1 2

二、思路概述

• 此题初始状态时，需要定义两个图，原图根据a–>b(a给b投了票），图中a点到b点就有一条边，反图就是原图中所有的边反向。
• 第一遍dfs遍历原图，得到逆后序序列。
• 第二遍 dfs 在反图中按照逆后序序列进行遍历，对联通分量进行标记，并记录每一个联通分量中的点数。
• 根据反图进行缩点，生成一个新的图，此图里的每个点都是一个联通图，点数代表联通图的标记。
• 对新图中入度为0的点进行dfs，计算能够到达的点数。属于第 i 个 SCC 的点来说，答案分为两部分，令 SCC[i] 表示第 i 个 SCC 中点的个数。当前 SCC 中的点，ans += SCC[i] – 1（去除自己）其它 SCC 中的点SUM ( SCC[j] )，其中 j 可到达 I找出最大值。
• 确定哪几个联通分支中的点的ans最大，就输出这些点。

三、细节

1. 没有初始化：此题中用到的数组、变量非常多，所以初始化一定要检查是否正确，完整。
2. 此题中有三个图，所以他们使用起来不要弄混了，也不要笔误。

四、完整代码

#include<iostream>
#include<algorithm>
#include<set>
#include<map>
#include<cstring>
using namespace std;
const int maxn=5010;//点的数量最多5000个
const int maxm=30020; 

int f[maxn],ivf[maxn],vis1[maxn],vis2[maxn],scc[maxn],scc2[maxn],degree[maxn];//f用来记录后序序列,ivf用来记录逆后序序列
int n,fcnt,acnt1,acnt2,acnt3,scnt,s;//共有n个点
int head1[maxn],head2[maxn],head3[maxn],ans[maxn],ano[maxn];
 
struct edge{
	int to,next;
}e1[maxm],e2[maxm],e3[maxn];//e1为原图，e2为反向图 


void init()
{
	fcnt=-1;
	acnt1=0;
	acnt2=0;
	acnt3=0;
	scnt=0;
	
	for(int i=0;i<n+3;i++)
	{
		ans[i]=0;
		ano[i]=0;
		f[i]=0;
		ivf[i]=0;
		vis1[i]=0;
		scc[i]=0;
		scc2[i]=0;
		degree[i]=0; 
		head1[i]=-1;
		head2[i]=-1;
		head3[i]=-1;
	}
} 

void addedge1(int from,int to)
{
	e1[acnt1].to =to;
	e1[acnt1].next =head1[from];
	head1[from] = acnt1;
    acnt1++;
}

void addedge2(int from,int to)
{
	e2[acnt2].to =to;
	e2[acnt2].next =head2[from];
	head2[from]=acnt2;
    acnt2++;
}

void addedge3(int from,int to)
{
	e3[acnt3].to =to;
	e3[acnt3].next =head3[from];
	head3[from]=acnt3;
    acnt3++;
      
    degree[to]++;//记录入度
}

void dfs1(int s)
{
	vis1[s]=1;
	for(int i=head1[s];i!=-1;i=e1[i].next)
	{
		int to=e1[i].to;
		if(!vis1[to])
		dfs1(to);
	}
	f[++fcnt]=s;	//f[i]的值从0开始，到n-1 
}
 
void dfs2(int s,int scnt)
{//s可到达的点均用scnt标记 
	scc[s]=scnt;
	scc2[scnt]++;
	for(int i=head2[s];i!=-1;i=e2[i].next)
	{
		int to=e2[i].to;
		if(!scc[to])
		{
		    dfs2(to,scnt);
		}		
	}	
}

void Kosaraju()//找到有向图中所有的 SCC
{
	//第一遍 dfs 确定原图的逆后序序列，
	dfs1(0);
	for(int i=1;i<n;i++)
	{
		if(!vis1[i])dfs1(i);
	}

	//第二遍 dfs 在反图中按照逆后序序列进行遍历
	for(int i=0;i<=fcnt;i++)
	{
		
		int x=f[fcnt-i];
		if(scc[x]==0)
		{scnt++;
		dfs2(x,scnt);	
		}
	}	
	
} 

void dfs3(int x)
{
	vis2[x]=1;
	//找出每个入度为0的点，进行dfs，计算SUM(SCC[j])
	for(int i=head3[x];i!=-1;i=e3[i].next)
	{	
		int to=e3[i].to;
		if(vis2[to]==0)
		{		
		    ans[s]=ans[s]+scc2[to];
		    dfs3(to);
		} 		 
	}
}

void answer()//最后答案的函数 
{//共n个点 
    
    //缩点，反向图
    
    int tol=0;
	for(int i=0;i<n;i++)
	{
		for(int j=head2[i];j!=-1;j=e2[j].next)
		{
			int y=e2[j].to; 
			if(scc[i]==scc[y])continue; 
			addedge3(scc[i],scc[y]);	
		}
    }  
    
    
	//此图里的每个点都是一个联通图，点数代表联通图的标记
	int maxs=0; 

	for(int i=1;i<=scnt;i++)
	{		
		if(degree[i]==0)
		{
			memset(vis2,0,sizeof(vis2));			
			ans[i]+=scc2[i]-1;
			
		    s=i;
			dfs3(i);
			if(maxs<ans[i])maxs=ans[i];			
		}		 
	}
	
	
	for(int i=1;i<=scnt;i++)
	{
		if(ans[i]==maxs) 
		{
			ano[i]=1;
		}
	}
	
	printf("%d\n",maxs);
	
	int anscnt=0;
	for(int i=0;i<n;i++)
	{			
	    if(ano[scc[i]]==1)
		{
			if(anscnt==0)
			{
				printf("%d",i);anscnt++;	
			}
			else printf(" %d",i);
		} 
	}
	
}
int main()
{
	int T;
	cin>>T;//T表示数据组数
	
	for(int i=1;i<=T;i++)
	{		
		int m;//n表示n个同学 ,编号为0~n-1 
		scanf("%d%d",&n,&m);//(2 <= n <= 5000, 0 <m <= 30000)
		init();
		while(m--)//接下来有 M 行包含两个整数 A 和 B(A != B) 表示 A 认为 B 合适
		{
			int a,b;
			scanf("%d%d",&a,&b);
			addedge1(a,b);//图， 
			addedge2(b,a);//反向图	
		}
		//cout<<acnt1<<" "<<acnt2<<endl;
		cout<<"Case "<<i<<": ";
		Kosaraju();
		answer();
		cout<<endl;
	}
 }