|
%I #8 Mar 02 2015 16:08:49
%S 2,3,5,6,25,30,36,38,39,42,60,78,84,90,106,114,120,122,126,130,150,
%T 152,156,171,178,180,183,186,187,194,198,216,217,218,219,221,222,228,
%U 230,240,244,252,255,258,259,260,262,264,270,287,294,297,299,300,303,321
%N Numbers n such that sum of digits in base 6 is a divisor of sum of prime divisors (A008472).
%p A076384 := proc(n) local i,j,t,t1, sod, sopd; t := NULL; for i from 2 to n do t1 := i; sod := 0; while t1 <> 0 do sod := sod + (t1 mod 6); t1 := floor(t1/6); od; sopd := 0; j := 1; while ithprime(j) <= i do if i mod ithprime(j) = 0 then sopd := sopd+ithprime(j); fi; j := j+1; od; if sopd mod sod = 0 then t := t,i; fi; od; t; end;
%t Select[Range[2,350],Divisible[Total[Transpose[FactorInteger[#]][[1]]], Total[ IntegerDigits[#,6]]]&] (* _Harvey P. Dale_, May 26 2013 *)
%Y Cf. A075657, A076380 - A076387.
%K nonn,base
%O 1,1
%A _Floor van Lamoen_, Oct 08 2002
|
