2018.09.26【BZOJ4602】【洛谷P4079】【SDOI2016】齒輪(搜尋)(圖的遍歷)
阿新 • • 發佈:2018-12-12
洛谷傳送門
解析:
這道題就過了啊。。。 然而官方標解是並查集。。。 我。。。
思路:
只要按照圖的遍歷一邊跑一遍驗證就行了。
比較好寫。
程式碼:
#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;
}