bzoj 1064 noi2008 假面舞会题解

莫名其妙的变成了我们的noip互测题...

其实这题思想还是比较简单的，只是分类不好分而已

其实就是一个dfs的事

首先，非常明显，原题目中的所有关系可以抽象成一个图（这是...显而易见的吧...）

接下来，我们仅需在图上讨论即可

当然，这个图有几个部分组成其实并没有那么重要，毕竟，这些部分基本是互不干扰的。

所以接下来我们只需要对每一个块分别处理即可

我们来分类：

首先，如果所有块都是树，我们只需求出每个树上的最长链即可

接下来，如果存在环（包括真实的环和类环，即1-2-3-1和1-2-4+1-3-4两种），那么种类数最多显然是所有

环大小的gcd（至于其他的树，可以完全不必考虑了）

于是问题就变成了怎么求环的大小

请大家注意一点，就是我们所说的环的大小是指的一个环中至多可以有几种面具

也就是说，对于一个这样的环：1-2-4和1-3-4，很显然2和3的编号应该是一样的，这样我们说这个类环的大小是3！

接下来我们讨论一下怎么求

其实求法很简单：化有向图为带权无向图！

即：如果这条边是正向的，我们把他的边权设为+1，反之设为-1

这样做的目的在于，还是以上面的图为例：1-2-4和1-3-4，我们双向建边就能求出环的大小（自己画一下，一下就出来）

剩下的部分就聊尽人事了

#include <cstdio>
#include <cmath>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <queue>
#include <stack>
using namespace std;
struct Edge
{
    int next;
    int to;
    int val;
}edge[];
int head[];
int dep[];
int cnt=;
int d;
int n,m;
int maxdep,mindep=0x3f3f3f3f;
bool used[];
void init()
{
    memset(head,-,sizeof(head));
    cnt=;
}
int gcd(int x,int y)
{
    if(y==)
    {
        return x;
    }
    return gcd(y,x%y);
}
void add(int l,int r,int w)
{
    edge[cnt].next=head[l];
    edge[cnt].to=r;
    edge[cnt].val=w;
    head[l]=cnt++;
}
void dfs(int x,int deep)
{
    used[x]=;
    dep[x]=deep;
    maxdep=max(maxdep,dep[x]);
    mindep=min(mindep,dep[x]);
    for(int i=head[x];i!=-;i=edge[i].next)
    {
        int to=edge[i].to;
        if(!used[to])
        {
            dfs(to,deep+edge[i].val);
        }else
        {
            d=gcd(d,abs(deep+edge[i].val-dep[to]));
        }
    }
}
inline int read()
{
    int f=,x=;char ch=getchar();
    while(ch<''||ch>''){if(ch=='-')f=-;ch=getchar();}
    while(ch>=''&&ch<=''){x=x*+ch-'';ch=getchar();}
    return x*f;
}
int main()
{
//    freopen("party.in","r",stdin);
//    freopen("party.out","w",stdout);
    n=read(),m=read();
    init();
    for(int i=;i<=m;i++)
    {
        int x=read(),y=read();
        add(x,y,);
        add(y,x,-);
    }
    int ret=;
    for(int i=;i<=n;i++)
    {
        if(!used[i])
        {
            maxdep=;
            mindep=0x3f3f3f3f;
            dfs(i,);
            ret+=maxdep-mindep+;
        }
    }
    if(!d)
    {
        if(ret>=)
        {
            printf("%d 3\n",ret);
            return ;
        }else
        {
            printf("-1 -1\n");
            return ;
        }
    }else
    {
        if(d<)
        {
            printf("-1 -1\n");
            return ;
        }
        for(int i=;i<=d;i++)
        {
            if(d%i==)
            {
                printf("%d %d\n",d,i);
                return ;
            }
        }
    }
}