|
%I #15 Nov 21 2021 07:34:28
%S 1,4,6,16,8,24,9,64,36,32,10,96,12,36,48,256,14,144,15,128,54,40,16,
%T 384,64,48,216,144,18,192,20,1024,60,56,72,576,21,60,72,512,22,216,24,
%U 160,288,64,25,1536,81,256,84,192,26,864,80,576,90,72,27,768,28,80,324
%N In the prime factorization of n replace the k-th prime with the k-th composite number, k > 0.
%H Harry J. Smith, <a href="/A066260/b066260.txt">Table of n, a(n) for n = 1..1000</a>
%H <a href="/index/Di#divseq">Index to divisibility sequences</a>
%F Completely multiplicative with a(p) = A002808(A049084(p)), p prime.
%e a(10) = a(2*5) = a(prime(1)*prime(3)) = a(prime(1))*a(prime(3)) = comp(1)*comp(3) = 4 * 8 = 32.
%t nmax = 100;
%t compos = Select[Range[FindRoot[n == nmax + PrimePi[n] + 1,
%t {n, nmax, 2 nmax}][[1, 2]] // Floor], CompositeQ];
%t a[n_] := If[n == 1, 1, Product[{p, e} = pe; compos[[PrimePi[p]]]^e,
%t {pe, FactorInteger[n]}]];
%t Array[a, nmax] (* _Jean-François Alcover_, Nov 21 2021 *)
%o (PARI) Composite(n) = local(k); k=n + primepi(n) + 1; while (k != n + primepi(k) + 1, k = n + primepi(k) + 1); return(k)
%o for (n=1, 1000, f=factor(n); a=1; for (i=1, matsize(f)[1], a*=Composite(primepi(f[i, 1]))^f[i, 2]); write("b066260.txt", n, " ", a) ) \\ _Harry J. Smith_, Feb 07 2010
%Y Cf. A002808, A049084, A066261, A066262.
%K nonn,mult
%O 1,2
%A _Reinhard Zumkeller_, Dec 10 2001
|
