You can use Solve
on an equation (==
instead of =
). It is essentially a quadratic equation, so there are two solutions. Either there is no inverse or there are two, depending whether you restrict the domain appropriately.
eqn = y == (x^(1/t) - a^(1/t))*(x^(1/t) - c*b^(1/t))/(x^(1/t) - a^(1/t))*(x^(1/t) -
c*b^(1/t)) + (x^(1/t) - b^(1/t))*x^(1/t) - c*a^(1/t);
soln = Solve[eqn, x]
(* {{x -> 4^-t (-b^((1/t)) (-1 - 2 c) -
Sqrt[b^(2/t) + 8 a^(1/t) c + 4 b^(2/t) c - 4 b^(2/t) c^2 + 8 y])^t},
{x -> 4^-t (-b^((1/t)) (-1 - 2 c) +
Sqrt[b^(2/t) + 8 a^(1/t) c + 4 b^(2/t) c - 4 b^(2/t) c^2 + 8 y])^t}} *)
Check both solutions:
eqn /. soln // PowerExpand // Expand
(* {True, True} *)