OFFSET
1,2
COMMENTS
a(n) < 3^n.
LINKS
Max Alekseyev, Table of n, a(n) for n = 1..796 (first 100 terms from Alois P. Heinz)
FORMULA
MAPLE
M:= proc(n) M(n):= numtheory[divisors](3^n-1) minus U(n-1) end:
U:= proc(n) U(n):= `if`(n=0, {}, M(n) union U(n-1)) end:
a:= n-> min(M(n)[]):
seq(a(n), n=1..60);
MATHEMATICA
M[n_] := M[n] = Divisors[3^n - 1]~Complement~U[n - 1];
U[n_] := U[n] = If[n == 0, {}, M[n]~Union~U[n - 1]];
a[n_] := Min[M[n]];
Table[a[n], {n, 1, 60}] (* Jean-François Alcover, Oct 24 2022, after Alois P. Heinz *)
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 27 2012
STATUS
approved