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 54 results. Next

A360015 Numbers whose exponent of 2 in their canonical prime factorization is equal to the maximal exponent.

Original entry on oeis.org

1, 2, 4, 6, 8, 10, 12, 14, 16, 20, 22, 24, 26, 28, 30, 32, 34, 36, 38, 40, 42, 44, 46, 48, 52, 56, 58, 60, 62, 64, 66, 68, 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 92, 94, 96, 100, 102, 104, 106, 110, 112, 114, 116, 118, 120, 122, 124, 128, 130, 132, 134, 136, 138
Offset: 1

Views

Author

Amiram Eldar, Jan 21 2023

Keywords

Comments

Numbers k such that A007814(k) = A051903(k).
The powers of 2 (A000079) are all terms.
The product of any two terms (not necessarily distinct) is also a term.
This sequence is a disjoint union of {1} and the subsequences of numbers m of the form 2^(k-1)*o where o = A000265(m), the odd part of m, is a k-free number, for k >= 2. These subsequences include, for k = 2, numbers of the form 2*o where o is an odd squarefree number (A056911); for k = 3, numbers of the form 4*o where o is an odd cubefree number; etc.
The asymptotic density of this sequence is Sum_{k>=2} 1/(zeta(k)*(2^k-1)) = 0.44541445377638761933... .
The asymptotic mean of the exponent of 2 in the prime factorization of the terms of this sequence is Sum_{k>=2} (k-1)/(zeta(k)*(2^k-1)) / Sum_{k>=2} 1/(zeta(k)*(2^k-1)) = 2.10346728882748723133... . [corrected by Amiram Eldar, Jul 10 2025]
Also numbers whose multiset of prime factors has low (i.e. least) mode 2. Here, a mode in a multiset is an element that appears at least as many times as each of the others; for example, the modes in {a,a,b,b,b,c,d,d,d} are {b,d}. - Gus Wiseman, Jul 14 2023

Examples

			From _Gus Wiseman_, Jul 14 2023: (Start)
108 = 2*2*3*3*3 is missing because its mode is not 2.
180 = 2*2*3*3*5 is present because it has low mode 2.
The terms together with their prime factorizations begin:
   1 =
   2 = 2
   4 = 2*2
   6 = 2*3
   8 = 2*2*2
  10 = 2*5
  12 = 2*2*3
  14 = 2*7
  16 = 2*2*2*2
  20 = 2*2*5
  22 = 2*11
  24 = 2*2*2*3
  26 = 2*13
  28 = 2*2*7
  30 = 2*3*5
  32 = 2*2*2*2*2
  34 = 2*17
  36 = 2*2*3*3
(End)
		

Crossrefs

Partitions of this type are counted by A241131.
The case of unique mode is A360013, complement here A360014.
For unique minimal prime exponent we have A364061, counted by A364062.
For minimal prime exponent we have A364158, counted by A364159.
A027746 lists prime factors (with multiplicity).
A112798 lists prime indices, length A001222, sum A056239.
A362611 counts modes in prime indices, triangle A362614.

Programs

  • Mathematica
    q[n_] := IntegerExponent[n, 2] == Max[FactorInteger[n][[;; , 2]]]; q[1] = True; Select[Range[150], q]
  • PARI
    is(n) = n == 1 || vecmax(factor(n)[,2]) == valuation(n, 2);

Formula

Disjoint union of A360013 and A360014.
a(n) = A360013(n)/2. - Gus Wiseman, Jul 14 2023

A362606 Numbers without a unique least prime exponent, or numbers whose prime factorization has more than one element of least multiplicity.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 26, 30, 33, 34, 35, 36, 38, 39, 42, 46, 51, 55, 57, 58, 60, 62, 65, 66, 69, 70, 74, 77, 78, 82, 84, 85, 86, 87, 90, 91, 93, 94, 95, 100, 102, 105, 106, 110, 111, 114, 115, 118, 119, 120, 122, 123, 126, 129, 130, 132, 133, 134, 138, 140
Offset: 1

