A060683 -id:A060683 - cp's OEIS frontend

A060681 Largest difference between consecutive divisors of n (ordered by size).

0, 1, 2, 2, 4, 3, 6, 4, 6, 5, 10, 6, 12, 7, 10, 8, 16, 9, 18, 10, 14, 11, 22, 12, 20, 13, 18, 14, 28, 15, 30, 16, 22, 17, 28, 18, 36, 19, 26, 20, 40, 21, 42, 22, 30, 23, 46, 24, 42, 25, 34, 26, 52, 27, 44, 28, 38, 29, 58, 30, 60, 31, 42, 32, 52, 33, 66, 34, 46, 35, 70, 36, 72, 37

Offset: 1

Labos Elemer, Apr 19 2001

Is a(n) the least m > 0 such that n - m divides n! + m? - Clark Kimberling, Jul 28 2012
Is a(n) the least m > 0 such that L(n-m) divides L(n+m), where L = A000032 (Lucas numbers)? - Clark Kimberling, Jul 30 2012
Records give A006093. - Omar E. Pol, Oct 26 2013
Divide n by its smallest prime factor p, then multiply with (p-1), with a(1) = 0 by convention. Compare also to A366387. - Antti Karttunen, Oct 23 2023
a(n) is also the smallest LCM of positive integers x and y where x + y = n. - Felix Huber, Aug 28 2024

			For n = 35, divisors are {1, 5, 7, 35}; differences are {4, 2, 28}; a(35) = largest difference = 28 = 35 - 35/5.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Antal Balog, Paul Erdős, and Gérald Tenenbaum,, On Arithmetic Functions Involving Consecutive Divisors, In: Analytical Number Theory, pp. 77-90, Birkhäuser, Basel, 1990.

Cf. A013661, A020639, A060680, A060682, A060683, A060685, A064097 (number of iterations needed to reach 1).

Cf. also A171462, A366387.

a060681 n = div n p * (p - 1) where p = a020639 n
-- Reinhard Zumkeller, Apr 06 2015

read("transforms") :
A060681 := proc(n)
    if n = 1 then
        0 ;
    else
        sort(convert(numtheory[divisors](n),list)) ;
        DIFF(%) ;
        max(op(%)) ;
    end if;
end proc:
seq(A060681(n),n=1..60) ; # R. J. Mathar, May 23 2018
# second Maple program:
A060681:=n->if(n=1,0,min(map(x->ilcm(x,n-x),[$1..1/2*n]))); seq(A060681(n),n=1..74); # Felix Huber, Aug 28 2024

