cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-20 of 36 results. Next

A344416 Heinz numbers of integer partitions whose sum is even and is at most twice the greatest part.

Original entry on oeis.org

3, 4, 7, 9, 10, 12, 13, 19, 21, 22, 25, 28, 29, 30, 34, 37, 39, 40, 43, 46, 49, 52, 53, 55, 57, 61, 62, 63, 66, 70, 71, 76, 79, 82, 84, 85, 87, 88, 89, 91, 94, 101, 102, 107, 111, 112, 113, 115, 116, 117, 118, 121, 129, 130, 131, 133, 134, 136, 138, 139, 146
Offset: 1

Views

Author

Gus Wiseman, May 20 2021

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.
Also numbers m whose sum of prime indices A056239(m) is even and is at most twice the greatest prime index A061395(m).

Examples

			The sequence of terms together with their prime indices begins:
      3: {2}         37: {12}          71: {20}
      4: {1,1}       39: {2,6}         76: {1,1,8}
      7: {4}         40: {1,1,1,3}     79: {22}
      9: {2,2}       43: {14}          82: {1,13}
     10: {1,3}       46: {1,9}         84: {1,1,2,4}
     12: {1,1,2}     49: {4,4}         85: {3,7}
     13: {6}         52: {1,1,6}       87: {2,10}
     19: {8}         53: {16}          88: {1,1,1,5}
     21: {2,4}       55: {3,5}         89: {24}
     22: {1,5}       57: {2,8}         91: {4,6}
     25: {3,3}       61: {18}          94: {1,15}
     28: {1,1,4}     62: {1,11}       101: {26}
     29: {10}        63: {2,2,4}      102: {1,2,7}
     30: {1,2,3}     66: {1,2,5}      107: {28}
     34: {1,7}       70: {1,3,4}      111: {2,12}
		

Crossrefs

