【BZOJ4602】齿轮，带权并查集

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

---

传送门
Ps
这段时间在狂补文化课，10天就碰过一次键盘，结果回来发现连快速读入都敲不对了233
接下来一段时间仍保持这种状态了……
思路：
R2唯一可做的一道题？
在学校弄了个30分做法
回来一交A了……
感觉分解质因数的方法高大上啊有木有
sum[x]表示x是其所在集合的代表元素fa[x]的几分之几（所以sum[fa[x]]=1,如果W[fa[x]]W[x]=AB，那么sum[x]=BA）
询问判断：如果u,v在同一集合里直接判断即可，反之将两者合并
合并思路：如果我们要合并u，v两个元素所在集合且uv=xy，那么我们需要解决的就是fa[u],fa[v]之间的数量关系，让sum[fa[u]]=sum[v]∗xsum[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();}