Dismiss Notice
Join Physics Forums Today!
The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

Partitions with restrictions

  1. Nov 4, 2011 #1
    After long and careful search on the web and in literature,
    I could not find the solution of the following problem.

    I need calculate p(N,K,L) - the number of partitions of N into
    no more than K parts not exceeding L.

    Example: N = 7, K = 4, L = 5

    1) 2+5
    2) 3+4
    3) 1+1+5
    4) 1+2+4
    5) 1+3+3
    6) 2+2+3
    7) 1+1+1+4
    8) 1+1+2+3
    9) 1+2+2+2

    So, here is p(7,4,5)=9

    I found the different formulas for the partitions with different restrictions,
    but not both for quantity and size of parts.
    May be I was bad looking, but might be the solution has not yet been found?

    Of course, it is not a problem to write search algorithms, to solve this problem,
    as in Pascal:

    *********************************************************************

    var K,L,M,i,i1,i2,itog:longint; b:real;
    Code:integer;

    procedure Box(pr:longint;ostatok:longint;nbox:longint);
    var j1,j2,j:longint;
    a: real;
    begin
    a:=ostatok/nbox;
    if frac(a)>0
    then j1:=trunc(a)+1
    else j1:=trunc(a);
    if pr>ostatok then j2:=ostatok else j2:= pr;
    for j:=j1 to j2 do
    if nbox>1 then Box(j,ostatok-j,nbox-1) else itog:=itog+1
    end;

    begin

    val(ParamStr(1),K,Code);
    val(ParamStr(2),L,Code);
    val(ParamStr(3),M,Code);
    {
    writeln(K);
    writeln(L);
    writeln(M);
    }
    itog:=0;
    If L*M<K Then writeln('It is impossible')
    Else
    begin
    If M>K Then i2:=K Else i2:=M;
    b:=K/L;
    If frac(b)>0 Then i1:=trunc(b)+1 Else i1:=trunc(b);

    For i:=i1 to i2 do Box(i,K-i,L-1)

    end;
    writeln (Itog)
    end.

    **********************************************************

    but it bogged down with a slight increase of parameters.
    Roughly speaking, if the value of parameters begins to run into the hundreds,
    then modern computer begins to squeak.

    Is there some "nice" formula to calculate p(N,K,L) ?
     
    Last edited: Nov 4, 2011
  2. jcsd
  3. Nov 5, 2011 #2
    Nobody knows, but it seems unlikely that a formula exists.
     
  4. Nov 5, 2011 #3
    Last edited by a moderator: Apr 26, 2017
  5. Nov 5, 2011 #4
    The question was if there is a "nice" formula, not a "nice" algorithm. The formula, for now, isn't discovered, maybe it doesn't exists, but the algorithm to calculate p(n,l,k) is easy to make.
     
  6. Nov 5, 2011 #5
    Easy-to-make algorithms will go into deep-deep recursion, so it is not easy to make really fast one.
    For example, shown code in Pascal will fall even on low parameters like p(200,20,20).
    After some thinking I wrote algorithm wich gives result in suitable time for p(125000,500,500).
    It uses arbitrary length integer arithmetics and gives result of near thousand bits.
    But I need about p(millions,thousands,thousands) job.
     
  7. Nov 5, 2011 #6
    In this case the problem is a P problem in information theory, in other words the time to calculate p(n,l,k) is Polynomial, so it will spend more time to calculate p(millions,thousands,thousands). Only a formula will give you a gift, but it is not discovered yet, and maybe it's likely that it doesn't exist.
     
Know someone interested in this topic? Share this thread via Reddit, Google+, Twitter, or Facebook




Similar Discussions: Partitions with restrictions
  1. Monotonic partitioning (Replies: 0)

  2. Partition Function (Replies: 0)

  3. Partitioning problem (Replies: 2)

  4. Goldbach Partitions (Replies: 12)

  5. Cosets / Partitions (Replies: 4)

Loading...