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-5 of 5 results.

A137921 Number of divisors d of n such that d+1 is not a divisor of n.

Original entry on oeis.org

1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 2, 3, 2, 3, 4, 4, 2, 4, 2, 4, 4, 3, 2, 5, 3, 3, 4, 5, 2, 5, 2, 5, 4, 3, 4, 6, 2, 3, 4, 6, 2, 5, 2, 5, 6, 3, 2, 7, 3, 5, 4, 5, 2, 6, 4, 6, 4, 3, 2, 7, 2, 3, 6, 6, 4, 6, 2, 5, 4, 7, 2, 8, 2, 3, 6, 5, 4, 6, 2, 8, 5, 3, 2, 8, 4, 3, 4, 7, 2, 8, 4, 5, 4, 3, 4, 9, 2, 5, 6, 7, 2, 6, 2, 7, 8
Offset: 1

Views

Author

Reinhard Zumkeller, Feb 23 2008

Keywords

Comments

a(n) = number of "divisor islands" of n. A divisor island is any set of consecutive divisors of a number where no pairs of consecutive divisors in the set are separated by 2 or more. - Leroy Quet, Feb 07 2010

Examples

			The divisors of 30 are 1,2,3,5,6,10,15,30. The divisor islands are (1,2,3), (5,6), (10), (15), (30). (Note that the differences between consecutive divisors 5-3, 10-6, 15-10 and 30-15 are all > 1.) There are 5 such islands, so a(30)=5.

Links

Crossrefs

Bisections: A099774, A174199.
First appearance of n is at position A173569(n).
Numbers whose divisors have no non-singleton runs are A005408.
The longest run of divisors of n has length A055874(n).
The number of successive pairs of divisors of n is A129308(n).
Cf. A000005, A001620, A027750, A060680, A088723, A088725, A181063, A199970, A328165, A328166, A328448, A328450.

Programs

  • Haskell
    a137921 n = length $ filter (> 0) $
   map ((mod n) . (+ 1)) [d | d <- [1..n], mod n d == 0]
