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");
else
printf("No\n");
}
return 0;
}