|
%I #13 Feb 06 2022 18:17:11
%S 0,1,2,1,4,5,6,7,8,3,2,11,12,13,14,1,16,17,18,19,4,7,22,23,24,1,26,9,
%T 28,29,30,31,32,11,34,5,36,37,38,13,8,41,6,43,44,3,46,47,48,7,2,17,52,
%U 53,54,11,8,19,58,59,60,61,62,3,64,65,66,67,68,23,2,71,72,73,74,1,76,11,78,79,16,27,82,83,12,17
%N Numerator of n / A276086(n).
%H Antti Karttunen, <a href="/A351250/b351250.txt">Table of n, a(n) for n = 0..30030</a>
%F a(n) = n / A324198(n) = n / gcd(n, A276086(n)).
%F a(n) = n / gcd(n, A328386(n)).
%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++]; Numerator[#/m]] &, 86, 0] (* _Michael De Vlieger_, Feb 06 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 A351250(n) = numerator(n/A276086(n));
%Y Cf. A276086, A324198, A328386, A328387 (positions of ones), A351251 (denominators).
%Y Cf. A324583 (the positions of fixed points after the zero).
%Y Cf. also A351230.
%K nonn,frac,easy
%O 0,3
%A _Antti Karttunen_, Feb 05 2022
|
