OFFSET
1,1
LINKS
FORMULA
a(n) = 4*A263941(n) + 1 for n>=2, a(n) = 4*floor((9*n)/2) + 1 for all n. - Hugo Pfoertner, Dec 08 2020
From Chai Wah Wu, Aug 18 2025: (Start)
a(n) = a(n-1) + a(n-2) - a(n-3) for n > 3.
G.f.: x*(-x^2 + 20*x + 17)/((x - 1)^2*(x + 1)). (End)
MATHEMATICA
Select[Range[2, 10^3], Function[b, AllTrue[{2, 3}, PowerMod[b, (# - 1), #^2] == 1 &]]] (* Michael De Vlieger, Dec 10 2020 *)
PROG
(PARI) is(n) = forprime(p=1, 3, if(Mod(n, p^2)^(p-1)!=1, return(0))); 1
CROSSREFS
KEYWORD
nonn
AUTHOR
Felix Fröhlich, Dec 08 2020
STATUS
approved
