OFFSET
1,1
COMMENTS
See A187950.
LINKS
Chai Wah Wu, Table of n, a(n) for n = 1..10000
MAPLE
r := sqrt(5) ;
for n from 1 to 1200 do
a := floor((n+4)*r)-floor(n*r)-floor(4*r) ;
if a = 0 then
printf("%d, ", n);
end if;
end do: # R. J. Mathar, Jun 03 2011
MATHEMATICA
Flatten[Position[Table[Floor[(n+4)Sqrt[5]]-Floor[n*Sqrt[5]]-Floor[ 4*Sqrt[5]], {n, 1100}], 0]] (* Harvey P. Dale, Sep 23 2014 *)
PROG
(Python)
from gmpy2 import isqrt
A188290_list = [n for n in range(1, 10**6) if isqrt(5*(n+4)**2) - isqrt(5*n**2) == 8] # Chai Wah Wu, Oct 08 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Clark Kimberling, Mar 26 2011
EXTENSIONS
All values corrected by R. J. Mathar, Jun 03 2011
STATUS
approved