A351251 - OEIS

%I #10 Feb 06 2022 19:05:50

%S 1,2,3,2,9,18,5,10,15,10,9,90,25,50,75,10,225,450,125,250,75,250,1125,

%T 2250,625,50,1875,1250,5625,11250,7,14,21,14,63,18,35,70,105,70,63,

%U 630,25,350,525,70,1575,3150,875,250,105,1750,7875,15750,4375,1750,1875,8750,39375,78750,49,98,147,14,441,882,245,490

%N Denominator of n / A276086(n).

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

%F a(n) = A276086(n) / gcd(n,A276086(n)) = A276086(n) / A324198(n).

%F a(n) = A276086(A351253(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++]; Denominator[#/m]] &, 68, 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 A351251(n) = denominator(n/A276086(n));

%Y Cf. A276086, A324198, A351250 (numerators), A351253.

%Y Cf. also A351231.

%K nonn,frac,easy

%O 0,2

%A _Antti Karttunen_, Feb 05 2022