Total circles area: Difference between revisions

The python code is now simpler to port to other languages
(Created the Total_Circles_Area task)
 
(The python code is now simpler to port to other languages)
Line 43:
=={{header|Python}}==
<lang python>from collections import namedtuple
 
standard_data = """\
1.6417233788 1.6121789534 0.0848270516
-1.4944608174 1.2077959613 1.1039549836
0.6110294452 -0.6907087527 0.9089162485
0.3844862411 0.2923344616 0.2375743054
-0.2495892950 -0.3832854473 1.0845181219
1.7813504266 1.6178237031 0.8162655711
-0.1985249206 -0.8343333301 0.0538864941
-1.7011985145 -0.1263820964 0.4776976918
-0.4319462812 1.4104420482 0.7886291537
0.2178372997 -0.9499557344 0.0357871187
-0.6294854565 -1.3078893852 0.7653357688
1.7952608455 0.6281269104 0.2727652452
1.4168575317 1.0683357171 1.1016025378
1.4637371396 0.9463877418 1.1846214562
-0.5263668798 1.7315156631 1.4428514068
-1.2197352481 0.9144146579 1.0727263474
-0.1389358881 0.1092805780 0.7350208828
1.5293954595 0.0030278255 1.2472867347
-0.5258728625 1.3782633069 1.3495508831
-0.1403562064 0.2437382535 1.3804956588
0.8055826339 -0.0482092025 0.3327165165
-0.6311979224 0.7184578971 0.2491045282
1.4685857879 -0.8347049536 1.3670667538
-0.6855727502 1.6465021616 1.0593087096
0.0152957411 0.0638919221 0.9771215985""".splitlines()
 
Circle = namedtuple("Circle", "x y r")
circles = [Circle(*map(float, r.split())) for r in standard_data]
 
circles = [
# compute the bounding box of the circles
Circle( 1.6417233788, 1.6121789534, 0.0848270516),
x_min = min(c.x - c.r for c in circles)
Circle(-1.4944608174, 1.2077959613, 1.1039549836),
x_max = max(c.x + c.r for c in circles)
Circle( 0.6110294452, -0.6907087527, 0.9089162485),
y_min = min(c.y - c.r for c in circles)
Circle( 0.3844862411, 0.2923344616, 0.2375743054),
y_max = max(c.y + c.r for c in circles)
Circle(-0.2495892950, -0.3832854473, 1.0845181219),
Circle( 1.7813504266, 1.6178237031, 0.8162655711),
Circle(-0.1985249206, -0.8343333301, 0.0538864941),
Circle(-1.7011985145, -0.1263820964, 0.4776976918),
Circle(-0.4319462812, 1.4104420482, 0.7886291537),
Circle( 0.2178372997, -0.9499557344, 0.0357871187),
Circle(-0.6294854565, -1.3078893852, 0.7653357688),
Circle( 1.7952608455, 0.6281269104, 0.2727652452),
Circle( 1.4168575317, 1.0683357171, 1.1016025378),
Circle( 1.4637371396, 0.9463877418, 1.1846214562),
Circle(-0.5263668798, 1.7315156631, 1.4428514068),
Circle(-1.2197352481, 0.9144146579, 1.0727263474),
Circle(-0.1389358881, 0.1092805780, 0.7350208828),
Circle( 1.5293954595, 0.0030278255, 1.2472867347),
Circle(-0.5258728625, 1.3782633069, 1.3495508831),
Circle(-0.1403562064, 0.2437382535, 1.3804956588),
Circle( 0.8055826339, -0.0482092025, 0.3327165165),
Circle(-0.6311979224, 0.7184578971, 0.2491045282),
Circle( 1.4685857879, -0.8347049536, 1.3670667538),
Circle(-0.6855727502, 1.6465021616, 1.0593087096),
Circle( 0.0152957411, 0.0638919221, 0.9771215985""".splitlines()]
 
def main():
# compute the bounding box of the circles
x_min = min(c.x - c.r for c in circles)
x_max = max(c.x + c.r for c in circles)
y_min = min(c.y - c.r for c in circles)
y_max = max(c.y + c.r for c in circles)
 
box_side = 500
 
dx = (x_max - x_min) / box_side
box_side = 500
dy = (y_max - y_min) / box_side
 
count = 0
dx = (x_max - x_min) / box_side
dy = (y_max - y_min) / box_side
 
for cr in xrange(box_side):
count = 0
y = y_min + r * dy
for rc in xrange(box_side):
x = x_min + c * dx
for circle in circles:
if (x - circle.x)**2 + (y - circle.y)**2 <= (circle.r ** 2):
count += 1
break
 
print "Approximated area:", count * dx * dy</lang>
for r in xrange(box_side):
y = y_min + r * dy
for c in xrange(box_side):
x = x_min + c * dx
for circle in circles:
if (x - circle.x)**2 + (y - circle.y)**2 <= (circle.r ** 2):
count += 1
break
 
main()</lang>
print "Approximated area:", count * dx * dy</lang>
{{out}}
<pre>Approximated area: 21.561559772</pre>