题目
1615: [Usaco2008 Mar]The Loathesome Hay Baler麻烦的干草打包机
Time Limit: 5 Sec Memory Limit: 64 MBDescription
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的。但是要注意转速和答案请用double类型。我Wa两次,一次没学好齿轮联动,一次脑补输出正数以为直接整形存储就行【没考虑小数进位啊,魂淡!
代码
/*Author:WNJXYK*/ #include<cstdio> #include<iostream> #include<cstring> #include<string> #include<algorithm> #include<queue> #include<set> #include<map> using namespace std; #define LL long long #define Inf 2147483647 #define InfL 10000000000LL inline void swap(int &x,int &y){int tmp=x;x=y;y=tmp;} inline void swap(LL &x,LL &y){LL tmp=x;x=y;y=tmp;} inline int remin(int a,int b){if (a<b) return a;return b;} inline int remax(int a,int b){if (a>b) return a;return b;} inline LL remin(LL a,LL b){if (a<b) return a;return b;} inline LL remax(LL a,LL b){if (a>b) return a;return b;} struct Edge{ int u,v,nxt; Edge(){} Edge(int a,int b,int d){u=a;v=b;nxt=d;} }; Edge e[20005]; int nume; int h[1151]; queue<int> que; double dist[1151]; int father[1151]; int lx[1151],ly[1151],lr[1151]; inline void addEdge(int u,int v){ e[++nume]=Edge(u,v,h[u]); h[u]=nume; e[++nume]=Edge(v,u,h[v]); h[v]=nume; } int n,fx,fy; double Ans=0; bool spfa(int x,int endss){ while(!que.empty()) que.pop(); que.push(x); memset(dist,-127,sizeof(dist)); dist[x]=10000; while(!que.empty()){ int p=que.front(); que.pop(); for (int i=h[p];i;i=e[i].nxt){ int v=e[i].v; if (dist[v]<0){ dist[v]=dist[p]*(double)lr[p]/(double)lr[v]; que.push(v); father[v]=p; if (v==endss) return true; } } } return true; } int main(){ scanf("%d%d%d",&n,&fx,&fy); for (int i=1;i<=n;i++) scanf("%d%d%d",&lx[i],&ly[i],&lr[i]); for (int i=1;i<=n;i++){ for (int j=i+1;j<=n;j++){ if ((lx[i]-lx[j])*(lx[i]-lx[j]) +(ly[i]-ly[j])*(ly[i]-ly[j]) ==(lr[i]+lr[j])*(lr[i]+lr[j])){ addEdge(i,j); } } } int starts=1,ends=1; for (int i=2;i<=n;i++){ if (lx[i]==0 && ly[i]==0) starts=i; if (lx[i]==fx && ly[i]==fy) ends=i; } father[starts]=0; spfa(starts,ends); for (int i=ends;i;i=father[i]) Ans+=dist[i]; printf("%d ",(int)Ans); return 0; }