OFFSET
1,2
COMMENTS
Terms from Robert G. Wilson v.
LINKS
Robert Israel, Table of n, a(n) for n = 1..1007
MAPLE
f:= proc(n) local x, y;
for x from 2 to 10 do
y:= x^n;
if floor(y/10^ilog10(y)) = 1 then return x^n fi
od
end proc:
f(1):= 1:
map(f, [$1..50]); # Robert Israel, Aug 13 2019
MATHEMATICA
a = {}; Do[k = 2; While[First[IntegerDigits[k^n]] != 1, k++ ]; a = Append[a, k^n], {n, 2, 25}]; a (* Robert G. Wilson v *)
PROG
(PARI) a(n) = {if (n==1, return (1)); k=2; while (! (ispower(k, n) && (digits(k)[1] == 1)), k++); k; } \\ Michel Marcus, Mar 18 2015
(Magma) m:=1; sol:=[1]; for n in [2..24] do k:=2; while Reverse(Intseq(k^n))[1] ne 1 do; k:=k+1; end while; sol[m+1]:=k^n; m:=m+1; end for; sol; // Marius A. Burtea, Aug 15 2019
(Python)
print(1, 1)
n = 1
while n < 20:
n, p = n+1, 2
s = str(p**n)
while s[0] != "1":
p = p+1
s = str(p**n)
print(n, p**n) # A.H.M. Smeets, Aug 16 2019
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Amarnath Murthy, Feb 05 2002
STATUS
approved