OFFSET
1,1
COMMENTS
The radical part is actually sqrt(A007913(A000037(n))) where A007913(m) is the squarefree part of m. - Michel Marcus, Jun 26 2020
How does this sequence differ from A048942? The definitions of both sequences are identical, but the second comment in A048942 states the terms differ from n = 14 onwards. - Felix Fröhlich, Jun 16 2022
LINKS
S. R. Finch, Class number theory
Steven R. Finch, Class number theory [Cached copy, with permission of the author]
Eric Weisstein's World of Mathematics, Fundamental Unit.
PROG
(PARI) f(n) = {if (issquare(n), return (0)); if (!issquarefree(n), m = core(n), m = n); my(u = abs(2*polcoeff(lift(bnfinit(x^2-m, 1).fu[1]), 0))); if (u^2==1, return (1)); if (u^2==4, return (sqrtint((u^2+4)/m)); ); if (u^2 < 4, return((u^2+4)/n)); my(v2 = [(u^2-4)/m, (u^2+4)/m]); sqrtint(vecmin(select(x->denominator(x)==1, v2))); }
lista(nn) = apply(f, select(x->!issquare(x), [1..nn])); \\ Michel Marcus, Jun 25 2020; corrected Jun 16 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved