cp's OEIS Frontend

A308317 Multiplicative with a(prime(k)^e) = A005117(e+1)^(2^(k-1)).

%I A308317 #11 May 20 2019 14:15:31
%S A308317 1,2,4,3,16,8,256,5,9,32,65536,12,4294967296,512,64,6,
%T A308317 18446744073709551616,18,340282366920938463463374607431768211456,48,
%U A308317 1024,131072,115792089237316195423570985008687907853269984665640564039457584007913129639936,20,81,8589934592,25
%N A308317 Multiplicative with a(prime(k)^e) = A005117(e+1)^(2^(k-1)).
%C A308317 This sequence is a permutation of the natural numbers (with inverse A308328).
%C A308317 The property of being a bijection is easily deduced from the Fermi-Dirac representation of a number.
%C A308317 The first known fixed points are: 1, 2, 9, 12, 18, 3584, 32256.
%H A308317 OEIS Wiki, <a href="/wiki/%22Fermi-Dirac_representation%22_of_n">"Fermi-Dirac representation" of n</a>
%H A308317 <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>
%F A308317 a(2^e) = A005117(e+1) for any e >= 0.
%F A308317 a(prime(k)) = A001146(k-1) for any k > 0.
%F A308317 A000120(A267116(a(n))) = A001221(n).
%e A308317 a(3) = a(prime(2)) = A001146(1) = 2^(2^1) = 4.
%e A308317 a(2^5) = A005117(6) = 7.
%e A308317 a(96) = a(2^5 * 3) = a(2^5) * a(3) = 7 * 4 = 28.
%o A308317 (PARI) A005117(n) = for (k=1, oo, if (issquarefree(k), if (n--==0, return (k))))
%o A308317 a(n) = my (f=factor(n)); prod (i=1, #f~, A005117(1+f[i,2])^(2^(primepi(f[i,1])-1)))
%Y A308317 Cf. A000120, A001146, A001221, A005117, A267116, A308328 (inverse), A318363.
%K A308317 nonn,mult
%O A308317 1,2
%A A308317 _Rémy Sigrist_, May 19 2019