login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A145714
a(n) = ceiling(sqrt(2*A145701(n))).
2
3, 5, 6, 11, 12, 19, 25, 29, 42, 48, 65, 72, 79, 84, 90, 101, 174, 204, 205, 209, 221, 245, 264, 275, 289, 299, 306, 349, 371, 372, 408, 409, 415, 449, 456, 474, 521, 535, 546, 569, 571, 575, 594
OFFSET
1,1
COMMENTS
If p_n is the n-th prime, then ceiling(sqrt(2p_n)) is in the sequence iff A145236(n) = A145236(n+1).
MAPLE
A145701 :=proc(n) local p; if n = 1 then return 3; end if; p := nextprime(procname(n-1)) ; while true do if isprime(p) and isprime(p+2) then if floor(sqrt(2*p)) <> floor( sqrt(2*p+4)) then return p; end if; end if; p := nextprime(p) ; end do; end proc:
A145714 := proc(n) ceil(sqrt(2*A145701(n))) ; end proc:
for n from 1 to 80 do printf("%d, \n", A145714(n)) ; end do; # R. J. Mathar, Aug 02 2010
CROSSREFS
Sequence in context: A001838 A285785 A080759 * A326310 A047443 A173593
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, Oct 16 2008
EXTENSIONS
More terms from R. J. Mathar, Aug 02 2010
STATUS
approved