第八周作业-C

题目描述：
大学班级选班长，N 个同学均可以发表意见 若意见为 A B 则表示 A 认为 B 合适，意见具有传递性，即 A 认为 B 合适，B 认为 C 合适，则 A 也认为 C 合适 勤劳的 TT 收集了M条意见，想要知道最高票数，并给出一份候选人名单，即所有得票最多的同学，你能帮帮他吗？
输入：
本题有多组数据。第一行 T 表示数据组数。每组数据开始有两个整数 N 和 M (2 <= n <= 5000, 0 <m <= 30000)，接下来有 M 行包含两个整数 A 和 B(A != B) 表示 A 认为 B 合适。
输出：
对于每组数据，第一行输出 “Case x: ”，x 表示数据的编号，从1开始，紧跟着是最高的票数。 接下来一行输出得票最多的同学的编号，用空格隔开，不忽略行末空格！
sample：
input：
2
4 3
3 2
2 0
2 1

3 3
1 0
2 1
0 2
output：
Case 1: 2
0 1
Case 2: 2
0 1 2
题目分析：
本题非常复杂，这涉及到强连通分量的问题，首先是常见的构件图，不过本题需要dfs3次，所以就需要建立三个图。

struct Edge
{
	int to;
	int next;
}e1[30010],e2[30010],e3[30010];
int head1[5010];
int head2[5010];
int head3[5010];
int tot;
void add3(int u,int v)
{
	e3[++tot].to=v;
	e3[tot].next=head3[u];
	head3[u]=tot;
}
void add(int u,int v)
{
	e1[++tot].to=v;
	e1[tot].next=head1[u];
	head1[u]=tot;
	e2[tot].to=u;
	e2[tot].next=head2[v];
	head2[v]=tot;
}

首先我们要找强连通分量，第一遍dfs要找逆后续，这样从上面我们建立边的时候就能看到顺序是反过来的了。然后找逆后续序列进行遍历，找到所有的SCC。

void dfs1(int s)
{
	vis[s]=true;
	for(int i=head1[s];i!=-1;i=e1[i].next)
	{
		int u=e1[i].to;
		if(vis[u]==false)
		    dfs1(u);
	}
	dfn[++dfnt]=s;
}
void dfs2(int s)
{
	c[s]=scc;
	for(int i=head2[s];i!=-1;i=e2[i].next)
	{
		int u=e2[i].to;
		if(!c[u])
		    dfs2(u);
	}
}
void kosaraju(int n)
{
	dfnt=0;
	scc=0;
	for(int i=0;i<=n;i++)
	{
		vis[i]=false;
		c[i]=0;
	}
	for(int i=0;i<n;i++)
	{
		if(vis[i]==false)
		    dfs1(i);
	}
	for(int i=n;i>=1;i--)
	{
		if(!c[dfn[i]])
		{
			++scc;
			dfs2(dfn[i]);
		}
	}
}

这样我们就完成了缩点，把互相可达的和单项联通的分开考虑。每个强连通分量x中的点的支持数为sum(cnt(j))+cnt(x)-1，其中由j可达x，cnt(j)为联通分量j内的人数，sum为求和。因此我们对强连通分量进行缩点（即每个强连通分量看作一个点），支持人数多的联通分量必然在出度为0的连通分量中。我们对缩点后的图求反图，则出度为0的点变为入度为0点，由该点进行搜索，搜到的点，即为原图中可达该点的点，同时对搜到的点所代表的连通分量的人数求和，可求得该点的支持人数。

		bool vvis[scc+1];
		while(!q.empty())
		{
			int temp=q.front();
			q.pop();
			queue<int>qq;
			qq.push(temp);
			for(int i=1;i<=scc;i++)
				vvis[i]=false;
			vvis[temp]=true; 
			while(!qq.empty())
			{
				int now=qq.front();
				qq.pop();
				cnt[temp]=cnt[temp]+sum[now];
				for(int i=head3[now];i!=-1;i=e3[i].next)
				{
					int ii=e3[i].to;
					if(vvis[ii]==false)
					{
						qq.push(ii);
						vvis[ii]=true;
					}
				}
			}
			Max=max(Max,cnt[temp]);
		}

最终输出最大支持人数的联通分量内的点即可，最终答案的联通分量数可能不止一个。和上一个题一样，既然不能最后有空格，那就先输出一个数字，后面的输出是一个数字+一个空格。

		cout<<"Case "<<x<<": "<<Max-1<<endl;
		auto it=last.begin();
		cout<<*it;
		it++;
		while(it!=last.end())
		{
			cout<<" "<<*it;
			it++;
		}
		cout<<endl;

