|
|
A177966
|
|
Indices m for which A177961(m) = 2 + m.
|
|
3
|
|
|
2, 5, 8, 11, 12, 14, 20, 23, 26, 27, 29, 35, 41, 42, 44, 50, 53, 56, 57, 65, 68, 74, 83, 86, 87, 89, 95, 98, 113, 116, 117, 119, 125, 128, 131, 132, 134, 140, 146, 147, 155, 158, 173, 176, 177, 179, 191, 192, 194, 200, 209, 215, 221, 222, 224, 230, 233, 239, 245, 251, 252, 254
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
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
|
|
|
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:
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));
|
|
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;
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
Corrected (11, 23, 27, etc. inserted) and extended by R. J. Mathar, Oct 25 2010
|
|
STATUS
|
approved
|
|
|
|