A038041 -id:A038041 - cp's OEIS frontend

A371734 Maximal length of a factorization of n into factors > 1 all having different sums of prime indices.

0, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 2, 3, 1, 2, 2, 3, 1, 3, 1, 2, 2, 2, 1, 3, 1, 2, 2, 2, 1, 3, 2, 3, 2, 2, 1, 3, 1, 2, 2, 3, 2, 3, 1, 2, 2, 3, 1, 3, 1, 2, 2, 2, 2, 3, 1, 3, 2, 2, 1, 3, 2, 2, 2, 3, 1, 3, 2, 2, 2, 2, 2, 3, 1, 2, 2, 3, 1, 3, 1, 3, 3

Offset: 1

Gus Wiseman, Apr 13 2024

nonn

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798. Sum of prime indices is given by A056239.

Factorizations into factors > 1 all having different sums of prime indices are counted by A321469.

			The factorizations of 90 of this type are (2*3*15), (2*5*9), (2*45), (3*30), (5*18), (6*15), (90), so a(90) = 3.

For set partitions of binary indices we have A000120, same sums A371735.
Positions of 1's are A000430.
Positions of terms > 1 are A080257.
Factorizations of this type are counted by A321469, same sums A321455.
For same instead of different sums we have A371733.
A001055 counts factorizations.
A002219 (aerated) counts biquanimous partitions, ranks A357976.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A321451 counts non-quanimous partitions, ranks A321453.
A321452 counts quanimous partitions, ranks A321454.
Cf. A035470, A279787, A305551, A321142, A322794, A326515, A326518, A326534, A336137, A371783, A371789, A371791.

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&, Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
Table[Max[Length/@Select[facs[n],UnsameQ@@hwt/@#&]],{n,100}]

A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
all_have_different_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == #facs));
A371734(n, m=n, facs=List([])) = if(1==n, if(all_have_different_sum_of_pis(facs),#facs,0), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s = max(s,A371734(n/d, d, newfacs)))); (s)); \\ Antti Karttunen, Jan 20 2025

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A371735 Maximal length of a set partition of the binary indices of n into blocks all having the same sum.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 2, 1, 1, 2, 3, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 2, 1, 1, 1, 1, 3, 2, 1, 2, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 1, 1, 2, 1

Offset: 0

Gus Wiseman, Apr 14 2024

nonn

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.

If a(n) = k then the binary indices of n (row n of A048793) are k-quanimous (counted by A371783).

			The binary indices of 119 are {1,2,3,5,6,7}, and the set partitions into blocks with the same sum are:
  {{1,7},{2,6},{3,5}}
  {{1,5,6},{2,3,7}}
  {{1,2,3,6},{5,7}}
  {{1,2,3,5,6,7}}
So a(119) = 3.

Set partitions of this type are counted by A035470, A336137.
A version for factorizations is A371733.
Positions of 1's are A371738.
Positions of terms > 1 are A371784.
A001055 counts factorizations.
A002219 (aerated) counts biquanimous partitions, ranks A357976.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A070939 gives length of binary expansion.
A321452 counts quanimous partitions, ranks A321454.
A326031 gives weight of the set-system with BII-number n.
A371783 counts k-quanimous partitions.
A371789 counts non-quanimous sets, differences A371790.
A371796 counts quanimous sets, differences A371797.
Cf. A006827, A038041, A096111, A279787, A305551, A321451, A321455, A322794, A326534, A371731, A371734.

sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]& /@ sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
Table[Max[Length/@Select[sps[bix[n]],SameQ@@Total/@#&]],{n,0,100}]

A087905 a(n) = n! * Sum_{d|n} (d/n)^d.

Original entry on oeis.org

1, 3, 8, 36, 144, 1010, 5760, 50400, 416640, 4250232, 43545600, 553106400, 6706022400, 95865541200, 1410695430144, 22720842144000, 376610217984000, 6888030445296000, 128047474114560000, 2587520533615041024

Offset: 1

Vladeta Jovovic, Oct 14 2003

nonn

G. C. Greubel, Table of n, a(n) for n = 1..445

Cf. A061095, A038041, A057625, A038048, A005225.

a[n_]:= n!*DivisorSum[n, (#/n)^# &]; Array[a, 50] (* G. C. Greubel, May 16 2018 *)

{a(n)= n!*sumdiv(n, d, (d/n)^d)};
for(n=1, 30, print1(a(n), ", ")) \\ G. C. Greubel, May 16 2018

E.g.f.: Sum_{k>0} x^k/(k-x^k).

A200473 Irregular triangle read by rows: T(n,k) = number of ways to assign n people to d_k unlabeled groups of equal size (where d_k is the k-th divisor of n).

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 10, 15, 1, 1, 1, 1, 35, 105, 1, 1, 280, 1, 1, 126, 945, 1, 1, 1, 1, 462, 5775, 15400, 10395, 1, 1, 1, 1, 1716, 135135, 1, 1, 126126, 1401400, 1, 1, 6435, 2627625, 2027025, 1, 1, 1, 1, 24310, 2858856, 190590400, 34459425, 1, 1

Offset: 1

Dennis P. Walsh, Nov 18 2011

This sequence is A200472 with zeros removed.

			T(n,k) begins:
1;
1,      1;
1,      1;
1,      3,       1;
1,      1;
1,     10,      15,       1;
1,      1;
1,     35,     105,       1;
1,    280,       1;
1,    126,     945,       1;
1,      1;
1,    462,    5775,   15400, 10395,   1;
1,      1;
1,   1716,  135135,       1;
1, 126126, 1401400,       1;
1,   6435, 2627625, 2027025,     1;

Alois P. Heinz, Rows n = 1..250, flattened
Dennis P. Walsh, Note on assigning n people to k unlabeled groups of equal size

Cf. A200472, A000005 (row lengths).

Cf. A038041 (row sums).

with(numtheory):
S:= n-> sort([divisors(n)[]]):
T:= (n,k)-> n!/(S(n)[k])!/((n/(S(n)[k]))!)^(S(n)[k]):
seq(seq(T(n, k), k=1..tau(n)), n=1..10);

row[n_] := (n!/#!)/(n/#)!^#& /@ Divisors[n];
Table[row[n], {n, 1, 20}] // Flatten (* Jean-François Alcover, Mar 24 2017 *)

T(n,k) = (n!/d_k!)/(n/d_k)!^d_k, n>=1, 1<=k<=tau(n), d_k = k-th divisor of n.

Sum_{k=1..tau(k)} T(n,k) = A038041(n). - Alois P. Heinz, Jul 22 2016

A262280 Number of ways to select a nonempty subset s from an n-set and then partition s into blocks of equal size.

Original entry on oeis.org

0, 1, 4, 11, 29, 72, 190, 527, 1552, 5031, 18087, 66904, 266381, 1164516, 5215644, 23868103, 117740143, 609872350, 3268548406, 18110463455, 102867877414, 620476915965, 4005216028161, 25747549921338, 166978155172420, 1168774024335203, 8556355097320141

Offset: 0

Alois P. Heinz, Sep 17 2015

nonn

			a(3) = 11: 1, 2, 3, 12, 1|2, 13, 1|3, 23, 2|3, 123, 1|2|3.

Alois P. Heinz, Table of n, a(n) for n = 0..616

Cf. A038041, A262320.

b:= proc(n) option remember;
      add(1/(d!*(n/d)!^d), d=numtheory[divisors](n))
    end:
a:= n-> n! * add(b(k)/(n-k)!, k=1..n):
seq(a(n), n=0..30);

b[n_] := b[n] = DivisorSum[n, 1/(#!*(n/#)!^#)&]; a[n_] := n!*Sum[b[k]/(n-k)!, {k, 1, n}]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 15 2017, translated from Maple *)

E.g.f.: exp(x) * Sum_{k>=1} (exp(x^k/k!)-1).
a(n) = Sum_{k=1..n} C(n,k) * A038041(k).
a(n) = A262320(n) - 1.

A275389 Number of set partitions of [n] with a strongly unimodal block size list.

Original entry on oeis.org

1, 1, 1, 4, 7, 19, 71, 219, 759, 2697, 12395, 47477, 231950, 1040116, 4851742, 26690821, 131478031, 736418510, 4262619682, 24680045903, 145629814329, 935900941506, 5778263418232, 37626913475878, 257550263109475, 1782180357952449, 12526035635331581

Offset: 0

Alois P. Heinz, Jul 26 2016

nonn

Strongly unimodal means strictly increasing then strictly decreasing.

			a(3) = 4: 123, 12|3, 13|2, 1|23.
a(4) = 7: 1234, 123|4, 124|3, 134|2, 1|234, 1|23|4, 1|24|3.
a(5) = 19: 12345, 1234|5, 1235|4, 123|45, 1245|3, 124|35, 125|34, 12|345, 1345|2, 134|25, 135|24, 13|245, 145|23, 14|235, 15|234, 1|2345, 1|234|5, 1|235|4, 1|245|3.

Alois P. Heinz, Table of n, a(n) for n = 0..677
Wikipedia, Partition of a set

Cf. A007837, A038041, A059618, A275309, A275310, A275311, A275312, A275313, A286077.

b:= proc(n, i, t) option remember; `if`(t=0 and n>i*(i-1)/2, 0,
     `if`(n=0, 1, add(b(n-j, j, 0)*binomial(n-1, j-1), j=1..min(n, i-1))
     +`if`(t=1, add(b(n-j, j, 1)*binomial(n-1, j-1), j=i+1..n), 0)))
    end:
a:= n-> b(n, 0, 1):
seq(a(n), n=0..30);

b[n_, i_, t_] := b[n, i, t] = If[t==0 && n > i*(i-1)/2, 0, If[n==0, 1, Sum[b[n-j, j, 0]*Binomial[n-1, j-1], {j, 1, Min[n, i-1]}] + If[t==1, Sum[b[n-j, j, 1]*Binomial[n-1, j-1], {j, i+1, n}], 0]]]; a[n_] := b[n, 0, 1]; Table[a[n], {n, 0, 30}] (* Jean-François Alcover, Feb 07 2017, translated from Maple *)

A301924 Regular triangle where T(n,k) is the number of unlabeled k-uniform connected hypergraphs spanning n vertices.

Original entry on oeis.org

1, 0, 1, 0, 2, 1, 0, 6, 3, 1, 0, 21, 29, 4, 1, 0, 112, 2101, 150, 5, 1, 0, 853, 7011181, 7013164, 1037, 6, 1, 0, 11117, 1788775603301, 29281354507753847, 1788782615612, 12338, 7, 1, 0, 261080, 53304526022885278403, 234431745534048893449761040648508, 234431745534048922729326772799024, 53304527811667884902, 274659, 8, 1

Offset: 1

Gus Wiseman, Jun 19 2018

			Triangle begins:
   1
   0    1
   0    2       1
   0    6       3       1
   0   21      29       4    1
   0  112    2101     150    5 1
   0  853 7011181 7013164 1037 6 1
   ...
The T(4,2) = 6 hypergraphs:
  {{1,3},{2,4},{3,4}}
  {{1,4},{2,4},{3,4}}
  {{1,2},{1,3},{2,4},{3,4}}
  {{1,4},{2,3},{2,4},{3,4}}
  {{1,3},{1,4},{2,3},{2,4},{3,4}}
  {{1,2},{1,3},{1,4},{2,3},{2,4},{3,4}}

Row sums are A301920.
Columns k=2..3 are A001349(n > 1), A003190(n > 1).
Cf. A003465, A006129, A038041, A055621, A298422, A299353, A299354, A299471, A301481, A301922, A306017-A306021, A309858.

InvEulerT(v)={my(p=log(1+x*Ser(v))); dirdiv(vector(#v,n,polcoeff(p,n)), vector(#v,n,1/n))}
permcount(v)={my(m=1, s=0, k=0, t); for(i=1, #v, t=v[i]; k=if(i>1&&t==v[i-1], k+1, 1); m*=t*k; s+=t); s!/m}
rep(typ)={my(L=List(), k=0); for(i=1, #typ, k+=typ[i]; listput(L, k); while(#L0, u=vecsort(apply(f, u)); d=lex(u, v)); !d}
Q(n, k, perm)={my(t=0); forsubset([n, k], v, t += can(Vec(v), t->perm[t])); t}
U(n, k)={my(s=0); forpart(p=n, s += permcount(p)*2^Q(n, k, rep(p))); s/n!}
A(n)={Mat(vector(n, k, InvEulerT(vector(n,i,U(i,k)-U(i-1,k)))~))}
{ my(T=A(8)); for(n=1, #T, print(T[n,1..n])) } \\ Andrew Howroyd, Aug 26 2019

Column k is the inverse Euler transform of column k of A301922. - Andrew Howroyd, Aug 26 2019

Terms a(16) and beyond from Andrew Howroyd, Aug 26 2019

A305552 Number of uniform normal multiset partitions of weight n.

Original entry on oeis.org

1, 1, 3, 5, 12, 17, 47, 65, 170, 277, 655, 1025, 2739, 4097, 10281, 17257, 41364, 65537, 170047, 262145, 660296, 1094457, 2621965, 4194305, 10898799, 16792721, 41945103, 69938141, 168546184, 268435457, 694029255, 1073741825, 2696094037, 4474449261, 10737451027

Offset: 0

Gus Wiseman, Jun 20 2018

nonn

A multiset is normal if it spans an initial interval of positive integers. A multiset partition m is uniform if all parts have the same size, and normal if all parts are normal. The weight of m is the sum of sizes of its parts.

			The a(4) = 12 uniform normal multiset partitions:
{1111}, {1222}, {1122}, {1112}, {1233}, {1223}, {1123}, {1234},
{11,11}, {11,12}, {12,12},
{1,1,1,1}.

Andrew Howroyd, Table of n, a(n) for n = 0..1000

Cf. A000005, A001315, A007716, A034691, A038041, A074854, A289078, A305552, A306017.

Table[Sum[Binomial[2^(n/k-1)+k-1,k],{k,Divisors[n]}],{n,35}]

a(n)={if(n<1, n==0, sumdiv(n, d, binomial(2^(n/d - 1) + d - 1, d)))} \\ Andrew Howroyd, Jun 22 2018

a(n) = Sum_{d|n} binomial(2^(n/d - 1) + d - 1, d).

A306437 Regular triangle read by rows where T(n,k) is the number of non-crossing set partitions of {1, ..., n} in which all blocks have size k.

Original entry on oeis.org

1, 1, 1, 1, 0, 1, 1, 2, 0, 1, 1, 0, 0, 0, 1, 1, 5, 3, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 14, 0, 4, 0, 0, 0, 1, 1, 0, 12, 0, 0, 0, 0, 0, 1, 1, 42, 0, 0, 5, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 132, 55, 22, 0, 6, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 429, 0, 0, 0, 0, 7, 0, 0, 0, 0, 0, 0, 1

Offset: 1

Gus Wiseman, Feb 15 2019

			Triangle begins:
  1
  1   1
  1   0   1
  1   2   0   1
  1   0   0   0   1
  1   5   3   0   0   1
  1   0   0   0   0   0   1
  1  14   0   4   0   0   0   1
  1   0  12   0   0   0   0   0   1
  1  42   0   0   5   0   0   0   0   1
  1   0   0   0   0   0   0   0   0   0   1
  1 132  55  22   0   6   0   0   0   0   0   1
Row 6 counts the following non-crossing set partitions (empty columns not shown):
  {{1}{2}{3}{4}{5}{6}}  {{12}{34}{56}}  {{123}{456}}  {{123456}}
                        {{12}{36}{45}}  {{126}{345}}
                        {{14}{23}{56}}  {{156}{234}}
                        {{16}{23}{45}}
                        {{16}{25}{34}}

Germain Kreweras, Sur les partitions non croisées d'un cycle, Discrete Math. 1 333-350 (1972).
Wikipedia, Noncrossing partition.

Row sums are A194560. Column k=2 is A126120. Trisection of column k=3 is A001764.

Cf. A000108, A000110, A000296, A001006, A001263, A001610, A016098, A038041, A061095, A125181, A134264.

T:= (n, k)-> `if`(irem(n, k)=0, binomial(n, n/k)/(n-n/k+1), 0):
seq(seq(T(n,k), k=1..n), n=1..14);  # Alois P. Heinz, Feb 16 2019

Table[Table[If[Divisible[n,d],d/n*Binomial[n,n/d-1],0],{d,n}],{n,15}]

If d|n, then T(n, d) = binomial(n, n/d)/(n - n/d + 1); otherwise T(n, k) = 0 [Theorem 1 of Kreweras].

A322531 Heinz numbers of integer partitions whose parts all have the same number of prime factors (counted with or without multiplicity) and whose product of parts is a squarefree number.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 11, 13, 15, 16, 17, 29, 31, 32, 33, 41, 43, 47, 51, 55, 59, 64, 67, 73, 79, 83, 85, 93, 101, 109, 113, 123, 127, 128, 137, 139, 149, 155, 157, 163, 165, 167, 177, 179, 181, 187, 191, 199, 201, 205, 211, 233, 241, 249, 255, 256, 257, 269, 271

Offset: 1

Gus Wiseman, Dec 14 2018

nonn

The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).
All entries are themselves squarefree numbers (except the powers of 2).
The first odd term not in this sequence but in A302521 is 141, which is the MM-number (see A302242) of {{1},{2,3}}.

			The sequence of all integer partitions whose parts all have the same number of prime factors and whose product of parts is a squarefree number begins: (), (1), (2), (1,1), (3), (1,1,1), (5), (6), (3,2), (1,1,1,1), (7), (10), (11), (1,1,1,1,1), (5,2), (13), (14), (15), (7,2), (5,3), (17), (1,1,1,1,1,1).

Cf. A003963, A005117, A038041, A056239, A073576, A112798, A302242, A302505, A306017, A319056, A319169, A320324, A321717, A321718, A322526, A322528.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Select[Range[100],And[SameQ@@PrimeOmega/@primeMS[#],SquareFreeQ[Times@@primeMS[#]]]&]

cp's OEIS Frontend

A371734 Maximal length of a factorization of n into factors > 1 all having different sums of prime indices.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A371735 Maximal length of a set partition of the binary indices of n into blocks all having the same sum.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A087905 a(n) = n! * Sum_{d|n} (d/n)^d.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A200473 Irregular triangle read by rows: T(n,k) = number of ways to assign n people to d_k unlabeled groups of equal size (where d_k is the k-th divisor of n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A262280 Number of ways to select a nonempty subset s from an n-set and then partition s into blocks of equal size.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A275389 Number of set partitions of [n] with a strongly unimodal block size list.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

A301924 Regular triangle where T(n,k) is the number of unlabeled k-uniform connected hypergraphs spanning n vertices.

Views

Author

Keywords

Examples

Crossrefs

Programs

PARI

Formula

Extensions

A305552 Number of uniform normal multiset partitions of weight n.

Views

Author

Keywords