A107819 Slowest increasing sequence where a(n)+n is prime.

1, 3, 4, 7, 8, 11, 12, 15, 20, 21, 26, 29, 30, 33, 38, 43, 44, 49, 52, 53, 58, 61, 66, 73, 76, 77, 80, 81, 84, 97, 100, 105, 106, 115, 116, 121, 126, 129, 134, 139, 140, 149, 150, 153, 154, 165, 176, 179, 180, 183

Offset

1,2

Comments

Appears to be closely related to A014692.

Links

Table of n, a(n) for n=1..50.

Example

(1)+1=2, which is prime; (3)+2=5=prime; (4)+3=7=prime; (7)+4=11=prime; (8)+5=13=prime; etc.

Maple

A107819 := proc(n)
        option remember;
        if n = 1 then
                return 1;
        else
                for a from procname(n-1)+1 do
                        if isprime(a+n) then
                                return a;
                        end if;
                end do:
        end if;
end proc: # R. J. Mathar, Sep 05 2011

Prog

(Sage)
@cached_function
def A107819(n):
    prev = (-1 if n == 1 else A107819(n-1))
    return next_prime(prev+n)-n  # D. S. McNeil, Sep 05 2011

Crossrefs

Sequence in context: A309388 A188259 A058235 * A300950 A368083 A067186

Adjacent sequences: A107816 A107817 A107818 * A107820 A107821 A107822

Keyword

easy,nonn

Author

Eric Angelini, Jun 11 2005

Extensions

Corrected by R. J. Mathar, Sep 05 2011

Status

approved