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-10 of 12 results. Next

A064174 Number of partitions of n with nonnegative rank.

Original entry on oeis.org

1, 1, 2, 3, 4, 6, 9, 12, 17, 23, 31, 42, 56, 73, 96, 125, 161, 207, 265, 336, 426, 536, 672, 840, 1046, 1296, 1603, 1975, 2425, 2970, 3628, 4417, 5367, 6503, 7861, 9482, 11412, 13702, 16423, 19642, 23447, 27938, 33231, 39453, 46767, 55342, 65386, 77135
Offset: 1

Views

Author

Vladeta Jovovic, Sep 20 2001

Keywords

Comments

The rank of a partition is the largest summand minus the number of summands.
This sequence (up to proof) equals "partitions of 2n with even number of parts, ending in 1, with max descent of 1, where the number of odd parts in odd places equals the number of odd parts in even places. (See link and 2nd Mathematica line.) - Wouter Meeussen, Mar 29 2013
Number of partitions p of n such that max(max(p), number of parts of p) is a part of p. - Clark Kimberling, Feb 28 2014
From Gus Wiseman, Mar 09 2019: (Start)
Also the number of integer partitions of n with maximum part greater than or equal to the number of parts. The Heinz numbers of these integer partitions are given by A324521. For example, the a(1) = 1 through a(8) = 12 partitions are:
(1) (2) (3) (4) (5) (6) (7) (8)
(21) (22) (32) (33) (43) (44)
(31) (41) (42) (52) (53)
(311) (51) (61) (62)
(321) (322) (71)
(411) (331) (332)
(421) (422)
(511) (431)
(4111) (521)
(611)
(4211)
(5111)
Also the number of integer partitions of n with maximum part less than or equal to the number of parts. The Heinz numbers of these integer partitions are given by A324562. For example, the a(1) = 1 through a(8) = 12 partitions are:
(1) (11) (21) (22) (221) (222) (322) (332)
(111) (211) (311) (321) (331) (2222)
(1111) (2111) (2211) (2221) (3221)
(11111) (3111) (3211) (3311)
(21111) (4111) (4211)
(111111) (22111) (22211)
(31111) (32111)
(211111) (41111)
(1111111) (221111)
(311111)
(2111111)
(11111111)
(End)

Examples

			a(20) = p(19) - p(15) + p(8) = 490 - 176 + 22 = 336.
		

Crossrefs

Programs

  • Maple
    f:= n -> add((-1)^(k+1)*combinat:-numbpart(n-(3*k^2-k)/2),k=1..floor((1+sqrt(24*n+1))/6)):
    map(f, [$1..100]); # Robert Israel, Aug 03 2015
  • Mathematica
    Table[Count[IntegerPartitions[n], q_ /; First[q] >= Length[q]], {n, 16}]
    (* also *)
    Table[Count[IntegerPartitions[2n],q_/;Last[q]===1 && Max[q-PadRight[Rest[q],Length[q]]]<=1 && Count[First/@Partition[q,2],?OddQ]==Count[Last/@Partition[q,2],?OddQ]],{n,16}]
    (* also *)
    Table[Count[IntegerPartitions[n], p_ /; MemberQ[p, Max[Max[p], Length[p]]]], {n, 50}] (* Clark Kimberling, Feb 28 2014 *)
  • PARI
    {a(n) = my(A=1); A = sum(m=0,n,x^m*prod(k=1,m,(1-x^(m+k-1))/(1-x^k +x*O(x^n)))); polcoeff(A,n)}
    for(n=1,60,print1(a(n),", ")) \\ Paul D. Hanna, Aug 03 2015
    
  • PARI
    my(N=50, x='x+O('x^N)); Vec(1/prod(k=1, N, 1-x^k)*sum(k=1, N, (-1)^(k-1)*x^(k*(3*k-1)/2))) \\ Seiichi Manyama, May 21 2023

Formula

a(n) = (A000041(n) + A047993(n))/2.
a(n) = p(n-1) - p(n-5) + p(n-12) - ... -(-1)^k*p(n-(3*k^2-k)/2) + ..., where p() is A000041(). - Vladeta Jovovic, Aug 04 2004
G.f.: Sum_{n>=1} x^n * Product_{k=1..n} (1 - x^(n+k-1))/(1 - x^k). - Paul D. Hanna, Aug 03 2015
A064173(n) + a(n) = A000041(n). - R. J. Mathar, Feb 22 2023
G.f.: (1/Product_{k>=1} (1-x^k)) * Sum_{k>=1} (-1)^(k-1) * x^(k*(3*k-1)/2). - Seiichi Manyama, May 21 2023

