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

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

Time Limit: 5 Sec Memory Limit: 64 MB Submit: 980 Solved: 403 [Submit][Status][Discuss]

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

bfs

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
char buf[10000000], *ptr = buf - 1;
inline int readint(){
    int f = 1, n = 0;
    char ch = *++ptr;
    while(ch < '0' || ch > '9'){
        if(ch == '-') f = -1;
        ch = *++ptr;
    }
    while(ch <= '9' && ch >= '0'){
        n = (n << 1) + (n << 3) + ch - '0';
        ch = *++ptr;
    }
    return f * n;
}
const int maxn = 1000 + 100;
inline int sqr(const int &x){
    return x * x;
}
int x[maxn], y[maxn], r[maxn];
double v[maxn];
int q[maxn], h, t, pre[maxn] = {0};
bool vis[maxn] = {false};
int main(){
    fread(buf, sizeof(char), sizeof(buf), stdin);
    int N, ex, ey, s, e;
    N = readint();
    ex = readint();
    ey = readint();
    for(int i = 1; i <= N; i++){
        x[i] = readint();
        y[i] = readint();
        r[i] = readint();
        if(!x[i] && !y[i]) s = i;
        if(x[i] == ex && y[i] == ey) e = i;
    }
    v[s] = 10000;
    h = t = 0;
    q[t++] = s;
    vis[s] = true;
    int now;
    while(h != t){
        now = q[h++];
        for(int i = 1; i <= N; i++)
            if(!vis[i] && sqr(x[now] - x[i]) + sqr(y[now] - y[i]) == sqr(r[now] + r[i])){
                vis[i] = true;
                pre[i] = now;
                v[i] = v[now] * r[now] / r[i];
                q[t++] = i;
            }
    }
    double ans = 0;
    for(int i = e; i; i = pre[i]) ans += v[i];
    printf("%d
", (int)ans);
    return 0;
}