「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
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
HINT
输出说明:
齿轮 位置 半径 转速
1 (0,0) 10 10,000
2 (0,30) 20 -5,000
3 (32,54) 20 5,000
——
齿轮转速绝对值之和:20,000
题解
直接bfs，但是注意齿轮带动判的是相切T T
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 |
#include<iostream> #include<cstdio> #include<algorithm> using namespace std; 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; } int n,ex,ey,first,last; double ans,v[1105]; int q[1105],fa[1105]; bool vis[1105]; struct data{int x,y,r;}c[1105]; inline int sqr(int x) { return x*x; } inline int dis(int x,int y) { return sqr(c[x].x-c[y].x)+sqr(c[x].y-c[y].y); } inline bool tangent(int x,int y) { if(sqr(c[x].r+c[y].r)==dis(x,y))return 1; return 0; } void bfs() { int head=0,tail=1; q[head]=first;v[first]=10000;vis[first]=1; while(head!=tail) { int now=q[head];head++; for(int i=1;i<=n;i++) { if(vis[i])continue; if(tangent(now,i)) { vis[i]=1;fa[i]=now;q[tail++]=i; v[i]=v[now]*c[now].r/c[i].r; if(i==last)return; } } } } int main() { n=read();ex=read();ey=read(); for(int i=1;i<=n;i++) c[i].x=read(),c[i].y=read(),c[i].r=read(); for(int i=1;i<=n;i++) { if(!c[i].x&&!c[i].y)first=i; if(c[i].x==ex&&c[i].y==ey)last=i; } bfs(); for(int i=last;i;i=fa[i]) ans+=v[i]; printf("%d\n",(int)ans); return 0; } |
OI 课件=>
