A306479 Squarefree composite numbers m such that rad(p-1) = rad(m-1) for every prime p dividing m.
1729, 46657, 1525781251
Offset: 1
Links
- Carlos Rivera, Puzzle 969. Rad(m - 1) = Rad(phi(m)), The Prime Puzzles & Problems Connection.
Programs
-
PARI
rad(n) = factorback(factorint(n)[, 1]); \\ A007947 isok(m) = {if ((m==1) || isprime(m) || !issquarefree(m), return(0)); my(f = factor(m)[,1], r = rad(m-1)); for (i=1, #f, if (rad(f[i]-1) != r, return (0));); return (1);} \\ Michel Marcus, Feb 18 2019
Comments