OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
MAPLE
tri:= n -> n*(n+1)/2:
istri:= x -> issqr(1+8*x):
select(t -> istri(tri(t+2) &^ tri(t+1) mod tri(t)), [$1..10000]);
PROG
(PARI) tri(n) = n*(n+1)/2; \\ A000217
isok(n) = ispolygonal(lift(Mod(tri(n+2), tri(n))^tri(n+1)), 3); \\ Michel Marcus, Jan 25 2021
CROSSREFS
KEYWORD
nonn
AUTHOR
J. M. Bergot and Robert Israel, Jan 24 2021
STATUS
approved