OFFSET
1,2
COMMENTS
This is the 7th row of the table A138180.
REFERENCES
See A002071.
MATHEMATICA
pMax = 17; smoothMax = 10^12; smoothNumbers[p_?PrimeQ, max_] := Module[{a, aa, k, pp, iter}, k = PrimePi[p]; aa = Array[a, k]; pp = Prime[Range[k]]; iter = Table[{a[j], 0, PowerExpand@Log[pp[[j]], max/Times @@ (Take[pp, j - 1]^Take[aa, j - 1])]}, {j, 1, k}]; Table[Times @@ (pp^aa), Sequence @@ iter // Evaluate] // Flatten // Sort]; Select[(Sqrt[1 + 4*smoothNumbers[pMax, smoothMax]] - 1)/2, IntegerQ]
PROG
(PARI) is(n)=my(t=510510); n*=n+1; while((t=gcd(n, t))>1, n/=t); n==1 \\ Charles R Greathouse IV, Nov 13 2016
CROSSREFS
KEYWORD
nonn,fini,full
AUTHOR
Jean-François Alcover, Nov 13 2016
STATUS
approved