1. 程式人生 > >2018.09.26【BZOJ4602】【洛谷P4079】【SDOI2016】齒輪(搜尋)(圖的遍歷)

2018.09.26【BZOJ4602】【洛谷P4079】【SDOI2016】齒輪(搜尋)(圖的遍歷)

洛谷傳送門

解析:

這道題O(V+E)O(V+E)就過了啊。。。 然而官方標解是並查集。。。 我。。。

思路:

只要按照圖的遍歷一邊跑一遍驗證就行了。

DFSDFS比較好寫。

程式碼:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
#define re register
#define gc getchar
#define pc putchar
#define cs const

inline
ll getint(){
    re ll num;
    re char c;
re bool f=0; while(!isdigit(c=gc()))f^=c=='-';num=c^48; while(isdigit(c=gc()))num=(num<<1)+(num<<3)+(c^48); return f?-num:num; } inline void outint(ll a){ static char ch[23]; if(a==0)pc('0'); if(a<0)pc('-'),a=-a; while(a)ch[++ch[0]]=a-a/10*10,a/=10; while(ch[
0])pc(ch[ch[0]--]^48); } cs int N=1002,M=100002; cs double begin=1.0*2*3*5*7*9*11*13; bitset<N> vis; double sp[N]; int n,m; int last[N],nxt[M<<1],to[M<<1],ecnt=1; double w[M<<1]; inline void addedge(int u,int v,int x,int y){ nxt[++ecnt]=last[u],last[u]=ecnt,to[ecnt]=v,w[ecnt]=(double
)y/x; nxt[++ecnt]=last[v],last[v]=ecnt,to[ecnt]=u,w[ecnt]=(double)x/y; } inline bool dfs(int u){ vis[u]=true; for(int re e=last[u],v=to[e];e;v=to[e=nxt[e]]){ if(vis[v]){ if(fabs(sp[u]*w[e]-sp[v])>1e-4)return false; } else { sp[v]=sp[u]*w[e]; if(!dfs(v))return false; } } return true; } signed main(){ int T=getint(),tt=0; while(T--){ ++tt; printf("Case #%d: ",tt); vis.reset();ecnt=1; memset(last,0,sizeof last); n=getint(),m=getint(); for(int re i=1;i<=m;++i){ int u=getint(),v=getint(),x=getint(),y=getint(); addedge(u,v,x,y); } bool flag=true; for(int re i=1;i<=n;++i){ if(!vis[i]){ sp[i]=begin; if(!dfs(i))flag=false; } } puts(flag?"Yes":"No"); } return 0; }