Во-первых, вы упустили число

;-)
Ну а далее, пока не пришёл уважаемый
gris с изящными теоретическими построениями, попробуем тупо в лоб.
Код:
goodNumberQ[n_] := AllTrue[{n, Total[IntegerDigits[n]], Times @@ IntegerDigits[n], DivisorSigma[0, n]}, IntegerQ[Sqrt[#]] &]
Select[Range[10^6], goodNumberQ]
Среди первого миллиона таких чисел нашлось 55. Начинается последовательность так:

И по понятным причинам их бесконечно много.