|
%I #8 Mar 30 2012 18:40:41
%S 29,53,97,113,161,159,145,269,244,232,231,247,261,373,399,386,328,350,
%T 375,371,395,547,559,572,537,541,577,635,679,663,607,621,687,673,658,
%U 769,871,853,839,856,832,881,947,939,1003,1007,955,915,907,889,941,989
%N Prime(4almostprime(n))-4almostprime(prime(n)). Commutator [A000040,A014613] at n.
%H Robert G. Wilson v, <a href="/A124284/b124284.txt">Table of n, a(n) for n=1..1000</a>
%F a(n) = prime(4almostprime(n)) - 4almostprime(prime(n)) = A000040(A014613(n)) -A014613(A000040(n)).
%e a(1) = prime(4almostprime(1)) - 4almostprime(prime(1)) = 53 - 24 = 29.
%e a(2) = prime(4almostprime(2)) - 4almostprime(prime(2)) = 89 - 36 = 53.
%e a(3) = prime(4almostprime(3)) - 4almostprime(prime(3)) = 151 - 54 = 97.
%e It is mere coincidence that the first 4 values are all primes.
%t FourAlmostPrimePi[n_] := Sum[PrimePi[n/(Prime@i*Prime@j*Prime@k)] - k + 1, {i, PrimePi[n^(1/4)]}, {j, i, PrimePi[(n/Prime@i)^(1/3)]}, {k, j, PrimePi@ Sqrt[n/(Prime@i*Prime@j)]}];
%t FourAlmostPrime[n_] := Block[{e = Floor[Log[2, n] + 1], a, b}, a = 2^e; Do[b = 2^p; While[ FourAlmostPrimePi@a < n, a = a + b]; a = a - b/2, {p, e, 0, -1}]; a + b/2];
%t Table[ Prime@ FourAlmostPrime@ n - FourAlmostPrime@ Prime@ n, {n, 52}]
%Y Cf. Primes indexed by 4-almost primes = A124282. 4-almost primes indexed by primes = A124283. Primes indexed by 3-almost primes = A124268. 3-almost primes indexed by primes = A124269. prime(3almostprime(n)) - 3almostprime(prime(n)) = A124270. See also A106349 Primes indexed by semiprimes. See also A106350 Semiprimes indexed by primes. See also A122824 Prime(semiprime(n)) - semiprime(prime(n)).
%Y Cf. A000040, A014613, A114414, A122824, A124269, A124270, A124282, A124283.
%K easy,nonn
%O 1,1
%A _Jonathan Vos Post_, Oct 24 2006
%E More terms from _Robert G. Wilson v_, Aug 31 2007
|
