|
%I #14 Jul 28 2020 11:07:44
%S 1,10,26,41,44,47,62,68,89,100,104,119,125,128,131,134,143,146,155,
%T 161,167,182,188,215,224,233,242,257,260,269,275,293,311,317,323,359,
%U 368,371,377,386,401,404,410,422,431,437,440,449,452,470,494,518,521,527,545
%N Numbers n such that the sum of the digits of n plus the sum of the digits of n^2 plus the sum of the digits of n^3 is prime.
%C Contains A011557. Except for those, all terms == 2 (mod 3). - _Robert Israel_, Jul 28 2020
%H Robert Israel, <a href="/A186279/b186279.txt">Table of n, a(n) for n = 1..10000</a>
%p filter:= n -> isprime(convert(convert(n,base,10),`+`)+convert(convert(n^2,base,10),`+`)+convert(convert(n^3,base,10),`+`)):
%p select(filter,[$1..1000]); # _Robert Israel_, Jul 28 2020
%t Select[Range[1000],PrimeQ[Total[Flatten[IntegerDigits[#^{1,2,3}]]]]&] (* _Harvey P. Dale_, Jul 18 2015 *)
%Y Cf. A011557, A186278.
%K nonn,easy,base
%O 1,2
%A _Harvey P. Dale_, Feb 16 2011
|
