A382661 The unitary Jordan totient function applied to the exponentially odd numbers (A268335).
1, 3, 8, 24, 24, 48, 63, 72, 120, 168, 144, 192, 288, 360, 384, 360, 528, 504, 504, 728, 840, 576, 960, 1023, 960, 864, 1152, 1368, 1080, 1344, 1512, 1680, 1152, 1848, 1584, 2208, 2304, 2808, 2184, 2880, 3024, 2880, 2520, 3480, 3720, 2880, 4032, 2880, 4488, 4224
Offset: 1
Links
- Amiram Eldar, Table of n, a(n) for n = 1..10000
Crossrefs
Programs
-
Mathematica
f[p_, e_] := p^(2*e)-1; uj2[1] = 1; uj2[n_] := Times @@ f @@@ FactorInteger[n]; expOddQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], OddQ]; uj2 /@ Select[Range[100], expOddQ]
-
PARI
uj2(n) = {my(f = factor(n)); prod(i = 1, #f~, f[i, 1]^(2*f[i, 2])-1);} isexpodd(n) = {my(f = factor(n)); for(i=1, #f~, if(!(f[i, 2] % 2), return (0))); 1;} list(lim) = apply(uj2, select(isexpodd, vector(lim, i, i)));