A073491 -id:A073491 - cp's OEIS frontend

A268193 Triangle read by rows: T(n,k) (n>=1, k>=0) is the number of partitions of n which have k distinct parts i such that i+1 is also a part.

1, 2, 2, 1, 4, 1, 4, 3, 8, 2, 1, 8, 6, 1, 13, 7, 2, 15, 11, 4, 22, 15, 4, 1, 24, 24, 7, 1, 37, 26, 12, 2, 40, 42, 16, 3, 57, 50, 22, 6, 64, 72, 33, 6, 1, 89, 84, 46, 11, 1, 98, 122, 60, 15, 2, 135, 141, 82, 24, 3, 149, 198, 106, 32, 5, 199, 231, 144, 45, 8, 224, 309, 187, 61, 10, 1

Offset: 1

Emeric Deutsch, Feb 13 2016

T(n,k) = number of partitions of n having k singleton parts other than the largest part. Example: T(5,1) = 3 because we have [4,1'], [3,2'], [2,2,1'] (the counted singletons are marked). These partitions are connected by conjugation to those in the definition.
From Gus Wiseman, Jul 10 2025: (Start)
Also the number of integer partitions of n with k maximal subsequences of consecutive parts not decreasing by 1 (anti-runs). For example, row n = 8 counts partitions with the following anti-runs:
  ((8))                ((3,3),(2))          ((3),(2,2),(1))
  ((4,4))              ((4),(3,1))          ((3),(2),(1,1,1))
  ((5,3))              ((5,2),(1))
  ((6,2))              ((4,2),(1,1))
  ((7,1))              ((2,2,2),(1,1))
  ((4,2,2))            ((2,2),(1,1,1,1))
  ((6,1,1))            ((2),(1,1,1,1,1,1))
  ((2,2,2,2))
  ((3,3,1,1))
  ((5,1,1,1))
  ((4,1,1,1,1))
  ((3,1,1,1,1,1))
  ((1,1,1,1,1,1,1,1))
(End)

			T(5,1) = 3 because we have [3,2], [2,2,1], and [2,1,1,1].
T(9,2) = 4 because we have [3,2',1,1,1,1'], [3,2,2',1,1'], [3,3,2',1'], and [4,3',2'] (the i's are marked).
Triangle starts:
  1;
  2;
  2,1;
  4,1;
  4,3;
  8,2,1;
  8,6,1;
From _Gus Wiseman_, Jul 11 2025: (Start)
Row n = 8 counts the following partitions by number of singleton parts other than the largest part:
  (8)                (5,3)        (4,3,1)
  (4,4)              (6,2)        (5,2,1)
  (4,2,2)            (7,1)
  (6,1,1)            (3,3,2)
  (2,2,2,2)          (3,2,2,1)
  (3,3,1,1)          (4,2,1,1)
  (5,1,1,1)          (3,2,1,1,1)
  (2,2,2,1,1)
  (4,1,1,1,1)
  (2,2,1,1,1,1)
  (3,1,1,1,1,1)
  (2,1,1,1,1,1,1)
  (1,1,1,1,1,1,1,1)
(End)

Alois P. Heinz, Rows n = 1..800, flattened

Row sums are A000041.
Row lengths are A003056.
For distinct parts instead of anti-runs we have A116608.
Column k = 1 is A116931.
For runs instead of anti-runs we have A384881.
The strict case is A384905.
The corresponding rank statistic is A356228, non-strict version A384906.
The proper case is A385814, runs A385815.
A007690 counts partitions with no singletons, complement A183558.
A034296 counts flat or gapless partitions, ranks A066311 or A073491.
Cf. A000009, A008284, A024786, A047993, A098859, A116674, A287170, A325325, A384882, A385213.

g := add(x^j*mul(1+t*x^i+x^(2*i)/(1-x^i), i = 1 .. j-1)/(1-x^j), j = 1 .. 80): gser := simplify(series(g, x = 0, 27)): for n from 0 to 25 do P[n] := sort(coeff(gser, x, n)) end do: for n to 25 do seq(coeff(P[n], t, k), k = 0 .. degree(P[n])) end do; # yields sequence in triangular form
# second Maple program:
b:= proc(n, i, t) option remember; expand(`if`(n=0, 1,
      `if`(i<1, 0, add(b(n-i*j, i-1, t or j>0)*
      `if`(t and j=1, x, 1), j=0..n/i))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=0..degree(p)))(b(n$2, false)):
seq(T(n), n=1..20);  # Alois P. Heinz, Feb 13 2016

b[n_, i_, t_] := b[n, i, t] = Expand[If[n == 0, 1, If[i < 1, 0, Sum[b[n - i*j, i - 1, t || j > 0]*If[t && j == 1, x, 1], {j, 0, n/i}]]]]; T[n_] := Function[p, Table[Coefficient[p, x, i], {i, 0, Exponent[p, x]}]][b[n, n, False]]; Table[T[n], {n, 1, 20}] // Flatten (* Jean-François Alcover, Dec 21 2016, after Alois P. Heinz *)
Table[Length[Select[IntegerPartitions[n],Length[Split[#,#1!=#2+1&]]==k&]],{n,0,10},{k,0,n}] (* Delete zeros for A268193. Gus Wiseman, Jul 10 2025 *)

T(n,0) = A116931(n).
Sum_{k>=1} T(n, k) = A000041(n) (the partition numbers).
Sum_{k>=1} k*T(n,k) = A024786(n-1).
G.f.: G(t,x) = Sum_{j>=1} ((x^j/(1-x^j))*Product_{i=1..j-1} (1 + tx^i + x^{2i}/(1-x^i))).

A356224 Number of divisors of n whose prime indices cover an initial interval of positive integers.

Original entry on oeis.org

1, 2, 1, 3, 1, 3, 1, 4, 1, 2, 1, 5, 1, 2, 1, 5, 1, 4, 1, 3, 1, 2, 1, 7, 1, 2, 1, 3, 1, 4, 1, 6, 1, 2, 1, 7, 1, 2, 1, 4, 1, 3, 1, 3, 1, 2, 1, 9, 1, 2, 1, 3, 1, 5, 1, 4, 1, 2, 1, 7, 1, 2, 1, 7, 1, 3, 1, 3, 1, 2, 1, 10, 1, 2, 1, 3, 1, 3, 1, 5, 1, 2, 1, 5, 1, 2, 1

Offset: 1

Gus Wiseman, Aug 04 2022

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.

			The a(n) gapless divisors of n = 1..24:
  1  2  1  4  1  6  1  8  1  2  1  12  1  2  1  16  1  18  1  4  1  2  1  24
     1     2     2     4     1     6      1     8      6      2     1     12
           1     1     2           4            4      2      1           8
                       1           2            2      1                  6
                                   1            1                         4
                                                                          2
                                                                          1
For example, the divisors of 12 are {1,2,3,4,6,12}, of which {1,2,4,6,12} belong to A055932, so a(12) = 5.

These divisors belong to A055932, a subset of A073491 (complement A073492).
The complement is A356225.
A001223 lists the prime gaps.
A328338 has third-largest divisor prime.
A356226 gives the lengths of maximal gapless intervals of prime indices.
Cf. A000005, A001222, A028334, A029709, A055874, A056239, A070824, A112798, A119313, A137921, A287170, A289509, A356223.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
normQ[m_]:=m=={}||Union[m]==Range[Max[m]];
Table[Length[Select[Divisors[n],normQ[primeMS[#]]&]],{n,100}]

A116674 Triangle read by rows: T(n,k) is the number of partitions of n into odd parts and having exactly k distinct parts (n>=1, k>=1).

Original entry on oeis.org

1, 1, 2, 1, 1, 2, 1, 2, 2, 2, 3, 1, 5, 3, 4, 1, 2, 7, 1, 2, 8, 2, 2, 10, 3, 2, 11, 5, 2, 13, 7, 4, 12, 11, 1, 19, 11, 1, 2, 18, 17, 1, 3, 20, 21, 2, 2, 22, 27, 3, 2, 25, 32, 5, 4, 24, 41, 7, 2, 30, 46, 11, 2, 31, 56, 15, 2, 36, 62, 22, 3, 33, 80, 25, 1, 2, 39, 87, 36, 1, 4, 38, 103, 45, 2, 2, 45

Offset: 1

Emeric Deutsch, Feb 22 2006

Row n has floor(sqrt(n)) terms. Row sums yield A000009. T(n,1)=A001227(n) (n>=1). Sum(k*T(n,k),k>=1)=A038348(n-1) (n>=1).

Conjecture: Also the number of strict integer partitions of n with k maximal runs of consecutive parts decreasing by 1. - Gus Wiseman, Jun 24 2025

			From _Gus Wiseman_, Jun 24 2025: (Start)
Triangle begins:
   1:  1
   2:  1
   3:  2
   4:  1  1
   5:  2  1
   6:  2  2
   7:  2  3
   8:  1  5
   9:  3  4  1
  10:  2  7  1
  11:  2  8  2
  12:  2 10  3
  13:  2 11  5
  14:  2 13  7
  15:  4 12 11
  16:  1 19 11  1
  17:  2 18 17  1
  18:  3 20 21  2
  19:  2 22 27  3
  20:  2 25 32  5
Row n = 9 counts the following partitions into odd parts by number of distinct parts:
  (9)                  (7,1,1)          (5,3,1)
  (3,3,3)              (3,3,1,1,1)
  (1,1,1,1,1,1,1,1,1)  (5,1,1,1,1)
                       (3,1,1,1,1,1,1)
Row n = 9 counts the following strict partitions by number of maximal runs:
  (9)      (6,3)    (5,3,1)
  (5,4)    (7,2)
  (4,3,2)  (8,1)
           (6,2,1)
(End)

Alois P. Heinz, Rows n = 1..1000, flattened

Row sums are A000009, strict case of A000041.
Row lengths are A000196.
Leading terms are A001227.
A007690 counts partitions with no singletons, complement A183558.
A034296 counts flat partitions, ranks A066311 or A073491.
A047993 counts partitions with max part = length.
A152140 counts partitions into odd parts by length.
A268193 counts partitions by number of maximal anti-runs, strict A384905.
A384881 counts partitions by number of maximal runs.
Cf. A008284, A038348, A047966, A089259, A325325, A384178, A384886, A384887.

g:=product(1+t*x^(2*j-1)/(1-x^(2*j-1)),j=1..35): gser:=simplify(series(g,x=0,34)): for n from 1 to 29 do P[n]:=coeff(gser,x^n) od: for n from 1 to 29 do seq(coeff(P[n],t,j),j=1..floor(sqrt(n))) od; # yields sequence in triangular form
# second Maple program:
with(numtheory):
b:= proc(n, i) option remember; expand(`if`(n=0, 1,
      `if`(i<1, 0, add(b(n-i*j, i-2)*`if`(j=0, 1, x), j=0..n/i))))
    end:
T:= n-> (p-> seq(coeff(p, x, i), i=1..degree(p)))(
         b(n, iquo(n+1, 2)*2-1)):
seq(T(n), n=1..30);  # Alois P. Heinz, Mar 08 2015

b[n_, i_] := b[n, i] = Expand[If[n == 0, 1, If[i<1, 0, Sum[b[n-i*j, i-2]*If[j == 0, 1, x], {j, 0, n/i}]]]]; T[n_] := Function[{p}, Table[Coefficient[p, x, i], {i, 1, Exponent[p, x]}]][b[n, Quotient[n+1, 2]*2-1]]; Table[T[n], {n, 1, 30}] // Flatten (* Jean-François Alcover, May 22 2015, after Alois P. Heinz *)
Table[Length[Select[IntegerPartitions[n],OddQ[Times@@#]&&Length[Union[#]]==k&]],{n,1,12},{k,1,Floor[Sqrt[n]]}] (*  Gus Wiseman, Jun 24 2025 *)

G.f.: product(1+tx^(2j-1)/(1-x^(2j-1)), j=1..infinity).

A356233 Number of integer factorizations of n into gapless numbers (A066311).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 1, 1, 4, 1, 1, 2, 5, 1, 4, 1, 2, 1, 1, 1, 7, 2, 1, 3, 2, 1, 4, 1, 7, 1, 1, 2, 9, 1, 1, 1, 3, 1, 2, 1, 2, 4, 1, 1, 12, 2, 2, 1, 2, 1, 7, 1, 3, 1, 1, 1, 8, 1, 1, 2, 11, 1, 2, 1, 2, 1, 2, 1, 16, 1, 1, 4, 2, 2, 2, 1, 5, 5, 1, 1, 4, 1, 1

Offset: 1

Gus Wiseman, Aug 28 2022

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. We define a number to be gapless (listed by A066311) iff its prime indices cover an interval of positive integers.

			The counted factorizations of n = 2, 4, 8, 12, 24, 36, 48:
  (2)  (4)    (8)      (12)     (24)       (36)       (48)
       (2*2)  (2*4)    (2*6)    (3*8)      (4*9)      (6*8)
              (2*2*2)  (3*4)    (4*6)      (6*6)      (2*24)
                       (2*2*3)  (2*12)     (2*18)     (3*16)
                                (2*2*6)    (3*12)     (4*12)
                                (2*3*4)    (2*2*9)    (2*3*8)
                                (2*2*2*3)  (2*3*6)    (2*4*6)
                                           (3*3*4)    (3*4*4)
                                           (2*2*3*3)  (2*2*12)
                                                      (2*2*2*6)
                                                      (2*2*3*4)
                                                      (2*2*2*2*3)

The shortest of these factorizations is listed at A356234, length A287170.
A000005 counts divisors.
A001055 counts factorizations.
A001221 counts distinct prime factors, sum A001414.
A003963 multiplies together the prime indices.
A132747 counts non-isolated divisors, complement A132881.
A356069 counts gapless divisors, initial A356224 (complement A356225).
A356226 lists the lengths of maximal gapless submultisets of prime indices:
- length: A287170
- minimum: A356227
- maximum: A356228
- bisected length: A356229
- standard composition: A356230
- Heinz number: A356231
- positions of first appearances: A356232
Cf. A001222, A060680-A060683, A073491-A073495, A193829, A328195, A328335-A328458.

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

A107428 Number of gap-free compositions of n.

Original entry on oeis.org

1, 2, 4, 6, 11, 21, 39, 71, 141, 276, 542, 1070, 2110, 4189, 8351, 16618, 33134, 66129, 131937, 263483, 526453, 1051984, 2102582, 4203177, 8403116, 16800894, 33593742, 67174863, 134328816, 268624026, 537192064, 1074288649, 2148414285, 4296543181, 8592585289

Offset: 1

N. J. A. Sloane, May 26 2005

nonn

A gap-free composition contains all the parts between its smallest and largest part. a(5)=11 because we have: 5, 3+2, 2+3, 2+2+1, 2+1+2, 1+2+2, 2+1+1+1, 1+2+1+1, 1+1+2+1, 1+1+1+2, 1+1+1+1+1. - Geoffrey Critzer, Apr 13 2014

			From _Gus Wiseman_, Oct 04 2022: (Start)
The a(0) = 1 through a(5) = 11 gap-free compositions:
  ()  (1)  (2)   (3)    (4)     (5)
           (11)  (12)   (22)    (23)
                 (21)   (112)   (32)
                 (111)  (121)   (122)
                        (211)   (212)
                        (1111)  (221)
                                (1112)
                                (1121)
                                (1211)
                                (2111)
                                (11111)
(End)

Alois P. Heinz, Table of n, a(n) for n = 1..1000
Alois P. Heinz, Plot of (a(n)-2^(n-2))/2^(n-2) for n = 42..1000
P. Hitczenko and A. Knopfmacher, Gap-free compositions and gap-free samples of geometric random variables, Discrete Math., 294 (2005), 225-239.

The unordered version (partitions) is A034296, ranked by A073491.
The initial case is A107429, unordered A000009, ranked by A333217.
The unordered complement is counted by A239955, ranked by A073492.
These compositions are ranked by A356841.
The complement is counted by A356846, ranked by A356842
A356230 ranks gapless factorization lengths, firsts A356603.
A356233 counts factorizations into gapless numbers.
Cf. A055932, A073493, A137921, A251729, A356224, A356843, A356845.

b:= proc(n, i, t) option remember; `if`(n=0, t!,
      `if`(i<1 or n add(b(n, i, 0), i=1..n):
seq(a(n), n=1..40);  # Alois P. Heinz, Apr 14 2014

Table[Length[Select[Level[Map[Permutations,IntegerPartitions[n]],{2}],Length[Union[#]]==Max[#]-Min[#]+1&]],{n,1,20}] (* Geoffrey Critzer, Apr 13 2014 *)
b[n_, i_, t_] := b[n, i, t] = If[n == 0, t!, If[i < 1 || n < i, 0, Sum[b[n - i*j, i - 1, t + j]/j!, {j, 1, n/i}]]]; a[n_] := Sum[b[n, i, 0], {i, 1, n}]; Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Aug 30 2016, after Alois P. Heinz *)

a(n) ~ 2^(n-2). - Alois P. Heinz, Dec 07 2014

G.f.: Sum_{j>0} Sum_{k>=j} C({j..k},x) where C({s},x) = Sum_{i in {s}} (C({s}-{i},x)*x^i)/(1 - Sum_{i in {s}} (x^i)) is the g.f. for compositions such that the set of parts equals {s} with C({},x) = 1. - John Tyler Rascoe, Jun 01 2024

More terms from Vladeta Jovovic, May 26 2005

A133808 Numbers that are primally tight, have 2 as first prime and weakly ascending powers.

Original entry on oeis.org

1, 2, 4, 6, 8, 16, 18, 30, 32, 36, 54, 64, 108, 128, 150, 162, 210, 216, 256, 324, 450, 486, 512, 648, 750, 900, 972, 1024, 1296, 1458, 1470, 1944, 2048, 2250, 2310, 2916, 3750, 3888, 4096, 4374, 4500, 5832, 6750, 7350, 7776, 8192, 8748, 10290, 11250

Offset: 1

Olivier Gérard, Sep 23 2007

nonn

All numbers of the form 2^k1*p_2^k2*...*p_n^k_n, where k1 <= k2 <= ... <= k_n and the p_i are the n first primes.

Subset of A073491, A133810.

			10 = 2*5 with missing prime factor 3 between 2 and 5 is not in the sequence.
12 = 2^2*3 with 2's exponent > 3's exponent is not in the sequence.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A025487, A087980, A073491, A133809-A133813.

Cf. A027748, A124010, A151800, A000040.

import Data.Set (singleton, deleteFindMin, insert)
a133808 n = a133808_list !! (n-1)
a133808_list = 1 : f (singleton (2, 2, 1)) where
   f s = y : f (insert (y * p, p, e + 1) $ insert (y * q^e, q, e) s')
             where q = a151800 p
                   ((y, p, e), s') = deleteFindMin s
-- Reinhard Zumkeller, Apr 13 2015

isok(n) = {my(f = factor(n)); my(nbf = #f~); if (prod(i=1, nbf, prime(i)) ! = prod(i=1, nbf, f[i, 1]), return (0)); for (j=2, nbf, if (f[j,2] < f[j-1,2], return (0));); return (1);} \\ Michel Marcus, Jun 04 2014

A356226 Irregular triangle giving the lengths of maximal gapless submultisets of the prime indices of n.

Original entry on oeis.org

1, 1, 2, 1, 2, 1, 3, 2, 1, 1, 1, 3, 1, 1, 1, 2, 4, 1, 3, 1, 2, 1, 1, 1, 1, 1, 1, 4, 2, 1, 1, 3, 2, 1, 1, 3, 1, 5, 1, 1, 1, 1, 2, 4, 1, 1, 1, 1, 1, 3, 1, 1, 2, 1, 1, 2, 1, 3, 1, 1, 1, 5, 2, 1, 2, 1, 1, 2, 1, 1, 4, 1, 1, 3, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 2, 1, 6

Offset: 1

Gus Wiseman, Aug 10 2022

A sequence is gapless if it covers an unbroken interval of positive integers. For example, the multiset {2,3,5,5,6,9} has three maximal gapless submultisets: {2,3}, {5,5,6}, {9}.

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.

			Triangle  begins: {}, {1}, {1}, {2}, {1}, {2}, {1}, {3}, {2}, {1,1}, {1}, {3}, {1}, {1,1}, {2}, {4}, {1}, {3}, {1}, {2,1}, ... For example, the prime indices of 20 are {1,1,3}, which separates into maximal gapless submultisets {{1,1},{3}}, so row 20 is (2,1).
The prime indices of 18564 are {1,1,2,4,6,7}, which separates into {1,1,2}, {4}, {6,7}, so row 18564 is (3,1,2). This corresponds to the factorization 18564 = 12 * 7 * 221.

Row sums are A001222.
Singleton row positions are A073491, complement A073492.
Length-2,3,4 row positions are A073493-A073495.
Row lengths are A287170, firsts A066205.
Row minima are A356227.
Row maxima are A356228.
Bisected run-lengths are A356229.
Standard composition numbers of rows are A356230.
Heinz numbers of rows are A356231.
Positions of first appearances are A356232.
A001221 counts distinct prime factors, with sum A001414.
A001223 lists the prime gaps, reduced A028334.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A132747 counts non-isolated divisors, complement A132881.
A356069 counts gapless divisors, initial A356224 (complement A356225).
Cf. A000005, A055874, A060680-A060683, A137921, A193829, A286470, A328166.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Length/@Split[primeMS[n],#1>=#2-1&],{n,100}]

A325160 Products of distinct, non-consecutive primes. Squarefree numbers not divisible by any two consecutive primes.

Original entry on oeis.org

1, 2, 3, 5, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 31, 33, 34, 37, 38, 39, 41, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 67, 69, 71, 73, 74, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109, 110, 111, 113, 115, 118, 119

Offset: 1

Gus Wiseman, Apr 05 2019

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions into distinct non-consecutive parts (counted by A003114). The nonsquarefree case is A319630, which gives the Heinz numbers of integer partitions with no consecutive parts (counted by A116931).

The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 6, 52, 515, 5146, 51435, 514416, 5144232, 51442384, ... . Apparently, the asymptotic density of this sequence exists and equals 0.51442... . - Amiram Eldar, Sep 24 2022

			The sequence of terms together with their prime indices begins:
   1: {}
   2: {1}
   3: {2}
   5: {3}
   7: {4}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  17: {7}
  19: {8}
  21: {2,4}
  22: {1,5}
  23: {9}
  26: {1,6}
  29: {10}
  31: {11}
  33: {2,5}
  34: {1,7}
  37: {12}

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001227, A003114, A005117, A025157, A034296, A056239, A073485, A073491, A089995, A112798, A116931, A319630, A325161, A325162.

Select[Range[100],Min@@Differences[Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]>1&]

isok(k) = {if (issquarefree(k), my(v = apply(primepi, factor(k)[,1])); ! #select(x->(v[x+1]-v[x] == 1), [1..#v-1]));} \\ Michel Marcus, Jan 09 2021

A356225 Number of divisors of n whose prime indices do not cover an initial interval of positive integers.

Original entry on oeis.org

0, 0, 1, 0, 1, 1, 1, 0, 2, 2, 1, 1, 1, 2, 3, 0, 1, 2, 1, 3, 3, 2, 1, 1, 2, 2, 3, 3, 1, 4, 1, 0, 3, 2, 3, 2, 1, 2, 3, 4, 1, 5, 1, 3, 5, 2, 1, 1, 2, 4, 3, 3, 1, 3, 3, 4, 3, 2, 1, 5, 1, 2, 5, 0, 3, 5, 1, 3, 3, 6, 1, 2, 1, 2, 5, 3, 3, 5, 1, 5, 4, 2, 1, 7, 3, 2, 3

Offset: 1

Gus Wiseman, Aug 13 2022

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.

			The a(70) = 6 divisors: 5, 7, 10, 14, 35, 70.

These divisors belong to the complement of A055932, a subset of A073491.
These divisors belong to A080259, a superset of A073492.
The complement is counted by A356224.
A001223 lists the prime gaps.
A328338 has third-largest divisor prime, smallest A119313.
A356226 gives the lengths of maximal gapless intervals of prime indices.
Cf. A000005, A001222, A028334, A055874, A056239, A070824, A112798, A137921, A287170, A356233-A356237.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
normQ[m_]:=m=={}||Union[m]==Range[Max[m]];
Table[Length[Select[Divisors[n],!normQ[primeMS[#]]&]],{n,100}]

a(n) = A000005(n) - A356224(n).

A365830 Heinz numbers of incomplete integer partitions, meaning not every number from 0 to A056239(n) is the sum of some submultiset.

Original entry on oeis.org

3, 5, 7, 9, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 27, 28, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77, 78, 79, 81, 82, 83, 85, 86, 87, 88, 89

Offset: 1

Gus Wiseman, Sep 26 2023

nonn

First differs from A325798 in lacking 156.
The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
The complement (complete partitions) is A325781.

			The terms together with their prime indices begin:
   3: {2}
   5: {3}
   7: {4}
   9: {2,2}
  10: {1,3}
  11: {5}
  13: {6}
  14: {1,4}
  15: {2,3}
  17: {7}
  19: {8}
  21: {2,4}
  22: {1,5}
  23: {9}
  25: {3,3}
  26: {1,6}
  27: {2,2,2}
  28: {1,1,4}
For example, the submultisets of (1,1,2,6) (right column) and their sums (left column) are:
   0: ()
   1: (1)
   2: (2)  or (11)
   3: (12)
   4: (112)
   6: (6)
   7: (16)
   8: (26) or (116)
   9: (126)
  10: (1126)
But 5 is missing, so 156 is in the sequence.

For prime indices instead of sums we have A080259, complement of A055932.
The complement is A325781, counted by A126796, strict A188431.
Positions of nonzero terms in A325799, complement A304793.
These partitions are counted by A365924, strict A365831.
A056239 adds up prime indices, row sums of A112798.
A276024 counts positive subset-sums of partitions, strict A284640
A299701 counts distinct subset-sums of prime indices.
A365918 counts distinct non-subset-sums of partitions, strict A365922.
A365923 counts partitions by distinct non-subset-sums, strict A365545.
Cf. A001223, A046663, A073491, A098743, A304792, A365658, A365919.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
nmz[y_]:=Complement[Range[Total[y]],Total/@Subsets[y]];
Select[Range[100],Length[nmz[prix[#]]]>0&]

cp's OEIS Frontend

A268193 Triangle read by rows: T(n,k) (n>=1, k>=0) is the number of partitions of n which have k distinct parts i such that i+1 is also a part.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A356224 Number of divisors of n whose prime indices cover an initial interval of positive integers.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A116674 Triangle read by rows: T(n,k) is the number of partitions of n into odd parts and having exactly k distinct parts (n>=1, k>=1).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A356233 Number of integer factorizations of n into gapless numbers (A066311).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A107428 Number of gap-free compositions of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A133808 Numbers that are primally tight, have 2 as first prime and weakly ascending powers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

PARI

A356226 Irregular triangle giving the lengths of maximal gapless submultisets of the prime indices of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A325160 Products of distinct, non-consecutive primes. Squarefree numbers not divisible by any two consecutive primes.

Views

Author

Keywords