NOI2008 假面舞会party

 

假面舞会

【问题描述】

一年一度的假面舞会又开始了，栋栋也兴致勃勃的参加了今年的舞会。今年的面具都是主办方特别定制的。每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具。每个面具都有一个编号，主办方会把此编号告诉拿该面具的人。为了使舞会更有神秘感，主办方把面具分为k (k≥3)类，并使用特殊的技术将每个面具的编号标在了面具上，只有戴第i 类面具的人才能看到戴第i+1 类面具的人的编号，戴第k 类面具的人能看到戴第1 类面具的人的编号。 参加舞会的人并不知道有多少类面具，但是栋栋对此却特别好奇，他想自己算出有多少类面具，于是他开始在人群中收集信息。

栋栋收集的信息都是戴第几号面具的人看到了第几号面具的编号。如戴第2号面具的人看到了第5 号面具的编号。栋栋自己也会看到一些编号，他也会根据自己的面具编号把信息补充进去。由于并不是每个人都能记住自己所看到的全部编号，因此，栋栋收集的信 息不能保证其完整性。现在请你计算，按照栋栋目前得到的信息，至多和至少有多少类面具。由于主办方已经声明了k≥3，所以你必须将这条信息也考虑进去。

【输入格式】

输入文件party.in 第一行包含两个整数n, m，用一个空格分隔，n 表示主办方总共准备了多少个面具，m 表示栋栋收集了多少条信息。

接下来m 行，每行为两个用空格分开的整数a, b，表示戴第a 号面具的人看到了第b 号面具的编号。相同的数对a, b 在输入文件中可能出现多次。

【输出格式】

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

【输入样例一】

6 5

1 2

2 3

3 4

4 1

3 5

【输出样例一】

4 4

【输入样例二】

3 3

1 2

2 1

2 3

【输出样例二】

-1 -1

【数据规模和约定】

• 50%的数据，满足n ≤ 300, m ≤ 1000；
• 100%的数据，满足n ≤ 100000, m ≤ 1000000。

 

相当经典的一道题。

首先易知如果是一棵树，那么最大的可能答案即为最长链，最小的为3。

如果有环又会如何呢？

初始想法是答案为所有环长度的gcd。

但是看看下面这个例子：

边为

1 2

2 3

3 4

1 4,

看起来是没有环的约束的。但是答案如果为4，就不能满足1 4边的要求。其实解为2。说明上面那个想法是有问题的。

解决方法是把原图中的有向边拆成两种边：正边（跟原边方向相同，权值为1），反边（跟原边方向相反，权值为-1）。然后进行一次边的遍历dfs，如果当前点无标号，就把当前长度设为它的标号，否则把abs（当前长度-它的标号）作为一个环的长度，压入表中。

最后对所有表中元素求一个gcd。

这个想法非常巧妙！

其中树的情况有一个小技巧，不必再写一个过程遍历。以某个点为根进行遍历的过程中，需要记录最长路径和最短路径（即为最长的指向根的路径）。多个联通块的最长路径-最短路径加起来即可。

 

突然想到：NOI不可以开无限栈怎么办？

{$M 10000000} program xqz; const maxn=200000; maxm=2000000; oo=1 shl 30; var i,j,k,m,n,ss,tt,m1,m2,e,ans:longint; edge,z:array[0..maxn] of longint; point,next,a:array[-maxm..maxm] of longint; function gcd(a,b:longint):longint; var k:longint; begin if a<b then begin k:=a; a:=b; b:=k; end; while b<>0 do begin k:=a; a:=b; b:=k mod b; end; gcd:=a; end; procedure dfs(i,k:longint); var b,p:longint; begin if z[i]>0 then begin if z[i]=k then exit; inc(a[0]); a[a[0]]:=abs(z[i]-k); exit; end; z[i]:=k; b:=edge[i]; if k>m1 then m1:=k; if k<m2 then m2:=k; while b<>0 do begin p:=point[b]; if b>0 then dfs(p,k+1) else dfs(p,k-1); b:=next[b]; end; end; begin assign(input,'party.in'); reset(input); assign(output,'party.out'); rewrite(output); read(n,m); for i:=1 to m do begin read(ss,tt); inc(e); point[e]:=tt; next[e]:=edge[ss]; edge[ss]:=e; point[-e]:=ss; next[-e]:=edge[tt]; edge[tt]:=-e; end; ans:=0; for i:=1 to n do if z[i]=0 then begin m1:=0; m2:=oo; //必须从n开始，因为防止z变成0 dfs(i,n); ans:=ans+m1-m2+1; end; if a[0]>0 then begin ans:=0; for i:=1 to a[0] do ans:=gcd(ans,a[i]); for j:=3 to ans do if ans mod j=0 then break; if ans>=3 then writeln(ans,' ',j) else writeln('-1 -1'); end else if ans>=3 then writeln(ans,' ',3) else writeln('-1 -1'); close(input); close(output); end.