OFFSET
1,2
PROG
(Magma) [Min([k: k in [1..5*10^5] | 1 + #[d: d in Divisors(k) | -Modexp(d, k div d, k) mod k eq d] eq n]): n in [1..6]];
(PARI) f(k) = sumdiv(k, d, -Mod(d, k)^(k/d) == d);
list(len) = {my(v = vector(len), c = 0, k = 1, i); while(c < len, i = f(k); if(i <= len && v[i] == 0, c++; v[i] = k); k++); v; } \\ Amiram Eldar, Nov 06 2025
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Juri-Stepan Gerasimov, Nov 03 2025
EXTENSIONS
a(10)-a(11) from Amiram Eldar, Nov 06 2025
STATUS
approved
