A055874 -id:A055874 - cp's OEIS frontend

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

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

A088725 Numbers having no divisors d>1 such that also d+1 is a divisor.

Original entry on oeis.org

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

Offset: 1

Reinhard Zumkeller, Oct 12 2003

nonn

Complement of A088723.
Union of A132895 and A005408, the odd numbers. - Ray Chandler, May 29 2008
The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 9, 79, 778, 7782, 77813, 778055, 7780548, 77805234, 778052138, 7780519314, ... . Apparently, the asymptotic density of this sequence exists and equals 0.77805... . - Amiram Eldar, Jun 14 2022

			From _Gus Wiseman_, Oct 16 2019: (Start)
The sequence of terms together with their divisors > 1 begins:
   1: {}
   2: {2}
   3: {3}
   4: {2,4}
   5: {5}
   7: {7}
   8: {2,4,8}
   9: {3,9}
  10: {2,5,10}
  11: {11}
  13: {13}
  14: {2,7,14}
  15: {3,5,15}
  16: {2,4,8,16}
  17: {17}
  19: {19}
  21: {3,7,21}
  22: {2,11,22}
  23: {23}
  25: {5,25}
(End)

Positions of 0's and 1's in A129308.
Positions of 0's and 1's in A328457 (also).
Numbers whose divisors (including 1) have no non-singleton runs are A005408.
The number of runs of divisors of n is A137921(n).
The longest run of divisors of n has length A055874(n).
Cf. A000005, A027750, A060680, A088722, A088723 (complement), A088724, A088726, A328166, A328450.

