OFFSET
1,4
COMMENTS
a(n) = a(n-1) iff n is prime. Thus a(1)=a(2)=a(3)=1 is the only triple in this sequence. - Franz Vrabec, Sep 10 2005
a(k) = a(k+1) for k in A006093. - Lekraj Beedassy, Aug 03 2006
Partial products of A048671. - Peter Luschny, Sep 09 2009
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..500
Liam Solus, Simplices for Numeral Systems, arXiv:1706.00480 [math.CO], 2017. Mentions this sequence.
FORMULA
log a(n) = n log n - 2n + O(n/log^4 n). (The error term can be improved. On the Riemann Hypothesis it is O(n^k) for any k > 1/2.) - Charles R Greathouse IV, Oct 16 2012
a(n) = A205957(n), 1 <= n <= 11. - Daniel Forgues, Apr 22 2014
EXAMPLE
a(5) = 2 as 5!/lcm(1..5) = 120/60 = 2.
MAPLE
seq(n!/lcm($1..n), n=1..30);
A025527 := proc(n) option remember; `if`(n < 3, 1, ilcm(op(numtheory[divisors](n) minus{1, n}))*A025527(n-1)) end:
seq(A025527(i), i=1..26); # Peter Luschny, Mar 23 2011
MATHEMATICA
Table[n!/Apply[LCM, Range[n]], {n, 1, 26}] (* Geoffrey Critzer, Jun 17 2013 *)
PROG
(Sage)
def A025527(n) :
if n < 2 : return 1
else :
D = divisors(n); D.pop()
return lcm(D)*A025527(n-1)
[A025527(i) for i in (1..26)] # Peter Luschny, Feb 03 2012
(PARI) a(n)=n!/lcm([2..n]) \\ Charles R Greathouse IV, Mar 06 2014
(GAP) List([1..30], n->Factorial(n)/Lcm([1..n])); # Muniru A Asiru, Apr 02 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Dec 11 1999
STATUS
approved