a[n_ ] := n - n/FactorInteger[n][[1, 1]]
Array[Max[Differences[Divisors[#]]] &, 80, 2] (* Harvey P. Dale, Oct 26 2013 *)

diff(v)=vector(#v-1,i,v[i+1]-v[i])
a(n)=vecmax(diff(divisors(n))) \\ Charles R Greathouse IV, Sep 02 2015

a(n) = if (n==1, 0, n - n/factor(n)[1,1]); \\ Michel Marcus, Oct 24 2015

first(n) = n = max(n, 1); my(res = vector(n)); res[1] = 0; forprime(p = 2, n, for(i = 1, n \ p, if(res[p * i] == 0, res[p * i] = i*(p-1)))); res \\ David A. Corneth, Jan 08 2019

from sympy import primefactors
def A060681(n): return n-n//min(primefactors(n),default=1) # Chai Wah Wu, Jun 21 2023

a(n) = n - n/A020639(n).
a(n) = n - A032742(n). - Omar E. Pol, Aug 31 2011
a(2n) = n, a(3*(2n+1)) = 2*(2n+1) = 4n + 2. - Antti Karttunen, Oct 23 2023
Sum_{k=1..n} a(k) ~ (1/2 - c) * n^2, where c is defined in the corresponding formula in A032742. - Amiram Eldar, Dec 21 2024

Edited by Dean Hickerson, Jan 22 2002

a(1)=0 added by N. J. A. Sloane, Oct 01 2015 at the suggestion of Antti Karttunen

A060680 Smallest difference between consecutive divisors of n.

Original entry on oeis.org

1, 2, 1, 4, 1, 6, 1, 2, 1, 10, 1, 12, 1, 2, 1, 16, 1, 18, 1, 2, 1, 22, 1, 4, 1, 2, 1, 28, 1, 30, 1, 2, 1, 2, 1, 36, 1, 2, 1, 40, 1, 42, 1, 2, 1, 46, 1, 6, 1, 2, 1, 52, 1, 4, 1, 2, 1, 58, 1, 60, 1, 2, 1, 4, 1, 66, 1, 2, 1, 70, 1, 72, 1, 2, 1, 4, 1, 78, 1, 2, 1, 82, 1, 4, 1, 2, 1, 88, 1, 6, 1, 2, 1, 4

Offset: 2

Labos Elemer, Apr 19 2001

nonn

a(n) = 1 if n is even and a(n) is even if n is odd.

a(n) = least m>0 such that n!+1+m and n-m are not relatively prime. - Clark Kimberling, Jul 21 2012

			For n = 35, divisors = {1,5,7,35}; differences = {4,2,28}; a(35) = smallest difference = 2.

Reinhard Zumkeller, Table of n, a(n) for n = 2..10000
Antal Balog, Paul Erdős and Gérald Tenenbaum, On Arithmetic Functions Involving Consecutive Divisors, In: Analytical Number Theory, pp. 77-90, Birkhäuser, Basel, 1990.

Cf. A060681 (largest difference), A060682, A060683, A060684.

Cf. A193829, A027750.

a060680 = minimum . a193829_row  -- Reinhard Zumkeller, Jun 25 2015

read("transforms") :
A060680 := proc(n)
    sort(convert(numtheory[divisors](n),list)) ;
    DIFF(%) ;
    min(op(%)) ;
end proc:
seq(A060680(n),n=2..60) ; # R. J. Mathar, May 23 2018

a[n_] := Min@@(Drop[d=Divisors[n], 1]-Drop[d, -1]);
(* Second program: *)
a[n_] := Min[Differences[Divisors[n]]];
Table[a[n], {n, 2, 100}] (* Jean-François Alcover, Oct 16 2024 *)

a(n) = {my(m = n, d1); fordiv(n, d, if(d > 1 && d - d1 < m, m = d - d1); d1 = d); m;} \\ Amiram Eldar, Mar 17 2025

a(2n+1) = A060684(n).

Corrected by David W. Wilson, May 04 2001

Edited by Dean Hickerson, Jan 22 2002

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}]

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}]

A356230 The a(n)-th composition in standard order is the sequence of lengths of maximal gapless submultisets of the prime indices of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 1, 4, 2, 3, 1, 4, 1, 3, 2, 8, 1, 4, 1, 5, 3, 3, 1, 8, 2, 3, 4, 5, 1, 4, 1, 16, 3, 3, 2, 8, 1, 3, 3, 9, 1, 5, 1, 5, 4, 3, 1, 16, 2, 6, 3, 5, 1, 8, 3, 9, 3, 3, 1, 8, 1, 3, 5, 32, 3, 5, 1, 5, 3, 6, 1, 16, 1, 3, 4, 5, 2, 5, 1, 17, 8, 3, 1, 9, 3

Offset: 1

Gus Wiseman, Aug 16 2022

nonn

The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.
A multiset 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}. These have lengths (3,1,2), which is the 38th composition in standard order, so a(18564) = 38.

