cp's OEIS Frontend

A351251 Denominator of n / A276086(n).

%I A351251 #10 Feb 06 2022 19:05:50
%S A351251 1,2,3,2,9,18,5,10,15,10,9,90,25,50,75,10,225,450,125,250,75,250,1125,
%T A351251 2250,625,50,1875,1250,5625,11250,7,14,21,14,63,18,35,70,105,70,63,
%U A351251 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 A351251 Denominator of n / A276086(n).
%H A351251 Antti Karttunen, <a href="/A351251/b351251.txt">Table of n, a(n) for n = 0..11550</a>
%F A351251 a(n) = A276086(n) / gcd(n,A276086(n)) = A276086(n) / A324198(n).
%F A351251 a(n) = A276086(A351253(n)).
%t A351251 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 A351251 (PARI)
%o A351251 A276086(n) = { my(m=1, p=2); while(n, m *= (p^(n%p)); n = n\p; p = nextprime(1+p)); (m); };
%o A351251 A351251(n) = denominator(n/A276086(n));
%Y A351251 Cf. A276086, A324198, A351250 (numerators), A351253.
%Y A351251 Cf. also A351231.
%K A351251 nonn,frac,easy
%O A351251 0,2
%A A351251 _Antti Karttunen_, Feb 05 2022