OFFSET
1,1
COMMENTS
Conjecture (verified for k <= 10^6 by M. F. Hasler): A024812(k) = (A024813(k)+k-1)/2, k=1,2,.... - L. Edson Jeffery, Mar 21 2013
FORMULA
a(k) = floor(cot(Pi/(2*A024813(k)))). - L. Edson Jeffery, Mar 21 2013
Conjecture: a(n) = a(n-1) + a(n-3) - a(n-4); g.f.: x*(x^15-x^14+3*x^2+2*x+2) / ((x-1)^2*(x^2+x+1)). - Colin Barker, Jan 03 2014
MATHEMATICA
f[n_] := Floor[Cot[Pi/(2 n)]]; f[ Select[ Range[2, 245], f[# - 1] < f[#] < f[# + 1] &]] (* Robert G. Wilson v, Mar 27 2013 *)
PROG
(PARI) {my(f(m)=floor(cotan(Pi/(2*m)))); for(m=2, 999, f(m-1)<f(m) & f(m)<f(m+1) & print1(f(m)", "))} \\ Note: Depending on default(realprecision), e.g. when this is set to 99, floor(cotan(Pi/4)) may yield 0 and erroneous output of f(3)=1. [M. F. Hasler, Mar 20 2013]
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Definition corrected, following posts to the SeqFan list from Harvey P. Dale and Don Reble, by M. F. Hasler, Mar 20 2013
STATUS
approved