bzoj1615 麻烦的干草打包机

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

HINT

 

输出说明:


齿轮 位置  半径     转速

1 (0,0)     10     10,000

2 (0,30)    20     -5,000

3 (32,54)   20      5,000

                   ------

齿轮转速绝对值之和:20,000

 
此题看懂了以后模拟即可,注意判断圆相切,然后整个是一棵树。  
 1 #include<bits/stdc++.h>
 2 using namespace std;  
 3 #define R register int
 4 #define rep(i,a,b) for(R i=a;i<=b;i++) 
 5 #define Rep(i,a,b) for(R i=a;i>=b;i--)
 6 #define ms(i,a)    memset(a,i,sizeof(a)) 
 7 #define LL         long long 
 8 #define gc()       getchar() 
 9 #define sqr(x)     ((x)*(x)) 
10 template<class T>void read(T &x){
11   x=0; char c=0; int w=0;  
12   while (!isdigit(c)) w+=c=='-',c=gc();  
13   while (isdigit(c)) x=x*10+(c^48),c=gc();
14   if(w) x=-x; 
15 }
16 int const maxn=1050+3;  
17 int n,x[maxn],y[maxn],xt,yt,r[maxn],a[maxn][maxn],vis[maxn],tin[maxn],tout[maxn],sum,id;  
18 double ans;  
19 int dist(int i,int j){return sqr(x[i]-x[j])+sqr(y[i]-y[j]);}
20 void dfs(int o,int fa){
21   tin[o]=++sum;  
22   rep(i,1,n) if(a[o][i] && i!=fa)  dfs(i,o);  
23   tout[o]=++sum;  
24 }
25 int ancestor(int x,int y){ 
26   return tin[x]<=tin[y] && tout[x]>=tout[y];   
27 }
28 void dfs2(int o,int fa,double v){
29   ans+=fabs(v);  
30   rep(i,1,n) if(a[o][i] && i!=fa && ancestor(i,id))  dfs2(i,o,1.0*r[o]/r[i]*v);  
31 }
32 int main(){
33   read(n); read(xt); read(yt);  
34   rep(i,1,n) read(x[i]),read(y[i]),read(r[i]);  
35   rep(i,1,n) rep(j,i+1,n) if(dist(i,j)==sqr(r[i]+r[j])) a[i][j]=a[j][i]=1;    
36   rep(i,1,n) if(x[i]==xt && y[i]==yt) id=i;  
37   rep(i,1,n) if(x[i]==0 && y[i]==0)  dfs(i,i);  
38   rep(i,1,n) if(x[i]==0 && y[i]==0)  dfs2(i,i,10000);  
39   cout<<int(ans)<<endl; 
40   return 0;  
41 }
View Code

 

posted @ 2019-02-08 11:44  zjxxcn  阅读(135)  评论(0编辑  收藏  举报