login
Brilliant numbers (A078972) whose digital sum is also brilliant.
1

%I #12 Nov 24 2020 16:38:43

%S 4,6,9,15,121,253,361,437,451,473,527,671,703,799,851,901,923,979,

%T 1003,1027,1067,1081,1139,1157,1207,1247,1261,1333,1513,1517,1643,

%U 1711,2183,2291,2701,2813,3007,3403,3551,3589,3713,4307,4399,4489,5041,5063,5767

%N Brilliant numbers (A078972) whose digital sum is also brilliant.

%H Robert Israel, <a href="/A085648/b085648.txt">Table of n, a(n) for n = 1..10000</a>

%e 7470269 is a term because 7470269 = 1249*5981 and 7+4+7+0+2+6+9 = 35 = 5*7.

%p isbril:= proc(n) local F;

%p F:= ifactors(n)[2];

%p (nops(F)=1 and F[1][2]=2) or (nops(F)=2 and F[1][2]=1 and F[2][2]=1 and ilog10(F[1][1])=ilog10(F[2][1]))

%p end proc:

%p R:= NULL:

%p for d from 1 to 2 do

%p if d = 1 then P:= [2,3,5,7]

%p else P:= select(isprime, [seq(i,i=10^(d-1)+1..10^d-1,2)])

%p fi;

%p R:= R, op(sort(select(t -> isbril(convert(convert(t,base,10),`+`)), [seq(seq(P[i]*P[j],j=1..i),i=1..nops(P))])))

%p od:

%p R; # _Robert Israel_, Nov 24 2020

%Y Cf. A078972.

%K base,nonn

%O 1,1

%A _Amarnath Murthy_, _Jason Earls_, Jul 11 2003

%E Offset changed by _Robert Israel_, Nov 24 2020