login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A067889 Primes sandwiched between two numbers having same number of divisors. 14

%I #30 Jun 28 2022 02:57:11

%S 7,19,41,103,137,199,307,349,491,739,823,919,1013,1061,1193,1277,1289,

%T 1409,1433,1447,1481,1543,1609,1667,1721,1747,2153,2357,2441,2617,

%U 2683,2777,3259,3319,3463,3581,3593,3769,3797,3911,3943,4013,4217,4423,4457

%N Primes sandwiched between two numbers having same number of divisors.

%C Primes p such that tau(p+1) = tau(p-1) where tau(k) = A000005(k).

%C These are the primes in sequence A067888 of numbers n such that tau(n+1) = tau(n-1). - _M. F. Hasler_, Aug 06 2015

%H T. D. Noe, <a href="/A067889/b067889.txt">Table of n, a(n) for n=1..10000</a>

%F a(n) seems curiously to be asymptotic to 25*n*log(n). [From the number of terms up to 10^8, 10^9, 10^10 and 10^11, i.e., 306147, 2616930, 22835324 and 202105198, this constant can be estimated by 25.858..., 25.858..., 25.845... and 25.872..., respectively. - _Amiram Eldar_, Jun 28 2022]

%e 7 is a member as 6 and 8 both have 4 divisors; 19 is a member as 18 and 20 both have 6 divisors each.

%p with(numtheory):j := 0:for i from 1 to 10000 do b := ithprime(i): if nops(divisors(b-1))=nops(divisors(b+1)) then j := j+1:a[j] := b:fi:od:seq(a[k],k=1..j);

%t Prime[ Select[ Range[ 700 ], Length[ Divisors[ Prime[ #1 ] - 1 ]] == Length[ Divisors[ Prime[ #1 ] + 1 ]] & ]]

%t Select[Prime[Range[1000]],DivisorSigma[0,#-1]==DivisorSigma[0,#+1]&] (* _Harvey P. Dale_, Jun 08 2018 *)

%o (PARI) is_A067889(p)=numdiv(p-1)==numdiv(p+1)&&isprime(p) \\ _M. F. Hasler_, Jul 31 2015

%Y Cf. A067888, A067891 (analog with sigma).

%K easy,nonn

%O 1,1

%A _Benoit Cloitre_, Mar 02 2002

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 29 06:15 EDT 2024. Contains 371265 sequences. (Running on oeis4.)