BZOJ 1064: [Noi2008]假面舞会 图论， BFS, 脑洞

Description

一年一度的假面舞会又开始了，栋栋也兴致勃勃的参加了今年的舞会。今年的面具都是主办方特别定制的。每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具。每个面具都有一个编号，主办方会把此编号告诉拿该面具的人。为了使舞会更有神秘感，主办方把面具分为k (k≥3)类，并使用特殊的技术将每个面具的编号标在了面具上，只有戴第i 类面具的人才能看到戴第i+1 类面具的人的编号，戴第k 类面具的人能看到戴第1 类面具的人的编号。 参加舞会的人并不知道有多少类面具，但是栋栋对此却特别好奇，他想自己算出有多少类面具，于是他开始在人群中收集信息。 栋栋收集的信息都是戴第几号面具的人看到了第几号面具的编号。如戴第2号面具的人看到了第5 号面具的编号。栋栋自己也会看到一些编号，他也会根据自己的面具编号把信息补充进去。由于并不是每个人都能记住自己所看到的全部编号，因此，栋栋收集的信 息不能保证其完整性。现在请你计算，按照栋栋目前得到的信息，至多和至少有多少类面具。由于主办方已经声明了k≥3，所以你必须将这条信息也考虑进去。
Input

第一行包含两个整数n, m，用一个空格分隔，n 表示主办方总共准备了多少个面具，m 表示栋栋收集了多少条信息。接下来m 行，每行为两个用空格分开的整数a, b，表示戴第a 号面具的人看到了第b 号面具的编号。相同的数对a, b 在输入文件中可能出现多次。
Output

包含两个数，第一个数为最大可能的面具类数，第二个数为最小可能的面具类数。如果无法将所有的面具分为至少3 类，使得这些信息都满足，则认为栋栋收集的信息有错误，输出两个-1。
Sample Input
【输入样例一】

6 5

1 2

2 3

3 4

4 1

3 5

【输入样例二】

3 3

1 2

2 1

2 3
Sample Output
【输出样例一】

4 4

【输出样例二】

-1 -1
HINT

100%的数据，满足n ≤ 100000, m ≤ 1000000。

解题思路： 下面分析过程来自 qzqzgfy ！

分联通块考虑，如果每个联通块都是没有有向环的话，那么各个联通块中，最长链就是最大答案，3就是最小答案。

只要有一个联通块有环，那么答案一定是这个环长度的因数，最大答案，就是这些环长度的gcd

不过，要是有这个非正常的环怎么办?

我们可以看到,4->3和2->3都指向了3，这怎么办？那么我们只要在一开始建图的时候，原来的有向边权值为1，再同时建一个反向边权值为-1，把有向图变成无向图。

为什么?，因为如图，4可以到3,2也可以到3，说明2的编号和4相同，所以2->3->4的路径实际上是”走出去一步，又走回来一步”，也就是我常说的”有来有回”,按照我们刚才的建图方式，这个环的长度就是:1+1+1+1-1=3，事实上，答案也是如此，3和5编号相同，2和4编号相同，这样图中实际上是只有3种面具。

复杂度: O(m) , 通过时间 : 180ms

代码如下：

#include <bits/stdc++.h>
using namespace std;
const int maxn = 110000;
const int maxm = 2200000;
int n, m, cnt, head[maxn];
bool vis[maxn];
int dist[maxn], comgcd; //每个环上正向边数-反向边数的最大公约数，也就是面具种类数最大值
struct edge{
    int v, w, nxt;
    edge(){}
    edge(int v, int w, int nxt) : v(v), w(w), nxt(nxt) {}
}E[maxm];
void init(){
    memset(head, -1, sizeof(head)); cnt = 0;
}
void addedge(int u, int v, int w){
    E[cnt].v = v, E[cnt].w = w, E[cnt].nxt = head[u], head[u] = cnt++;
}
int bfs(int s){
    int maxdis = 0, mindis = 0; //maxdist是这个联通块中S出发的最长链距离,mindist反之
    queue <int> q;
    q.push(s);
    dist[s] = 0;
    vis[s] = 1;
    q.push(s);
    while(!q.empty()){
        int u = q.front(); q.pop();
        for(int i = head[u]; ~i; i = E[i].nxt){
            int v = E[i].v;
            if(vis[v]){//找到环了
                comgcd = __gcd(comgcd, dist[u] + E[i].w - dist[v]); //找到了一个环，那么这个环上的点的个数是最终答案
                continue;
            }
            vis[v] = true;
            dist[v] = dist[u] + E[i].w;
            maxdis = max(maxdis, dist[v]);
            mindis = min(mindis, dist[v]);
            q.push(v);
        }
    }
    return maxdis - mindis + 1;
}

int main(){
    int sum = 0;
    init();
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= m; i++){
        int u, v;
        scanf("%d%d", &u, &v);
        addedge(u, v, 1);
        addedge(v, u, -1);//反向边边权为-1
    }
    for(int i = 1; i <= n; i++){
        if(!vis[i]){
            sum += bfs(i);
        }
    }
    comgcd = abs(comgcd);//注意最大公约数要取绝对值
    if(comgcd){ //有环，就不需要考虑树
        if(comgcd < 3){//面具种类数最大不到3，那么输入数据有问题
            printf("-1 -1\n");
            return 0;
        }
        int j;
        for(j = 3; j < comgcd && comgcd % j; j++);
        printf("%d %d\n", comgcd, j);
        return 0;
    }
    else if(sum < 3){
        printf("-1 -1\n");
        return 0;
    }
    printf("%d 3\n", sum);
    return 0;
}