def kanta(r,a,b):
return a*b/(2*r) - 2*b*r/a + a
def kylki(r,a,b):
return (a^2+4*r^2)*(a^2*b+2*a^2*r-4*b*r^2) / (4*a*r*(a^2-4*r^2))
r = 5
for a in range(2*r+1, 100):
for b in range(1, 100):
v1 = kanta(r,a,b)
v2 = kylki(r,a,b)
if all(int(v)==v for v in [v1,v2]):
print (a,b)