We study the equation -Delta u = |x|(alpha)u(p alpha & lowast;+epsilon)+lambda(epsilon)|x|(beta)u in Omega, under the condition u = 0 on partial derivative Omega, where Omega is a smooth bounded domain in R-N, N >= 5, which is symmetric respect to x(1), x(2), 2026;, x(N) and contains the origin, alpha > -2, -2 < beta < N - 4, p(alpha)(& lowast; )= N+2 alpha+2/N-2, epsilon > 0 is a small parameter and lambda(epsilon) > 0 depends on epsilon, with lambda(epsilon) -> 0 as epsilon -> 0. Our main focus lies in finding positive solutions that take the form of a tower of bubbles of order alpha, exhibiting concentration at the origin as epsilon tends to zero. Furthermore, we extend our study to the equation -Delta u = |x|(alpha)up(alpha & lowast;-epsilon)-lambda(epsilon)|x|beta u in R-N\B-1, where B-1 is the unit ball centered at the origin, under Dirichlet zero boundary condition and an additional vanishing condition at infinity. In this context, we discover positive solutions that take the form of a tower of bubbles of order alpha, progressively flattening as epsilon tends to zero.