清华OJ旅行商(TSP)

描述

Shrek是一个大山里的邮递员，每天负责给所在地区的n个村庄派发信件。但杯具的是，由于道路狭窄，年久失修，村庄间的道路都只能单向通过，甚至有些村庄无法从任意一个村庄到达。这样我们只能希望尽可能多的村庄可以收到投递的信件。

Shrek希望知道如何选定一个村庄A作为起点（我们将他空投到该村庄），依次经过尽可能多的村庄，路途中的每个村庄都经过仅一次，最终到达终点村庄B，完成整个送信过程。这个任务交给你来完成。

输入

第一行包括两个整数n，m，分别表示村庄的个数以及可以通行的道路的数目。

以下共m行，每行用两个整数v1和v2表示一条道路，两个整数分别为道路连接的村庄号，道路的方向为从v1至v2，n个村庄编号为[1, n]。

输出

输出一个数字，表示符合条件的最长道路经过的村庄数。

样例

见英文题面

限制

1 ≤ n ≤ 1,000,000

0 ≤ m ≤ 1,000,000

输入保证道路之间没有形成环

时间：2 sec

空间：256 MB

Example

Input

4 3
1 4
2 4
4 3

Output

3

思路：

        因为只能单向通行，实际上是求图中最长的一条路径，包含的结点个数最多。最长路径有个特点就是，起点入度为0，终点出度为0，因为如果不是的话，我们向外扩展一个结点，都会得到更长的一条路径。故可以按照拓扑排序的思想，用队列（或者其他任何一个集合）维护当前图中入度为0的结点，每次取出一个结点，根据当前结点更新所有邻接点的路径长度，而后将当前结点删除。删除后，更新其每一个邻接点的入度，如果邻接点入度变成0，邻接点入队，直到队列为空，得到最长的路径所包含的结点数量。

代码如下：

#include<cstring>
#include<stdio.h>
#include<iostream>
using namespace std;
const int N = 1e6 + 5;
int h[N],e[N],ne[N],idx;    //邻接表方式存储图
int q[N], n, m;             //数组模拟队列
int d[N], maxVilleges = 1;  //d[i]表示顶点i的入度, maxVilleges为最后的答案
int villeges[N];            //村子个数

void add(int a, int b){     //添加一条a到b的边
    e[idx] = b, ne[idx] = h[a], h[a] = idx ++;
}

void topu(){
    int hh = 0, tt = -1;

    /*0入度边入队*/
    for(int i = 1; i <= n; i++){
        if(!d[i]) q[++tt] = i; 
    }

    /*初始化每个顶点目前的村庄个数*/
    for(int i = 1; i <= n; i++){
        villeges[i] = 1;
    }

    //遍历入度为0的点，直到队列为空
    while(hh <= tt){
        int t = q[hh++];

        //遍历顶点t 的所有邻接点j
        for(int i = h[t]; i != - 1; i = ne[i]){
            int j = e[i];
            if(villeges[t] + 1 > villeges[j]){
                villeges[j] = villeges[t] + 1;
            }
            maxVilleges = max(maxVilleges, villeges[j]);
            d[j] --;
            if(!d[j]) q[++tt] = j;
        }
    }
}

int main(){
    scanf("%d%d", &n, &m);
    memset(h, -1, sizeof(h));

    while(m --){
        int a, b;
        scanf("%d%d", &a, &b);
        add(a, b);
        d[b]++;
    }

    topu();

    printf("%d\n", maxVilleges);

    return 0;
}

结果：