Select[Range[100],FreeQ[Differences[Rest[Divisors[#]]],1]&] (* Harvey P. Dale, Sep 16 2017 *)

isok(n) = {my(d=setminus(divisors(n), [1])); #setintersect(d, apply(x->x+1, d)) == 0;} \\ Michel Marcus, Oct 28 2019

A088722(a(n)) = 0.

Extended by Ray Chandler, May 29 2008

A328166 Heinz number of the run-lengths of the divisors of n.

Original entry on oeis.org

2, 3, 4, 6, 4, 10, 4, 12, 8, 12, 4, 28, 4, 12, 16, 24, 4, 40, 4, 36, 16, 12, 4, 112, 8, 12, 16, 48, 4, 120, 4, 48, 16, 12, 16, 224, 4, 12, 16, 144, 4, 120, 4, 48, 64, 12, 4, 448, 8, 48, 16, 48, 4, 160, 16, 144, 16, 12, 4, 832, 4, 12, 64, 96, 16, 160, 4, 48, 16

Offset: 1

Gus Wiseman, Oct 07 2019

nonn

The Heinz number of an integer partition or multiset {y_1,...,y_k} is prime(y_1)*...*prime(y_k).

			Splitting the divisors of 30 into runs gives {{1, 2, 3}, {5, 6}, {10}, {15}, {30}}, and the Heinz number of {1, 1, 1, 2, 3} is 120, so a(30) = 120.
More examples from _Antti Karttunen_, Dec 09 2021: (Start)
Splitting the divisors of 1 into runs gives {1}, and the Heinz number of that is 2.
Splitting the divisors of 2 into runs gives {1, 2}, and the Heinz number of that is 3. [one run of length 2, therefore a(2) = prime(2)^1].
Splitting the divisors of 3 into runs gives {1} and {3}, and the Heinz number of that is 4. [two runs of length 1, therefore a(3) = prime(1)^2].
Splitting the divisors of 4 into runs gives {1, 2} and {4}, and the Heinz number of that is 6. [one run of length 1, and other run of length 2, therefore a(4) = prime(1)*prime(2)].
Splitting the divisors of 5 into runs gives {1} and {5}, and the Heinz number of that is 4. [two runs of length 1, therefore a(5) = prime(1)^2].
(End)

Antti Karttunen, Table of n, a(n) for n = 1..20000
Wikipedia, Run (cards)
Index entries for sequences related to Heinz numbers

The longest run of divisors of n has length A055874(n).
Numbers whose divisors > 1 have no non-singleton runs are A088725.
The number of successive pairs of divisors of n is A129308(n).
The Heinz number of the set of divisors of n is A275700(n).
Numbers whose divisors do not have weakly decreasing run-lengths are A328165.
Cf. A000005, A056239, A060680, A060775, A119313, A137921, A181063, A199970.

Table[Times@@Prime/@Length/@Split[Divisors[n],#2==#1+1&],{n,30}]

A328166(n) = { my(rl=0,pd=0,v=vector(numdiv(n)),m=1); fordiv(n, d, if(d>(1+pd), v[rl]++; rl=0); pd=d; rl++); v[rl]++; for(i=1,#v, m *= prime(i)^v[i]); (m); }; \\ Antti Karttunen, Dec 09 2021

A001222(a(n)) = A137921(n).

A056239(a(n)) = A000005(n).

A076934 Smallest integer of the form n/k!.

Original entry on oeis.org

1, 1, 3, 2, 5, 1, 7, 4, 9, 5, 11, 2, 13, 7, 15, 8, 17, 3, 19, 10, 21, 11, 23, 1, 25, 13, 27, 14, 29, 5, 31, 16, 33, 17, 35, 6, 37, 19, 39, 20, 41, 7, 43, 22, 45, 23, 47, 2, 49, 25, 51, 26, 53, 9, 55, 28, 57, 29, 59, 10, 61, 31, 63, 32, 65, 11, 67, 34, 69, 35, 71

Offset: 1

Amarnath Murthy, Oct 19 2002

Equivalently, n divided by the largest factorial divisor of n.
Also, the smallest r such that n/r is a factorial number.
Positions of 1's are the factorial numbers A000142. Is every positive integer in this sequence? - Gus Wiseman, May 15 2019
Let m = A055874(n), the largest integer such that 1,2,...,m divides n. Then a(n*m!) = n since m+1 does not divide n, showing that every integer is part of the sequence. - Etienne Dupuis, Sep 19 2020

David A. Corneth, Table of n, a(n) for n = 1..10000

Cf. A000142, A011371, A022559, A055874, A055881, A071626, A111701, A115627, A135291, A229020.

Cf. A325272, A325273, A325508, A325509, A325543, A325544.

Table[n/Max@@Intersection[Divisors[n],Array[Factorial,n]],{n,100}] (* Gus Wiseman, May 15 2019 *)
a[n_] := Module[{k=1}, While[Divisible[n, k!], k++]; n/(k-1)!]; Array[a, 100] (* Amiram Eldar, Dec 25 2023 *)

first(n) = {my(res = [1..n]); for(i = 2, oo, k = i!; if(k <= n, for(j = 1, n\k, res[j*k] = j ) , return(res) ) ) } \\ David A. Corneth, Sep 19 2020

From Amiram Eldar, Dec 25 2023: (Start)
a(n) = n/A055881(n)!.
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = BesselI(2, 2) = 0.688948... (A229020). (End)

More terms from David A. Corneth, Sep 19 2020

A181063 Smallest positive integer with a discrete string of exactly n consecutive divisors, or 0 if no such integer exists.

Original entry on oeis.org

1, 2, 6, 12, 3960, 60, 420, 840, 17907120, 2520, 411863760, 27720, 68502634200, 447069823200, 360360, 720720, 7600186994400, 12252240, 9524356075634400, 81909462250455840, 1149071006394511200, 232792560, 35621201198229847200, 5354228880, 91351145008363640400

Offset: 1

Matthew Vandermast, Oct 07 2010

nonn

The word "discrete" is used to describe a string of consecutive divisors that is not part of a longer such string.
Does a(n) ever equal 0?
a(n) = A003418(n) iff n belongs to A181062; otherwise, a(n) > A003418(n). a(A181062(n)) = A051451(n).

			a(5) = 3960 is divisible by 8, 9, 10, 11, and 12, but not 7 or 13. It is the smallest positive integer with a string of 5 consecutive divisors that is not part of a longer string.
From _Gus Wiseman_, Oct 16 2019: (Start)
The sequence of terms together with their divisors begins:
     1: {1}
     2: {1,2}
     6: {1,2,3,6}
    12: {1,2,3,4,6,12}
  3960: {1,2,...,8,9,10,11,12,...,1980,3960}
    60: {1,2,3,4,5,6,...,30,60}
   420: {1,2,3,4,5,6,7,...,210,420}
   840: {1,2,3,4,5,6,7,8,...,420,840}
(End)

David W. Wilson, Table of n, a(n) for n = 1..1000

The version taking only the longest run is A328449.
The longest run of divisors of n has length A055874(n).
Numbers whose divisors > 1 have no non-singleton runs are A088725.
The number of successive pairs of divisors of n is A129308(n).
Cf. A000005, A003418, A027750, A060775, A181064, A199970, A328166, A328448.

tav=Table[Length/@Split[Divisors[n],#2==#1+1&],{n,10000}];
Table[Position[tav,i][[1,1]],{i,Split[Union@@tav,#2==#1+1&][[1]]}] (* Assumes there are no zeros. - Gus Wiseman, Oct 16 2019 *)

A181062 Prime powers minus 1.

Original entry on oeis.org

0, 1, 2, 3, 4, 6, 7, 8, 10, 12, 15, 16, 18, 22, 24, 26, 28, 30, 31, 36, 40, 42, 46, 48, 52, 58, 60, 63, 66, 70, 72, 78, 80, 82, 88, 96, 100, 102, 106, 108, 112, 120, 124, 126, 127, 130, 136, 138, 148, 150, 156, 162, 166, 168, 172, 178, 180, 190, 192, 196, 198, 210, 222, 226

Offset: 1

Matthew Vandermast, Oct 07 2010

If 0 is excluded, a(n) gives the possible lengths of the longest string of consecutive divisors of a positive integer: range of values of A055874.
a(n) is the largest number m such that A051451(n) = A003418(m).
From Jianing Song, Nov 01 2023: (Start)
Let q = A000961(n) for n > 1. Then:
 - a(n) is the number of units in the finite field F_q.
 - a(n) is the number of solutions to x*y = t for any t != 0 in F_q.
 - If q is odd, then a(n) is also the number of solutions to x^2 - y^2 = t for any t != 0 in F_q. (End)

			Any integer that is divisible by 5 consecutive integers will be divisible by at least 6 consecutive integers. Hence 5 is not in the sequence.

Matthew House, Table of n, a(n) for n = 1..10000

Cf. A000961, A181063, A181064, A181121.
Includes A006093 as a subsequence.
Cf. A003418, A051451, A055874.

Join[{0}, Select[Range@225, PrimePowerQ] - 1] (* Ivan Neretin, Aug 04 2016 *)

isA181062(n)=ispower(n++,,&n);isprime(n) \\ Charles R Greathouse IV, Jun 29 2011

isok(n) = (n==0) || isprimepower(n++); \\ Michel Marcus, Aug 05 2016

a(n) = A000961(n)-1.

Entry revised by N. J. A. Sloane, Jan 06 2013

A356069 Number of divisors of n whose prime indices cover an interval of positive integers (A073491).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 28 2022

nonn

First differs from A000005 at 10, 14, 20, 21, 22, ... = A307516.

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) counted divisors of n = 1, 2, 4, 6, 12, 16, 24, 30, 36, 48, 72, 90:
  1   2   4   6  12  16  24  30  36  48  72  90
      1   2   3   6   8  12  15  18  24  36  45
          1   2   4   4   8   6  12  16  24  30
              1   3   2   6   5   9  12  18  18
                  2   1   4   3   6   8  12  15
                  1       3   2   4   6   9   9
                          2   1   3   4   8   6
                          1       2   3   6   5
                                  1   2   4   3
                                      1   3   2
                                          2   1
                                          1

These divisors belong to A073491, a superset of A055932, complement A073492.
The initial case is A356224.
The complement in the initial case is counted by A356225.
A000005 counts divisors.
A001223 lists the prime gaps.
A056239 adds up prime indices, row sums of A112798, lengths A001222.
A328338 has third-largest divisor prime.
A356226 gives the lengths of maximal gapless intervals of prime indices.
Cf. A028334, A029709, A055874, A070824, A119313, A137921, A287170, A307516, A356223, A356233-A356237.

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

A356228 Greatest size of a gapless submultiset of the prime indices of n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 13 2022

nonn

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 intervals: {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.

			The prime indices of 700 are {1,1,3,3,4}, with maximal gapless submultisets {1,1}, {3,3,4}, so a(700) = 3.
The prime indices of 18564 are {1,1,2,4,6,7}, with maximal gapless submultisets {1,1,2}, {4}, {6,7}, so a(18564) = 3.

Positions of first appearances are A000079.
The maximal gapless submultisets are counted by A287170, firsts A066205.
These are the row-maxima of A356226, firsts A356232.
The smallest instead of greatest size is A356227.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A001223 lists the prime gaps, reduced A028334.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gapless prime indices, cf. A073492-A073495.
A356069 counts gapless divisors.
A356224 counts even gapless divisors, complement A356225.
Cf. A000005, A055874, A060680-A060683, A132747, A132881, A137921, A193829, A286470, A328162, A328457, A356229.

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

a(n) = A333766(A356230(n)).

a(n) = A061395(A356231(n)).

A055926 Numbers k such that {largest m such that 1, 2, ..., m divide k} is different from {largest m such that m! divides k}; numbers k which are either odd multiples of 12 or the largest m such that (m-1)! divides k is a composite number > 5.

Original entry on oeis.org

12, 36, 60, 84, 108, 120, 132, 156, 180, 204, 228, 240, 252, 276, 300, 324, 348, 360, 372, 396, 420, 444, 468, 480, 492, 516, 540, 564, 588, 600, 612, 636, 660, 684, 708, 732, 756, 780, 804, 828, 840, 852, 876, 900, 924, 948, 960, 972, 996, 1020, 1044, 1068

Offset: 1

Leroy Quet, Jul 16 2000

From Antti Karttunen, Nov 20 - Dec 06 2013: (Start)
This sequence has several interpretations:
Numbers k such that A055874(k) differs from A055881(k). [Leroy Quet's original definition of the sequence. Note that A055874(k) >= A055881(k) for all k.]
Numbers k such that {largest m such that m! divides k^2} is different from {largest m such that m! divides k}, i.e., numbers k for which A232098(k) > A055881(k).
Numbers k which are either 12 times an odd number (A073762) or the largest m such that (m-1)! divides k is a composite number > 5 (A232743).
Please see my attached notes for the proof of the equivalence of these interpretations.
Additional implications based on that proof:
A232099 is a subset of this sequence.
A055881(a(n))+1 is always composite. In the range n = 1..17712, only values 4, 6, 8, 9 and 10 occur.
The new definition can be also rephrased by saying that the sequence contains all the positive integers k whose factorial base representation of (A007623(k)) either ends as '...200' (in which case k is an odd multiple of 12, 12 = '200', 36 = '1200', 60 = '2200', ...) or the number of trailing zeros + 2 in that representation is a composite number greater than or equal to 6, e.g. 120 = '10000' (in other words, A055881(k) is one of the terms of A072668 after the initial 3). Together these conditions also imply that all the terms are divisible by 12.
(End)

			12 is included because 3! is the largest factorial to divide 12, but 1, 2, 3 and 4 all divide 12. Equally, 12 is included because it is one of the terms of A073762, or equally, because its factorial base representation ends with digits '...200': A007623(12) = 200.
840 (= 3*5*7*8) is included because the largest factorial which divides 840 is 5! (840 = 7*120), but all positive integers up to 8 divide 840. Equally, 840 is included because it is one of the terms of A232743 as 5+1 = 6 is a composite number larger than 5. Note that A007623(840) = 110000.

Antti Karttunen, Table of n, a(n) for n = 1..17712
Antti Karttunen, A proof for the equivalence of three alternative definitions of A055926
Wikipedia, Wilson's theorem (See the section "Composite modulus")

Union of A073762 and A232743. Equivalently, setwise difference of A232742 and A017593. Subset: A232099.

Cf. A055874, A055881, A072668, A232098, A232100, A232741, A232744, A232745.

More terms from Antti Karttunen, Dec 01 2013

A356227 Smallest size of a maximal gapless submultiset of the prime indices of n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 13 2022

nonn

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.

			The prime indices of 18564 are {1,1,2,4,6,7}, with maximal gapless submultisets {1,1,2}, {4}, {6,7}, so a(18564) = 1.

Positions of first appearances are A000079.
The maximal gapless submultisets are counted by A287170, firsts A066205.
These are the row-minima of A356226, firsts A356232.
The greatest instead of smallest size is A356228.
A001221 counts distinct prime factors, with sum A001414.
A001222 counts prime factors with multiplicity.
A001223 lists the prime gaps, reduced A028334.
A003963 multiplies together the prime indices of n.
A056239 adds up prime indices, row sums of A112798.
A073491 lists numbers with gapless prime indices, cf. A073492-A073495.
A356224 counts even gapless divisors, complement A356225.
Cf. A000005, A055874, A060680-A060683, A132747, A132881, A137921, A193829, A286470, A356229.

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

a(n) = A333768(A356230(n)).

a(n) = A055396(A356231(n)).

cp's OEIS Frontend

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A088725 Numbers having no divisors d>1 such that also d+1 is a divisor.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A328166 Heinz number of the run-lengths of the divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A076934 Smallest integer of the form n/k!.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A181063 Smallest positive integer with a discrete string of exactly n consecutive divisors, or 0 if no such integer exists.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A181062 Prime powers minus 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A356069 Number of divisors of n whose prime indices cover an interval of positive integers (A073491).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica