2016河南省第九届ACM程序设计竞赛，问题 E: 机器设备

题目链接：http://acm.codedream.ren/JudgeOnline/

问题 E: 机器设备

时间限制: 1 Sec   内存限制: 64 MB
提交: 11   解决: 5
[ 提交][ 状态][ 讨论版]

题目描述

Alpha 公司设计出一种节能的机器设备。它的内部结构是由 N 个齿轮组成。整个机器设备有 一个驱动齿轮，当启动它时，它立即按 10,000 圈/小时转速顺时针转动，然后它又带动与它相切 的齿轮反方向，即逆时针转动。齿轮之间互相作用，每个齿轮都可能驱动着多个齿轮，最终带动 一个工作齿轮完成相应的任务。 在这套设备中，记录了每个齿轮的圆心坐标和齿轮半径。已知驱动齿轮位于（0,0），最终的 工作齿轮位于（Xt, Yt）。 Alpha 公司想知道传动序列中所有齿轮的转速。所谓传动序列，即能量由驱动齿轮传送，最 后到达工作齿轮的过程中用到的所有齿轮。能量传送过程是，在一个半径为 R，转速为 S 圈/每小 时的齿轮的带动下，与它相切的半径为 R’的齿轮的转速为-S*R/R’ 转/小时。负号的意思是， 表 示按反方向转动。

 

已知，机器设备中除了驱动齿轮以外，所有齿轮都可能被另外某个齿轮带动，并且不会出现 2 个不同的齿轮带动同一个齿轮的情况。 你的任务是计算整个传动序列中所有齿轮的能量之和。即所有齿轮转速的绝对值之和。 

输入

第一行： T 表示以下有 T 组测试数据（1≤T ≤8） 对每组测试数据： 第 1 行： N Xt Yt （2≤N ≤1100） 接下来有 N 行， Xi Yi Ri 表示 N 个齿轮的坐标和半径 i=1,2,….,N （ -5000 ≤Xi ，Yi ≤ 5000 3 ≤ Ri ≤ 1000 ） 坐标以及半径是整数

输出

每组测试数据，输出占一行，即所有齿轮转速的绝对值之和 在double范围内，输出整数部分

样例输入

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

样例输出

20000

简单的搜索，深搜广搜都可以，这题做了好长时间，怕超时所以用邻接表写了一下午，结果这题好像不用担心时间。

最后广搜过了。

AC代码：代码不复杂，只是结构体比较多，看着有点长。（最后的答案要的是整数部分，%.0lf输出是四舍五入）

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<vector>
#include<map>
#include<string>
using namespace std;
typedef long long int ll;
const int mod = 1e7+7;
const int INF = 0x3f3f3f3f;
const int maxn = 5e4 +10;
struct node
{
    int x,y,r,k;
    double z,sum;
}point[1150];
int n,gx,gy,yuandian,zhongdian,t;
int v[1150],start,endd;
double ans;
queue<node>que;
void init()
{
    memset(v,0,sizeof(v));
    ans=0.0;
    while(!que.empty())
        que.pop();
    yuandian=1;
    zhongdian=1;
}
bool xiangqie(int x,int y,int r,int xx,int yy,int rr)
{
      if((xx-x)*(xx-x)+(yy-y)*(yy-y)==(rr+r)*(rr+r))
        return 1;
      return 0;
}
int main()
{
    scanf("%d",&t);
    while(t--)
    {
        init();
        scanf("%d%d%d",&n,&gx,&gy);
        for(int i=1; i<=n; i++)
        {
            scanf("%d%d%d",&point[i].x,&point[i].y,&point[i].r);
            point[i].k=i;
            if(yuandian&&point[i].x==0&&point[i].y==0)
            {
                start=i;
                yuandian=0;
                point[i].z=10000.00;
            }
            if(zhongdian&&point[i].x==gx&&point[i].y==gy)
            {
                zhongdian=0;
                endd=i;
            }
        }
        if(start==endd)
        {
            printf("10000\n");
            continue;
        }
        node ka,kb;
        v[start]=1;
        ka.x=point[start].x,ka.y=point[start].y,ka.r=point[start].r,ka.z=10000.00,ka.k=start,ka.sum=10000.0;
        que.push(ka);
        while(!que.empty())
        {
            kb=que.front();
            que.pop();
            for(int i=1; i<=n; i++)
            {
                if(!v[i]&&xiangqie(kb.x,kb.y,kb.r,point[i].x,point[i].y,point[i].r))
                {
                    ka.x=point[i].x;
                    ka.y=point[i].y;
                    ka.r=point[i].r;
                    ka.k=i;
                    ka.z=kb.z*((double)kb.r)/((double)ka.r);
                    ka.sum=kb.sum+ka.z;
                    v[ka.k]=1;
                    que.push(ka);
                    if(i==endd)
                    {
                        ans=ka.sum;
                        break;
                    }
                }
            }
            if(ans>0.0)
                break;
        }
        printf("%d\n",(int)(ans+0.000001));//取整，处理误差
    }
}