login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A335324
Square part of 4th-power-free part of n.
6
1, 1, 1, 4, 1, 1, 1, 4, 9, 1, 1, 4, 1, 1, 1, 1, 1, 9, 1, 4, 1, 1, 1, 4, 25, 1, 9, 4, 1, 1, 1, 1, 1, 1, 1, 36, 1, 1, 1, 4, 1, 1, 1, 4, 9, 1, 1, 1, 49, 25, 1, 4, 1, 9, 1, 4, 1, 1, 1, 4, 1, 1, 9, 4, 1, 1, 1, 4, 1, 1, 1, 36, 1, 1, 25, 4, 1, 1, 1, 1, 1, 1, 1, 4, 1
OFFSET
1,4
COMMENTS
Equivalently, biquadratefree (4th-power-free) part of square part of n.
Multiplicative. The terms are squares of squarefree numbers (A062503).
Every positive integer n is the product of a unique subset S_n of the terms of A050376 (sometimes called Fermi-Dirac primes). a(n) is the product of the members of S_n that are squares of prime numbers (A001248).
LINKS
Eric Weisstein's World of Mathematics, Biquadratefree.
Eric Weisstein's World of Mathematics, Square part.
FORMULA
a(n) = A053165(A008833(n)) = A008833(A053165(n)).
a(n) = A053165(n) / A007913(n).
a(n) = A008833(n) / A008835(n).
n = A007913(n) * a(n) * A008835(n).
a(n) = A225546(A038500(A225546(n))).
a(n^2) = A007913(n)^2.
a(A003961(n)) = A003961(a(n)).
a(A331590(n, k)) = A331590(a(n), a(k)).
a(p^e) = p^(2*floor(e/2) - 4*floor(e/4)). - Amiram Eldar, Jun 01 2020
From Amiram Eldar, Sep 21 2023: (Start)
Dirichlet g.f.: zeta(s) * zeta(2*s-2) * zeta(4*s)/(zeta(2*s) * zeta(4*s-4)).
Sum_{k=1..n} a(k) ~ (4*zeta(3/2)*zeta(4))/(21*zeta(3)) * n^(3/2). (End)
EXAMPLE
Removing the 4th powers from 192 = 2^6 * 3^1 gives 2^(6 - 4) * 3^1 = 2^2 * 3 = 12. So the 4th-power-free part of 192 is 12. The square part of 12 (largest square dividing 12) is 4. So a(192) = 4.
MATHEMATICA
f[p_, e_] := p^(2*Floor[e/2] - 4*Floor[e/4]); a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Jun 01 2020 *)
PROG
(PARI) A053165(n)=my(f=factor(n)); f[, 2]=f[, 2]%4; factorback(f);
a(n) = my(m=A053165(n)); m/core(m); \\ Michel Marcus, Jun 01 2020
(Python)
from math import prod
from sympy import factorint
def A335324(n): return prod(p**(e&2) for p, e in factorint(n).items()) # Chai Wah Wu, Aug 07 2024
CROSSREFS
A007913, A008833, A008835, A053165 are used in formulas defining the sequence.
Column 1 of A352780.
Range of values is A062503.
Positions of 1's: A252895.
Related to A038500 by A225546.
The formula section details how the sequence maps the terms of A003961, A331590.
Sequence in context: A085731 A131301 A350698 * A366245 A083730 A373440
KEYWORD
nonn,easy,mult
AUTHOR
Peter Munn, May 31 2020
STATUS
approved