|
|
A000792
|
|
a(n) = max{(n - i)*a(i) : i < n}; a(0) = 1.
(Formerly M0568 N0205)
|
|
73
|
|
|
1, 1, 2, 3, 4, 6, 9, 12, 18, 27, 36, 54, 81, 108, 162, 243, 324, 486, 729, 972, 1458, 2187, 2916, 4374, 6561, 8748, 13122, 19683, 26244, 39366, 59049, 78732, 118098, 177147, 236196, 354294, 531441, 708588, 1062882, 1594323, 2125764, 3188646, 4782969, 6377292
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
COMMENTS
|
Numbers of the form 3^k, 2*3^k, 4*3^k with a(0) = 1 prepended.
If a set of positive numbers has sum n, this is the largest value of their product.
In other words, maximum of products of partitions of n: maximal value of Product k_i for any way of writing n = Sum k_i. To find the answer, take as many of the k_i's as possible to be 3 and then use one or two 2's (see formula lines below).
a(n) is also the maximal size of an Abelian subgroup of the symmetric group S_n. For example when n = 6, one of the Abelian subgroups with maximal size is the subgroup generated by (123) and (456), which has order 9. [Bercov and Moser] - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 19 2001
Also the maximum number of maximal cliques possible in a graph with n vertices (cf. Capobianco and Molluzzo). - Felix Goldberg (felixg(AT)tx.technion.ac.il), Jul 15 2001. [Corrected by Jim Nastos and Tanya Khovanova, Mar 11 2009]
Every triple of alternate terms {3*k, 3*k+2, 3*k+4} in the sequence forms a G.P. with first term 3^k and common ratio 2. - Lekraj Beedassy, Mar 28 2002
For n > 4, a(n) is the least multiple m of 3 not divisible by 8, for which omega(m) <= 2 and sopfr(m) = n. - Lekraj Beedassy, Apr 24 2003
Maximal number of divisors that are possible amongst numbers m such that A080256(m) = n. - Lekraj Beedassy, Oct 13 2003
Or, numbers of form 2^p*3^q with p <= 2, q >= 0 and 2p + 3q = n. Largest number obtained using only the operations +,* and () on the parts 1 and 2 of any partition of n into these two summands where the former exceeds the latter. - Lekraj Beedassy, Jan 07 2005
a(n) corresponds also to the ultimate occurrence of n in A001414 and thus stands for the highest number m such that sopfr(m) = n, for n >= 2. - Lekraj Beedassy, Apr 29 2002
a(n) for n >= 1 is a paradigm shift sequence with procedural length p = 0, in the sense of A193455. - Jonathan T. Rowell, Jul 26 2011
For n >= 2, a(n) is the largest number whose prime divisors (with multiplicity) add to n, whereas the smallest such number (resp. smallest composite number) is A056240(n) (resp. A288814(n)). - David James Sycamore, Nov 23 2017
For n >= 3, a(n+1) = a(n)*(1 + 1/s), where s is the smallest prime factor of a(n). - David James Sycamore, Apr 10 2018
|
|
REFERENCES
|
B. R. Barwell, Cutting String and Arranging Counters, J. Rec. Math., 4 (1971), 164-168.
B. R. Barwell, Journal of Recreational Mathematics, "Maximum Product": Solution to Prob. 2004;25(4) 1993 Baywood NY.
M. Capobianco and J. C. Molluzzo, Examples and Counterexamples in Graph Theory, p. 207. North-Holland: 1978.
S. L. Greitzer, International Mathematical Olympiads 1959-1977, Prob. 1976/4 pp. 18;182-3 NML vol. 27 MAA 1978
J. L. Gross and J. Yellen, eds., Handbook of Graph Theory, CRC Press, 2004; p. 396.
P. R. Halmos, Problems for Mathematicians Young and Old, Math. Assoc. Amer., 1991, pp. 30-31 and 188.
L. C. Larson, Problem-Solving Through Problems. Problem 1.1.4 pp. 7. Springer-Verlag 1983.
D. J. Newman, A Problem Seminar. Problem 15 pp. 5;15. Springer-Verlag 1982.
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1 + x + 2*x^2 + x^4)/(1 - 3*x^3). - Simon Plouffe in his 1992 dissertation.
a(3n) = 3^n; a(3*n+1) = 4*3^(n-1) for n > 0; a(3*n+2) = 2*3^n.
a(n) = n if n <= 2, otherwise a(n-1) + Max{gcd(a(i), a(j)) | 0 < i < j < n}. - Reinhard Zumkeller, Feb 08 2002
A007600(a(n)) = n; Andrew Chi-Chih Yao attributes this observation to D. E. Muller. - Vincent Vatter, Apr 24 2006
a(n) = 3^(n - 2 - 2*floor((n - 1)/3))*2^(2 - (n - 1) mod 3) for n > 1. - Hieronymus Fischer, Nov 11 2007
From Kiyoshi Akima (k_akima(AT)hotmail.com), Aug 31 2009: (Start)
a(n) = 3^floor(n/3)/(1 - (n mod 3)/4), n > 1.
a(n) = 3^(floor((n - 2)/3))*(2 + ((n - 2) mod 3)), n > 1. (End)
a(n) = (2^b)*3^(C - (b + d))*(4^d), n > 1, where C = floor((n + 1)/3), b = max(0, ((n + 1) mod 3) - 1), d = max(0, 1 - ((n + 1) mod 3)). - Jonathan T. Rowell, Jul 26 2011
G.f.: 1 / (1 - x / (1 - x / (1 + x / (1 - x / (1 + x / (1 + x^2 / (1 + x))))))). - Michael Somos, May 12 2012
3*a(n) = 2*a(n+1) if n > 1 and n is not divisible by 3. - Michael Somos, Jan 23 2014
a(n) = a(n-1) + largest proper divisor of a(n-1), n > 2. - Ivan Neretin, Apr 13 2015
a(n) = max{a(i)*a(n-i) : 0 < i < n} for n >= 4. - Jianing Song, Feb 15 2020
|
|
EXAMPLE
|
a{8} = 18 because we have 18 = (8-5)*a(5) = 3*6 and one can verify that this is the maximum.
a(5) = 6: the 7 partitions of 5 are (5), (4, 1), (3, 2), (3, 1, 1), (2, 2, 1), (2, 1, 1, 1), (1, 1, 1, 1, 1) and the corresponding products are 5, 4, 6, 3, 4, 2 and 1; 6 is the largest.
G.f. = 1 + x + 2*x^2 + 3*x^3 + 4*x^4 + 6*x^5 + 9*x^6 + 12*x^7 + 18*x^8 + ...
|
|
MAPLE
|
m := floor(n/3) ;
if n mod 3 = 0 then
3^m ;
elif n mod 3 = 1 then
4*3^(m-1) ;
else
2*3^m ;
end if;
floor(%) ;
|
|
MATHEMATICA
|
a[1] = 1; a[n_] := 4* 3^(1/3 *(n - 1) - 1) /; (Mod[n, 3] == 1 && n > 1); a[n_] := 2*3^(1/3*(n - 2)) /; Mod[n, 3] == 2; a[n_] := 3^(n/3) /; Mod[n, 3] == 0; Table[a[n], {n, 0, 40}]
CoefficientList[Series[(1 + x + 2x^2 + x^4)/(1 - 3x^3), {x, 0, 50}], x] (* Harvey P. Dale, May 01 2011 *)
f[n_] := Max[ Times @@@ IntegerPartitions[n, All, Prime@ Range@ PrimePi@ n]]; f[1] = 1; Array[f, 43, 0] (* Robert G. Wilson v, Jul 31 2012 *)
a[ n_] := If[ n < 2, Boole[ n > -1], 2^Mod[-n, 3] 3^(Quotient[ n - 1, 3] + Mod[n - 1, 3] - 1)]; (* Michael Somos, Jan 23 2014 *)
|
|
PROG
|
(PARI) {a(n) = floor( 3^(n - 4 - (n - 4) \ 3 * 2) * 2^( -n%3))}; /* Michael Somos, Jul 23 2002 */
(PARI) lista(nn) = {print1("1, 1, "); print1(a=2, ", "); for (n=1, nn, a += a/divisors(a)[2]; print1(a, ", "); ); } \\ Michel Marcus, Apr 14 2015
(Haskell)
a000792 n = a000792_list !! n
a000792_list = 1 : f [1] where
f xs = y : f (y:xs) where y = maximum $ zipWith (*) [1..] xs
(Magma) I:=[1, 1, 2, 3, 4]; [n le 5 select I[n] else 3*Self(n-3): n in [1..45]]; // Vincenzo Librandi, Apr 14 2015
|
|
CROSSREFS
|
Cf. array A064364, rightmost (nonvanishing) numbers in row n >= 2.
See A056240 and A288814 for the minimal numbers whose prime factors sums up to n.
|
|
KEYWORD
|
nonn,easy,nice
|
|
AUTHOR
|
|
|
EXTENSIONS
|
More terms and better description from Therese Biedl (biedl(AT)uwaterloo.ca), Jan 19 2000
|
|
STATUS
|
approved
|
|
|
|