cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A128247 a(n) = A001783(n)/A038610(n).

This page as a plain text file.
%I A128247 #15 Mar 17 2018 23:10:29
%S A128247 1,1,1,1,2,1,12,1,8,3,1440,1,17280,3,112,45,29030400,1,522547200,3,
%T A128247 12800,945,4828336128000,1,38626689024,4725,512512000,2025,
%U A128247 3796230997278720000,1,113886929918361600000,42525,10436608000,1403325,551910745571328,175
%N A128247 a(n) = A001783(n)/A038610(n).
%H A128247 T. D. Noe, <a href="/A128247/b128247.txt">Table of n, a(n) for n = 1..200</a>
%e A128247 The positive integers which are <= 10 and are coprime to 10 are 1,3,7,9. So a(10) = 1*3*7*9/lcm(1,3,7,9) = 189/63 = 3.
%t A128247 Table[s = Select[Range[1, n], GCD[#, n] == 1 &]; (Times @@ s)/(LCM @@ s), {n, 30}] (* _T. D. Noe_, Oct 04 2012 *)
%o A128247 (Sage)
%o A128247 def Gauss_factorial(N, n): return mul(j for j in (1..N) if gcd(j, n) == 1)
%o A128247 def L(N, n): return lcm([j for j in (1..N) if gcd(j, n) == 1])
%o A128247 def A128247(n): return Gauss_factorial(n, n)/L(n, n)
%o A128247 [A128247(n) for n in (1..34)] # _Peter Luschny_, Oct 02 2012
%Y A128247 Cf. A038610, A001783.
%K A128247 nonn
%O A128247 1,5
%A A128247 _Leroy Quet_, May 03 2007
%E A128247 More terms from _Sean A. Irvine_, Jun 21 2011