Views

Author

Gus Wiseman, May 05 2023

Keywords

Comments

First differs from A130092 in lacking 180.
First differs from A351295 in lacking 180 and having 216.
First differs from A362605 in having 60.

Examples

			The prime factorization of 1800 is {2,2,2,3,3,5,5}, and the parts of least multiplicity are {3,5}, so 1800 is in the sequence.
The terms together with their prime indices begin:
    6: {1,2}
   10: {1,3}
   14: {1,4}
   15: {2,3}
   21: {2,4}
   22: {1,5}
   26: {1,6}
   30: {1,2,3}
   33: {2,5}
   34: {1,7}
   35: {3,4}
   36: {1,1,2,2}
   38: {1,8}
   39: {2,6}
   42: {1,2,4}
		

Crossrefs

The complement is A359178, counted by A362610.
For mode we have A362605, counted by A362607.
Partitions of this type are counted by A362609.
These are the positions of terms > 1 in A362613.
A112798 lists prime indices, length A001222, sum A056239.
A362614 counts partitions by number of modes.
A362615 counts partitions by number of co-modes.

Programs

  • Mathematica
    Select[Range[100],Count[Last/@FactorInteger[#],Min@@Last/@FactorInteger[#]]>1&]

A363486 Low mode in the multiset of prime indices of n.

Original entry on oeis.org

0, 1, 2, 1, 3, 1, 4, 1, 2, 1, 5, 1, 6, 1, 2, 1, 7, 2, 8, 1, 2, 1, 9, 1, 3, 1, 2, 1, 10, 1, 11, 1, 2, 1, 3, 1, 12, 1, 2, 1, 13, 1, 14, 1, 2, 1, 15, 1, 4, 3, 2, 1, 16, 2, 3, 1, 2, 1, 17, 1, 18, 1, 2, 1, 3, 1, 19, 1, 2, 1, 20, 1, 21, 1, 3, 1, 4, 1, 22, 1, 2, 1
Offset: 1

Views

Author

Gus Wiseman, Jun 23 2023

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.
A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes in {a,a,b,b,b,c,d,d,d} are {b,d}.
Extending the terminology of A124943, the "low mode" in a multiset is its least mode.

Crossrefs

Positions of first appearances are 1 and A000040.
Positions of 1's are A360013, counted by A241131.
For greatest instead of least we have A363487.
The version for median is A363941, triangle A124943.
The high version for median is A363942, triangle A124944.
The version for mean instead of mode is A363943, high A363944.
A112798 lists prime indices, length A001222, sum A056239.
A326567/A326568 gives mean of prime indices.
A359178 ranks partitions with a unique co-mode, counted by A362610.
A356862 ranks partitions with a unique mode, counted by A362608.
A362605 ranks partitions with more than one mode, counted by A362607.
A362606 ranks partitions with more than one co-mode, counted by A362609.
A362611 counts modes in prime indices, triangle A362614.
A362613 counts co-modes in prime indices, triangle A362615.
A362616 ranks partitions (max part) = (unique mode), counted by A362612.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    modes[ms_]:=Select[Union[ms],Count[ms,#]>=Max@@Length/@Split[ms]&];
    Table[If[n==1,0,First[modes[prix[n]]]],{n,30}]

A363727 Numbers whose prime indices satisfy (mean) = (median) = (mode), assuming there is a unique mode.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 19, 23, 25, 27, 29, 31, 32, 37, 41, 43, 47, 49, 53, 59, 61, 64, 67, 71, 73, 79, 81, 83, 89, 90, 97, 101, 103, 107, 109, 113, 121, 125, 127, 128, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199
Offset: 1

Views

Author

Gus Wiseman, Jun 23 2023

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.
A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes in {a,a,b,b,b,c,d,d,d} are {b,d}.
The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

Examples

			The terms together with their prime indices begin:
     2: {1}          29: {10}              79: {22}
     3: {2}          31: {11}              81: {2,2,2,2}
     4: {1,1}        32: {1,1,1,1,1}       83: {23}
     5: {3}          37: {12}              89: {24}
     7: {4}          41: {13}              90: {1,2,2,3}
     8: {1,1,1}      43: {14}              97: {25}
     9: {2,2}        47: {15}             101: {26}
    11: {5}          49: {4,4}            103: {27}
    13: {6}          53: {16}             107: {28}
    16: {1,1,1,1}    59: {17}             109: {29}
    17: {7}          61: {18}             113: {30}
    19: {8}          64: {1,1,1,1,1,1}    121: {5,5}
    23: {9}          67: {19}             125: {3,3,3}
    25: {3,3}        71: {20}             127: {31}
    27: {2,2,2}      73: {21}             128: {1,1,1,1,1,1,1}
		

Crossrefs

These partitions are counted by A363719, factorizations A363741.
For unequal instead of equal we have A363730, counted by A363720.
Excluding primes gives A363722.
Excluding prime-powers gives A363729, counted by A363728.
A112798 lists prime indices, length A001222, sum A056239.
A326567/A326568 gives mean of prime indices.
A356862 ranks partitions with a unique mode, counted by A362608.
A359178 ranks partitions with multiple modes, counted by A362610.
A360005 gives twice the median of prime indices.
A362611 counts modes in prime indices, triangle A362614.
A362613 counts co-modes in prime indices, triangle A362615.
A363486 gives least mode in prime indices, A363487 greatest.
Just two statistics:
- (mean) = (median): A359889, counted by A240219.
- (mean) != (median): A359890, counted by A359894.
- (mean) = (mode): counted by A363723, see A363724, A363731.
- (median) = (mode): counted by A363740.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    modes[ms_]:=Select[Union[ms],Count[ms,#]>=Max@@Length/@Split[ms]&];
    Select[Range[100],{Mean[prix[#]]}=={Median[prix[#]]}==modes[prix[#]]&]

Formula

Assuming there is a unique mode, we have A326567(a(n))/A326568(a(n)) = A360005(a(n))/2 = A363486(a(n)) = A363487(a(n)).

A360013 Numbers whose exponent of 2 in their canonical prime factorization is larger than all the other exponents.

Original entry on oeis.org

2, 4, 8, 12, 16, 20, 24, 28, 32, 40, 44, 48, 52, 56, 60, 64, 68, 72, 76, 80, 84, 88, 92, 96, 104, 112, 116, 120, 124, 128, 132, 136, 140, 144, 148, 152, 156, 160, 164, 168, 172, 176, 184, 188, 192, 200, 204, 208, 212, 220, 224, 228, 232, 236, 240, 244, 248, 256
Offset: 1

Views

Author

Amiram Eldar, Jan 21 2023

Keywords

Comments

Numbers k such that A007814(k) > A051903(A000265(k)).
The powers of 2 (A000079), except for 1, are all terms.
The product of any two terms (not necessarily distinct) is also a term.
This sequence is a disjoint union of {2} and the subsequences of numbers m of the form 2^k*o where o = A000265(m), the odd part of m, is a k-free number, for k >= 2. These subsequences include, for k = 2, numbers of the form 4*o where o is an odd squarefree number (A056911); for k = 3, numbers of the form 8*o where o is an odd cubefree number; etc.
The asymptotic density of this sequence is Sum_{k>=2} 1/(zeta(k)*2*(2^k-1)) = 0.222707226888193809... .
The asymptotic mean of the exponent of 2 in the prime factorization of the terms of this sequence is Sum_{k>=2} k/(zeta(k)*2*(2^k-1)) / Sum_{k>=2} 1/(zeta(k)*2*(2^k-1)) = 3.10346728882748723133... . [corrected by Amiram Eldar, Jul 10 2025]
This sequence is a subsequence of A360015 and the asymptotic density of this sequence within A360015 is exactly 1/2.
Also even numbers whose multiset of prime factors has unique mode 2. - Gus Wiseman, Jul 10 2023

Examples

			From _Gus Wiseman_, Jul 09 2023: (Start)
108 = 2*2*3*3*3 is missing because its mode is not 2.
180 = 2*2*3*3*5 is missing because 2 is not the unique mode.
120 = 2*2*2*3*5 is present because its unique mode is 2.
The terms together with their prime factorizations begin:
   2 = 2
   4 = 2*2
   8 = 2*2*2
  12 = 2*2*3
  16 = 2*2*2*2
  20 = 2*2*5
  24 = 2*2*2*3
  28 = 2*2*7
  32 = 2*2*2*2*2
  40 = 2*2*2*5
  44 = 2*2*11
  48 = 2*2*2*2*3
  52 = 2*2*13
  56 = 2*2*2*7
  60 = 2*2*3*5
  64 = 2*2*2*2*2*2
(End)
		

Crossrefs

Equals A360015 \ A360014.
Partitions of this type are counted by A241131.
Allowing any unique mode gives A356862, complement A362605.
Allowing any unique co-mode gives A359178, complement A362606.
Not requiring the mode to be unique gives A360015.
The opposite version is A362616, counted by A362612.
For co-mode instead of mode we have A364061, counted by A364062.
With least prime factor instead of 2, we have A364160, counted by A364193.
With a different factorization, we have the subsequence A335738.
A124010 gives prime signature, ordered A118914.
A362611 counts modes in prime factorization, triangle A362614.
A362613 counts co-modes in prime factorization, triangle A362615.
A363486 gives least mode in prime indices, A363487 greatest.

Programs

  • Mathematica
    q[n_] := Module[{e = IntegerExponent[n, 2], m}, m = n/2^e; (m == 1 && e > 0) || AllTrue[FactorInteger[m][[;; , 2]], # < e &]]; Select[Range[256], q]
  • PARI
    is(n) = {my(e = valuation(n, 2), m = n >> e); (m == 1 && e > 0) || (m > 1 && vecmax(factor(m)[,2]) < e)};

Formula

a(n) = 2*A360015(n). - Gus Wiseman, Jul 10 2023

A363487 High mode in the multiset of prime indices of n.

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 4, 1, 2, 3, 5, 1, 6, 4, 3, 1, 7, 2, 8, 1, 4, 5, 9, 1, 3, 6, 2, 1, 10, 3, 11, 1, 5, 7, 4, 2, 12, 8, 6, 1, 13, 4, 14, 1, 2, 9, 15, 1, 4, 3, 7, 1, 16, 2, 5, 1, 8, 10, 17, 1, 18, 11, 2, 1, 6, 5, 19, 1, 9, 4, 20, 1, 21, 12, 3, 1, 5, 6, 22, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Jul 04 2023

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.
A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes in {a,a,b,b,b,c,d,d,d} are {b,d}.
Extending the terminology of A124944, the "high mode" in a multiset is its greatest mode.

Crossrefs

Positions of first appearances are 1 and A000040.
Positions of 1's are A360015, counted by A241131.
For low instead of high mode we have A363486.
The version for low median is A363941, triangle A124943.
The version for high median is A363942, triangle A124944.
The version for mean instead of mode is A363944, low A363943.
A112798 lists prime indices, length A001222, sum A056239.
A326567/A326568 gives mean of prime indices.
A359178 ranks partitions with a unique co-mode, counted by A362610.
A356862 ranks partitions with a unique mode, counted by A362608.
A362605 ranks partitions with more than one mode, counted by A362607.
A362606 ranks partitions with more than one co-mode, counted by A362609.
A362611 counts modes in prime indices, triangle A362614.
A362613 counts co-modes in prime indices, triangle A362615.
A362616 ranks partitions (max part) = (unique mode), counted by A362612.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    modes[ms_]:=Select[Union[ms],Count[ms,#]>=Max@@Length/@Split[ms]&];
    Table[If[n==1,0,Last[modes[prix[n]]]],{n,30}]

A362616 Numbers in whose prime factorization the greatest factor is the unique mode.

Original entry on oeis.org

2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 18, 19, 23, 25, 27, 29, 31, 32, 37, 41, 43, 47, 49, 50, 53, 54, 59, 61, 64, 67, 71, 73, 75, 79, 81, 83, 89, 97, 98, 101, 103, 107, 108, 109, 113, 121, 125, 127, 128, 131, 137, 139, 147, 149, 150, 151, 157, 162, 163, 167
Offset: 1

Views

Author

Gus Wiseman, May 05 2023

Keywords

Comments

First differs from A329131 in lacking 450 and having 1500.
A mode in a multiset is an element that appears at least as many times as each of the others. For example, the modes of {a,a,b,b,b,c,d,d,d} are {b,d}.

Examples

			The factorization of 90 is 2*3*3*5, modes {3}, so 90 is missing.
The factorization of 450 is 2*3*3*5*5, modes {3,5}, so 450 is missing.
The factorization of 900 is 2*2*3*3*5*5, modes {2,3,5}, so 900 is missing.
The factorization of 1500 is 2*2*3*5*5*5, modes {5}, so 1500 is present.
The terms together with their prime indices begin:
     2: {1}          27: {2,2,2}           67: {19}
     3: {2}          29: {10}              71: {20}
     4: {1,1}        31: {11}              73: {21}
     5: {3}          32: {1,1,1,1,1}       75: {2,3,3}
     7: {4}          37: {12}              79: {22}
     8: {1,1,1}      41: {13}              81: {2,2,2,2}
     9: {2,2}        43: {14}              83: {23}
    11: {5}          47: {15}              89: {24}
    13: {6}          49: {4,4}             97: {25}
    16: {1,1,1,1}    50: {1,3,3}           98: {1,4,4}
    17: {7}          53: {16}             101: {26}
    18: {1,2,2}      54: {1,2,2,2}        103: {27}
    19: {8}          59: {17}             107: {28}
    23: {9}          61: {18}             108: {1,1,2,2,2}
    25: {3,3}        64: {1,1,1,1,1,1}    109: {29}
		

Crossrefs

First term with given bigomega is A000079.
For median instead of mode we have A053263.
Partitions of this type are counted by A362612.
A112798 lists prime indices, length A001222, sum A056239.
A356862 ranks partitions with a unique mode, counted by A362608.
A359178 ranks partitions with a unique co-mode, counted by A362610.
A362605 ranks partitions with more than one mode, counted by A362607.
A362606 ranks partitions with more than one co-mode, counted by A362609.
A362614 counts partitions by number of modes, ranked by A362611.
A362615 counts partitions by number of co-modes, ranked by A362613.

Programs

  • Mathematica
    prifacs[n_]:=If[n==1,{},Flatten[ConstantArray@@@FactorInteger[n]]];
    Select[Range[100],Commonest[prifacs[#]]=={Max[prifacs[#]]}&]

A363942 High median in the multiset of prime indices of n.

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 4, 1, 2, 3, 5, 1, 6, 4, 3, 1, 7, 2, 8, 1, 4, 5, 9, 1, 3, 6, 2, 1, 10, 2, 11, 1, 5, 7, 4, 2, 12, 8, 6, 1, 13, 2, 14, 1, 2, 9, 15, 1, 4, 3, 7, 1, 16, 2, 5, 1, 8, 10, 17, 2, 18, 11, 2, 1, 6, 2, 19, 1, 9, 3, 20, 1, 21, 12, 3, 1, 5, 2, 22, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Jul 01 2023

Keywords

Comments

The high median (see A124944) in a multiset is either the middle part (for odd length), or the greatest of the two middle parts (for even length).
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.

Examples

			The prime indices of 90 are {1,2,2,3}, with high median 2, so a(90) = 2.
The prime indices of 150 are {1,2,3,3}, with high median 3, so a(150) = 3.
		

Crossrefs

Positions of first appearances are 1 and A000040.
The triangle for this statistic (high median) is A124944, low A124943.
Regular median of prime indices is A360005(n)/2.
For mode instead of median we have A363487, low A363486.
The low version is A363941.
For mean instead of median we have A363944, triangle A363946, low A363943.
A061395 give maximum prime index, A055396 minimum.
A112798 lists prime indices, length A001222, sum A056239.
A362611 counts modes in prime indices, triangle A362614.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    merr[y_]:=If[Length[y]==0,0, If[OddQ[Length[y]],y[[(Length[y]+1)/2]],y[[1+Length[y]/2]]]];
    Table[merr[prix[n]],{n,100}]

A171979 Number of partitions of n such that smaller parts do not occur more frequently than greater parts.

Original entry on oeis.org

1, 1, 2, 3, 4, 5, 8, 8, 12, 14, 19, 21, 30, 31, 42, 50, 62, 69, 91, 99, 126, 144, 175, 198, 246, 275, 331, 379, 452, 509, 612, 686, 811, 922, 1076, 1219, 1428, 1604, 1863, 2108, 2434, 2739, 3162, 3551, 4075, 4593, 5240, 5885, 6721, 7527, 8556, 9597, 10870
Offset: 0

Views

Author

Reinhard Zumkeller, Jan 20 2010

Keywords

Comments

A000009(n) <= a(n) <= A000041(n).
Equivalently, the number of partitions of n such that (maximal multiplicity of parts) = (multiplicity of the maximal part), as in the Mathematica program. - Clark Kimberling, Apr 04 2014
Also the number of integer partitions of n whose greatest part is a mode, meaning it appears at least as many times as each of the others. The name "Number of partitions of n such that smaller parts do not occur more frequently than greater parts" seems to describe A100882 = "Number of partitions of n in which the sequence of frequencies of the summands is nonincreasing," which first differs from this at n = 10 due to the partition (3,3,2,1,1). - Gus Wiseman, May 07 2023

Examples

			a(5) = #{5, 4+1, 3+2, 2+2+1, 5x1} = 5;
a(6) = #{6, 5+1, 4+2, 3+3, 3+2+1, 2+2+2, 2+2+1+1, 6x1} = 8;
a(7) = #{7, 6+1, 5+2, 4+3, 4+2+1, 3+3+1, 2+2+2+1, 7x1} = 8;
a(8) = #{8, 7+1, 6+2, 5+3, 5+2+1, 4+4, 4+3+1, 3+3+2, 3+3+1+1, 2+2+2+2, 2+2+2+1+1, 8x1} = 12.
		

Crossrefs

For median instead of mode we have A053263.
The complement is counted by A240302.
The case where the maximum is the only mode is A362612.
A000041 counts integer partitions, strict A000009.
A362608 counts partitions with a unique mode, complement A362607.
A362611 counts modes in prime factorization.
A362614 counts partitions by number of modes.

Programs

  • Mathematica
    z = 60; f[n_] := f[n] = IntegerPartitions[n]; m[p_] := Max[Map[Length, Split[p]]]  (* maximal multiplicity *)
    Table[Count[f[n], p_ /; m[p] == Count[p, Max[p]]], {n, 0, z}] (* this sequence *)
    Table[Count[f[n], p_ /; m[p] > Count[p, Max[p]]], {n, 0, z}]  (* A240302 *)
    (* Clark Kimberling, Apr 04 2014 *)
    (* Second program: *)
    b[n_, i_, k_] := b[n, i, k] = If[n == 0, If[k == 0, 1, 0],
         If[i < 1, 0, b[n, i - 1, k] + Sum[b[n - i*j, i - 1,
         If[k == -1, j, If[k == 0, 0, If[j > k, 0, k]]]], {j, 1, n/i}]]];
    a[n_] := PartitionsP[n] - b[n, n, -1];
    a /@ Range[0, 70] (* Jean-François Alcover, Jun 05 2021, after Alois P. Heinz in A240302 *)
    Table[Length[Select[IntegerPartitions[n],MemberQ[Commonest[#],Max[#]]&]],{n,0,30}] (* Gus Wiseman, May 07 2023 *)
  • PARI
    { my(N=53, x='x+O('x^N));
    my(gf=1+sum(i=1,N,sum(j=1,floor(N/i),x^(i*j)*prod(k=1,i-1,(1-x^(k*(j+1)))/(1-x^k)))));
    Vec(gf) } \\ John Tyler Rascoe, Mar 09 2024

Formula

a(n) = p(n,0,1,1) with p(n,i,j,k) = if k<=n then p(n-k,i,j+1,k) +p(n,max(i,j),1,k+1) else (if j0 then 0 else 1).
a(n) + A240302(n) = A000041(n). - Clark Kimberling, Apr 04 2014.
G.f.: 1 + Sum_{i, j>0} x^(i*j) * Product_{k=1..i-1} ((1 - x^(k*(j+1)))/(1 - x^k)). - John Tyler Rascoe, Mar 09 2024

A362621 One and numbers whose multiset of prime factors (with multiplicity) has the same median as maximum.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 16, 17, 18, 19, 23, 25, 27, 29, 31, 32, 37, 41, 43, 47, 49, 50, 53, 54, 59, 61, 64, 67, 71, 73, 75, 79, 81, 83, 89, 97, 98, 101, 103, 107, 108, 109, 113, 121, 125, 127, 128, 131, 137, 139, 147, 149, 151, 157, 162, 163, 167, 169
Offset: 1

Author

Gus Wiseman, May 12 2023

Keywords

Comments

First differs from A334965 in having 750 and lacking 2250.
The median of a multiset is either the middle part (for odd length), or the average of the two middle parts (for even length).

Examples

			The prime factorization of 108 is 2*2*3*3*3, and the multiset {2,2,3,3,3} has median 3 and maximum 3, so 108 is in the sequence.
The prime factorization of 2250 is 2*3*3*5*5*5, and the multiset {2,3,3,5,5,5} has median 4 and maximum 5, so 2250 is not in the sequence.
The terms together with their prime indices begin:
     1: {}           25: {3,3}           64: {1,1,1,1,1,1}
     2: {1}          27: {2,2,2}         67: {19}
     3: {2}          29: {10}            71: {20}
     4: {1,1}        31: {11}            73: {21}
     5: {3}          32: {1,1,1,1,1}     75: {2,3,3}
     7: {4}          37: {12}            79: {22}
     8: {1,1,1}      41: {13}            81: {2,2,2,2}
     9: {2,2}        43: {14}            83: {23}
    11: {5}          47: {15}            89: {24}
    13: {6}          49: {4,4}           97: {25}
    16: {1,1,1,1}    50: {1,3,3}         98: {1,4,4}
    17: {7}          53: {16}           101: {26}
    18: {1,2,2}      54: {1,2,2,2}      103: {27}
    19: {8}          59: {17}           107: {28}
    23: {9}          61: {18}           108: {1,1,2,2,2}
		

Crossrefs

Partitions of this type are counted by A053263.
For mode instead of median we have A362619, counted by A171979.
For parts at middle position (instead of median) we have A362622.
The complement is A362980, counted by A237821.
A027746 lists prime factors, A112798 indices, length A001222, sum A056239.
A362611 counts modes in prime factorization, triangle version A362614.
A362613 counts co-modes in prime factorization, triangle version A362615.

Programs

  • Mathematica
    Select[Range[100],(y=Flatten[Apply[ConstantArray,FactorInteger[#],{1}]];Max@@y==Median[y])&]
Previous Showing 11-20 of 54 results. Next