

A319813


a(n) is the smallest a such that n is divisible by a^n + 1, or 0 if no such a exists.


1



1, 1, 2, 0, 4, 0, 6, 0, 2, 3, 10, 0, 12, 0, 14, 0, 16, 0, 18, 0, 5, 0, 22, 0, 4, 5, 2, 0, 28, 0, 30, 0, 32, 13, 34, 0, 36, 0, 17, 0, 40, 0, 42, 0, 14, 0, 46, 0, 6, 3, 50, 0, 52, 0, 19, 0, 8, 17, 58, 0, 60, 0, 5, 0, 64, 0, 66, 0, 68, 0, 70, 0, 72, 31, 14, 0, 76
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,3


COMMENTS

a(n) = 0 iff n is even and 1 is not a square modulo n, that is, n is even and not in A008784. For other n > 2, 2 <= a(n) <= n  1.
a(p) = p  1 for primes p. For composite n, a(n) = n  1 iff gcd(n, phi(n)) = 1, that is, n is in A050384.


LINKS



FORMULA

For n = 9, 9 is divisible by a^9 + 1 implies a == 2 (mod 3), so a(9) = 2.
For n = 10, 10 is divisible by a^10 + 1 implies a == 3, 7 (mod 10), so a(10) = 3.
For n = 34, 34 is divisible by a^34 + 1 implies a == 13, 21 (mod 34), so a(34) = 13.


PROG

(PARI) a(n) = if(!(n%2)&&!issquare(Mod(1, n)), 0, my(i=1); while(Mod(i, n)^n!=n1, i++); i)


CROSSREFS

Cf. A074792 (a^n  1 instead of a^n + 1).


KEYWORD

nonn


AUTHOR



STATUS

approved



