OFFSET
1,1
COMMENTS
MATHEMATICA
Select[Range[9000], FactorInteger[#^2-1][[-1, 1]]==7&]
PROG
(Magma) [ n: n in [2..9000] | m eq 7 where m is D[#D] where D is PrimeDivisors(n^2-1) ]; // Klaus Brockhaus, Feb 17 2011
(PARI) is(n)=n=n^2-1; forprime(p=2, 5, n/=p^valuation(n, p)); n>1 && 7^valuation(n, 7)==n \\ Charles R Greathouse IV, Jul 01 2013
CROSSREFS
KEYWORD
fini,full,nonn
AUTHOR
Artur Jasinski, Oct 21 2010
STATUS
approved