
COMMENTS

Notice: Not all a(n) are 1 or primes, the first example is a(11) = 50521, it equals 19*2659.
a(2n) is a product of powers of Bernoulli irregular primes (A000928), with the exception of n = 0,1,2,3,4,5,7.
a(2n+1) is a product of powers of Euler irregular primes (A120337), with the exception of n = 0,1,2.
Conjectures: All terms are squarefree, and there are infinitely many n such that a(n) is prime.
a(n) = 1 iff n is in the set {0, 1, 2, 3, 4, 5, 6, 8, 10, 14}.
a(n) is prime for n = {7, 9, 12, 16, 17, 18, 26, 34, 36, 38, 39, 42, 49, 74, 114, 118, ...}.
All prime factors of a(n) are irregular primes (Bernoulli or Euler) and with an irregular pair to n: (61, 7), (277, 9), (19, 11), (2659, 11), (691, 12), (43, 13), (967, 13), (47, 15), (4241723, 15), (3617, 16), (228135437, 17), (43867, 18), (79, 19), (349, 19), (84224971, 19), ...
Number of ns such that a prime p divides a(n) is the irregular index of p, for example, 67 divides both a(27) and a(58), so it has irregular index two.
a(149) is the first a(n) which is not completely factored (with a 202digit composite remaining).
