python3 <<'PY'
import math
a=0.25; b=1.0/12
def G(t,y):
    y1=3*t-2.0/3
    al=2*a+2*b-y1-t; mu1=a+b-y1; mu2=b+a-t; nu1=a+b
    r1=mu1/al; r2=mu2/al; s1=nu1/al; s3=1-s1
    A=al*(math.sqrt(max(0,a*r1*s1))+math.sqrt(max(0,b*r1*s3))+math.sqrt(max(0,b*r2*s1))+math.sqrt(max(0,a*r2*s3)))
    beta=a+b-y; B=math.sqrt(beta)*(b+math.sqrt(max(0,a*(a-y))))
    gam=y1+t+y; C=math.sqrt(gam)*(math.sqrt(max(0,b*y1))+math.sqrt(max(0,a*t))+math.sqrt(max(0,a*y)))
    return A+B+C
def opty(t):
    lo=0.0; hi=a; gr=(math.sqrt(5)-1)/2
    for _ in range(200):
      c=hi-gr*(hi-lo); d=lo+gr*(hi-lo)
      if G(t,c)<G(t,d): lo=c
      else: hi=d
    return (lo+hi)/2
lo=2.0/9; hi=0.25
for _ in range(200):
  m1=lo+(hi-lo)/3; m2=hi-(hi-lo)/3
  if G(m1,opty(m1))<G(m2,opty(m2)): lo=m1
  else: hi=m2
t=(lo+hi)/2; y=opty(t)
val=G(t,y)
print("t",repr(t))
print("y",repr(y))
print("y1",repr(3*t-2/3))
print("aff",repr(val))
print("cost",repr(1-val))
print("floor100",math.floor(100*(1-val)))
print("endpoint t",G(2/9,opty(2/9)),G(0.25,opty(0.25)))
PY