OFFSET
1,4
LINKS
PROG
(PARI)
up_to = 2310;
DirSqrt(v) = {my(n=#v, u=vector(n)); u[1]=1; for(n=2, n, u[n]=(v[n]/v[1] - sumdiv(n, d, if(d>1&&d<n, u[d]*u[n/d], 0)))/2); u}; \\ From A317937.
A034386(n) = prod(i=1, primepi(n), prime(i));
A342002(n) = { my(s=0, m=1, p=2, e); while(n, e = (n%p); m *= p^(e>0); s += (e/p); n = n\p; p = nextprime(1+p)); (s*m); };
vA346103aux = DirSqrt(vector(up_to, n, A342920(n)));
A346103(n) = numerator(vA346103aux[n]);
CROSSREFS
KEYWORD
sign,frac
AUTHOR
Antti Karttunen, Jul 09 2021
STATUS
approved