BZOJ1064：假面舞会（DFS & 思维）

1064: [Noi2008]假面舞会

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 2264   Solved: 1094
[ Submit][ Status][ Discuss]

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。

题意：略。

思路：比较烧脑的一道图论，参考了别人的思路。要处理的有三种结构：单方向的环，混杂方向的环，树。对于单方向的环，显然K最大是取这些环长度的GCD；对于混杂方向的环，可以看作前进和后退；若只有树结构，K最大就是取每个连通分量的最长链的和。那么我们每条边都建一条反向边权为-1的边，dfs时就可以顺带处理第二种和第三种结构了。那么dfs可能会没有搜到一些环，比如A环和B环有公共边C，DFS可能没法比较到三个环，事实上比较其中两个就够了，因为有性质GCD(X,Y) = GCD(X,Y+X)。

# include <iostream>
# include <cstdio>
# include <cstring>
# include <vector>
# include <algorithm>
# define mp make_pair
# define pb push_back
# define INF 0x3f3f3f3f
using namespace std;
const int maxn = 1e5+30;
vector<pair<int,int> >g[maxn*20];
int vis[maxn], val[maxn], fa[maxn], imax[maxn], imin[maxn];
int anc, ans=0;
int gcd(int x, int y)
{
    return y==0?x:gcd(y, x%y);
}
int find(int x)
{
    if(x != fa[x])
        fa[x] = find(fa[x]);
    return fa[x];
}
void dfs(int cur)
{
    imax[anc] = max(imax[anc], val[cur]);
    imin[anc] = min(imin[anc], val[cur]);
    vis[cur] = 1;
    for(int i=0; i<g[cur].size(); ++i)
    {
        int to = g[cur][i].first;
        int num = g[cur][i].second;
        if(!vis[to])
        {
            val[to] = val[cur] + num;
            dfs(to);
        }
        else
            ans = gcd(ans, abs(val[cur]+num-val[to]));
    }
}
int main()
{
    int n, m, u, v;
    scanf("%d%d",&n,&m);
    for(int i=1; i<=n; ++i) fa[i]=i, imax[i] = -INF, imin[i] = INF;
    while(m--)
    {
        scanf("%d%d",&u,&v);
        g[u].pb(mp(v, 1));
        g[v].pb(mp(u, -1));
        fa[find(u)] = fa[find(v)];
    }
    for(int i=1; i<=n; ++i)
    {
        if(!vis[i])
        {
            anc = find(i);
            dfs(anc);
        }
    }
    if(ans == 0)
    {
        int sum = 0;
        for(int i=1; i<=n; ++i)
            if(i == fa[i])
                sum += imax[i]-imin[i]+1;
        if(sum < 3) printf("-1 -1");
        else printf("%d %d\n",sum, 3);
    }
    else
    {
        if(ans < 3) printf("-1 -1\n");
        else
        {
            int ans2;
            for(int i=3; i<=ans; ++i) if(ans%i==0){ans2=i;break;}
            printf("%d %d\n",ans,ans2);
        }
    }
    return 0;
}