【BZOJ】1615: [Usaco2008 Mar]The Loathesome Hay Baler麻烦的干草打包机（模拟+bfs）

http://www.lydsy.com/JudgeOnline/problem.php?id=1615

这种题。。。。。

#include <cstdio>

#include <cstring>

#include <cmath>

#include <string>

#include <iostream>

#include <algorithm>

#include <queue>

using namespace std;

#define rep(i, n) for(int i=0; i<(n); ++i)

#define for1(i,a,n) for(int i=(a);i<=(n);++i)

#define for2(i,a,n) for(int i=(a);i<(n);++i)

#define for3(i,a,n) for(int i=(a);i>=(n);--i)

#define for4(i,a,n) for(int i=(a);i>(n);--i)

#define CC(i,a) memset(i,a,sizeof(i))

#define read(a) a=getint()

#define print(a) printf("%d", a)

#define dbg(x) cout << #x << " = " << x << endl

#define printarr(a, n, m) rep(aaa, n) { rep(bbb, m) cout << a[aaa][bbb]; cout << endl; }

inline const int getint() { int r=0, k=1; char c=getchar(); for(; c<'0'||c>'9'; c=getchar()) if(c=='-') k=-1; for(; c>='0'&&c<='9'; c=getchar()) r=r*10+c-'0'; return k*r; }

inline const int max(const int &a, const int &b) { return a>b?a:b; }

inline const int min(const int &a, const int &b) { return a<b?a:b; }

const int N=1105;

struct dat { int x, y, r; double w; }a[N];

bool vis[N];

int q[N], front, tail, X, Y, n, id;

double d[N], ans[N];

double dis(const dat &x, const dat &y) { return sqrt((x.x-y.x)*(x.x-y.x)+(x.y-y.y)*(x.y-y.y)); }

int main() {

	read(n); read(X); read(Y);

	for1(i, 1, n) {

		read(a[i].x), read(a[i].y), read(a[i].r);

		if(a[i].x==0 && a[i].y==0) id=i;

	}

	vis[id]=1; q[tail++]=id; d[id]=ans[id]=10000;

	while(front!=tail) {

		int x=q[front++]; if(front==N) front=0;

		if(a[x].x==X && a[x].y==Y) { printf("%d", abs((int)ans[x])); break; }

		for1(i, 1, n) if(!vis[i] && abs(dis(a[x], a[i])-(a[x].r+a[i].r))<1e-8) {

			d[i]=-d[x]*a[x].r/a[i].r;

			ans[i]+=abs(d[i])+ans[x];

			vis[i]=1;

			q[tail++]=i; if(tail==N) tail=0;

		}

	}

	return 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

Source

Silver