傳送門

題解：

總結幾種做法，並指出其中的優劣；

主要是兩種演算法，一種是暴力dfs，另一種是並查集。

事實上這種“每條邊都考慮”大概都可以用上述兩種辦法處理，並查集複雜度略高，但是其實近乎線性。

暴力dfs就是，建一張無向圖（不能是有向的），然後對於每個聯通塊，第一個元素設為1，然後通過第一個元素，

算出其它元素的值，然後對於點x，以及邊(x,y)，如果y也訪問過了，就判斷一下通過x算的y和已經算出來的y是否相等即可。

並查集是類似的，因為是無向圖，先加邊加成樹（或者森林），然後之後在加邊就判斷即可。

其實這題最主要的是精度問題。

但是出題人特別良心沒有卡精度，double都水過去了。

精度控制的幾種辦法：

第一種，最好想的，直接上double(或者long double)。

問題：資料理論最大是100^1000的，理論上會爆（事實沒有）。

第二種，比較容易想到的，取對數，變成加減運算。

問題：可以卡精度，例如1e18和(1e9-1)*(1e9+1)。但是也沒有卡。

第三種，在模質數意義下進行運算，除法變成乘逆元，判斷類似雜湊。

問題：理論上會衝突，實際上不會（幾乎不會）。

第四種，注意到x,y<=100可以質因數分解。

問題：這個應該能夠保證正確，但是常數略大，（儘管對於這個題是足夠的）。

程式碼：用的是第一種SB方法竟然沒有被卡QwQ

//BZOJ 4602
//SDOI 2016
#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#define MAXN 1100
#define MAXM 21000
using namespace std;
struct edges{
	int to,pre;
	double wgt;
}e[MAXM];
int h[MAXN],etop;
int add_edge(int u,int v,double wgt)
{
	etop++;
	e[etop].to=v;
	e[etop].pre=h[u];
	e[etop].wgt=wgt;
	h[u]=etop;
	return 0;
}
double val[MAXN];
const double ep=0.000001;
bool vis[MAXN];
bool dcmp(double x)
{
	if(fabs(x)<ep) return false;
	else return true;
}
bool dfs(int x)
{
	vis[x]=true;
	for(int i=h[x];i;i=e[i].pre)
		if(!vis[e[i].to])
		{
			val[e[i].to]=val[x]*e[i].wgt;
			if(!dfs(e[i].to)) return false;
		}
		else if(dcmp(val[e[i].to]-val[x]*e[i].wgt)) return false;
	return true;
}
int main()
{
	int T;scanf("%d",&T);
	for(int t=1;t<=T;t++)
	{
		int n,m;scanf("%d%d",&n,&m);
		memset(h,0,sizeof(h));etop=0;
		memset(vis,false,sizeof(vis));
		for(int i=1;i<=m;i++)
		{
			int u,v,x,y;scanf("%d%d%d%d",&u,&v,&x,&y);
			add_edge(u,v,(double)y/x);
			add_edge(v,u,(double)x/y);
		}
		for(int i=1;i<=n;i++)
			if(!vis[i])
			{
				val[i]=1.0;
				if(!dfs(i))
				{
					printf("Case #%d: No\n",t);
					goto loop;
				}
			}
		printf("Case #%d: Yes\n",t);
		loop:;
	}
	return 0;
}