2018.09.26【BZOJ4602】【洛谷P4079】【SDOI2016】齿轮（搜索）（图的遍历）

洛谷传送门

BZOJ传送门

---

解析：

这道题$O(V+E)$就过了啊。。。
然而官方标解是并查集。。。
我。。。

思路：

只要按照图的遍历一边跑一遍验证就行了。

$DFS$比较好写。

---

代码：

#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;
}