A351225 - OEIS

%I #14 May 02 2022 11:12:49

%S 1,1,1,3,5,13,-1,3,7,21,35,79,13,37,61,135,209,433,107,231,355,729,

%T 1103,2227,601,1225,1849,3723,5597,11221,-23,-17,-11,9,29,91,-1,33,67,

%U 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

%N a(n) = A276086(n) - n, where A276086 is the primorial base exp-function.

%C Sequence contains no zeros because the parity of A276086(n) is opposite to that of n.

%H Antti Karttunen, <a href="/A351225/b351225.txt">Table of n, a(n) for n = 0..11550</a>

%H <a href="/index/Pri#primorialbase">Index entries for sequences related to primorial base</a>

%t 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 *)

%o (PARI)

%o A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };

%o A351225(n) = (A276086(n)-n);

%Y Cf. A276086, A351226 (positions of negative terms), A351227 (of positive terms).

%Y Cf. also A350074.

%K sign,base

%O 0,4

%A _Antti Karttunen_, Feb 05 2022