login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A225652 a(n) = (1/n) * lcm(n,p1,p2,...,pk) for that partition of n which maximizes this value among all partitions [p1,p2,...,pk] of n. 6
1, 1, 2, 3, 6, 5, 12, 15, 20, 21, 30, 35, 60, 45, 56, 105, 210, 77, 420, 99, 220, 315, 840, 385, 924, 1155, 1540, 585, 2520, 364, 4620, 3465, 3640, 4620, 3432, 5005, 13860, 8190, 6160, 9009, 30030, 4290, 60060, 9945, 12376, 45045, 120120, 17017, 51480, 36036 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,3

LINKS

Alois P. Heinz, Table of n, a(n) for n = 1..100

FORMULA

a(n) = A225646(n)/n.

MAPLE

b:= proc(n, i) option remember; `if`(n=0, {1},

      `if`(i<1, {}, {seq(map(x->ilcm(x, `if`(j=0, 1, i)),

       b(n-i*j, i-1))[], j=0..n/i)}))

    end:

a:= n-> max(seq(ilcm(n, i), i=b(n$2)))/n:

seq(a(n), n=1..50);  # Alois P. Heinz, May 25 2013

MATHEMATICA

b[n_, i_] := b[n, i] = If[n==0, {1}, If[i<1, {}, Table[Map[Function[{x}, LCM[x, If[j==0, 1, i]]], b[n-i*j, i-1]], {j, 0, n/i}]]]; a[n_] := Max[Table[LCM[n, i], {i, b[n, n]}]]/n; Table[Print["a(", n, ") = ", an = a[n]]; an, {n, 1, 50}] (* Jean-Fran├žois Alcover, Jul 29 2015, after Alois P. Heinz *)

PROG

(Scheme): (define (A225652 n) (/ (A225646 n) (max 1 n)))

CROSSREFS

Cf. A225656, A225636, A225637, A225558.

Sequence in context: A253568 A053570 A129647 * A136183 A100211 A071257

Adjacent sequences:  A225649 A225650 A225651 * A225653 A225654 A225655

KEYWORD

nonn

AUTHOR

Antti Karttunen, May 15 2013

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 16 04:09 EDT 2019. Contains 325064 sequences. (Running on oeis4.)