OFFSET
1,3
LINKS
Robert Israel, Table of n, a(n) for n = 1..216
EXAMPLE
a(3) = 37 because 37*3^3 + 1 = 1000 = 10^3.
MAPLE
f:= proc(n) local S, x;
S:= subs(1=n^3+1, map(t -> rhs(op(t)), [msolve(x^n-1, n^3)]));
x:= min(S);
(x^n-1)/n^3
end proc:
f(1):= 1:
map(f, [$1..20]); # Robert Israel, Aug 26 2020
MATHEMATICA
a={}; Do[k = 2; While[ !IntegerQ[(k^n - 1)/n^3], k++ ]; AppendTo[a, (k^n-1)/n^3], {n, 1, 20}]; a
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Sep 10 2012
STATUS
approved