%I #12 Jul 10 2014 03:40:47
%S 4,9,10,22,121,169,178,187,202,259,295,301,358,394,466,493,529,538,
%T 565,583,655,718,745,763,781,799,817,835,862,871,889,898,934,943,961,
%U 979,1003,1111,1159,1177,1186,1195,1267,1285,1294,1339,1357,1366,1393,1438,1465
%N Semiprimes whose digit sum is a perfect square.
%C Subsequence of A028839.
%H K. D. Bajpai, <a href="/A244863/b244863.txt">Table of n, a(n) for n = 1..10000</a>
%e 178 is in the sequence because 178 = 2*89 (semiprime) and 1+7+8 = 16 (square).
%e 187 is in the sequence because 187 = 11*17 (semiprime) and 1+8+7 = 16 (square).
%p select(n -> numtheory:-bigomega(n)=2 and issqr(convert(convert(n,base,10),`+`)),
%p [$1..3000]); # _Robert Israel_, Jul 09 2014
%t Select[Range[3000], PrimeOmega[#] == 2 && IntegerQ[Sqrt[Apply [Plus, IntegerDigits[#]]]] &]
%Y Cf. A007953, A028839, A048519.
%Y Cf. A107288 (Primes whose digit sum is square).
%K nonn,base
%O 1,1
%A _K. D. Bajpai_, Jul 07 2014
|