OFFSET
1,1
COMMENTS
All m for which 2*m+1 is in A003627 are in the sequence:
This concerns m=2, 5, 8, 11, 14, 20, 23, 26, 29, 35,...
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
isA177966 := proc(n) A177961(m) = m+2 ; end proc:
for m from 1 to 800 do if isA177966(m) then printf("%d, ", m) ; end if; end do:
# R. J. Mathar, Oct 25 2010
N:= 1000: # to get all terms <= N
A1:= map(t -> (t-1)/2, select(isprime, {seq(6*i-1, i=1..(N+1)/3)})):
A2:= map(t -> (t+1)/2, select(isprime, {seq(23+30*i, i=0..(N-12)/15)})):
sort(convert(A1 union A2, list));
# Robert Israel, Jul 31 2015
MATHEMATICA
M = 1000; (* to get all terms <= M *)
A1 = (Select[Table[6 i - 1, {i, 1, (M + 1)/3}], PrimeQ] - 1)/2;
A2 = (Select[Table[23 + 30 i, {i, 0, (M - 12)/15}], PrimeQ] + 1)/2;
Union[A1, A2] (* Jean-François Alcover, Jul 17 2020, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Vladimir Shevelev, May 16 2010
EXTENSIONS
Corrected (11, 23, 27, etc. inserted) and extended by R. J. Mathar, Oct 25 2010
STATUS
approved