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

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

唉……这水题做了老半天……主要是题意没看清楚

题意是告诉你两个点，起点是（0,0），终点是输入第一行给定的。要找出从起点到终点的那条链，求速度之和，连通性的判定在题目中已经讲的很清楚了

就是爆搜啦

#include<cstdio>
struct point{
	int x,y,r;
	double v;
}p[10001];
int n,cnt,t,w=1,start,end;
int head[10001],q[10001],fa[10001];
bool mark[10001];
double tot;
inline double abs(double x)
{if (x<0)x=-x;return x;}
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
inline bool mrk(int x,int y)
{
	int xx=p[x].x-p[y].x,yy=p[x].y-p[y].y,rr=p[x].r+p[y].r;
	int dist=xx*xx+yy*yy;
	if (dist!=rr*rr)return 0;
	return 1;
}
inline void bfs()
{
	q[1]=start;mark[start]=1;
	while (t<w)
	{
		int now=q[++t];
		for (int i=1;i<=n;i++)
		{
			if(mark[i])continue;
			if (!mrk(now,i))continue;
			q[++w]=i;
			fa[i]=now;
			p[i].v=p[now].v*p[now].r/p[i].r;
			mark[i]=1;
			if (i==end) return;
		}
	}
}
int main()
{
	n=read();p[0].x=read();p[0].y=read();
	for (int i=1;i<=n;i++)
	  {
	  	p[i].x=read();
	  	p[i].y=read();
	  	p[i].r=read();
	  }
	for (int i=1;i<=n;i++)
	  {
	  	if(!p[i].x&&!p[i].y) start=i;
	  	if(p[i].x==p[0].x&&p[i].y==p[0].y) end=i;
	  }
	p[start].v=10000;
	bfs();
	for (int i=end;i;i=fa[i])
	  tot+=abs(p[i].v);
	printf("%d",(int)tot);
	return 0;
}