Numbers grouped by number of gaps in prime indices are A073491-A073495.
These are the standard composition numbers of rows of A356226.
Using Heinz numbers instead of standard compositions gives A356231.
Positions of first appearances are A356603, sorted A356232.
A001221 counts distinct prime factors, with sum A001414.
A003963 multiplies together the prime indices.
A056239 adds up the prime indices, row sums of A112798.
A066099 lists compositions in standard order.
A132747 counts non-isolated divisors, complement A132881.
A333627 represents the run-lengths of standard compositions.
A356069 counts gapless divisors, initial A356224 (complement A356225).
Cf. A000120, A001222, A060680-A060683, A066205, A286470, A287170, A328166, A356227, A356228, A356229.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
stcinv[q_]:=Total[2^(Accumulate[Reverse[q]])]/2;
Table[stcinv[Length/@Split[primeMS[n],#1>=#2-1&]],{n,100}]

A000120(a(n)) = A287170(n).
A333766(a(n)) = A356228(n).
A333768(a(n)) = A356227(n).

A258409 Greatest common divisor of all (d-1)'s, where the d's are the positive divisors of n.

Original entry on oeis.org

1, 2, 1, 4, 1, 6, 1, 2, 1, 10, 1, 12, 1, 2, 1, 16, 1, 18, 1, 2, 1, 22, 1, 4, 1, 2, 1, 28, 1, 30, 1, 2, 1, 2, 1, 36, 1, 2, 1, 40, 1, 42, 1, 2, 1, 46, 1, 6, 1, 2, 1, 52, 1, 2, 1, 2, 1, 58, 1, 60, 1, 2, 1, 4, 1, 66, 1, 2, 1, 70, 1, 72, 1, 2, 1, 2, 1, 78, 1

Offset: 2

Ivan Neretin, May 29 2015

nonn

a(n) = 1 for even n; a(p) = p-1 for prime p.
a(n) is even for odd n (since all divisors of n are odd).
It appears that a(n) = A052409(A005179(n)), i.e., it is the largest integer power of the smallest number with exactly n divisors. - Michel Marcus, Nov 10 2015
Conjecture: GCD of all (p-1) for prime p|n. - Thomas Ordowski, Sep 14 2016
Conjecture is true, because the set of numbers == 1 (mod g) is closed under multiplication. - Robert Israel, Sep 14 2016
Conjecture: a(n) = A289508(A328023(n)) = GCD of the differences between consecutive divisors of n. See A328163 and A328164. - Gus Wiseman, Oct 16 2019

			65 has divisors 1, 5, 13, and 65, hence a(65) = gcd(1-1,5-1,13-1,65-1) = gcd(0,4,12,64) = 4.

Ivan Neretin, Table of n, a(n) for n = 2..10000

Cf. A049559, A057237, A060680, A063994, A187730, A027750.
Cf. A084190 (similar but with LCM).
Looking at prime indices instead of divisors gives A328167.
Partitions whose parts minus 1 are relatively prime are A328170.
Cf. A000005, A060681, A060683, A193829, A289508.

a258409 n = foldl1 gcd $ map (subtract 1) $ tail $ a027750_row' n
-- Reinhard Zumkeller, Jun 25 2015

f:= n -> igcd(op(map(`-`,numtheory:-factorset(n),-1))):
map(f, [$2..100]); # Robert Israel, Sep 14 2016

Table[GCD @@ (Divisors[n] - 1), {n, 2, 100}]

a(n) = my(g=0); fordiv(n, d, g = gcd(g, d-1)); g; \\ Michel Marcus, May 29 2015

a(n) = gcd(apply(x->x-1, divisors(n))); \\ Michel Marcus, Nov 10 2015

a(n)=if(n%2==0, return(1)); if(n%3==0, return(2)); if(n%5==0 && n%4 != 1, return(2)); gcd(apply(p->p-1, factor(n)[,1])) \\ Charles R Greathouse IV, Sep 19 2016

A060682 Number of distinct differences between consecutive divisors of n (ordered by size).

Original entry on oeis.org

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

Offset: 2

Labos Elemer, Apr 19 2001

nonn

Number of all differences for n is d(n)-1 = A000005(n)-1. Increments are not necessarily different, so a(n)<=d(n)-1.

			For n=70, divisors={1,2,5,7,10,14,35,70}; differences={1,3,2,3,4,21,35}; a(70) = number of distinct differences = 6.

Reinhard Zumkeller, Table of n, a(n) for n = 2..10000
A. Balog, P. Erdős and G. Tenenbaum, On Arithmetic Functions Involving Consecutive Divisors, In: Analytical Number Theory, pp. 77-90, Birkhäuser, Basel, 1990.
Jason Earls, Smarandache iterations of the first kind on functions involving divisors and prime factors, in Smarandache Notions Journal (2004), Vol. 14.1, page 264.

Cf. A000005, A060680, A060681, A060683.

Cf. A193829.

import Data.List (nub, genericLength)
a060682 = genericLength . nub . a193829_row
-- Reinhard Zumkeller, Jun 25 2015

a[n_ ] := Length[Union[Drop[d=Divisors[n], 1]-Drop[d, -1]]]

a(n) = my(d=divisors(n)); #vecsort(vector(#d-1, k, d[k+1] - d[k]),,8); \\ Michel Marcus, Jul 04 2017

Edited by Dean Hickerson, Jan 22 2002

A356231 Heinz number of the sequence (A356226) of lengths of maximal gapless submultisets of the prime indices of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 2, 5, 3, 4, 2, 5, 2, 4, 3, 7, 2, 5, 2, 6, 4, 4, 2, 7, 3, 4, 5, 6, 2, 5, 2, 11, 4, 4, 3, 7, 2, 4, 4, 10, 2, 6, 2, 6, 5, 4, 2, 11, 3, 6, 4, 6, 2, 7, 4, 10, 4, 4, 2, 7, 2, 4, 6, 13, 4, 6, 2, 6, 4, 6, 2, 11, 2, 4, 5, 6, 3, 6, 2, 14, 7, 4, 2, 10

Offset: 1

Gus Wiseman, Aug 18 2022

nonn

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.
A multiset 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}. These have lengths (3,1,2), with Heinz number 30, so a(18564) = 30.

Positions of prime terms are A073491, complement A073492.
Positions of terms with bigomega 2-4 are A073493-A073495.
Applying bigomega gives A287170, firsts A066205, even bisection A356229.
These are the Heinz numbers of the rows of A356226.
Minimal/maximal prime indices are A356227/A356228.
A version for standard compositions is A356230, firsts A356232/A356603.
A001221 counts distinct prime factors, with sum A001414.
A003963 multiplies together the prime indices.
A056239 adds up the prime indices, row sums of A112798.
A132747 counts non-isolated divisors, complement A132881.
A356069 counts gapless divisors, initial A356224 (complement A356225).
Cf. A000005, A001222, A055932, A060680-A060683, A193829, A286470, A328166, A356233-A356237.

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

A001222(a(n)) = A287170(n).
A055396(a(n)) = A356227(n).
A061395(a(n)) = A356228(n).

A356234 Irregular triangle read by rows where row n is the ordered factorization of n into maximal gapless divisors.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 8, 9, 2, 5, 11, 12, 13, 2, 7, 15, 16, 17, 18, 19, 4, 5, 3, 7, 2, 11, 23, 24, 25, 2, 13, 27, 4, 7, 29, 30, 31, 32, 3, 11, 2, 17, 35, 36, 37, 2, 19, 3, 13, 8, 5, 41, 6, 7, 43, 4, 11, 45, 2, 23, 47, 48, 49, 2, 25, 3, 17, 4, 13, 53, 54, 5, 11, 8

Offset: 1

Gus Wiseman, Aug 28 2022

Row-products are the positive integers 1, 2, 3, ...

			The first 16 rows:
   1 =
   2 = 2
   3 = 3
   4 = 4
   5 = 5
   6 = 6
   7 = 7
   8 = 8
   9 = 9
  10 = 2 * 5
  11 = 11
  12 = 12
  13 = 13
  14 = 2 * 7
  15 = 15
  16 = 16
The factorization of 18564 is 18564 = 12*7*221, so row 18564 is {12,7,221}.

Row-lengths are A287170, firsts A066205, even bisection A356229.
Applying bigomega to all parts gives A356226, statistics A356227-A356232.
A001055 counts factorizations.
A001221 counts distinct prime factors, sum A001414.
A003963 multiplies together the prime indices.
A056239 adds up the prime indices, row sums of A112798.
A132747 counts non-isolated divisors, complement A132881.
A356069 counts gapless divisors, initial A356224 (complement A356225).
Cf. A000005, A001222, A060680-A060683, A073491-A073495, A193829, A330103, A356233-A356237.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Times@@Prime/@#&/@Split[primeMS[n],#1>=#2-1&],{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)).

cp's OEIS Frontend

A060681 Largest difference between consecutive divisors of n (ordered by size).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Python

Formula

Extensions

A060680 Smallest difference between consecutive divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

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

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A356230 The a(n)-th composition in standard order is the sequence of lengths of maximal gapless submultisets of the prime indices of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A258409 Greatest common divisor of all (d-1)'s, where the d's are the positive divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

A060682 Number of distinct differences between consecutive divisors of n (ordered by size).

Views

Author