bzoj 4602: [Sdoi2016]齿轮

Description

现有一个传动系统，包含了N个组合齿轮和M个链条。每一个链条连接了两个组合齿轮u和v，并提供了一个传动比x 

: y。即如果只考虑这两个组合齿轮，编号为u的齿轮转动x圈，编号为v的齿轮会转动y圈。传动比为正表示若编号

为u的齿轮顺时针转动，则编号为v的齿轮也顺时针转动。传动比为负表示若编号为u的齿轮顺时针转动，则编号为v

的齿轮会逆时针转动。若不同链条的传动比不相容，则有些齿轮无法转动。我们希望知道，系统中的这Ｎ个组合齿

轮能否同时转动。

Input

有多组数据，第一行给定整数Ｔ，表示总的数据组数，之后依次给出Ｔ组数据。每一组数据的第一行给定整数Ｎ和

Ｍ，表示齿轮总数和链条总数。之后有Ｍ行，依次描述了每一个链条，其中每一行给定四个整数u，v，x和y，表示

只考虑这一组联动关系的情况下，编号为u的齿轮转动x圈，编号为v的齿轮会转动y圈。请注意，x为正整数，而y为

非零整数，但是y有可能为负数。

T<=32，N<=1000，M<=10000且x与y的绝对值均不超过100

Output

输出T行，对应每一组数据。首先应该输出标识这是第几组数据，参见样例输出。之后输出判定结果，如果N个组合

齿轮可以同时正常运行，则输出Yes，否则输出No。

Sample Input

2
3 3
1 2 3 5
2 3 5 -7
1 3 3 -7
3 3
1 2 3 5
2 3 5 -7
1 3 3 7

Sample Output

Case #1: Yes

Case #2: No

似乎精度不会爆炸啊。。直接用double连边，边权为x/y

dfs看是否矛盾即可

#include<cstdio>#include<string>#include<cstring>#include<algorithm>using namespace std;struct line{int s,t;double x;int next;}a[100001];int head[1001];int edge;inline void add(int s,int t,double x){a[edge].next=head[s];head[s]=edge;a[edge].s=s;a[edge].t=t;a[edge].x=x;}double eps=1e-8;inline double absx(double x){if(x<0)x=-x;return x;}double sx[1001];bool v[1001];inline bool dfs(int d){bool flag=true;int i;for(i=head[d];i!=0;i=a[i].next){int t=a[i].t;if(!v[t]){v[t]=true;sx[t]=sx[d]*a[i].x;if(!dfs(t))flag=false;}else if(absx(sx[t]-sx[d]*a[i].x)>eps)flag=false;}return flag;}int main(){int T;int k=0;scanf("%d",&T);while(T>0){k++;T--;memset(head,0,sizeof(head));memset(a,0,sizeof(a));edge=0;int n,m;scanf("%d%d",&n,&m);int i;int s,t;double x,y;for(i=1;i<=m;i++){scanf("%d%d%lf%lf",&s,&t,&x,&y);edge++;add(s,t,x/y);edge++;add(t,s,y/x);}memset(v,false,sizeof(v));memset(sx,0,sizeof(sx));bool flag=true;for(i=1;i<=n;i++){if(!v[i]){sx[i]=1.0;v[i]=true;if(!dfs(i))flag=false;}}printf("Case #%d: ",k);if(flag)printf("Yes\n");elseprintf("No\n");}return 0;}