Considering a semilinear elliptic equation {-Delta u + lambda u = mu g(x, u) + b(x) in Omega, u = 0 on partial derivative Omega, in a bounded domain Omega subset of R-n with a smooth boundary, we apply a new variational principle introduced in Momeni (2011, 2017) to show the existence of a strong solution, where g can have critical growth. To be more accurate, assuming G(x, .) is the primitive of g(x, .) and G*(x, .) is the Fenchel dual of G(x, .), we shall find a minimum of the functional I[.] defined by I[u] = integral(Omega) mu G*(x, -Delta u + lambda u - b(x)/mu) dx - integral(Omega) mu G(x, u) + b(x) u dx, over a convex set K, consisting of bounded functions in an appropriate Sobolev space. The symmetric nature of the functional I[.], provided by existence of a function G and its Fenchel dual G*, alleviate the difficulty and shorten the process of showing the existence of solutions for problems with supercritical nonlinearity. It also makes it an ideal choice among the other energy functionals including Euler-Lagrange functional. (C) 2017 Elsevier Ltd. All rights reserved.