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.

Showing 1-3 of 3 results.

A321469 Number of factorizations of n into factors > 1 with different sums of prime indices. Number of multiset partitions of the multiset of prime indices of n with distinct block-sums.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 2, 2, 1, 3, 1, 3, 2, 2, 1, 4, 1, 2, 2, 3, 1, 4, 1, 3, 2, 2, 2, 5, 1, 2, 2, 4, 1, 5, 1, 3, 3, 2, 1, 6, 1, 3, 2, 3, 1, 5, 2, 5, 2, 2, 1, 7, 1, 2, 2, 4, 2, 5, 1, 3, 2, 4, 1, 8, 1, 2, 3, 3, 2, 5, 1, 6, 2, 2, 1, 7, 2, 2, 2, 5, 1, 7, 2, 3, 2, 2, 2, 8, 1, 3, 3, 5, 1, 5, 1, 5, 5
Offset: 1

Views

Author

Gus Wiseman, Nov 11 2018

Keywords

Comments

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. The sum of prime indices of n is A056239(n).

Examples

			The a(72) = 8 multiset partitions with distinct block-sums:
    {{1,1,1,2,2}}
   {{1},{1,1,2,2}}
   {{2},{1,1,1,2}}
   {{1,1},{1,2,2}}
   {{1,2},{1,1,2}}
   {{2,2},{1,1,1}}
  {{1},{2},{1,1,2}}
  {{1},{1,1},{2,2}}
Missing from this list are:
    {{1},{1},{1,2,2}}
    {{1},{1,2},{1,2}}
    {{2},{2},{1,1,1}}
    {{2},{1,1},{1,2}}
   {{1},{1},{1},{2,2}}
   {{1},{1},{2},{1,2}}
   {{1},{2},{2},{1,1}}
  {{1},{1},{1},{2},{2}}
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    sps[{}]:={{}};sps[set:{i_,_}]:=Join@@Function[s,Prepend[#,s]&/@sps[Complement[set,s]]]/@Cases[Subsets[set],{i,_}];
    mps[set_]:=Union[Sort[Sort/@(#/.x_Integer:>set[[x]])]&/@sps[Range[Length[set]]]];
    Table[Length[Select[mps[primeMS[n]],UnsameQ@@Sort[Total/@#]&]],{n,100}]
  • PARI
    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));
    A321469(n, m=n, facs=List([])) = if(1==n, all_have_different_sum_of_pis(facs), my(s=0, newfacs); fordiv(n, d, if((d>1)&&(d<=m), newfacs = List(facs); listput(newfacs,d); s += A321469(n/d, d, newfacs))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

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

A371733 Maximal length of a factorization of n into factors > 1 all having the same sum of prime indices.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Apr 13 2024

Keywords

Comments

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 all having the same sum of prime indices are counted by A321455.

Examples

			The factorizations of 588 of this type are (7*7*12), (21*28), (588), so a(588) = 3.
The factorizations of 900 of this type are (5*5*6*6), (9*10*10), (25*36), (30*30), (900), so a(900) = 4.
		

Crossrefs

Positions of 1's are A321453, counted by A321451.
Positions of terms > 1 are A321454, counted by A321452.
Factorizations of this type are counted by A321455, different sums A321469.
For different sums instead of same sums we have A371734.
For set partitions of binary indices we have A371735.
A001055 counts factorizations.
A002219 (aerated) counts biquanimous partitions, ranks A357976.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A321142 and A371794 count non-biquanimous strict partitions.
A371789 counts non-quanimous sets, differences A371790.
A371796 counts quanimous sets, differences A371797.

Programs

  • Mathematica
    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],SameQ@@hwt/@#&]],{n,100}]
  • PARI
    A056239(n) = if(1==n, 0, my(f=factor(n)); sum(i=1, #f~, f[i, 2] * primepi(f[i, 1])));
    all_have_same_sum_of_pis(facs) = if(!#facs, 1, (#Set(apply(A056239,facs)) == 1));
    A371733(n, m=n, facs=List([])) = if(1==n, if(all_have_same_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, A371733(n/d, d, newfacs)))); (s)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(108) 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

Views

Author

Gus Wiseman, Apr 14 2024

Keywords

Comments

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).

Examples

			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.
		

Crossrefs

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.

Programs

  • Mathematica
    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}]
Showing 1-3 of 3 results.