%I A120736
%S A120736 1,2,6,8,9,10,12,14,18,22,24,26,30,34,36,38,40,42,46,54,56,58,60,62,66,
%T A120736 70,72,74,78,80,82,84,86,88,90,94,96,102,104,106,108,110,114,118,120,
%U A120736 122,126,128,130,132,134,136,138,142,146,150,152,154,156,158,166,168
%N A120736 A number n is included if it satisfies: the prime p divides n for all 
               p's where p divides d(n), d(n) = number of positive divisors of n.
%H A120736 Leroy Quet, <a href="http://www.prism-of-spirals.net/">Home Page</a> 
               (listed in lieu of email address)
%e A120736 d(26) = 4. 2 is the only prime dividing 4. 2 divides 26, so 26 is in 
               the sequence.
%p A120736 isA120736 := proc(n) local d,p; d := numtheory[tau](n) ; p := 2 ; while 
               p <= n do if ( d mod p ) = 0 then if ( n mod p ) <> 0 then RETURN(false) 
               ; fi ; fi ; p := nextprime(p) ; od ; RETURN(true) ; end: for n from 
               1 to 200 do if isA120736(n) then printf("%d,",n) ; fi ; od ; - R. 
               J. Mathar (mathar(AT)strw.leidenuniv.nl), Aug 17 2006
%Y A120736 Cf. A120737.
%Y A120736 Sequence in context: A003610 A077477 A095879 this_sequence A130099 A047397 
               A043341
%Y A120736 Adjacent sequences: A120733 A120734 A120735 this_sequence A120737 A120738 
               A120739
%K A120736 nonn
%O A120736 1,2
%A A120736 Leroy Quet Jun 29 2006
%E A120736 More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Aug 17 2006