These partitions are counted by A000070 = even-indexed terms of A025065.
The opposite version appears to be A320924, counted by A209816.
The opposite version with odd weights allowed appears to be A322109.
The conjugate opposite version allowing odds is A344291, counted by A110618.
The conjugate version is A344296, also counted by A025065.
The conjugate opposite version is A344413, counted by A209816.
Allowing odd weight gives A344414.
The case of equality is A344415, counted by A035363.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A058696 counts partitions of even numbers, ranked by A300061.
A265640 lists Heinz numbers of palindromic partitions.
A301987 lists numbers whose sum of prime indices equals their product.
A334201 adds up all prime indices except the greatest.
A340387 lists Heinz numbers of partitions whose sum is twice their length.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],EvenQ[Total[primeMS[#]]]&&Max[primeMS[#]]>=Total[primeMS[#]]/2&]

Formula

Intersection of A300061 and A344414.

A290104 a(n) = A003963(n) / A290103(n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 3, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 1, 1, 4, 3, 1, 1, 1, 4, 1, 1, 2, 1, 1, 1, 1, 1, 4, 1, 3, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 2, 1, 1, 8, 1, 1, 2, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 4, 2, 3, 1, 1, 1, 1, 2, 1, 1, 4, 1, 1, 2, 1, 1, 2, 3, 1, 4, 1, 1, 1
Offset: 1

Views

Author

Antti Karttunen, Aug 13 2017

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k). Then a(n) is the product divided by the LCM of the integer partition with Heinz number n. - Gus Wiseman, Aug 01 2018

Examples

			n = 21 = 3 * 7 = prime(2) * prime(4), thus A003963(21) = 2*4 = 8, while A290103(21) = lcm(2,4) = 4, so a(21) = 8/4 = 2.
		

Crossrefs

Differs from A290106 for the first time at n=21.

Programs

  • Mathematica
    Table[If[n == 1, 1, Apply[Times, Map[PrimePi[#1]^#2 & @@ # &, #]] / Apply[LCM, PrimePi[#[[All, 1]] ]]] &@ FactorInteger@ n, {n, 120}] (* Michael De Vlieger, Aug 14 2017 *)
  • Scheme
    (define (A290104 n) (/ (A003963 n) (A290103 n)))

Formula

a(n) = A003963(n) / A290103(n).
Other identities. For all n >= 1:
a(A181819(n)) = A005361(n)/A072411(n).

A316433 Number of integer partitions of n whose length is equal to the LCM of all parts.

Original entry on oeis.org

1, 0, 1, 1, 1, 0, 2, 1, 4, 3, 4, 4, 8, 5, 7, 8, 10, 8, 13, 13, 20, 18, 25, 25, 36, 34, 48, 52, 64, 64, 85, 85, 108, 106, 129, 133, 160, 158, 189, 194, 229, 228, 276, 279, 332, 336, 394, 402, 476, 489, 572, 599, 699, 728, 845, 889, 1032, 1094, 1251, 1332, 1523
Offset: 1

Views

Author

Gus Wiseman, Jul 02 2018

Keywords

Examples

			The a(13) = 8 partitions are (4441), (55111), (322222), (332221), (333211), (622111), (631111), (7111111).
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],LCM@@#==Length[#]&]],{n,30}]
  • PARI
    a(n) = {my(nb = 0); forpart(p=n, if (lcm(Vec(p))==#p, nb++);); nb;} \\ Michel Marcus, Jul 03 2018

A319333 Heinz numbers of integer partitions whose sum is equal to their LCM.

Original entry on oeis.org

2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 30, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 179, 181, 191, 193, 197, 198, 199, 211, 223, 227, 229, 233, 239, 241, 251, 257, 263
Offset: 1

Views

Author

Gus Wiseman, Sep 17 2018

Keywords

Comments

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of partitions whose Heinz numbers are in the sequence begins: (1), (2), (3), (4), (5), (6), (7), (8), (9), (10), (3,2,1), (11), (12).
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[2,100],LCM@@primeMS[#]==Total[primeMS[#]]&]

A340851 Number of factorizations of n such that every factor is a divisor of the number of factors.

Original entry on oeis.org

1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Feb 04 2021

Keywords

Comments

Also factorizations whose number of factors is divisible by their least common multiple.

Examples

			The a(n) factorizations for n = 8192, 46656, 73728:
  2*2*2*2*2*4*8*8          6*6*6*6*6*6              2*2*2*2*2*2*2*2*2*4*6*6
  2*2*2*2*4*4*4*8          2*2*2*2*2*2*3*3*3*3*3*3  2*2*2*2*2*2*2*2*3*4*4*6
  2*2*2*4*4*4*4*4                                   2*2*2*2*2*2*2*3*3*4*4*4
  2*2*2*2*2*2*2*2*2*2*2*4                           2*2*2*2*2*2*2*2*2*2*6*12
                                                    2*2*2*2*2*2*2*2*2*3*4*12
		

Crossrefs

The version for partitions is A340693, with reciprocal version A143773.
Positions of nonzero terms are A340852.
The reciprocal version is A340853.
A320911 can be factored into squarefree semiprimes.
A340597 have an alt-balanced factorization.
A340656 lack a twice-balanced factorization, complement A340657.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A316439 counts factorizations by product and length.
A339846 counts factorizations of even length.
A339890 counts factorizations of odd length.
A340101 counts factorizations into odd factors, odd-length case A340102.
A340653 counts balanced factorizations.
A340785 counts factorizations into even numbers, even-length case A340786.
A340831/A340832 count factorizations with odd maximum/minimum.
A340854 cannot be factored with odd least factor, complement A340855.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],And@@IntegerQ/@(Length[#]/#)&]],{n,100}]

A316556 Number of distinct LCMs of nonempty submultisets of the integer partition with Heinz number n.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 3, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 4, 1, 1, 3, 2, 3, 2, 1, 2, 2, 2, 1, 3, 1, 2, 3, 2, 1, 2, 1, 2, 3, 2, 1, 2, 3, 2, 2, 2, 1, 4, 1, 2, 2, 1, 2, 4, 1, 2, 3, 4, 1, 2, 1, 2, 3, 2, 3, 3, 1, 2, 1, 2, 1, 3, 3, 2, 2, 2, 1, 4, 3, 2, 3, 2, 3, 2, 1, 2, 3, 2, 1, 4, 1, 2, 5
Offset: 1

Views

Author

Gus Wiseman, Jul 06 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Number of distinct values obtained when A290103 is applied to all divisors of n larger than one. - Antti Karttunen, Sep 25 2018

Examples

			462 is the Heinz number of (5,4,2,1) which has possible LCMs of nonempty submultisets {1,2,4,5,10,20} so a(462) = 6.
		

Crossrefs

Cf. also A304793, A305611, A319685, A319695 for other similarly constructed sequences.

Programs

  • Mathematica
    Table[Length[Union[LCM@@@Rest[Subsets[If[n==1,{},Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]]]]]]],{n,100}]
  • PARI
    A290103(n) = lcm(apply(p->primepi(p),factor(n)[,1]));
    A316556(n) = { my(m=Map(),s,k=0); fordiv(n,d,if((d>1)&&!mapisdefined(m,s=A290103(d)), mapput(m,s,s); k++)); (k); }; \\ Antti Karttunen, Sep 25 2018

Extensions

More terms from Antti Karttunen, Sep 25 2018

A327778 Number of integer partitions of n whose LCM is a multiple of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 2, 1, 1, 1, 5, 1, 11, 1, 11, 23, 1, 1, 23, 1, 85, 85, 45, 1, 152, 1, 84, 1, 451, 1, 1787, 1, 1, 735, 260, 1925, 1908, 1, 437, 1877, 4623, 1, 14630, 1, 6934, 10519, 1152, 1, 6791, 1, 1817, 10159, 22556, 1, 2819, 47927, 69333, 22010, 4310, 1
Offset: 0

Views

Author

Gus Wiseman, Sep 25 2019

Keywords

Examples

			The partitions of n = 6, 10, 12, and 15 whose LCM is a multiple of n:
  (6)      (10)         (12)             (15)
  (3,2,1)  (5,3,2)      (5,4,3)          (6,5,4)
           (5,4,1)      (6,4,2)          (7,5,3)
           (5,2,2,1)    (8,3,1)          (9,5,1)
           (5,2,1,1,1)  (4,3,3,2)        (10,3,2)
                        (4,4,3,1)        (5,4,3,3)
                        (6,4,1,1)        (5,5,3,2)
                        (4,3,2,2,1)      (6,5,2,2)
                        (4,3,3,1,1)      (6,5,3,1)
                        (4,3,2,1,1,1)    (10,3,1,1)
                        (4,3,1,1,1,1,1)  (5,3,3,2,2)
                                         (5,3,3,3,1)
                                         (5,4,3,2,1)
                                         (5,5,3,1,1)
                                         (6,5,2,1,1)
                                         (5,3,2,2,2,1)
                                         (5,3,3,2,1,1)
                                         (5,4,3,1,1,1)
                                         (6,5,1,1,1,1)
                                         (5,3,2,2,1,1,1)
                                         (5,3,3,1,1,1,1)
                                         (5,3,2,1,1,1,1,1)
                                         (5,3,1,1,1,1,1,1,1)
		

Crossrefs

The Heinz numbers of these partitions are given by A327783.
Partitions whose LCM is equal to their sum are A074761.
Partitions whose LCM is greater than their sum are A327779.
Partitions whose LCM is less than their sum are A327781.

Programs

  • Maple
    a:= proc(m) option remember; local b; b:=
          proc(n, i, l) option remember; `if`(n=0 or i=1,
            `if`(l=m, 1, 0), `if`(i<2, 0, b(n, i-1, l))+
             b(n-i, min(n-i, i), igcd(m, ilcm(l, i))))
          end; `if`(isprime(m), 1, b(m$2, 1))
        end:
    seq(a(n), n=0..60);  # Alois P. Heinz, Sep 26 2019
  • Mathematica
    Table[Length[Select[IntegerPartitions[n],Divisible[LCM@@#,n]&]],{n,30}]
    (* Second program: *)
    a[m_] := a[m] = Module[{b}, b[n_, i_, l_] := b[n, i, l] = If[n == 0 || i == 1, If[l == m, 1, 0], If[i<2, 0, b[n, i - 1, l]] + b[n - i, Min[n - i, i], GCD[m, LCM[l, i]]]]; If[PrimeQ[m], 1, b[m, m, 1]]];
    a /@ Range[0, 60] (* Jean-François Alcover, May 18 2021, after Alois P. Heinz *)

Formula

a(n) = 1 <=> n in { A000961 }. - Alois P. Heinz, Sep 26 2019

A340693 Number of integer partitions of n where each part is a divisor of the number of parts.

Original entry on oeis.org

1, 1, 1, 2, 2, 3, 2, 5, 5, 7, 7, 10, 10, 14, 14, 17, 19, 24, 24, 32, 33, 42, 43, 58, 59, 75, 79, 98, 104, 124, 128, 156, 166, 196, 204, 239, 251, 292, 306, 352, 372, 426, 445, 514, 543, 616, 652, 745, 790, 896, 960, 1080, 1162, 1311, 1400, 1574, 1692, 1892
Offset: 0

Views

Author

Gus Wiseman, Jan 23 2021

Keywords

Comments

The only strict partitions counted are (), (1), and (2,1).
Is there a simple generating function?

Examples

			The a(1) = 1 through a(9) = 7 partitions:
  1  11  21   22    311    2211    331      2222      333
         111  1111  2111   111111  2221     4211      4221
                    11111          4111     221111    51111
                                   211111   311111    222111
                                   1111111  11111111  321111
                                                      21111111
                                                      111111111
		

Crossrefs

Note: Heinz numbers are given in parentheses below.
The reciprocal version is A143773 (A316428), with strict case A340830.
The case where length also divides n is A326842 (A326847).
The Heinz numbers of these partitions are A340606.
The version for factorizations is A340851, with reciprocal version A340853.
A018818 counts partitions of n into divisors of n (A326841).
A047993 counts balanced partitions (A106529).
A067538 counts partitions of n whose length/max divides n (A316413/A326836).
A067539 counts partitions with integer geometric mean (A326623).
A072233 counts partitions by sum and length.
A168659 = partitions whose greatest part divides their length (A340609).
A168659 = partitions whose length divides their greatest part (A340610).
A326843 = partitions of n whose length and maximum both divide n (A326837).
A330950 = partitions of n whose Heinz number is divisible by n (A324851).

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And@@IntegerQ/@(Length[#]/#)&]],{n,0,30}]

A256554 Number T(n,k) of cycle types of degree-n permutations having the k-th smallest possible order; triangle T(n,k), n>=0, 1<=k<=A009490(n), read by rows.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 3, 2, 2, 1, 2, 1, 3, 2, 2, 1, 3, 1, 1, 1, 1, 4, 2, 4, 1, 5, 1, 1, 1, 1, 1, 1, 4, 3, 4, 1, 7, 1, 1, 1, 2, 2, 1, 1, 1, 1, 5, 3, 6, 2, 9, 1, 2, 1, 3, 4, 1, 1, 1, 1, 1, 1, 5, 3, 6, 2, 12, 1, 2, 1, 4, 1, 6, 2, 2, 1, 2, 1, 1, 1, 2
Offset: 0

Views

Author

Alois P. Heinz, Apr 01 2015

Keywords

Comments

Sum_{k>=0} A256553(n,k)*T(n,k) = A181844(n).

Examples

			Triangle T(n,k) begins:
  1;
  1;
  1, 1;
  1, 1, 1;
  1, 2, 1, 1;
  1, 2, 1, 1, 1, 1;
  1, 3, 2, 2, 1, 2;
  1, 3, 2, 2, 1, 3, 1, 1, 1;
  1, 4, 2, 4, 1, 5, 1, 1, 1, 1, 1;
  1, 4, 3, 4, 1, 7, 1, 1, 1, 2, 2, 1, 1, 1;
  1, 5, 3, 6, 2, 9, 1, 2, 1, 3, 4, 1, 1, 1, 1, 1;
		

Crossrefs

Row sums give A000041.
Row lengths give A009490.
Columns k=1-9 give: A000012, A004526, A002264, A008642(n-4), A002266, A074752, A132270, A008643(n-8) for n>7, A008649(n-9) for n>8.
Last elements of rows give A074064.
Main diagonal gives A074761.

Programs

  • Maple
    b:= proc(n, i) option remember; `if`(n=0 or i=1, x,
          b(n, i-1)+(p-> add(coeff(p, x, t)*x^ilcm(t, i),
          t=1..degree(p)))(add(b(n-i*j, i-1), j=1..n/i)))
        end:
    T:= n->(p->seq((h->`if`(h=0, [][], h))(coeff(p, x, i))
         , i=1..degree(p)))(b(n$2)):
    seq(T(n), n=0..12);
  • Mathematica
    b[n_, i_] := b[n, i] = If[n == 0 || i == 1, x, b[n, i - 1] + Function[p, Sum[Coefficient[p, x, t]*x^LCM[t, i], {t, 1, Exponent[p, x]}]][Sum[b[n - i*j, i - 1], {j, 1, n/i}]]]; T[n_] := Function[p, Table[Function[h, If[h == 0, {{}, {}}, h]][Coefficient[p, x, i]], {i, 1, Exponent[p, x]}]][b[n, n]]; Table[T[n], {n, 0, 12}] // Flatten (* Jean-François Alcover, Jan 23 2017, translated from Maple *)

A317624 Number of integer partitions of n where all parts are > 1 and whose LCM is n.

Original entry on oeis.org

0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 3, 1, 5, 1, 1, 1, 17, 1, 1, 1, 7, 1, 60, 1, 1, 1, 1, 1, 76, 1, 1, 1, 55, 1, 105, 1, 11, 10, 1, 1, 187, 1, 6, 1, 13, 1, 30, 1, 111, 1, 1, 1, 5043, 1, 1, 15, 1, 1, 230, 1, 17, 1, 242, 1, 4173, 1, 1, 12, 19, 1
Offset: 0

Views

Author

Gus Wiseman, Aug 01 2018

Keywords

Examples

			The a(20) = 5 partitions are (20), (10,4,4,2), (10,4,2,2,2), (5,5,4,4,2), (5,5,4,2,2,2).
The a(45) = 10 partitions:
  (45),
  (15,15,9,3,3), (15,9,9,9,3),
  (15,9,9,3,3,3,3), (15,9,5,5,5,3,3), (9,9,9,5,5,5,3),
  (15,9,3,3,3,3,3,3,3), (9,9,5,5,5,3,3,3,3), (9,5,5,5,5,5,5,3,3),
  (9,5,5,5,3,3,3,3,3,3,3).
From _David A. Corneth_, Sep 08 2018: (Start)
Let sum(t) denote the sum of elements of a tuple t. The tuples t with distinct divisors of 45 that have lcm(t) = 45 and sum(t) <= 45 are {(45) and (3, 9, 15), (3, 5, 9, 15), (3, 5, 9), (5, 9), (9, 15), (5, 9, 15)}. For each such tuple t, find the number of partitions of 45 - s(t) into distinct parts of t.
For the tuple (45), there is 1 partition of 45 - 45 = 0 into parts with 45. That is: {()}.
For the tuple (3, 9, 15), there are 4 partitions of 45 - (3 + 9 + 15) = 18 into parts with 3, 9 and 15. They are {(3, 15), (9, 9), (3, 3, 3, 9), (3, 3, 3, 3, 3, 3)}.
For the tuple (3, 5, 9), there are 4 partitions of 45 - (3 + 5 + 9) = 28 into parts with 3, 5 and 9; they are {(5, 5, 9, 9), (3, 3, 3, 5, 5, 9), (3, 5, 5, 5, 5, 5), (3, 3, 3, 3, 3, 3, 5, 5)}.
For the tuple (3, 5, 9, 15), there is 1 partition of 45 - (3 + 5 + 9 + 15) = 13 into parts with 3, 5, 9 and 15. That is (3, 5, 5).
The other tuples, (5, 9), (9, 15), and (5, 9, 15); they give no extra tuples. That's because there is no solution to the Diophantine equation for 5x + 9y = 45 - (5 + 9), corresponding to the tuple (5, 9) with nonnegative x, y.
That also excludes (9, 15); if there is a solution for that, there would also be a solution for (5, 9). This could whittle down the number of seeds even further. Similarly, (5, 9, 15) gives no solution.
Therefore a(45) = 1 + 4 + 4 + 1 = 10.
(End)
In general, there are A318670(n) (<= A069626(n)) such seed sets of divisors where to start extending the partition from. (See the second PARI program which uses subroutine toplevel_starting_sets.) - _Antti Karttunen_, Sep 08 2018
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],And[Min@@#>=2,LCM@@#==n]&]],{n,30}]
  • PARI
    strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4);
    partitions_into_lcm(orgn,n,parts,from=1,m=1) = if(!n,(m==orgn),my(k = #parts, s=0); for(i=from,k,if(parts[i]<=n, s += partitions_into_lcm(orgn,n-parts[i],parts,i,lcm(m,parts[i])))); (s));
    A317624(n) = if(n<=1,0,partitions_into_lcm(n,n,strong_divisors_reversed(n))); \\ Antti Karttunen, Sep 07 2018
    
  • PARI
    strong_divisors_reversed(n) = vecsort(select(x -> (x>1), divisors(n)), , 4);
    partitions_into(n,parts,from=1) = if(!n,1, if(#parts==from, (0==(n%parts[from])), my(s=0); for(i=from,#parts,if(parts[i]<=n, s += partitions_into(n-parts[i],parts,i))); (s)));
    toplevel_starting_sets(orgn,n,parts,from=1,ss=List([])) = { my(k = #parts, s=0, newss); if(lcm(Vec(ss))==orgn,s += partitions_into(n,ss)); for(i=from,k,if(parts[i]<=n, newss = List(ss); listput(newss,parts[i]); s += toplevel_starting_sets(orgn,n-parts[i],parts,i+1,newss))); (s) };
    A317624(n) = if(n<=1,0,toplevel_starting_sets(n,n,strong_divisors_reversed(n))); \\ Antti Karttunen, Sep 08-10 2018
Previous Showing 11-20 of 36 results. Next