OFFSET
1,2
COMMENTS
LINKS
Arkadiusz Wesolowski, Table of n, a(n) for n = 1..102
PROG
(Magma) lst:=[]; pr:=1; for n in [0..203] do pd:=PrimeDivisors(Lucas(n)); d:=1; t:=0; for c in [1..#pd] do f:=pd[c]; if Gcd(pr, f) eq 1 then t+:=1; else d:=d*f; end if; end for; if t eq 1 then Append(~lst, n); end if; pr:=pr*Truncate(Lucas(n)/d); end for; lst;
CROSSREFS
KEYWORD
nonn
AUTHOR
Arkadiusz Wesolowski, Jun 22 2016
STATUS
approved