BZOJ 1615: [Usaco2008 Mar]The Loathesome Hay Baler麻烦的干草打包机

最新推荐文章于 2018-10-27 13:11:00 发布

jjs98

最新推荐文章于 2018-10-27 13:11:00 发布

阅读量847

点赞数

分类专栏：搜索文章标签： dfs

本文链接：https://blog.csdn.net/pzler/article/details/28410869

版权

搜索专栏收录该内容

41 篇文章 0 订阅

订阅专栏

Description

Farmer John新买的干草打包机的内部结构大概算世界上最混乱的了，它不象普通的机器一样有明确的内部传动装置，而是，N (2 <= N <= 1050)个齿轮互相作用，每个齿轮都可能驱动着多个齿轮。 FJ记录了对于每个齿轮i，记录了它的3个参数：X_i，Y_i表示齿轮中心的位置坐标(-5000 <= X_i <= 5000; -5000 <= Y_i <= 5000)；R_i表示该齿轮的半径(3 <= R_i <= 800)。驱动齿轮的位置为0,0，并且FJ也知道最终的工作齿轮位于X_t,Y_t。驱动齿轮顺时针转动，转速为10,000转/小时。你的任务是，确定传动序列中所有齿轮的转速。传动序列的定义为，能量由驱动齿轮传送到工作齿轮的过程中用到的所有齿轮的集合。对能量传送无意义的齿轮都应当被忽略。在一个半径为Rd，转速为S转/每小时的齿轮的带动下，与它相接的半径为Rx的齿轮的转速将为-S*Rd/Rx转/小时。S前的负号的意思是，一个齿轮带动的另一个齿轮的转向会与它的转向相反。 FJ只对整个传动序列中所有齿轮速度的绝对值之和感兴趣，你的任务也就相应转化成求这个值。机器中除了驱动齿轮以外的所有齿轮都被另外某个齿轮带动，并且不会出现2个不同的齿轮带动同一个齿轮的情况。相信你能轻易地写个程序来完成这些计算:)

Input

* 第1行: 3个用空格隔开的整数：N，X_t，Y_t

* 第2..N+1行: 第i+1描述了齿轮i的位置及半径：X_i，Y_i，以及R_i

Output

* 第1行: 输出所有在传动中起到作用的齿轮转速的绝对值，包括驱动齿轮和工作齿轮。只需要输出答案的整数部分

Sample Input

4 32 54
0 0 10
0 30 20
32 54 20
-40 30 20

机器里一共有4个齿轮，位于0,0的是半径为10的驱动齿轮，它带动了位于
0,30的，半径为20的某个齿轮。这个齿轮又间接带动了位于32,54，半径为20的
工作齿轮，以及一个位于-40,30，半径同样为20的冗余的齿轮。

Sample Output

20000

HINT

输出说明:

齿轮位置半径     转速

1 (0,0)     10     10,000

2 (0,30)    20     -5,000

3 (32,54)   20      5,000

                   ------

齿轮转速绝对值之和:20,000

题解

dfs，回溯。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cstdlib>
#include<cmath>
using namespace std;
int n,xt,yt,zz,head[1052];
int st,usd;
struct chilun {int x,y; double r;} a[1052];
struct bian {int frm,to,nx;} e[500002];
bool pd[1052];
double s=10000.0,q[1052];
bool check(int i,int j)
{
    double l=sqrt((a[i].x-a[j].x)*(a[i].x-a[j].x)+(a[i].y-a[j].y)*(a[i].y-a[j].y));
    if(l==a[i].r+a[j].r) return true;
    else return false;
}
void init()
{
    scanf("%d%d%d",&n,&xt,&yt);
    for(int i=1;i<=n;i++)
       {scanf("%d%d%lf",&a[i].x,&a[i].y,&a[i].r);
        if(a[i].x==0&&a[i].y==0) st=i;
       }
    for(int i=1;i<=n;i++)
    for(int j=i+1;j<=n;j++)
       {if(check(i,j))
           {zz++; e[zz].frm=i; e[zz].to=j; e[zz].nx=head[i]; head[i]=zz;
            zz++; e[zz].frm=j; e[zz].to=i; e[zz].nx=head[j]; head[j]=zz;
           }
       }
}
void dfs(int w)
{
    if(a[w].x==xt&&a[w].y==yt)
       {double ans=s;
        for(int i=1;i<=usd;i++) ans+=q[i];
        printf("%d",(int)ans);
        return ;
       }
    for(int i=head[w]; i;i=e[i].nx)
       {int p=e[i].to;
        if(pd[p]) continue;
        pd[p]=1;
        double yl=s;
        q[++usd]=s;s=s*a[w].r/a[p].r;
        dfs(p);
        --usd; s=yl;
       }
}
int main()
{
    init();
    pd[st]=1; dfs(st);
    return 0;
}