題目
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\n",(int)Ans);
return 0;
}