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.

A195150 Number of divisors d of n such that d-1 does not divide n.

Original entry on oeis.org

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

Views

Author

Omar E. Pol, Sep 19 2011

Keywords

Comments

Define "subdivisor" of n to be the positive integer b such that b = d - 1, if d divides n and b does not divide n. For the list of subdivisors of n see A195153.
First occurrence of k is given in A173569. - Robert G. Wilson v, Sep 23 2011

Examples

			a(24) = 4 since the divisors of 24 are 1,2,3,4,6,8,12,24, so the subdivisors of 24 are 5,7,11,23 because 6-1 = 5, 8-1 = 7, 12-1 = 11 and 24-1 = 23. Note that the positive integers 1,2,3 are not subdivisors of 24 because they are divisors of 24.
		

Crossrefs

Programs

  • Haskell
    a195150 n = length [d | d <- [3..n], mod n d == 0, mod n (d-1) /= 0]
    -- Reinhard Zumkeller, Sep 23 2011
  • Mathematica
    f[n_] := Module[{d = Divisors[n]}, Length[Select[Rest[d-1], Mod[n, #] > 0 &]]]; Table[f[n], {n, 100}] (* T. D. Noe, Sep 22 2011 *)

Formula

a(n) = A137921(n) - 1. - Robert G. Wilson v, Sep 23 2001
Sum_{k=1..n} a(k) ~ n * (log(n) + 2*gamma - 3), where gamma is Euler's constant (A001620). - Amiram Eldar, Jan 18 2024