【BZOJ4602】齿轮，带权并查集

Time:2016.06.09
Author:xiaoyimi
转载注明出处谢谢

---

传送门
$Ps$
这段时间在狂补文化课，10天就碰过一次键盘，结果回来发现连快速读入都敲不对了233
接下来一段时间仍保持这种状态了……
思路：
R2唯一可做的一道题？
在学校弄了个30分做法
回来一交A了……
感觉分解质因数的方法高大上啊有木有
$sum[x]$表示x是其所在集合的代表元素fa[x]的几分之几（所以$sum[fa[x]]$=1,如果$\frac {W[fa[x]]}{W[x]}=\frac AB$，那么sum[x]=$\frac BA$）
询问判断：如果u,v在同一集合里直接判断即可，反之将两者合并
合并思路：如果我们要合并u，v两个元素所在集合且$\frac uv=\frac xy$，那么我们需要解决的就是fa[u],fa[v]之间的数量关系，让$sum[fa[u]]=\frac{sum[v]*x}{sum[u]*y}$即可
find操作：先让$fa[x]$合并到一个新集合里，这时$sum[fa[x]]$已经更新过了，只要让$sum[x]=sum[x]*sum[fa[x]]$就行了
（以上一些式子的推导很简单，我就不列过程了）
注意：
暂无
代码：

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cmath>
using namespace std;
int T,n,m;
int fa[1003];
double sum[1003];
int in()
{
    int t=0;char ch=getchar();bool f=0;
    while (ch>'9'||ch<'0') f=(ch=='-'),ch=getchar();
    while (ch>='0'&&ch<='9') t=(t<<1)+(t<<3)+ch-48,ch=getchar();
    return f?-t:t;
}
int find(int x)
{
    if (x==fa[x]) return x;
    int t=fa[x];
    fa[x]=find(t);
    sum[x]*=sum[t];
    return fa[x];
}
void work()
{
    n=in();m=in();
    for (int i=1;i<=n;i++) fa[i]=i,sum[i]=1;
    int x,y,u,v;bool flag=0;
    while (m--)
    {
        u=in();v=in();x=in();y=in();
        int f1=find(u),f2=find(v);
        if (f1==f2)
            if (fabs(sum[u]/sum[v]-1.0*x/y)>1e-9) flag=1;
            else;
        else
            fa[f1]=f2,
            sum[f1]=sum[v]*x/(sum[u]*y);
    }
    if (flag) puts("No");
    else puts("Yes");
}
main()
{
    T=in();
    for (int i=1;i<=T;i++)
        printf("Case #%d: ",i),work();
}