Extensions

Mathematica programs modified by Clark Kimberling, Feb 12 2014

A324522 Numbers > 1 where the minimum prime index is equal to the number of prime factors counted with multiplicity.

Original entry on oeis.org

2, 9, 15, 21, 33, 39, 51, 57, 69, 87, 93, 111, 123, 125, 129, 141, 159, 175, 177, 183, 201, 213, 219, 237, 245, 249, 267, 275, 291, 303, 309, 321, 325, 327, 339, 381, 385, 393, 411, 417, 425, 447, 453, 455, 471, 475, 489, 501, 519, 537, 543, 573, 575, 579, 591
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2019

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.
Also Heinz numbers of integer partitions where the minimum part is equal to the number of parts (A006141). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
    2: {1}
    9: {2,2}
   15: {2,3}
   21: {2,4}
   33: {2,5}
   39: {2,6}
   51: {2,7}
   57: {2,8}
   69: {2,9}
   87: {2,10}
   93: {2,11}
  111: {2,12}
  123: {2,13}
  125: {3,3,3}
  129: {2,14}
  141: {2,15}
  159: {2,16}
  175: {3,3,4}
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    q:= n-> is(pi(min(factorset(n)))=bigomega(n)):
    select(q, [$2..600])[];  # Alois P. Heinz, Mar 07 2019
  • Mathematica
    Select[Range[2,100],PrimePi[FactorInteger[#][[1,1]]]==PrimeOmega[#]&]

Formula

A055396(a(n)) = A001222(a(n)).

A324521 Numbers > 1 where the maximum prime index is less than or equal to the number of prime factors counted with multiplicity.

Original entry on oeis.org

2, 4, 6, 8, 9, 12, 16, 18, 20, 24, 27, 30, 32, 36, 40, 45, 48, 50, 54, 56, 60, 64, 72, 75, 80, 81, 84, 90, 96, 100, 108, 112, 120, 125, 126, 128, 135, 140, 144, 150, 160, 162, 168, 176, 180, 189, 192, 196, 200, 210, 216, 224, 225, 240, 243, 250, 252, 256
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2019

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.
Also Heinz numbers of integer partitions with nonnegative rank (A064174). The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   2: {1}
   4: {1,1}
   6: {1,2}
   8: {1,1,1}
   9: {2,2}
  12: {1,1,2}
  16: {1,1,1,1}
  18: {1,2,2}
  20: {1,1,3}
  24: {1,1,1,2}
  27: {2,2,2}
  30: {1,2,3}
  32: {1,1,1,1,1}
  36: {1,1,2,2}
  40: {1,1,1,3}
  45: {2,2,3}
  48: {1,1,1,1,2}
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    q:= n-> is(pi(max(factorset(n)))<=bigomega(n)):
    select(q, [$2..300])[];  # Alois P. Heinz, Mar 07 2019
  • Mathematica
    Select[Range[2,100],PrimePi[FactorInteger[#][[-1,1]]]<=PrimeOmega[#]&]
  • PARI
    isok(m) = (m>1) && (primepi(vecmax(factor(m)[, 1])) <= bigomega(m)); \\ Michel Marcus, Nov 14 2022
    
  • Python
    from sympy import factorint, primepi
    def ok(n):
        f = factorint(n)
        return primepi(max(f)) <= sum(f.values())
    print([k for k in range(2, 257) if ok(k)]) # Michael S. Branicky, Nov 15 2022

Formula

A061395(a(n)) <= A001222(a(n)).

A039900 Number of partitions satisfying 0 < cn(0,5) + cn(1,5) + cn(4,5).

Original entry on oeis.org

0, 1, 1, 2, 4, 6, 9, 13, 19, 27, 38, 52, 71, 95, 127, 167, 220, 285, 370, 474, 607, 770, 976, 1226, 1540, 1920, 2391, 2960, 3660, 4501, 5529, 6760, 8254, 10038, 12190, 14750, 17825, 21470, 25825, 30975, 37101, 44322, 52879, 62937, 74811, 88733, 105110, 124261
Offset: 0

Views

Author

Keywords

Comments

For a given partition cn(i,n) means the number of its parts equal to i modulo n.
Short: o < 0 + 1 + 4 (OMZAAp).
Number of partitions of n such that (greatest part) >= (multiplicity of greatest part), for n >= 1. For example, a(6) counts these 9 partitions: 6, 51, 42, 411, 33, 321, 3111, 22111, 21111. See the Mathematica program at A240057 for the sequence as a count of these partitions, along with counts of related partitions. - Clark Kimberling, Apr 02 2014
The Heinz numbers of these integer partitions are given by A324561. - Gus Wiseman, Mar 09 2019
From Gus Wiseman, Mar 09 2019: (Start)
Also the number of integer partitions of n whose minimum part is less than or equal to the number of parts. The Heinz numbers of these integer partitions are given by A324560. For example, the a(1) = 1 through a(7) = 13 integer partitions are:
(1) (11) (21) (22) (32) (42) (52)
(111) (31) (41) (51) (61)
(211) (221) (222) (322)
(1111) (311) (321) (331)
(2111) (411) (421)
(11111) (2211) (511)
(3111) (2221)
(21111) (3211)
(111111) (4111)
(22111)
(31111)
(211111)
(1111111)
(End)

Examples

			From _Gus Wiseman_, Mar 09 2019: (Start)
The a(1) = 1 through a(7) = 13 integer partitions with at least one part equal to 0, 1, or 4 modulo 5:
  (1)  (11)  (21)   (4)     (5)      (6)       (43)
             (111)  (31)    (41)     (42)      (52)
                    (211)   (221)    (51)      (61)
                    (1111)  (311)    (321)     (331)
                            (2111)   (411)     (421)
                            (11111)  (2211)    (511)
                                     (3111)    (2221)
                                     (21111)   (3211)
                                     (111111)  (4111)
                                               (22111)
                                               (31111)
                                               (211111)
                                               (1111111)
(End)
		

Crossrefs

Programs

  • Maple
    b:= proc(n, i, t) option remember; `if`(n=0, t,
          `if`(i<1, 0, b(n, i-1, t)+ `if`(i>n, 0, b(n-i, i,
          `if`(irem(i, 5) in {2, 3}, t, 1)))))
        end:
    a:= n-> b(n$2, 0):
    seq(a(n), n=0..50);  # Alois P. Heinz, Apr 03 2014
  • Mathematica
    Table[Count[IntegerPartitions[n], p_ /; Min[p] <= Length[p]], {n, 40}] (* Clark Kimberling, Feb 13 2014 *)
    b[n_, i_, t_] := b[n, i, t] = If[n==0, t, If[i<1, 0, b[n, i-1, t] + If[i > n, 0, b[n-i, i, If[MemberQ[{2, 3}, Mod[i, 5]], t, 1]]]]]; a[n_] := b[n, n, 0]; Table[a[n], {n, 0, 50}] (* Jean-François Alcover, Nov 16 2015, after Alois P. Heinz *)
  • PARI
    my(N=66, x='x+O('x^N)); concat(0, Vec(sum(k=0, N, x^k*(1-x^k^2)/prod(j=1, k, 1-x^j)))) \\ Seiichi Manyama, Jan 13 2022

Formula

G.f.: Sum_{k>=0} x^k * (1-x^(k^2)) / Product_{j=1..k} (1-x^j). - Seiichi Manyama, Jan 13 2022
a(n) = A000041(n) - A003106(n). - Vaclav Kotesovec, Oct 20 2024

A324562 Numbers > 1 where the maximum prime index is greater than or equal to the number of prime factors counted with multiplicity.

Original entry on oeis.org

2, 3, 5, 6, 7, 9, 10, 11, 13, 14, 15, 17, 19, 20, 21, 22, 23, 25, 26, 28, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 56, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 82, 83, 84, 85
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2019

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.
Also Heinz numbers of the integer partitions enumerated by A064174. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   2: {1}
   3: {2}
   5: {3}
   6: {1,2}
   7: {4}
   9: {2,2}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  15: {2,3}
  17: {7}
  19: {8}
  20: {1,1,3}
  21: {2,4}
  22: {1,5}
  23: {9}
  25: {3,3}
  26: {1,6}
  28: {1,1,4}
		

Crossrefs

Programs

  • Maple
    with(numtheory):
    q:= n-> is(pi(max(factorset(n)))>=bigomega(n)):
    select(q, [$2..100])[];  # Alois P. Heinz, Mar 07 2019
  • Mathematica
    Select[Range[2,100],PrimePi[FactorInteger[#][[-1,1]]]>=PrimeOmega[#]&]

Formula

A061395(a(n)) >= A001222(a(n)).

A324517 Numbers > 1 where the maximum prime index equals the number of prime factors minus the number of distinct prime factors.

Original entry on oeis.org

4, 24, 27, 36, 54, 80, 200, 224, 240, 360, 405, 500, 540, 600, 625, 672, 675, 704, 784, 810, 900, 1008, 1120, 1125, 1250, 1350, 1500, 1512, 1664, 1701, 1875, 2112, 2250, 2268, 2352, 2744, 2800, 3168, 3360, 3402, 3520, 3528, 3750, 3872, 3920, 3969, 4352, 4752
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2019

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.
Also Heinz numbers of the integer partitions enumerated by A324518. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
    4: {1,1}
   24: {1,1,1,2}
   27: {2,2,2}
   36: {1,1,2,2}
   54: {1,2,2,2}
   80: {1,1,1,1,3}
  200: {1,1,1,3,3}
  224: {1,1,1,1,1,4}
  240: {1,1,1,1,2,3}
  360: {1,1,1,2,2,3}
  405: {2,2,2,2,3}
  500: {1,1,3,3,3}
  540: {1,1,2,2,2,3}
  600: {1,1,1,2,3,3}
  625: {3,3,3,3}
  672: {1,1,1,1,1,2,4}
  675: {2,2,2,3,3}
  704: {1,1,1,1,1,1,5}
  784: {1,1,1,1,4,4}
  810: {1,2,2,2,2,3}
  900: {1,1,2,2,3,3}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[2,1000],With[{f=FactorInteger[#]},PrimePi[f[[-1,1]]]==Total[Last/@f]-Length[f]]&]

Formula

A061395(a(n)) = A001222(a(n)) - A001221(a(n)) = A046660(a(n)).

A324519 Numbers > 1 where the minimum prime index equals the number of prime factors minus the number of distinct prime factors.

Original entry on oeis.org

4, 12, 18, 20, 27, 28, 44, 50, 52, 60, 68, 76, 84, 90, 92, 98, 116, 124, 126, 132, 135, 140, 148, 150, 156, 164, 172, 188, 189, 198, 204, 212, 220, 225, 228, 234, 236, 242, 244, 260, 268, 276, 284, 292, 294, 297, 306, 308, 316, 332, 338, 340, 342, 348, 350
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2019

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.
Also Heinz numbers of the integer partitions enumerated by A324520. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   4: {1,1}
  12: {1,1,2}
  18: {1,2,2}
  20: {1,1,3}
  27: {2,2,2}
  28: {1,1,4}
  44: {1,1,5}
  50: {1,3,3}
  52: {1,1,6}
  60: {1,1,2,3}
  68: {1,1,7}
  76: {1,1,8}
  84: {1,1,2,4}
  90: {1,2,2,3}
  92: {1,1,9}
  98: {1,4,4}
		

Crossrefs

Programs

  • Mathematica
    Select[Range[2,100],With[{f=FactorInteger[#]},PrimePi[f[[1,1]]]==Total[Last/@f]-Length[f]]&]

Formula

A055396(a(n)) = A001222(a(n)) - A001221(a(n)) = A046660(a(n)).

A324515 Numbers > 1 where the maximum prime index minus the minimum prime index equals the number of prime factors minus the number of distinct prime factors.

Original entry on oeis.org

2, 3, 5, 7, 11, 12, 13, 17, 18, 19, 23, 29, 31, 37, 40, 41, 43, 45, 47, 53, 59, 61, 67, 71, 73, 75, 79, 83, 89, 97, 100, 101, 103, 107, 109, 112, 113, 120, 127, 131, 137, 139, 149, 151, 157, 163, 167, 173, 175, 179, 180, 181, 189, 191, 193, 197, 199, 211, 223
Offset: 1

Views

Author

Gus Wiseman, Mar 06 2019

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.
Also Heinz numbers of the integer partitions enumerated by A324516. The Heinz number of an integer partition (y_1, ..., y_k) is prime(y_1) * ... * prime(y_k).

Examples

			The sequence of terms together with their prime indices begins:
   2: {1}
   3: {2}
   5: {3}
   7: {4}
  11: {5}
  12: {1,1,2}
  13: {6}
  17: {7}
  18: {1,2,2}
  19: {8}
  23: {9}
  29: {10}
  31: {11}
  37: {12}
  40: {1,1,1,3}
  41: {13}
  43: {14}
  45: {2,2,3}
		

Crossrefs

Programs

  • Maple
    filter:= proc(n) local F, Inds, t;
      if isprime(n) then return true fi;
      F:= ifactors(n)[2];
      Inds:= map(numtheory:-pi, F[..,1]);
      max(Inds) - min(Inds) = add(t[2],t=F) - nops(F)
    end proc:
    select(filter, [$2..300]); # Robert Israel, Nov 19 2023
  • Mathematica
    Select[Range[2,100],With[{f=FactorInteger[#]},PrimePi[f[[-1,1]]]-PrimePi[f[[1,1]]]==Total[Last/@f]-Length[f]]&]

Formula

A243055(a(n)) = A061395(a(n)) - A055396(a(n)) = A001222(a(n)) - A001221(a(n)) = A046660(a(n)).

A361204 Positive integers k such that 2*omega(k) <= bigomega(k).

Original entry on oeis.org

1, 4, 8, 9, 16, 24, 25, 27, 32, 36, 40, 48, 49, 54, 56, 64, 72, 80, 81, 88, 96, 100, 104, 108, 112, 121, 125, 128, 135, 136, 144, 152, 160, 162, 169, 176, 184, 189, 192, 196, 200, 208, 216, 224, 225, 232, 240, 243, 248, 250, 256, 272, 288, 289, 296, 297, 304
Offset: 1

Views

Author

Gus Wiseman, Mar 14 2023

Keywords

Examples

			The terms together with their prime indices begin:
     1: {}
     4: {1,1}
     8: {1,1,1}
     9: {2,2}
    16: {1,1,1,1}
    24: {1,1,1,2}
    25: {3,3}
    27: {2,2,2}
    32: {1,1,1,1,1}
    36: {1,1,2,2}
    40: {1,1,1,3}
    48: {1,1,1,1,2}
    49: {4,4}
    54: {1,2,2,2}
    56: {1,1,1,4}
    64: {1,1,1,1,1,1}
		

Crossrefs

These partitions are counted by A237363.
The complement is A361393.
A001221 (omega) counts distinct prime factors.
A001222 (bigomega) counts prime factors.
A112798 lists prime indices, sum A056239.
A360005 gives median of prime indices (times 2), distinct A360457.
Comparing twice the number of distinct parts to the number of parts:
less: A360254, ranks A360558
equal: A239959, ranks A067801
greater: A237365, ranks A361393
less or equal: A237363, ranks A361204
greater or equal: A361394, ranks A361395

Programs

  • Maple
    filter:= proc(n) local F,t;
      F:= ifactors(n)[2];
      add(t[2],t=F) >= 2*nops(F)
    end proc:
    select(filter, [$1..1000]); # Robert Israel, Mar 22 2023
  • Mathematica
    Select[Range[100],2*PrimeNu[#]<=PrimeOmega[#]&]

Formula

A001222(a(n)) >= 2*A001221(a(n)).

A361395 Positive integers k such that 2*omega(k) >= bigomega(k).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74
Offset: 1

Views

Author

Gus Wiseman, Mar 16 2023

Keywords

Comments

Differs from A068938 in having 1 and 4 and lacking 80.
Includes all squarefree numbers.

Examples

			The prime indices of 80 are {1,1,1,1,3}, with 5 parts and 2 distinct parts, and 2*2 < 5, so 80 is not in the sequence.
		

Crossrefs

Complement of A360558.
Positions of nonnegative terms in A361205.
These partitions are counted by A361394.
A001222 (bigomega) counts prime factors, distinct A001221 (omega).
A112798 lists prime indices, sum A056239.
A360005 gives median of prime indices (times 2), distinct A360457.
Comparing twice the number of distinct parts to the number of parts:
less: A360254, ranks A360558
equal: A239959, ranks A067801
greater: A237365, ranks A361393
less or equal: A237363, ranks A361204
greater or equal: A361394, ranks A361395

Programs

  • Mathematica
    Select[Range[100],2*PrimeNu[#]>=PrimeOmega[#]&]

Formula

A001222(a(n)) <= 2*A001221(a(n)).
Showing 1-10 of 12 results. Next