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!)
A211678 Twin primes p, p+2 with unique values of sigma(p) and sigma(p+2); sigma(n) = A000203(n) = sum of divisors of n. 2

%I #7 Apr 26 2012 13:54:16

%S 3,5,7,197,199,281,283,347,349,461,463,641,643,821,823,857,859,1289,

%T 1291,1697,1699,1721,1723,1787,1789,1877,1879,2081,2083,2141,2143,

%U 2381,2383,2549,2551,2801,2803,3257,3259,3539,3541,3557,3559,3929,3931,4019,4021

%N Twin primes p, p+2 with unique values of sigma(p) and sigma(p+2); sigma(n) = A000203(n) = sum of divisors of n.

%C Subsequence of A211656 and A211660.

%e Twin primes 197 and 199 are in sequence because sigma(197) = 198, sigma(199) = 200 and there are no other numbers m, n with sigma(m) = 198 or sigma(n) = 200.

%t d = DivisorSigma[1, Range[4100]]; t = Transpose[Select[Tally[Sort[d]], #[[2]] == 1 && #[[1]] <= Length[d] &]][[1]]; t2 = Sort[Flatten[Table[Position[d, i], {i, t}]]]; t3 = Select[t2, PrimeQ]; tp = {}; Do[If[t3[[i + 1]] - t3[[i]] == 2 && DivisorSigma[1, t3[[i]]] != DivisorSigma[1, t3[[i + 1]]], AppendTo[tp, t3[[i]]]; AppendTo[tp, t3[[i]] + 2]], {i, Length[t3] - 1}]; Union[tp] (* _T. D. Noe_, Apr 26 2012 *)

%Y Cf. A211767 (lesser of twin primes p, p+2 with unique values of sigma(p) and sigma(p+2), A211769 (greater of twin primes p, p+2 with unique values of sigma(p) and sigma(p+2).

%K nonn

%O 1,1

%A _Jaroslav Krizek_, Apr 20 2012

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 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)