1. 程式人生 > >BZOJ1064 NOI2008假面舞會(dfs樹)

BZOJ1064 NOI2008假面舞會(dfs樹)

則無 如果 online lag += col using 數量 noi

  將圖中的環的長度定義為正向邊數量-反向邊數量,那麽答案一定是所有環的環長的gcd。dfs一下就能找到圖中的一些環,並且圖中的所有環的環長都可以由這些環長加加減減得到(好像不太會證)。如果有環長為1或2則無解。

  沒有環的話圖就是一個有向樹。類似定義鏈的長度,那麽一個連通塊內答案就是最長鏈,也即dfs樹上最大深度-最小深度+1,對所有連通塊累加即可。

#include<iostream> 
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include
<algorithm> using namespace std; int read() { int x=0,f=1;char c=getchar(); while (c<0||c>9) {if (c==-) f=-1;c=getchar();} while (c>=0&&c<=9) x=(x<<1)+(x<<3)+(c^48),c=getchar(); return x*f; } #define N 100010 #define M 1000010 int n,m,p[N],deep[N],t=0
,ans=0,mx,mn; bool flag[N]; struct data{int to,nxt,len; }edge[M]; void addedge(int x,int y,int z){t++;edge[t].to=y,edge[t].nxt=p[x],edge[t].len=z,p[x]=t;} int gcd(int n,int m){return m==0?n:gcd(m,n%m);} void dfs(int k,int from) { if (ans==-1) return; flag[k]=1; for (int i=p[k];i;i=edge[i].nxt)
if (edge[i].to!=from) { if (!flag[edge[i].to]) deep[edge[i].to]=deep[k]+edge[i].len,dfs(edge[i].to,k); else if (abs(deep[k]-deep[edge[i].to]+edge[i].len)==1||abs(deep[k]-deep[edge[i].to]+edge[i].len)==2) {ans=-1;return;} else ans=gcd(ans,abs(deep[k]-deep[edge[i].to]+edge[i].len)); } } void find(int k) { flag[k]=1;mn=min(mn,deep[k]),mx=max(mx,deep[k]); for (int i=p[k];i;i=edge[i].nxt) if (!flag[edge[i].to]) find(edge[i].to); } int main() { #ifndef ONLINE_JUDGE freopen("bzoj1064.in","r",stdin); freopen("bzoj1064.out","w",stdout); const char LL[]="%I64d\n"; #else const char LL[]="%lld\n"; #endif n=read(),m=read(); for (int i=1;i<=m;i++) { int x=read(),y=read(); addedge(x,y,1),addedge(y,x,-1); } for (int i=1;i<=n;i++) if (!flag[i]) dfs(i,i); if (ans==-1) {cout<<-1<< <<-1;return 0;} else if (ans==0) { memset(flag,0,sizeof(flag)); for (int i=1;i<=n;i++) if (!flag[i]) mx=-N,mn=N,find(i),ans+=mx-mn+1; if (ans>2) cout<<ans<< <<3;else cout<<-1<< <<-1; } else for (int i=3;i<=ans;i++) if (ans%i==0) {cout<<ans<< <<i;break;} return 0; }

BZOJ1064 NOI2008假面舞會(dfs樹)