OFFSET
1,1
COMMENTS
EXAMPLE
a(2)=5 because set of prime divisors of 5^2-1 =2^3*3 is {2,3}, the same as for example for 7^2-1 = 2^4*3.
MATHEMATICA
aa = {}; bb = {}; cc = {}; ff = {}; Do[k = n^2 - 1; kk = FactorInteger[k]; b = {}; Do[AppendTo[b, kk[[m]][[1]]], {m, 1, Length[kk]}]; dd = Position[aa, b]; If[dd == {}, AppendTo[cc, n]; AppendTo[aa, b], AppendTo[ff, n]; AppendTo[bb, cc[[dd[[1]][[1]]]]]], {n, 2, 1000000}]; Take[Union[bb, ff], 100] (* Artur Jasinski *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Artur Jasinski, Oct 12 2010, Oct 21 2010
EXTENSIONS
Name improved by T. D. Noe, Nov 15 2010
STATUS
approved