A351225 a(n) = A276086(n) - n, where A276086 is the primorial base exp-function.
1, 1, 1, 3, 5, 13, -1, 3, 7, 21, 35, 79, 13, 37, 61, 135, 209, 433, 107, 231, 355, 729, 1103, 2227, 601, 1225, 1849, 3723, 5597, 11221, -23, -17, -11, 9, 29, 91, -1, 33, 67, 171, 275, 589, 133, 307, 481, 1005, 1529, 3103, 827, 1701, 2575, 5199, 7823, 15697, 4321, 8695, 13069, 26193, 39317, 78691, -11, 37, 85, 231
Offset: 0
Links
Crossrefs
Programs
-
Mathematica
Array[Block[{i, m, n = #, p}, m = i = 1; While[n > 0, p = Prime[i]; m *= p^Mod[n, p]; n = Quotient[n, p]; i++]; m - #] &, 64, 0] (* Michael De Vlieger, Feb 05 2022 *)
-
PARI
A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); }; A351225(n) = (A276086(n)-n);
Comments