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;
}