OFFSET
1,1
COMMENTS
a(n) = prime(n) if prime(n) is in A019546.
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 100: # to get a(1)..a(N)
pn:= ithprime(N):
count:= 0:
digs:= [2, 3, 5, 7]:
for d from 1 while count < N do
for m from 4^d to 2*4^d-1 while count < N do
L:= convert(m, base, 4);
n:= add(digs[L[i]+1]*10^(i-1), i=1..nops(L)-1);
ps:= select(p -> p <= pn and not assigned(A[p]), numtheory:-factorset(n));
count:= count + nops(ps);
for p in ps do A[p]:= n od:
od od:
seq(A[ithprime(i)], i=1..N); # Robert Israel, Dec 27 2018
CROSSREFS
KEYWORD
base,easy,nonn
AUTHOR
Amarnath Murthy, Nov 07 2005
EXTENSIONS
More terms from Nick Woods (njw130(AT)psu.edu), Apr 25 2006
a(44) and a(50) corrected by Robert Israel, Dec 27 2018
STATUS
approved