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!)
A200778 Least k >0 such that k*p*(k*p-1)-1 and k*p*(k*p-1)+1 is a twin prime pair, where p=prime(n). 2
2, 1, 5, 1, 2, 3, 3, 13, 9, 8, 10, 43, 69, 15, 17, 50, 3, 42, 1, 2, 3, 3, 20, 33, 3, 44, 7, 35, 49, 9, 6, 189, 15, 1, 113, 21, 7, 154, 3, 3, 18, 12, 29, 33, 20, 6, 27, 3, 2, 3, 23, 11, 10, 12, 18, 137, 41, 12, 36, 29, 54, 17, 10, 59, 55, 3, 51, 36 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
lim_{N->infinity) (Sum_{n=1..N} k(n)) / (Sum_{n=1..N} log(p(n))^2) = 1.
LINKS
EXAMPLE
2*2*(2*2 - 1) - 1 = 11, twin prime of 13, so a(1)=2.
MAPLE
A200778 := proc(n)
p := ithprime(n) ;
for k from 1 do
if isprime(k*p*(k*p-1)-1) and isprime(k*p*(k*p-1)+1) then
return k;
end if;
end do:
end proc:
seq(A200778(n), n=1..80) ; # R. J. Mathar, Nov 26 2011
MATHEMATICA
lktpp[n_]:=Module[{k=1, p=Prime[n]}, While[AnyTrue[k*p(k*p-1)+{1, -1}, CompositeQ], k++]; k]; Array[lktpp, 70] (* Harvey P. Dale, May 03 2019 *)
CROSSREFS
Cf. A200654.
Sequence in context: A222481 A351954 A353577 * A345355 A132601 A047818
KEYWORD
nonn
AUTHOR
Pierre CAMI, Nov 22 2011
STATUS
approved

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 April 20 00:26 EDT 2024. Contains 371798 sequences. (Running on oeis4.)