代码如下：

#include <iostream>
#include <set>
#include <queue>
#include <vector>
#include <cstdio>
using namespace std;
struct Edge
{
	int to;
	int next;
}e1[30010],e2[30010],e3[30010];
int head1[5010];
int head2[5010];
int head3[5010];
int tot;
void add3(int u,int v)
{
	e3[++tot].to=v;
	e3[tot].next=head3[u];
	head3[u]=tot;
}
void add(int u,int v)
{
	e1[++tot].to=v;
	e1[tot].next=head1[u];
	head1[u]=tot;
	e2[tot].to=u;
	e2[tot].next=head2[v];
	head2[v]=tot;
}
int scc,dfnt,c[5010],dfn[5010];
bool vis[5010];
void dfs1(int s)
{
	vis[s]=true;
	for(int i=head1[s];i!=-1;i=e1[i].next)
	{
		int u=e1[i].to;
		if(vis[u]==false)
		    dfs1(u);
	}
	dfn[++dfnt]=s;
}
void dfs2(int s)
{
	c[s]=scc;
	for(int i=head2[s];i!=-1;i=e2[i].next)
	{
		int u=e2[i].to;
		if(!c[u])
		    dfs2(u);
	}
}
void kosaraju(int n)
{
	dfnt=0;
	scc=0;
	for(int i=0;i<=n;i++)
	{
		vis[i]=false;
		c[i]=0;
	}
	for(int i=0;i<n;i++)
	{
		if(vis[i]==false)
		    dfs1(i);
	}
	for(int i=n;i>=1;i--)
	{
		if(!c[dfn[i]])
		{
			++scc;
			dfs2(dfn[i]);
		}
	}
}
int main() 
{
	int T;
	scanf("%d",&T);
	for(int x=1;x<=T;x++)
	{
		int n,m;
		scanf("%d%d",&n,&m);
		for(int i=0;i<=n;i++)
		{
			head1[i]=-1;
			head2[i]=-1;
		}
		tot=0;
		for(int i=m;i>0;i--)
		{
			int a,b;
			scanf("%d%d",&a,&b);
			add(a,b);
		}
		kosaraju(n);
		set<int>ans;
		int sum[scc+1];
		int cnt[scc+1];
		for(int i=0;i<=scc;i++)
		{
			sum[i]=0;
			cnt[i]=0;
			head3[i]=-1;
		}
		set<pair<int,int>> ee;
		for(int k=0;k<n;k++)
		{
			sum[c[k]]++;
			for(int i=head1[k];i!=-1;i=e1[i].next)
			{
				int u=e1[i].to;
				if(c[u]!=c[k])
				{
					ee.insert(make_pair(c[u],c[k]));
				}
			}
		}
		tot=0;
		for(auto it=ee.begin();it!=ee.end();it++)
		{
			add3(it->first,it->second);
			cnt[it->second]++;
		}
		int Max=-1;
		queue<int>q;
		for(int i=1;i<=scc;i++)
		{
			if(cnt[i]==0)
			    q.push(i);	
			cnt[i]=0;		
		}
		bool vvis[scc+1];
		while(!q.empty())
		{
			int temp=q.front();
			q.pop();
			queue<int>qq;
			qq.push(temp);
			for(int i=1;i<=scc;i++)
				vvis[i]=false;
			vvis[temp]=true; 
			while(!qq.empty())
			{
				int now=qq.front();
				qq.pop();
				cnt[temp]=cnt[temp]+sum[now];
				for(int i=head3[now];i!=-1;i=e3[i].next)
				{
					int ii=e3[i].to;
					if(vvis[ii]==false)
					{
						qq.push(ii);
						vvis[ii]=true;
					}
				}
			}
			Max=max(Max,cnt[temp]);
		}
		for(int i=1;i<=scc;i++)
		{
			if(cnt[i]==Max)
			    ans.insert(i);
		}
		vector<int>last;
		for(int i=0;i<n;i++)
		{
			if(ans.find(c[i])!=ans.end())
				last.push_back(i);
		}
		cout<<"Case "<<x<<": "<<Max-1<<endl;
		auto it=last.begin();
		cout<<*it;
		it++;
		while(it!=last.end())
		{
			cout<<" "<<*it;
			it++;
		}
		cout<<endl;
	}
	return 0;
}