A175953 Let a(1)=1; for n>1 a(n)=nextprime(a(n-1)+(a(n-1)+1)/4).

1, 2, 3, 5, 7, 11, 17, 23, 29, 37, 47, 59, 79, 101, 127, 163, 211, 269, 337, 431, 541, 677, 853, 1069, 1361, 1709, 2137, 2677, 3347, 4201, 5261, 6577, 8231, 10289, 12889, 16127, 20161, 25219, 31531, 39419, 49277, 61603, 77017, 96281, 120371, 150473

Offset

1,2

Comments

The following definition of nextprime(q) is used: if q is integer and prime, nextprime(q)=q. If q is integer and composite or rational, nextprime(q) is the smallest prime >q. [From R. J. Mathar, Oct 30 2010]

Links

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

Maple

Contribution from R. J. Mathar, Oct 30 2010: (Start)
nprime := proc(n) if type(n, 'integer') then if isprime(n) then return n; else return nextprime(n) ; end if; else return nextprime(floor(n)) ; end if; end proc:
A175953 := proc(n) option remember; if n= 1 then 1; else p := procname(n-1)+(procname(n-1)+1)/4 ; return nprime(p) ; end if; end proc:
seq(A175953(n), n=1..120) ; (End)

Crossrefs

Cf. A008578, A055496, A089571, A162336.

Sequence in context: A075551 A070866 A237285 * A362017 A040089 A267944

Adjacent sequences: A175950 A175951 A175952 * A175954 A175955 A175956

Keyword

nonn,easy

Author

Juri-Stepan Gerasimov, Oct 29 2010

Extensions

More terms from R. J. Mathar, Oct 30 2010

Status

approved