|
%I #16 Jan 27 2019 21:07:41
%S 1,7,10,36,54,62,600,658,692,823,888,907,1200,1389,1898,2100,2171,
%T 2824,3177,3344,3711,3938,4680,5029,5771,6737,6787,6813,7517,9600,
%U 9955,10436,11356,12169,12867,12967,13105,13588,13879,14707,14867,15400,15513
%N Numbers such that the sum of the factorials of the digits of the cube is a square.
%H Vincenzo Librandi, <a href="/A126076/b126076.txt">Table of n, a(n) for n = 1..949</a>
%e a(4)=36, 36^3 = 46656, (4! + 6! + 6! + 5! + 6!) = 2304 = 48^2.
%p a:=proc(n) local nn: nn:=convert(n^3,base,10): if type(sqrt(add(nn[j]!,j=1..nops(nn))),integer)=true then n else fi end: seq(a(n),n=1..20000); # _Emeric Deutsch_, Mar 06 2007
%t Select[Range[20000], IntegerQ[Sqrt[Total[IntegerDigits[#^3]!]]]&] (* _Vincenzo Librandi_, Mar 29 2014 *)
%K nonn,base
%O 1,2
%A _Yalcin Aktar_, Mar 02 2007
%E More terms from _Emeric Deutsch_, Mar 06 2007
%E Name changed by _Franklin T. Adams-Watters_, Mar 29 2014
|