-- Reinhard Zumkeller, Nov 23 2011
  • Maple
    with(numtheory): disl := proc (b) local ct, j: ct := 1: for j to nops(b)-1 do if 2 <= b[j+1]-b[j] then ct := ct+1 else end if end do: ct end proc: seq(disl(divisors(n)), n = 1 .. 120); # Emeric Deutsch, Feb 12 2010
  • Mathematica
    f[n_] := Length@ Split[ Divisors@n, #2 - #1 == 1 &]; Array[f, 105] (* f(n) from Bobby R. Treat *) (* Robert G. Wilson v, Feb 22 2010 *)
Table[Count[Differences[Divisors[n]],?(#>1&)]+1,{n,110}] (* _Harvey P. Dale, Jun 05 2012 *)
a[n_] := DivisorSum[n, Boole[!Divisible[n, #+1]]&]; Array[a, 100] (* Jean-François Alcover, Dec 01 2015 *)
  • PARI
    a(n)=my(d,s=0);if(n%2,numdiv(n),d=divisors(n);for(i=1,#d,if(n%(d[i]+1),s++));s)
  • PARI
    a(n)=sumdiv(n,d,(n%(d+1)!=0)); \\ Joerg Arndt, Jan 06 2015
  • Python
    from sympy import divisors
def A137921(n):
    return len([d for d in divisors(n,generator=True) if n % (d+1)])
# Chai Wah Wu, Jan 05 2015

Formula

a(n) <= A000005(n), with equality iff n is odd; a(A137922(n)) = 2.
a(n) = A000005(n) - A129308(n). - Michel Marcus, Jan 06 2015
a(n) = A001222(A328166(n)). - Gus Wiseman, Oct 16 2019
Sum_{k=1..n} a(k) ~ n * (log(n) + 2*gamma - 2), where gamma is Euler's constant (A001620). - Amiram Eldar, Jan 18 2024

Extensions

Corrected and edited by Charles R Greathouse IV, Apr 19 2010
Edited by N. J. A. Sloane, Aug 10 2010

A129308 a(n) is the number of positive integers k such that k*(k+1) divides n.

Original entry on oeis.org

0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 3, 0, 1, 0, 1, 0, 2, 0, 2, 0, 1, 0, 3, 0, 1, 0, 1, 0, 3, 0, 1, 0, 1, 0, 3, 0, 1, 0, 2, 0, 3, 0, 1, 0, 1, 0, 3, 0, 1, 0, 1, 0, 2, 0, 2, 0, 1, 0, 5, 0, 1, 0, 1, 0, 2, 0, 1, 0, 1, 0, 4, 0, 1, 0, 1, 0, 2, 0, 2, 0, 1, 0, 4, 0, 1, 0, 1, 0, 4, 0, 1, 0, 1, 0, 3, 0, 1, 0, 2, 0, 2, 0, 1, 0
Offset: 1

Views

Author

Leroy Quet, May 26 2007

Keywords

Comments

The usual OEIS policy is not to include sequences like this where alternate terms are zero; this is an exception.
In other words, a(n) is the number of oblong numbers (A002378) dividing n. - Bernard Schott, Jul 29 2022

Examples

			The divisors of 20 are 1,2,4,5,10,20. Of these there are two that are of the form k(k+1): 2 = 1*2 and 20 = 4*5. So a(2) = 2.

Links

Crossrefs

Positions of 0's and 1's are A088725, whose characteristic function is A360128.
First appearance of n is A287142(n), with sorted version A328450.
The longest run of divisors of n has length A055874(n).
One less than A195155.
Cf. A000005, A002378, A003601, A007862, A027750, A033676, A060680, A060681, A072627, A088722, A181063, A199970, A328026, A328165, A328166.
Cf. A005369, A059841, A287142, A344005.

Programs

  • Mathematica
    a = {}; For[n = 1, n < 90, n++, k = 1; co = 0; While[k < Sqrt[n], If[IntegerQ[ n/(k*(k + 1))], co++ ]; k++ ]; AppendTo[a, co]]; a (* Stefan Steinerberger, May 27 2007 *)
Table[Count[Differences[Divisors[n]],1],{n,30}] (* Gus Wiseman, Oct 15 2019 *)
  • PARI
    a(n)=sumdiv(n, d, n%(d+1)==0); \\ Michel Marcus, Jan 06 2015
  • Python
    from itertools import pairwise
from sympy import divisors
def A129308(n): return 0 if n&1 else sum(1 for a, b in pairwise(divisors(n)) if a+1==b) # Chai Wah Wu, Jun 09 2025

Formula

a(2n-1) = 0; a(2n) = A007862(n). - Ray Chandler, Jun 24 2008
G.f.: Sum_{n>=1} x^(n*(n+1))/(1-x^(n*(n+1))). - Joerg Arndt, Jan 30 2011 [modified by Ilya Gutkovskiy, Apr 14 2021]
a(n) = A000005(n) - A137921(n), where A137921(n) is the number of maximal runs of successive divisors of n. - Gus Wiseman, Oct 15 2019
a(n) = Sum_{d|n} A005369(d). - Ridouane Oudra, Jan 22 2021
a(n) = A195155(n)-1. - Antti Karttunen, Feb 21 2023
From Amiram Eldar, Dec 31 2023: (Start)
a(n) = A088722(n) + A059841(n).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1. (End)

Extensions

More terms from Stefan Steinerberger, May 27 2007
Extended by Ray Chandler, Jun 24 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

Views

Author

Gus Wiseman, Oct 07 2019

Keywords

Comments

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

Examples

			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)

Links

Crossrefs

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.

Programs

  • Mathematica
    Table[Times@@Prime/@Length/@Split[Divisors[n],#2==#1+1&],{n,30}]
  • PARI
    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

Formula

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

A328450 Numbers that are a smallest number with k pairs of successive divisors, for some k.

Original entry on oeis.org

1, 2, 6, 12, 60, 72, 180, 360, 420, 840, 1260, 2520, 3780, 5040, 13860, 27720, 36960, 41580, 55440, 83160, 166320, 277200, 360360, 471240, 491400, 720720, 1081080, 1113840, 2162160, 2827440, 3341520, 4324320, 5405400, 6126120
Offset: 0

Views

Author

Gus Wiseman, Oct 15 2019

Keywords

Comments

A sorted version of A287142.

Examples

			The divisors of 72 are {1, 2, 3, 4, 6, 8, 9, 12, 18, 24, 36, 72}, with pairs of successive divisors {{1, 2}, {2, 3}, {3, 4}, {8, 9}}, and no smaller number has 4 successive pairs, so 72 belongs to the sequence.

Crossrefs

Sorted positions of first appearances in A129308.
The longest run of divisors of n has length A055874(n).
Numbers whose divisors > 1 have no non-singleton runs are A088725.
The Heinz number of the multiset of run-lengths of divisors of n is A328166(n).
The smallest number whose divisors have a longest run of length n is A328449(n).
Cf. A000005, A003601, A027750, A033676, A060680, A060681, A072627, A181063, A199970, A287142, A328165.

Programs

  • Mathematica
    dat=Table[Count[Differences[Divisors[n]],1],{n,10000}];
Sort[Table[Position[dat,i][[1,1]],{i,Union[dat]}]]

A328458 Maximum run-length of the nontrivial divisors (greater than 1 and less than n) of n.

Original entry on oeis.org

1, 0, 0, 1, 0, 2, 0, 1, 1, 1, 0, 3, 0, 1, 1, 1, 0, 2, 0, 2, 1, 1, 0, 3, 1, 1, 1, 1, 0, 2, 0, 1, 1, 1, 1, 3, 0, 1, 1, 2, 0, 2, 0, 1, 1, 1, 0, 3, 1, 1, 1, 1, 0, 2, 1, 2, 1, 1, 0, 5, 0, 1, 1, 1, 1, 2, 0, 1, 1, 1, 0, 3, 0, 1, 1, 1, 1, 2, 0, 2, 1, 1, 0, 3, 1, 1, 1, 1, 0, 2, 1, 1, 1, 1, 1, 3, 0, 1, 1, 2, 0, 2, 0, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Oct 17 2019

Keywords

Comments

By convention, a(1) = 1, and a(p) = 0 for p prime.

Examples

			The non-singleton runs of the nontrivial divisors of 1260 are: {2,3,4,5,6,7} {9,10} {14,15} {20,21} {35,36}, so a(1260) = 6.

Links

Crossrefs

Positions of first appearances are A328459.
Positions of 0's and 1's are A088723.
The version that looks at all divisors is A055874.
The number of successive pairs of divisors > 1 of n is A088722(n).
The Heinz number of the multiset of run-lengths of divisors of n is A328166(n).
Cf. A033676, A060681, A060775, A070824, A088725, A129308, A181063, A199970, A328165, A163870, A328194, A328448, A328449.

Programs

  • Mathematica
    Table[Switch[n,1,1,?PrimeQ,0,,Max@@Length/@Split[DeleteCases[Divisors[n],1|n],#2==#1+1&]],{n,100}]
  • PARI
    A328458(n) = if(1==n,n,my(rl=0,pd=0,m=0); fordiv(n, d, if(1(1+pd), m = max(m,rl); rl=0); pd=d; rl++)); max(m,rl)); \\ Antti Karttunen, Feb 23 2023

Extensions

Data section extended up to a(105) by Antti Karttunen, Feb 23 2023
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.