OFFSET
1,1
COMMENTS
The sequence contains A056220(n) for n>1. For a given x, solutions (x,y) to x^2 - D*y^2 = 1 do not all stand for the least nontrivial positive solution pair. For instance, x=244 solves the latter equation for D=15 and D=135 with corresponding y=63 and y=21. While (244,21) is indeed the fundamental solution to x^2 - 135*y^2 = 1, (244,63) comes only as the third one after the solution pairs (4,1) and (31,8) to x^2 - 15*y^2 = 1.
PROG
(PARI) { for(x=2, 10^3, f=core(x^2-1, 1)[2]; if(f>1 && !isprime(f), print1(x, ", ")) ) } \\ Max Alekseyev, Jul 25 2009
CROSSREFS
KEYWORD
nonn
AUTHOR
Lekraj Beedassy, Feb 25 2002
EXTENSIONS
Corrected and extended by Max Alekseyev, Jul 25 2009
STATUS
approved