Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A126796
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A126796 Number of complete partitions of n. +0
1
1, 1, 1, 2, 2, 4, 5, 8, 10, 16, 20, 31, 39, 55, 71, 100, 125, 173, 218, 291, 366, 483, 600, 784, 971, 1244, 1538, 1957, 2395, 3023, 3693, 4605, 5604, 6942, 8397, 10347, 12471, 15235, 18309, 22267, 26619, 32219, 38414, 46216, 54941, 65838 (list; graph; listen)
OFFSET

0,4

COMMENT

A partition of n is complete if every number 1 to n can be represented as a sum of parts of the partition. This generalizes perfect partitions, where the representation for each number must be unique.

A partition is complete iff each part is no more than 1 more than the sum of all smaller parts. (This includes the smallest part, which thus must be 1.) - Franklin T. Adams-Watters (FrankTAW(AT)Netscape.net), Mar 22 2007

REFERENCES

SeungKyung Park, Complete Partitions, Fibonacci Quarterly, Vol. 36 (1998), pp. 354-360.

LINKS

David W. Wilson, Table of n, a(n) for n = 0..300

EXAMPLE

a(5) = 4 because there are 4 complete partitions of 5: {1, 1, 1, 1, 1}, {1, 1, 1, 2}, {1, 2, 2} and {1, 1, 3}.

MAPLE

isCompl := proc(p, n) local m, pers, reps, f, lst, s; reps := {}; pers := combinat[permute](p); for m from 1 to nops(pers) do lst := op(m, pers); for f from 1 to nops(lst) do s := add( op(i, lst), i=1..f); reps := reps union {s}; od; od; for m from 1 to n do if not m in reps then RETURN(false); fi; od; RETURN(true); end: A126796 := proc(n) local prts, res, p; prts := combinat[partition](n); res :=0; for p from 1 to nops(prts) do if isCompl(op(p, prts), n) then res := res+1; fi; od; RETURN(res); end: for n from 1 to 40 do printf("%d %d ", n, A126796(n)); od; - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 27 2007

At the beginning of the 2nd Maple program replace the current 15 by any other positive integer n in order to obtain a(n). - Emeric Deutsch (deutsch(AT)duke.poly.edu), Mar 04 2007

with(combinat): a:=proc(n) local P, b, k, p, S, j: P:=partition(n): b:=0: for k from 1 to numbpart(n) do p:=powerset(P[k]): S:={}: for j from 1 to nops(p) do S:=S union {add(p[j][i], i=1..nops(p[j]))} od: if nops(S)=n+1 then b:=b+1 else b:=b: fi: od: end: seq(a(n), n=1..30); - Emeric Deutsch (deutsch(AT)duke.poly.edu), Mar 04 2007

with(combinat): n:=15: P:=partition(n): b:=0: for k from 1 to numbpart(n) do p:=powerset(P[k]): S:={}: for j from 1 to nops(p) do S:=S union {add(p[j][i], i=1..nops(p[j]))} od: if nops(S)=n+1 then b:=b+1 else b:=b: fi: od: b; - Emeric Deutsch (deutsch(AT)duke.poly.edu), Mar 04 2007

PROGRAM

(PARI) T(n, k)=if(k<=1, 1, if(n<2*k-1, T(n, floor((n+1)/2)), T(n, k-1)+T(n-k, k))) a(n)=T(n, floor((n+1)/2)) /* If modified to save earlier results, this would be efficient. */ - Franklin T. Adams-Watters (FrankTAW(AT)Netscape.net), Mar 22 2007

CROSSREFS

Cf. A002033.

Sequence in context: A104504 A027337 A027193 this_sequence A157162 A109434 A089299

Adjacent sequences: A126793 A126794 A126795 this_sequence A126797 A126798 A126799

KEYWORD

nonn

AUTHOR

Brian Hopkins (bhopkins(AT)spc.edu), Feb 20 2007

EXTENSIONS

More terms from R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Feb 27 2007

More terms from Emeric Deutsch (deutsch(AT)duke.poly.edu), Mar 04 2007

Further terms from Franklin T. Adams-Watters (FrankTAW(AT)Netscape.net) and David W. Wilson, Mar 22 2007

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified November 25 13:47 EST 2009. Contains 167481 sequences.


AT&T Labs Research