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

%I #49 Jan 18 2024 01:41:19

%S 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,

%T 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,

%U 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

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

%C 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

%H Charles R Greathouse IV, <a href="/A137921/b137921.txt">Table of n, a(n) for n = 1..10000</a>

%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/DivisorFunction.html">Divisor Function</a>.

%F a(n) <= A000005(n), with equality iff n is odd; a(A137922(n)) = 2.

%F a(n) = A000005(n) - A129308(n). - _Michel Marcus_, Jan 06 2015

%F a(n) = A001222(A328166(n)). - _Gus Wiseman_, Oct 16 2019

%F Sum_{k=1..n} a(k) ~ n * (log(n) + 2*gamma - 2), where gamma is Euler's constant (A001620). - _Amiram Eldar_, Jan 18 2024

%e 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.

%p 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

%t 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 *)

%t Table[Count[Differences[Divisors[n]],_?(#>1&)]+1,{n,110}] (* _Harvey P. Dale_, Jun 05 2012 *)

%t a[n_] := DivisorSum[n, Boole[!Divisible[n, #+1]]&]; Array[a, 100] (* _Jean-François Alcover_, Dec 01 2015 *)

%o (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)

%o (PARI) a(n)=sumdiv(n,d,(n%(d+1)!=0)); \\ _Joerg Arndt_, Jan 06 2015

%o (Haskell)

%o a137921 n = length $ filter (> 0) $

%o map ((mod n) . (+ 1)) [d | d <- [1..n], mod n d == 0]

%o -- _Reinhard Zumkeller_, Nov 23 2011

%o (Python)

%o from sympy import divisors

%o def A137921(n):

%o ....return len([d for d in divisors(n,generator=True) if n % (d+1)])

%o # _Chai Wah Wu_, Jan 05 2015

%Y Bisections: A099774, A174199.

%Y First appearance of n is at position A173569(n).

%Y Numbers whose divisors have no non-singleton runs are A005408.

%Y The longest run of divisors of n has length A055874(n).

%Y The number of successive pairs of divisors of n is A129308(n).

%Y Cf. A000005, A001620, A027750, A060680, A088723, A088725, A181063, A199970, A328165, A328166, A328448, A328450.

%K nonn,nice

%O 1,3

%A _Reinhard Zumkeller_, Feb 23 2008

%E Corrected and edited by _Charles R Greathouse IV_, Apr 19 2010

%E Edited by _N. J. A. Sloane_, Aug 10 2010