Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import sympy
- from tqdm import tqdm_notebook
- r = 5
- circle = sympy.geometry.Circle((0,0), r)
- ans = 0
- for a1 in tqdm_notebook(range(r)):
- for a2 in range(r):
- for b1 in range(r):
- for b2 in range(r):
- for c1 in range(r):
- for c2 in range(r):
- a = (a1, a2)
- b = (b1, b2)
- c = (c1, c2)
- if circle.encloses_point(a) and circle.encloses_point(b) and circle.encloses_point(c):
- try:
- t = sympy.geometry.Triangle(a, b, c)
- if abs(t.area) == 1/2:
- ans = max(float(t.perimeter), ans)
- print(a, b, c, float(t.perimeter))
- except:
- pass
- ans
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement