强连通分量题。题目神描述,要锻炼题意理解能力的话此题为上等。题意:在有向图G中,给这样一个新的定义:对于在G中任何一个点v可达的点w,w都可达v,那么点v是一个sink。G图的bottom子图是由G图中所有的sink点构成,请按照顺序输出G图对应的bottom子图中的所有点编号,如果没有sink点,那么输出一个空行。Orz,题意描述得我好辛苦。
我的解题思路:首先如果v可达w,那么必须w可达v,v才算是sink点。那么很明显,v可达的所有点都必须可达v,那么v才算是sink点。一个强连通分量中的任意点都是相互可达的,当我们把图缩点成一个DAG之后,只有出度为0的强连通分量里面的点才满足sink点的性质。在一个DAG中至少存在一个出度为0的点,所以题目要求没有sink点时输出空行这个情况是绝对不会发生的。最后,我们按顺序把出度为0的强连通分量里面的点全部按照升序输出就可以了。
我的解题代码:
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cctype>
#include <cmath>
#include <algorithm>
#include <vector>
#include <queue>
using namespace std;
const int N = 5001;
vector<int> e[N];
int dfn[N], low[N];
int stack[N], stop;
bool instack[N];
int belong[N];
int outdeg[N]; //存储强连通分量的出度
int timer, sccn;
int n, m;
void InitRead();
void DataProcess();
void Tarjan(int x);
int main()
{
while (~scanf("%d", &n))
{
if (n == 0) break;
InitRead();
DataProcess();
}
return 0;
}
void InitRead()
{
timer = sccn = stop = 0;
memset(dfn, 0, sizeof(dfn));
memset(instack, false, sizeof(instack));
memset(outdeg, 0, sizeof(outdeg));
scanf("%d", &m);
for (int i=1; i<=n; ++i)
{
e[i].clear();
}
int a, b;
for (int i=0; i<m; ++i)
{
scanf("%d %d", &a, &b);
e[a].push_back(b);
}
return;
}
void DataProcess()
{
for (int i=1; i<=n; ++i)
{
if (!dfn[i]) Tarjan(i);
}
for (int i=1; i<=n; ++i) //计算强连通分量出度
{
int size = e[i].size();
for (int j=0; j<size; ++j)
{
if (belong[i] != belong[e[i][j]])
{
outdeg[belong[i]]++;
}
}
}
for (int i=1; i<=n; ++i) //按升序输出出度为0的强连通分量里的点
{
if (outdeg[belong[i]] == 0)
{
printf("%d", i++);
for (; i<=n; ++i)
{
if (outdeg[belong[i]] == 0)
{
printf(" %d", i);
}
}
break;
}
}
putchar('\n');
return;
}
void Tarjan(int x)
{
dfn[x] = low[x] = ++timer;
stack[stop++] = x;
instack[x] = true;
int size = e[x].size();
int y;
for (int i=0; i<size; ++i)
{
y = e[x][i];
if (!dfn[y])
{
Tarjan(y);
low[x] = min(low[x], low[y]);
}
else if (instack[y])
{
low[x] = min(low[x], dfn[y]);
}
}
if (dfn[x] == low[x])
{
sccn++;
do
{
y = stack[--stop];
instack[y] = false;
belong[y] = sccn;
} while (x != y);
}
return;
}