login
Permutation of nonnegative integers: a(n) = A156552(A005940(1+n)-1).
5

%I #12 Feb 28 2022 07:56:15

%S 0,1,2,3,4,7,8,5,6,17,16,23,64,65,10,9,32,19,128,129,512,67,256,47,24,

%T 4097,20,4099,32768,79,1024,11,18,31,14,29,4096,2049,66,515,514,263,

%U 65536,2053,8388608,524289,16384,87,16777216,2097153,70,524291,17179869184,287,524288,1037,8388610,289,134217728,1071,520,97,38,15,12,257,8192,63

%N Permutation of nonnegative integers: a(n) = A156552(A005940(1+n)-1).

%C Note the indexing: although the domain starts from 1, the range includes also zero.

%H Antti Karttunen, <a href="/A297164/b297164.txt">Table of n, a(n) for n = 1..379</a>

%H <a href="/index/Bi#binary">Index entries for sequences related to binary expansion of n</a>

%H <a href="/index/Per#IntegerPermutation">Index entries for sequences that are permutations of the natural numbers</a>

%H <a href="/index/Pri#prime_indices">Index entries for sequences computed from indices in prime factorization</a>

%F a(n) = A156552(A005940(1+n)-1).

%o (Scheme) (define (A297164 n) (A156552 (+ -1 (A005940 (+ 1 n)))))

%o (PARI)

%o A005940(n) = { my(p=2,t=1); n--; until(!n\=2, if(n%2, t*=p, p=nextprime(p+1))); t };

%o A156552(n) = {my(f = factor(n), p, p2 = 1, res = 0); for(i = 1, #f~, p = 1 << (primepi(f[i, 1]) - 1); res += (p * p2 * (2^(f[i, 2]) - 1)); p2 <<= f[i, 2]); res}; \\ From A156552

%o A297164(n) = A156552(A005940(1+n)-1); \\ _Antti Karttunen_, Feb 27 2022

%Y Inverse: A297163.

%Y Cf. A005940, A156552.

%Y Cf. also A297166, A341520, A351960.

%K nonn

%O 1,3

%A _Antti Karttunen_, Jan 05 2018