[go: up one dir, main page]

login
A328626
Inverse permutation to A328625.
9
0, 1, 2, 5, 4, 3, 6, 7, 20, 17, 16, 21, 12, 13, 8, 29, 28, 9, 18, 19, 26, 11, 10, 27, 24, 25, 14, 23, 22, 15, 30, 31, 32, 35, 34, 33, 126, 127, 80, 167, 166, 81, 162, 163, 128, 119, 118, 129, 78, 79, 116, 131, 130, 117, 114, 115, 164, 83, 82, 165, 60, 61, 62, 65, 64, 63, 36, 37, 140, 107, 106, 141, 102, 103, 38, 209, 208, 39, 138, 139, 206, 41, 40
OFFSET
0,3
FORMULA
a(n) = A276085(A328627(n)).
For all n, A328620(a(n)) = A328620(n).
PROG
(PARI)
A002110(n) = prod(i=1, n, prime(i));
A276085(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*A002110(primepi(f[k, 1])-1)); };
A328627(n) = { my(m=1, p=2, d=0); while(n, d = lift(Mod(n, p)/(d+1)); m *= (p^d); n = n\p; p = nextprime(1+p)); (m); };
CROSSREFS
Cf. A328625 (inverse permutation).
Cf. also A289234, A328623.
Sequence in context: A328625 A060125 A265361 * A265357 A265358 A361832
KEYWORD
nonn
AUTHOR
Antti Karttunen, Oct 25 2019
STATUS
approved