BZOJ1064 NOI2008假面舞會(dfs樹)
阿新 • • 發佈:2018-09-13
則無 如果 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樹)