bzoj 1064: [Noi2008]假面舞会

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。

建图。a可以看见b就a到b连边权为1的，b到a连边权为-1的。然后找环，环大小的gcd就是maxans,maxans大于等于3的最小约数就是minans。没环minans=3,maxans=所有点集相加。

http://hi.baidu.com/lydrainbowcat/item/f3fdc53164770bd76d15e980

这里很清楚

#include<cstdio>
#include<string>
#include<cstring>
#include<algorithm>
using namespace std;
struct edge
{
     int s,t;
     int x;
     int next;
}a[2000001];
int head[100001];
int edge;
inline void add(int s,int t,int x)
{
	 a[edge].next=head[s];
     head[s]=edge;
     a[edge].s=s;
     a[edge].t=t;
     a[edge].x=x;
}
inline int abs(int x)
{
     if(x<0)
          x=-x;
     return x;
}
inline int gcd(int x,int y)
{
	 if(x==0)
	      return y;
	 if(y==0)
	      return x;
     int m=x%y;
     while(m!=0)
     {
          x=y;
          y=m;
          m=x%y;
     }
     return y;
}
bool v[100001],vis[2000001];
int nu[100001];
int ans;
inline void dfs(int d)
{
     v[d]=true;
     int i;
     for(i=head[d];i!=0;i=a[i].next)
     {
     	  int t=a[i].t;
          if(!v[t])
          {
               nu[t]=nu[d]+a[i].x;
               dfs(t);
          }
          else
               ans=gcd(ans,abs(nu[d]+a[i].x-nu[t]));
     }
}
int deepx,depx;
inline void trdp(int d)
{
	 v[d]=true;
     deepx=max(deepx,nu[d]);
     depx=min(depx,nu[d]);
     int i;
     for(i=head[d];i!=0;i=a[i].next)
     {
          int t=a[i].t;
          if(!vis[i])
          {
               vis[i]=true;
               if(i%2==1)
			        vis[i+1]=true;
			   else
			        vis[i-1]=true;
               nu[t]=nu[d]+a[i].x;
               trdp(t);
          }
     }
}
int main()
{
     int n,m;
     scanf("%d%d",&n,&m);
     int i;
     int s,t;
     for(i=1;i<=m;i++)
     {
          scanf("%d%d",&s,&t);
          edge++;
          add(s,t,1);
          edge++;
          add(t,s,-1);
     }
     memset(v,false,sizeof(v));
     for(i=1;i<=n;i++)
          if(!v[i])
               dfs(i);
     int exans;
     if(ans!=0)
     {
     	  if(ans>=3)
     	  {
               for(i=3;i<=ans;i++)
               {
                    if(ans%i==0)
                    {
                         exans=i;
                         break;
                    }
               }
          }
     }
     else
     {
          memset(v,false,sizeof(v));
          memset(vis,false,sizeof(vis));
          for(i=1;i<=n;i++)
          {
          	   if(!v[i])
          	   {
                    deepx=0;
                    depx=0;
                    nu[i]=0;
                    trdp(i);
                    ans+=deepx-depx+1;
               }
          }
          exans=3;
     } 
	 if(ans<3)
     {
          ans=-1;
          exans=-1;
     }
     printf("%d %d\n",ans,exans);
     return 0;
}