Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- from decimal import Decimal
- from math import sqrt
- q = int(input())
- okr = list()
- cil = list()
- for i in range(q):
- a = list(eval(input()))
- a[0],a[1],a[2] = a[0]-4,a[1]-4,a[2]*4
- a = list(map(Decimal,a))
- cil.append(a)#1=5
- r = Decimal(2*(Decimal(50)/Decimal((50-a[2]))))
- d = Decimal(sqrt(a[0]**2 + a[1]**2))
- y = Decimal(d * (Decimal(a[2])/Decimal(50-a[2])))
- try:
- k = Decimal(abs(a[1]/a[0]))
- x1 = Decimal(sqrt(y**2/(k**2+1)))
- y1 = k*x1
- except:
- y1 = y
- x1 = 0
- if a[0]==0 and a[1] == 0:
- okr.append([(abs(a[0]) +abs(x1))*a[0]/(abs(a[0]+1)),(abs(a[1]) +abs(y1))*a[1]/(abs(a[1]+1)),r,a[0],a[1]])
- elif a[0] == 0:
- okr.append([(abs(a[0]) +abs(x1))*a[0]/(abs(a[0]+1)),(abs(a[1]) +abs(y1))*a[1]/(abs(a[1])),r,a[0],a[1]])
- elif a[1] == 0:
- okr.append([(abs(a[0]) +abs(x1))*a[0]/(abs(a[0])),(abs(a[1]) +abs(y1))*a[1]/(abs(a[1]+1)),r,a[0],a[1]])
- else:
- okr.append([(abs(a[0]) +abs(x1))*a[0]/(abs(a[0])),(abs(a[1]) +abs(y1))*a[1]/(abs(a[1])),r,a[0],a[1]])
- res = 0
- for i in range(q):
- for c in range(q):
- if okr[i][3] == cil[c][0] and okr[i][4] == cil[c][1]:
- continue
- d = sqrt((okr[i][0] - cil[c][0])**2 + (okr[i][1] - cil[c][1])**2)
- if d*5<= okr[i][2] + 2:
- pass
- print(d*5, okr[i][2] + 2)
- print(res)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement