Disciplined convex programming error: Illegal operation: pow_p( {convex}, {-1} )

I think you have to use the geometric mean

since

geom_mean((x;x;y)) = x^2/3*y^1/3 >= s(1)