C班长竞选--强连通分量

kosaraju算法

这是一个求连通分量的算法，利用的是DFS序；
主要是两遍DFS，第一遍求DFS后序序列；第二遍利用DFS的逆后序标记属于同一个强连通分量的节点。
代码如下：

vector<int> G1[5005],G2[5005];//原图，反图，

void dfs1(int x){
	vis[x]=1;
	for(int i=0; i<G1[x].size(); i++)
		if(!vis[G1[x][i]])
			dfs1(G1[x][i]);
	dfn[dcnt]=x;dcnt++;
}

void dfs2(int x)//分成SCC
{
	SCC[scnt]++;
	c[x]=scnt;
	for(int i=0; i<G2[x].size(); i++)
		if(!c[G2[x][i]])
			dfs2(G2[x][i]);
}

void kosaraju(){
	for(int i=0;i<n;i++)
		if(!vis[i]) dfs1(i);
	for(int i=n-1;i>=0;i--)//逆后序
		if(!c[dfn[i]]){
			++scnt,dfs2(dfn[i]);
		}
}

题目

大学班级选班长，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开始，紧跟着是最高的票数。 接下来一行输出得票最多的同学的编号，用空格隔开，不忽略行末空格！

分析

1、首先利用 kosaraju算法求出全连通分量（SCC）
2、利用SCC缩点：（得到图G3）
不属于同一个SCC的边就加入缩点之后的图当中，因为两个SCC之间的边可能有非常多条，所以为了防止重复，用visit数组记录哪些是重复有的边。
3、利用缩点之后的图G3，投票数最多的一定在入度为0的SCC中，再用DFS算得路过的SCC[i]个数;在加上自己的SCC[ ]-1(自己除外）

为啥我WA了14次
1、看错题目，开始以为输出的是投票最多的人数，后来才发现是投票最多的票数…
2、在看错题目之前，我一直修改的是初始化，因为我在初始化里吃了不少亏；然后我发现，初始化可以放在一个函数ini（）里面，更好管理。
3、不知道为什么，我用C++WA了，然后我改成G++AC了（同一份代码~~）

代码

#include<stdio.h>
#include<vector>
#include<algorithm>
#include<string.h>
using namespace std;
int n,tsum;
bool vis[5005],visit[5005][5005];
int SCC[5005],c[5005],dfn[5005],in_deg[5005],dcnt,scnt;
vector<int> G1[5005],G2[5005],G3[5005];//原图，反图，缩点之后的tu
void ini(){
	for(int j=0;j<5005;j++){
		G1[j].clear();G2[j].clear();G3[j].clear();
	}
	for(int i=0;i<5005;i++)
		for(int j=0;j<5005;j++)
			visit[i][j]=0;
	dcnt=0;scnt=0;tsum=0;
	memset(c,0,sizeof(c));
	memset(dfn,0,sizeof(dfn));
	memset(vis,0,sizeof(vis));
	memset(SCC,0,sizeof(SCC));
	memset(in_deg,0,sizeof(in_deg));
}
void dfs1(int x){
	vis[x]=1;
	for(int i=0; i<G1[x].size(); i++)
		if(!vis[G1[x][i]])
			dfs1(G1[x][i]);
	dfn[dcnt]=x;dcnt++;
}

void dfs2(int x)//分成SCC
{
	SCC[scnt]++;//printf("scnt:%d SCC:%d\n",scnt,SCC[scnt]);
	c[x]=scnt;
	for(int i=0; i<G2[x].size(); i++)
		if(!c[G2[x][i]])
			dfs2(G2[x][i]);
}

void kosaraju(){
//	printf("n:%d\n",n);
	for(int i=0;i<n;i++)
		if(!vis[i]) dfs1(i);
/*	for(int i=0;i<n;i++){
		printf("dfn:%d ",dfn[i]);
	}*/
	for(int i=n-1;i>=0;i--)//逆后序
		if(!c[dfn[i]]){
			++scnt,dfs2(dfn[i]);
		}

/*	for(int i=0;i<n;i++)
		printf("%d ",c[i]);*/
}

void suodian(){//缩点
	for(int i=0;i<n;i++){
		for(auto y:G2[i]){
		//	printf("i:%d j:%d x:%d y:%d\n",i,y,c[i],c[y]);
			if(c[y]!=c[i]&&visit[c[y]][c[i]]==0){//SCC间的连接
				visit[c[y]][c[i]]=1;//防止重复
				G3[c[i]].push_back(c[y]);
				in_deg[c[y]]++;
			}
		}
	}
}
void dfs3(int x){//in_deg=0的点,点从1开始
	vis[x]=1;
	for(auto y:G3[x])
		if(!vis[y]){
			tsum=tsum+SCC[y];
			dfs3(y);
		}
}
struct tmax{//记录最大的投票数，以及SCC的序号
	int thesum,xuhao;
	bool operator < (const tmax &tm)const{
		if(thesum!=tm.thesum) return thesum>tm.thesum;
	}
}point[50005];
int main()
{
	int T,M,A,B;
	scanf("%d",&T);
	for(int cishu=1;cishu<=T;cishu++)
	{
		scanf("%d %d",&n,&M);
		ini();
		for(int j=0;j<M;j++){
			scanf("%d %d",&A,&B);
			G1[A].push_back(B);
			G2[B].push_back(A);
		}
		kosaraju();
		suodian();
	/*	for(int j=1;j<=scnt;j++)
			printf("%d ",SCC[j]);*/
		//反图里面入度为0的
		int cnt=0;
		for(int m=1;m<=scnt;m++){
			if(in_deg[m]==0){
				//printf("M:%d\n",m);
				tsum=SCC[m]-1;
				memset(vis,0,sizeof(vis));
				dfs3(m);
				point[cnt].thesum=tsum;point[cnt].xuhao=m;
				cnt++;
			}
		}
		sort(point,point+cnt);
	/*	printf("count: %d\n",cnt);
		for(int t=0;t<cnt;t++){
			printf("thesum: %d xuhao:%d\n",point[t].thesum,point[t].xuhao);
		}*/
		int themax=point[0].thesum;//最多的
		int ans[5005],coun=0;
		int t=0;
		while(point[t].thesum==themax){
			for(int j=0;j<n;j++){
				if(c[j]==point[t].xuhao){
					ans[coun]=j;coun++;
				}
			}t++;
		}
		sort(ans,ans+coun);
		printf("Case %d: %d\n",cishu,themax);
		int j=0;
		while(j<coun-1){
			printf("%d ",ans[j]);j++;
		}
		printf("%d\n",ans[coun-1]);
	}
	return 0;
}