[已解决]设计一个函数，半径为R的大圆内放置半径为1的小圆，计算小圆的数量

杨朱2

条件：1，大圆的半径可以是任意正数。不论浮点和整形
          2，小圆与大圆最多相切，不可以交叉、干涉
          3，小圆与小圆同样最多只能相切，不可以交叉，干涉
          4 ，小圆半径为1
          5，要求在任意大于0的大圆半径r的情况下有效，与实际相符（哪怕在1<r<4.1的范围内）
难点：问题没有想象中的简单，关键在于堆叠排列小圆的方法，以及边界的限定问题。

最佳答案

月排行榜 / 总排行榜

syf040916

2019-3-31 23:00:24

1. 
   
2. def Dcd(D,d):
   
3.         sum=1
   
4.         n=int((D/2-d/2)//d)
   
5.         sum_res=0
   
6.         lift=d*3**0.5/2
   
7.         h=n*d*3**0.5/2
   
8.         for i in range(n):sum+=6*(i+1)
   
9.         if n%2==0:
   
10.                 for j in range(int(n/2)):
    
11.                         r=(2*j+1)*d/2
    
12.                         if (r**2+(h+lift)**2)**0.5+d/2<=D/2:sum_res+=2
    
13.                         else:break
    
14.         else:
    
15.                 for k in range(int((n+1)/2)):
    
16.                         r=k*d
    
17.                         if r==0:
    
18.                                 if (r**2+(h+lift)**2)**0.5+d/2<=D/2:
    
19.                                         sum_res+=1
    
20.                                 else:break
    
21.                         else:
    
22.                                 if (r**2+(h+lift)**2)**0.5+d/2<=D/2:
    
23.                                         sum_res+=2
    
24.                                 else:break
    
25.         return sum+sum_res*6

复制代码

跳转到最佳答案楼层

syf040916

1. 
   
2. def Dcd(D,d):
   
3.         sum=1
   
4.         n=int((D/2-d/2)//d)
   
5.         sum_res=0
   
6.         lift=d*3**0.5/2
   
7.         h=n*d*3**0.5/2
   
8.         for i in range(n):sum+=6*(i+1)
   
9.         if n%2==0:
   
10.                 for j in range(int(n/2)):
    
11.                         r=(2*j+1)*d/2
    
12.                         if (r**2+(h+lift)**2)**0.5+d/2<=D/2:sum_res+=2
    
13.                         else:break
    
14.         else:
    
15.                 for k in range(int((n+1)/2)):
    
16.                         r=k*d
    
17.                         if r==0:
    
18.                                 if (r**2+(h+lift)**2)**0.5+d/2<=D/2:
    
19.                                         sum_res+=1
    
20.                                 else:break
    
21.                         else:
    
22.                                 if (r**2+(h+lift)**2)**0.5+d/2<=D/2:
    
23.                                         sum_res+=2
    
24.                                 else:break
    
25.         return sum+sum_res*6

复制代码

Stubborn

syf040916 发表于 2019-3-31 23:00

6666666

syf040916

重新更新了一下算法