7-3 病毒溯源 (20 分)

病毒容易发生变异。某种病毒可以通过突变产生若干变异的毒株，而这些变异的病毒又可能被诱发突变产生第二代变异，如此继续不断变化。

现给定一些病毒之间的变异关系，要求你找出其中最长的一条变异链。

在此假设给出的变异都是由突变引起的，不考虑复杂的基因重组变异问题 —— 即每一种病毒都是由唯一的一种病毒突变而来，并且不存在循环变异的情况。

输入格式：

输入在第一行中给出一个正整数 N（≤104），即病毒种类的总数。于是我们将所有病毒从 0 到 N−1 进行编号。

随后 N 行，每行按以下格式描述一种病毒的变异情况：

k 变异株1 …… 变异株k

其中 k 是该病毒产生的变异毒株的种类数，后面跟着每种变异株的编号。第 i 行对应编号为 i 的病毒（0≤i<N）。题目保证病毒源头有且仅有一个。

输出格式：

首先输出从源头开始最长变异链的长度。

在第二行中输出从源头开始最长的一条变异链，编号间以 1 个空格分隔，行首尾不得有多余空格。如果最长链不唯一，则输出最小序列。

注：我们称序列 { a1​,⋯,an​ } 比序列 { b1​,⋯,bn​ } “小”，如果存在 1≤k≤n 满足 ai​=bi​ 对所有 i<k 成立，且 ak​<bk​。

输入样例：

10
3 6 4 8
0
0
0
2 5 9
0
1 7
1 2
0
2 3 1

输出样例：

4
0 4 9 1

代码 ：

#include <iostream>
#include<cmath>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
using namespace std;

bool vis[10001];//标记是否经过
bool judge[10001];//标记源头
int n;//病毒种类

vector<int> ans;//最终答案
vector<int> a;//记录当前寻找的变异链
int main()
{
	cin >> n;
	vector<vector<int>> inf(n);
    //memset(vis,false,sizeof(vis));
    //memset(judge,true,sizeof(judge));
    //这里说明一下，我当初是用memset来填写的，但不知道原因就是报错，编译错误。
    //所以改成了循环填写。
	for (int i = 0; i < n; i++)
	{
		vis[i] = false;//还没有经过
		judge[i] = true;//认为都是源头
	}
	for (int i = 0; i < n; i++)
	{
		int k;
		cin >> k;
		for (int j = 0; j < k; j++)
		{
			int b;
			cin >> b;
			judge[b] = false;//变异出现的病毒都设为false
			inf[i].push_back(b);
		}
	}
	int frist;//变异链源头
	for (int i = 0; i < n; i++)
	{
		if (judge[i] == true)
		{
			frist = i;
			break;
		}
	}
	stack<int> zhan;
	zhan.push(frist);
	while (!zhan.empty())
	{
		int s;
		s = zhan.top();
		if (inf[s].size() == 0)
		{
			a.push_back(s);
			vis[s] = true;
			if (a.size()>ans.size())//当前链长>ans长，交换。
			{
				ans.clear();
				for (int i = 0; i < a.size(); i++)
				{
					ans.push_back(a[i]);
				}
			}
			else if (a.size() == ans.size())//长度相同，选序列小的
			{
				int i = 0;
				while (a[i] == ans[i])
					i++;
				if (a[i] < ans[i])//交换
				{
					ans.clear();
					for (int i = 0; i < a.size(); i++)
					{
						ans.push_back(a[i]);
					}
				}
			}
		}
		//说明没有到变异链的结尾
		int sq = 0;//
		for (; sq < inf[s].size(); sq++)
		{
			if (vis[inf[s][sq]] == false)
			{
				if (vis[s] == false)
				{
					a.push_back(s);
					vis[s] = true;
				}
				zhan.push(inf[s][sq]);
				break;
			}
		}
		if (sq == inf[s].size())//说明当前变异种类中都找到了，回到上一个。
		{
			zhan.pop();
			a.pop_back();
		}
	}
	cout << ans.size() << endl;
	for (int i = 0; i < ans.size()-1; i++)
	{
		cout << ans[i] << " ";
	}
	cout << ans[ans.size() - 1];
	return 0;
}

思路：深度优先搜索来寻找变异链，变量动态数组ans 和a，用来储存所找到的变异链，ans为最终的链，每次找到一条变异链时，就与ans的长度（当长度相同时，就比较序列大小）并进行交换操作，使ans始终储存当前所找到的所有变异链中最长的。当将所有的变异链条找完后，ans中储存的自然是最长，最小序列。

其中的关键是实现深度优先搜索算法。我原本是用递归来做的，结果似乎是递归栈溢出了，导致部分正确，还有代码段错误，后来就尝试用stack来实现了，会要容易想一些，因为递归的次数一旦多起来的话确实感觉到有点难以想象。以后遇到类似的，还是多用stack来实现。至于代码的细节，我就不写了，对于以后复习的自己，建议对照样例的病毒变异图来写。