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

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

Time Limit: 5 Sec   Memory Limit: 64 MB
Submit: 972   Solved: 400
[ 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

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<stdio.h>
#include<vector>
using namespace std;
vector<int> G[1055];
typedef struct
{
	int x, y;
	int r;
}Res;
Res s[1055];
int len[1055];
void Sech(int u, int step)
{
	int i, v;
	len[u] = step;
	for(i=0;i<G[u].size();i++)
	{
		v = G[u][i];
		if(len[v]!=0)
			continue;
		Sech(v, step+1);
	}
}
int main(void)
{
	double ans, now;
	int n, x, i, j, S, T, y, u, v;
	scanf("%d%d%d", &n, &x, &y);
	for(i=1;i<=n;i++)
	{
		scanf("%d%d%d", &s[i].x, &s[i].y, &s[i].r);
		if(s[i].x==x && s[i].y==y)
			T = i;
		if(s[i].x==0 && s[i].y==0)
			S = i;
	}
	for(i=1;i<=n-1;i++)
	{
		for(j=i+1;j<=n;j++)
		{
			if((s[i].x-s[j].x)*(s[i].x-s[j].x)+(s[i].y-s[j].y)*(s[i].y-s[j].y)==(s[i].r+s[j].r)*(s[i].r+s[j].r))
			{
				G[i].push_back(j);
				G[j].push_back(i);
			}
		}
	}
	Sech(T, 1);
	u = S, ans = now = 10000;
	while(u!=T)
	{
		for(i=0;i<G[u].size();i++)
		{
			v = G[u][i];
			if(len[v]<len[u])
				break;
		}
		now = now*s[u].r/s[v].r;
		ans += now;
		u = v;
	}
	printf("%d\n", (int)ans);
	return 0;
}