A167333 Totally multiplicative sequence with a(p) = 2*(5p-1) = 10p-2 for prime p.
1, 18, 28, 324, 48, 504, 68, 5832, 784, 864, 108, 9072, 128, 1224, 1344, 104976, 168, 14112, 188, 15552, 1904, 1944, 228, 163296, 2304, 2304, 21952, 22032, 288, 24192, 308, 1889568, 3024, 3024, 3264, 254016, 368, 3384, 3584, 279936, 408, 34272, 428
Offset: 1
Links
- G. C. Greubel, Table of n, a(n) for n = 1..1000
Programs
-
Mathematica
a[1] = 1; a[n_] := (fi = FactorInteger[n]; Times @@ ((5*fi[[All, 1]] - 1)^fi[[All, 2]])); Table[a[n]*2^PrimeOmega[n], {n, 1, 100}] (* G. C. Greubel, Jun 06 2016 *) f[p_, e_] := (10*p-2)^e; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Oct 18 2023 *) -
PARI
a(n) = {my(f=factor(n)); for (k=1, #f~, f[k,1] = 10*f[k,1]-2;); factorback(f);} \\ Michel Marcus, Jun 06 2016