OFFSET
1,2
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
FORMULA
A333235(a(n)) = n.
MAPLE
N:= 1000: # for a(1) to a(A025528(N))
R:= NULL: p:= 2:
while p < N do
R:= R, seq(p^k, k=1..ilog[p](N));
p:= nextprime(p);
od:
L:= sort([R]):
M:= nops(L):
f:= proc(n) local F, t;
F:= ifactors(n)[2];
mul(ListTools:-BinarySearch(L, t[1]^t[2]), t=F)
end proc:
V:= Vector(M): count:= 0:
for n from 1 while count < M do
v:= f(n);
if v <= M and V[v] = 0 then count:= count+1; V[v]:= n fi;
od:
convert(V, list);
CROSSREFS
KEYWORD
nonn,look
AUTHOR
Robert Israel, Feb 11 2021
STATUS
approved