【BZOJ 1064】 [Noi2008]假面舞会

最新推荐文章于 2021-12-21 13:21:43 发布

Regina8023

最新推荐文章于 2021-12-21 13:21:43 发布

阅读量1.2k

点赞数 3

分类专栏： ~图论~ 文章标签： BZOJ OI 图论

本文链接：https://blog.csdn.net/regina8023/article/details/43369885

版权

~图论~ 专栏收录该内容

7 篇文章 0 订阅

订阅专栏

1064: [Noi2008]假面舞会

Time Limit: 10 Sec Memory Limit: 162 MB
Submit: 988 Solved: 507
[ Submit][ Status]

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号面具的人能看到戴1号面具的人!

如果a能看到b，我们从a到b连一条有向边。

分三种情况考虑：

1.无环且每个点入度为1：

k最大值就是最长链的长度；最小值显然就是3了。

2.有环：

在环中一定存在k看到1的情况，因此k必然是环长的约数。

k最大值就是所有环长最大公约数，最小值就是最大公约数的大于3的最小约数

3.无环，但是某些点的入度>1

红点的入度为2，他在上面那条链上是第五个，下面的链上是第三个，可以发现k此时只能是1或2。

因为红点只能有一个标号，所以上面的链长与下面的链长模k一定是同余的。

因此k一定是链长差值的约数。

对于这三种情况，应该如何实现呢？

有一个重要的技巧，对于a看到b，连a-->b权值为1的边，b-->a权值为-1的边。

这样连边之后第三种情况就可以和有环的情况一起讨论了！

#include <iostream>
#include <cstring>
#include <algorithm>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#define M 100005
using namespace std;
int ma,tot=0,mi,cirma,du[M],n,m,h[M],d[M],v[M];
struct edge
{
	int y,ne,v;
}e[M*20];
void read(int &tmp)
{
	tmp=0;
	char ch=getchar();
	int fu=1;
	for (;ch<'0'||ch>'9';ch=getchar())
		if (ch=='-') fu=-1;
	for (;ch>='0'&&ch<='9';ch=getchar())
		tmp=tmp*10+ch-'0';
	tmp*=fu;
}
void Addedge(int x,int y,int v)
{
	e[++tot].y=y;
	e[tot].v=v;
	e[tot].ne=h[x];
	h[x]=tot;
}
int gcd(int a,int b)
{
	return b==0?a:gcd(b,a%b);
}
void dfs(int root,int x)
{
	v[x]=1;
	for (int i=h[x];i;i=e[i].ne)
	{
		int y=e[i].y;
		if (!v[y]) 
		{
			d[y]=d[x]+e[i].v;
			dfs(root,y);
		}
		else
		{
			int k=abs(d[x]+e[i].v-d[y]);
			cirma=gcd(cirma,k);
		}
	}
}
void dfs2(int x)
{
	v[x]=1;
	for (int i=h[x];i;i=e[i].ne)
	{
		int y=e[i].y;
		if (v[y]) continue;
		d[y]=d[x]+e[i].v;
		dfs2(y);
	}
	if (ma<d[x]) ma=d[x];
	if (mi>d[x]) mi=d[x];
}
int main()
{
	read(n),read(m);
	cirma=0;
	for (int i=1;i<=m;i++)
	{
		int x,y;
		read(x),read(y);
		du[x]++;
        Addedge(x,y,1),Addedge(y,x,-1);
	}
	for (int i=1;i<=n;i++)
		v[i]=0,d[i]=0;
	for (int i=1;i<=n;i++)
		if (!v[i])
			d[i]=1,dfs(i,i);
	if (!cirma)
	{
		int tot=0;
		for (int i=1;i<=n;i++)
			v[i]=0;
		for (int i=1;i<=n;i++)
		{
			if (v[i]) continue;
			ma=0,mi=0,d[i]=0,dfs2(i),tot=tot+(ma-mi+1);
		}
		if (tot>=3) printf("%d 3\n",tot);
		else printf("-1 -1\n");
	}
	else
	{
		if (cirma<3) printf("-1 -1\n");
		else
		{
			int cirmin=0;
			for (int i=3;i<=cirma;i++)
				if (cirma%i==0)
				{
					cirmin=i;
					break;
				}
			printf("%d %d\n",cirma,cirmin);
		}
	}
	return 0;
}