OFFSET
1,1
LINKS
Zak Seidov, Table of n, a(n) for n = 1..3000
MAPLE
A000217inv:=proc(n)
local t1;
t1:=floor(sqrt(2*n));
if n = t1*(t1+1)/2 then
return t1 ;
else
return -1;
end if;
end proc:
for n from 1 to 90000 do
s := A034963(n)/4 ;
if type(s, 'integer') then
tr := A000217inv(s) ;
if tr >= 0 then
printf("%d, ", s) ;
end if;
end if;
end do: # R. J. Mathar, Jun 06 2013
MATHEMATICA
Select[Mean/@Partition[Prime[Range[20000]], 4, 1], OddQ[Sqrt[8#+1]]&] (* Harvey P. Dale, Dec 18 2015 *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alex Ratushnyak, May 30 2013
STATUS
approved