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。
题解
byvoid神犇:https://www.byvoid.com/blog/noi-2008-party/
顺便提一下。如果最后gcd=0,答案是各连通块最大标号-最小标号+1的和,可以用个并查集
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 |
#include<map> #include<cmath> #include<queue> #include<cstdio> #include<vector> #include<cstring> #include<cstdlib> #include<iostream> #include<algorithm> #define inf 1000000000 #define ll long long using namespace std; ll read() { ll x=0,f=1;char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();} return x*f; } int T,n,m,cnt,ans,ans2,tot; int mx[100005],mn[100005],fa[100005]; int last[100005],ind[100005],vis[100005],len[100005]; struct edge{ int to,next,v; }e[400005]; int find(int x) { return x==fa[x]?x:fa[x]=find(fa[x]); } int gcd(int a,int b) { return b==0?a:gcd(b,a%b); } void insert(int u,int v,int w) { e[++cnt]=(edge){v,last[u],w};last[u]=cnt; } void dfs(int x) { vis[x]=1; mx[T]=max(mx[T],ind[x]); mn[T]=min(mn[T],ind[x]); for(int i=last[x];i;i=e[i].next) { int t=ind[x]+e[i].v; if(!vis[e[i].to]) { ind[e[i].to]=t; dfs(e[i].to); } else ans=gcd(ans,abs(t-ind[e[i].to])); } } int main() { memset(mn,127,sizeof(mn)); n=read();m=read(); for(int i=1;i<=n;i++)fa[i]=i; for(int i=1;i<=m;i++) { int u=read(),v=read(); insert(u,v,1); insert(v,u,-1); fa[find(u)]=find(v); } for(int i=1;i<=n;i++) if(!vis[i]) { T=i; dfs(i); } for(int i=3;i<=ans;i++) if(ans%i==0){ans2=i;break;} ans2=max(ans2,3); if(!ans) { for(int i=1;i<=n;i++) { int t=find(i); mn[t]=min(mn[t],mn[i]); mx[t]=max(mx[t],mx[i]); } for(int i=1;i<=n;i++)if(i==fa[i])tot+=mx[i]-mn[i]+1; ans=tot; } if(ans<3)puts("-1 -1\n"); else printf("%d %d\n",ans,ans2); return 0; } |
黄学长,请问对于链的情况(gcd==0)对每个连通块做一遍拓扑排序求最大深度行不行?