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

A175304 A positive integer n is included if d(n+d(n)) = d(n), where d(n) is the number of divisors of n.

Original entry on oeis.org

3, 5, 6, 10, 11, 12, 17, 22, 29, 34, 35, 41, 44, 51, 58, 59, 60, 65, 70, 71, 72, 82, 84, 87, 91, 92, 96, 101, 102, 107, 111, 115, 118, 119, 125, 128, 129, 130, 137, 141, 142, 147, 149, 155, 174, 179, 182, 183, 191, 197, 201, 202, 205, 209, 213, 214, 215, 217, 222
Offset: 1

Views

Author

Leroy Quet, Mar 24 2010

Keywords

Comments

The sequence contains the smaller member of every pair of twin primes (A001359) and all squarefree semiprimes m such that m+4 is also a squarefree semiprime (A255746). Can one prove that this is an infinite sequence? - Vladimir Shevelev, Jul 11 2015
The sequence does not contain perfect squares. Indeed, let a(m)=k^2. Then d(k^2+d(k^2)) = d(k^2). Note that d(k^2) is odd. On the other hand, it is known (A046522) that d(k^2)<2*k. Hence, (k+1)^2 - k^2 > d(k^2). Thus k^2Vladimir Shevelev, Feb 10 2017
If p is prime and t+1 is odd prime, then p^t is not in the sequence. Indeed, if d(p^t+t+1)=t+1, then p^t+t+1=q^t, where q is prime > p (if p^t+t+1= say q^l*r^m, then (l+1)*(m+1)=t+1 which is impossible by the condition). But q>=p+2 and p^t+t+1>=p^t+2*t*p^(t-1) or t+1>=2*t*p^(t-1) which trivially has only solution t=1; however, by the condition t>=2. - Vladimir Shevelev, Feb 18 2017
If an odd integer k is in this sequence, so is 2k. - Charlie Neder, Jan 14 2019

Examples

			10 is in the sequence because d(10)=4 and d(10+d(10))=d(14)=4. - _Emeric Deutsch_, Apr 08 2010
		

Crossrefs

Positions of zeros in A286530.

Programs

  • Maple
    with(numtheory): a := proc (n) if tau(n+tau(n)) = tau(n) then n else end if end proc: seq(a(n), n = 1 .. 230); # Emeric Deutsch, Apr 08 2010
  • Mathematica
    Select[Range@ 224, Function[n, DivisorSigma[0, n + #] == # &@ DivisorSigma[0, n]]](* Michael De Vlieger, Sep 27 2015 *)
    Position[#, 0][[All, 1]] &@ Table[DivisorSigma[0, n + DivisorSigma[0, n]] - DivisorSigma[0, n], {n, 222}] (* Michael De Vlieger, May 21 2017 *)
  • PARI
    is(n)=numdiv(n+n=numdiv(n))==n \\ M. F. Hasler, Sep 27 2015

Extensions

More terms from Emeric Deutsch, Apr 08 2010
Showing 1-1 of 1 results.