La prueba está bien,
pero la magia de los valores
de a b
me molestan.
Así,
Voy a probar
para encontrar todas las soluciones.
Spoiler:
Yo no,
pero no le dan algunas condiciones.
La suma de los divisores
de paqb es
σ(paqb)=(p+1−1)(pb+1−1)(p−1)(q−1)
y el número de divisores es
d(paqb)=(a+1)(b+1).
Por lo que el promedio es de
r(paqb)=σ(paqb)d(paqb)=(p+1−1)(pb+1−1)(p−1)(q−1)(a+1)(b+1).
En su caso,
este es
a=(q−1)/2,b=1,
así
a+1=(q+1)/2 y
r(paqb)=(pa+1−1)(qb+1−1)(p−1)(q−1)(a+1)(b+1)=(p(q+1)/2−1)(q2−1)(p−1)(q−1)((q+1)/2))2=(p(q+1)/2−1)(q2−1)(p−1)(q−1)(q+1)=(p(q+1)/2−1)(p−1)
y este es un número entero.
Si b=1,
r(paqb)=(pa+1−1)(qb+1−1)(p−1)(q−1)(a+1)(b+1)=(pa+1−1)(q2−1)(p−1)(q−1)(a+1)2=(pa+1−1)(q+1)2(a+1)(p−1)
Para que esto sea un número entero
es suficiente con que
2(a+1)|(q+1),
así
2(a+1)=q+1 obras
si q es impar.
En general,
si q+1=2m(a+1),
a=p+12m−1.
Si q=2uv−1
donde u≥1 v es impar,
entonces
a=2uv2m−1=2u−1vm−1.
En particular,
la elección de m=v da
a=2u−1−1
y
m=1 da
OP de la solución.
Otra posibilidad es tener
a+1=2m , de modo que
(p+1−1)(q+1)2(a+1)(p−1)=(p2m−1)(q+1)2(a+1)(p−1)=(pm−1)(pm+1)(q+1)4m(p−1)
así que
4m|(pm+1)
iba a funcionar.
Para m=1, esto requiere que
p≡3mod;
para m=2, esto requiere que
p^2 \equiv 7 \bmod 8
lo que no puede suceder ya que
p^2 \equiv 1 \bmod 8.
Si a+1=2m b+1 = 2n
\begin{array}\\
r(p^aq^b)
&=\dfrac{(p^{a+1}-1)(q^{b+1}-1)}{(p-1)(q-1)(a+1)(b+1)}\\
&=\dfrac{(p^{2m}-1)(q^{2n}-1)}{(p-1)(q-1)4mn}\\
&=\dfrac{(p^{m}-1)(q^{n}-1)(p^{m}+1)(q^{n}+1)}{(p-1)(q-1)4mn}\\
\end{array}
así es suficiente si
4mn | (p^{m}+1)(q^{n}+1)
.
En este punto,
Yo no veo ninguna
soluciones fáciles,
lo voy a dejar aquí.