|
%I #15 Nov 21 2013 12:48:47
%S 15,24,27,33,48,51,55,64,71,75,81,99,105,108,111,119,120,123,126,132,
%T 141,147,150,156,160,162,171,175,177,189,198,199,204,208,215,219,222,
%U 224,249,252,258,261,263,264,267,270,272,280,285,291,294,300,304,335
%N Numbers n such that n-th heptagonal number is 6-almost prime.
%C Hep(2) = 7 is the only prime heptagonal number.
%H Vincenzo Librandi, <a href="/A114558/b114558.txt">Table of n, a(n) for n = 1..1000</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/HeptagonalNumber.html">Heptagonal Number.</a>
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/AlmostPrime.html">Almost Prime.</a>
%F n such that Hep(n) = n*(5*n-3)/2 is 6-almost prime. n such that A000566(n) is an element of A046306. n such that A001222(A000566(n)) = 6. n such that A001222(n*(5*n-3)/2) = 6.
%e a(1) = 15 because Hep(15) = 15*(5*15-3)/2 = 540 = 2^2 * 3^3 * 5 is 6-almost prime.
%e a(2) = 24 because Hep(24) = 24*(5*24-3)/2 = 1404 = 2^2 * 3^3 * 13.
%e a(7) = 55 because Hep(55) = 55*(5*55-3)/2 = 7480 = 2^3 * 5 * 11 * 17 is 6-almost prime [also 7480 = Hep(55) = Hep(Hep(5)) is an iterated heptagonal number].
%e a(11) = 81 because Hep(81) = 81*(5*81-3)/2 = 16281 = 3^5 * 67 [also 16281 = Hep(81) = Hep(Hep(6)) is an iterated heptagonal number].
%e a(24) = 156 because Hep(156) = 156*(5*156-3)/2 = 60606 = 2 * 3^2 * 7 * 13 * 37 is 6-almost prime (and a palindrome).
%e a(30) = 189 because Hep(189) = 189*(5*189-3)/2 = 89019 = 3^4 * 7 * 157 is 6-almost prime [also 89019 = Hep(189) = Hep(Hep(9)) is an iterated heptagonal number].
%t Select[Range[400],Total[Transpose[FactorInteger[# (5#-3)/2]][[2]]]==6&] (* _Harvey P. Dale_, May 15 2011 *)
%Y Cf. A000040, A000566, A001222, A001358, A046306.
%K nonn
%O 1,1
%A _Jonathan Vos Post_, Feb 15 2006
%E More terms from Harvey P. Dale, May 15 2011.
|
