OFFSET
1,2
LINKS
Antti Karttunen, Table of n, a(n) for n = 1..285
PROG
(PARI)
upto_e = 101; \\ 101 --> 211 terms
v018819 = vector(upto_e, n, A018819(n));
A050377(n) = factorback(apply(e -> v018819[e], factor(n)[, 2]));
A283980(n) = {my(f=factor(n)); prod(i=1, #f~, my(p=f[i, 1], e=f[i, 2]); if(p==2, 6, nextprime(p+1))^e)}; \\ From A283980
A330688list(e) = { my(lista = List([1, 2]), i=2, u = 2^e, t, m=0, v025487); while(lista[i] != u, if(2*lista[i] <= u, listput(lista, 2*lista[i]); t = A283980(lista[i]); if(t <= u, listput(lista, t))); i++); v025487 = vecsort(Vec(lista)); lista = List([]); for(i=1, #v025487, if((t=A050377(v025487[i]))>m, listput(lista, t); m=t)); Vec(lista); };
v330688 = A330688list(upto_e);
A330688(n) = v330688[n];
CROSSREFS
KEYWORD
nonn
AUTHOR
Antti Karttunen, Dec 28 2019
STATUS
approved