bzoj 4602: [Sdoi2016]齿轮
Posted lqybzx
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj 4602: [Sdoi2016]齿轮相关的知识,希望对你有一定的参考价值。
Description
现有一个传动系统,包含了N个组合齿轮和M个链条。每一个链条连接了两个组合齿轮u和v,并提供了一个传动比x : y。即如果只考虑这两个组合齿轮,编号为u的齿轮转动x圈,编号为v的齿轮会转动y圈。传动比为正表示若编号 为u的齿轮顺时针转动,则编号为v的齿轮也顺时针转动。传动比为负表示若编号为u的齿轮顺时针转动,则编号为v 的齿轮会逆时针转动。若不同链条的传动比不相容,则有些齿轮无法转动。我们希望知道,系统中的这N个组合齿 轮能否同时转动。Input
有多组数据,第一行给定整数T,表示总的数据组数,之后依次给出T组数据。每一组数据的第一行给定整数N和 M,表示齿轮总数和链条总数。之后有M行,依次描述了每一个链条,其中每一行给定四个整数u,v,x和y,表示 只考虑这一组联动关系的情况下,编号为u的齿轮转动x圈,编号为v的齿轮会转动y圈。请注意,x为正整数,而y为 非零整数,但是y有可能为负数。 T<=32,N<=1000,M<=10000且x与y的绝对值均不超过100Output
输出T行,对应每一组数据。首先应该输出标识这是第几组数据,参见样例输出。之后输出判定结果,如果N个组合 齿轮可以同时正常运行,则输出Yes,否则输出No。Sample Input
23 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: YesCase #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;
以上是关于bzoj 4602: [Sdoi2016]齿轮的主要内容,如果未能解决你的问题,请参考以下文章