Esto me pareció un buen punto de partida para examinar la biblioteca GMP C con un poco más de detalle, así que decidí escribir un programa rápido para imprimir números que se ajusten a la descripción de este problema. No estoy seguro de que esto sea una "respuesta", pero pensé que el programa podría interesar a las partes interesadas.
Código abajo, enlace a un gist de github aquí .
// This program is a quick brute force test of the question described at
// http://math.stackexchange.com/questions/612346/find-all-positive-integers-n-s-t-3n-5n-is-divisible-by-n2-1
// Copyright (C) 2014 Victor Robertson
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
// You should have received a copy of the GNU General Public License
// along with this program. If not, see <http://www.gnu.org/licenses/>.
#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>
size_t gLower = 1;
size_t gUpper = (0 - 1); // force underflow = size_t max
size_t i;
void sig_handler(int signo)
{
printf("tested %zu to %zu\n", gLower, i - 1);
exit(0);
}
int main(int argc, char* argv[]) {
// all n s.t. 3^n+5^n is divisible by n^2−1
if (signal(SIGINT, sig_handler) == SIG_ERR) {
return 1;
}
// set args if provided
if (argc > 1) {
gLower = atol(argv[1]);
if (argc > 2) {
gUpper = atol(argv[2]);
}
}
i = gLower;
mpz_t i_sqr, three_n, five_n, res;
mpz_init(i_sqr);
mpz_init(three_n);
mpz_init(five_n);
mpz_init(res);
// start by calculating 3^n and 5^n where n is the first number to test
mpz_ui_pow_ui(three_n, 3, i);
mpz_ui_pow_ui(five_n, 5, i);
// 3 is a known solution, might as well skip straight to 5
if (gLower < 5) {
printf("3\n");
gLower = 5;
}
// adjust the lower bound to exclude odds
if (gLower % 2 == 0) {
gLower -= 1;
}
for(i = gLower; i < gUpper; i += 2) {
// I think this can be accomplished in a better, possibly more efficient means though I'm not positive.
mpz_ui_pow_ui(i_sqr, i, 2);
mpz_sub_ui(i_sqr, i_sqr, 1);
mpz_add(res, three_n, five_n);
if (mpz_divisible_p(res, i_sqr)) {
printf("%zu\n",i);
}
// multiple the current 3^n and 5^n by 3 and 5 respectively to acquire 3^(n+1) and 5^(n+1)
mpz_mul_ui(three_n, three_n, 3);
mpz_mul_ui(five_n, five_n, 5);
}
printf("Done!\n");
mpz_clear(i_sqr);
mpz_clear(three_n);
mpz_clear(five_n);
mpz_clear(res);
}
2 votos
Las comprobaciones numéricas lo confirman hasta $240002$
0 votos
@TimRatigan $240002\not\equiv 0\pmod 3$ . ¿Está utilizando el hecho de que $n\equiv 0\pmod 3$ para acelerar sus cálculos?
1 votos
Sí. Incluso para $n$ , $6\mid n$ y para impar $n$ , $n$ es de la forma $3(8k+1)$ . Lo confirmé hasta $240\,000$ y lo extendió trivialmente a $240\,002$ . @IanMateus
4 votos
Me voy a dormir, pero me voy con la idea de que tal vez la relación de recurrencia $$ 3^{n+2}+5^{n+2} = 8(3^{n+1}-5^{n+1})-15(3^n+5^n) $$ puede ayudar.
0 votos
@Tim: incluso $n$ no puede ser una solución.
0 votos
@Marek ¿por qué? ${}{}{}{}$
0 votos
@Marek ¿por qué? Tim, Marek Este problema es muy bonito, deberíamos abrir un chat sobre ello mañana.
0 votos
@Tim: $3^n + 5^n$ es siempre par mientras que $n^2 - 1$ sería impar.
0 votos
@Marek la última vez que lo comprobé los números Impares pueden dividir a los pares.
0 votos
@IanMateus hagamos eso.
0 votos
@Tim: vale, debería irme oficialmente a dormir :D Por alguna razón pensé que el problema pedía igualdad ahora. Duh...
0 votos
@TimRatigan Tienes razón, lo pasé por alto por completo.
3 votos
Lo he probado numéricamente para $n=1,\ldots,10^6$ y sólo $n=3$ funciona.
0 votos
De hecho, $n-1$ divide $3^n+5^n$ , para $n=2, 3, 5, 9, 18, 39, 153, 222, 378, 630, 1685, 1749, 3003, 8178, 10422, 41310, 70338, 70338, 103833, 141669, 151590, 285390, 385578, 542793, 804870, 816750, 950418,$ mientras que $n+1$ divide $3^n+5^n$ , para $n=3, 7, 75, 2355, 11475, 31995, 57075, 80311, 196185, 215325, 335115, 991875$ - Una vez más, he probado esta divisibilidad para $n\le 10^6$ .
0 votos
@TimRatigan en retrospectiva, no veo cómo podemos decir $m=n/3\equiv 1\pmod 8$ . Todo lo que puedo ver $\mod 8$ si $m$ es impar es que $m^2\equiv 1\pmod 8$ . ¿Puedes poner algo de luz aquí?
0 votos
@IanMateus Considerar $\pmod{16}$ . Para los casos de impar, $3^n+5^n\equiv 8\pmod{16}$ Así que $n^2-1\not\equiv 0\pmod{16}$ . Ahora, aquí hice un pequeño woopsie ayer (era tarde para mí también), y deduje $n\equiv 3\pmod 8$ mientras que $n\equiv \pm 3\pmod 8$ pero sigue siendo más fuerte que $n\equiv 1\pmod 2$ .
1 votos
He creado esta habitación Se invita a todos los que quieran discutir este problema.
0 votos
@TimRatigan Espero que no arruine los cálculos que has hecho hasta ahora. Esta sección es ilegible, podemos discutir el problema en el sala de chat He creado.
2 votos
@YiorgosS.Smyrlis: tu comentario con todos los números fue presa del problema descrito en esta pregunta . Esto terminó por desordenar el formato de todos los demás comentarios a esta pregunta. Recuerda no tener un bloque de 79 caracteres sin espacios en blanco en los comentarios o en el texto.
0 votos
Puede demostrar que si $n$ es de la forma mencionada, entonces $$x\equiv 0 \mod 3 \\ x\equiv 3 \mod 5 \\ x\equiv 3,5 \mod 8$$ . Sólo pensé que alguien podría encontrar esto útil para dar una prueba de la conjetura de TI82.
0 votos
@benh podría añadir una prueba del segundo hecho en esta habitación ? No lo sabía, ¿cómo se deshizo de los casos $n\equiv 2\pmod 5$ y $n\equiv 0 \pmod 5$ ?
0 votos
No he pensado demasiado en esto sino que simplemente he comprobado a ciegas que no hay soluciones con $3 < n < 10^{11}$ . A efectos de comprobación: el $n$ entre $10^{10}$ y $10^{11}$ que satisface las condiciones enumeradas en el post de Benh y que satisface $n-1 \mid 5^n + 3^n$ son $n \in \{{11438893323, 12494973573, 17840109573, 20088721653, 38691236493, 39957787323, 55781538813, 72693298293\}}$ .
0 votos
@TI82 ¿puedo preguntarte dónde has encontrado esta cuestión y si sabes que existe una solución?
0 votos
Podría ampliar la secuencia que dio Tapio para $n < 10^{12}$ . Los siguientes términos son $\{107039949093,124041632013,132870562083,142742866893 ,413695605213,458667117363,540389902923,542759619213,604352071323,624162295053,669529883253,737854668243 ,932496818253\}$