Al principio explicaba por que d dividía a (2n+1)
Cuando suponemos que d es 401 entonces 401 divide a 2n+1 luego tendremos
2n+1=401m con m entero
Y es entero positivo porque también lo es n por hipótesis y la cantidad de la izquierda resulta ser positiva entonces.
Es verdad, el paréntesis es par si m es impar.
Tomemos ese número m impar
m = 2r+1 con r € NU{0}
2n+1=401(2r+1)
n = (802r + 401 - 1)/2 = 401r + 200
An = 100 + (401r+200)^2 = 100+(401^2)r^2+401·400r + 40000 =
Calculamos aparte 40000 + 100 = 40100 = 401·100
An=401(401r^2 + 400r + 100)
Recuerdo que estamos calculando el mcd(An, 2n+1) que es el mismo que mcd(An, An+1)
Y ese mcd es 401 porque es factor común de ambos y no puede ser del tipo 401 por algo porque antes habiamos demostrado que d dividía a 401.
Lo he comprobado con el ordenador mediante este programa en Pascal:
program mcd100mascuadradosseguidos;
function mcd(a,b:integer):integer;
begin
if a<b then
if (b mod a) = 0 then
mcd := a
else
mcd := mcd(b mod a, a)
else
if (a mod b) = 0 then
mcd := b
else
mcd := mcd(a mod b, b)
end;
var i,m:integer;
begin
for i:=1 to 46000 do
begin
m := mcd(100+i*i, 100+(i+1)*(i+1));
if (m <> 1) then
begin
writeln(i,' ',m);
if (m <> 401) then
begin
writeln('CASO ESPECIAL'); readln;
end;
end;
end;
readln;
end.
Y el resultado para esos 46000 primeros mcd ( no me deja probar con más más por que los cuadrados desbordaban el rango positivo de las variable integer de 16bits) es el previsto en la teoria que dice:
El mcd es 401 cuando n = 401r + 200 con r € NU{0} o sea: 200, 601, 1002, 1403, 1804, etc. y es 1 en los demás casos
Y eso es todo, si tienes algunba otra duda o hay algo mal me lo comunicas.