|
%I #60 Feb 03 2024 10:15:48
%S 1,1,3,1,5,3,7,1,3,5,11,3,13,7,15,1,17,3,19,5,21,11,23,3,5,13,3,7,29,
%T 15,31,1,33,17,35,3,37,19,39,5,41,21,43,11,15,23,47,3,7,5,51,13,53,3,
%U 55,7,57,29,59,15,61,31,21,1,65,33,67,17,69,35,71,3
%N Squarefree product of all odd primes dividing n, and 1 if n is a power of 2: A099985/2.
%C There are no odd primes dividing n iff n is a power of 2.
%C This sequence coincides with the bisection of A007947 (even indices), which is A099985, dividing out the even prime 2 in the squarefree kernel.
%C a(n) divides A106609(n) for n>=1. - _Alexander R. Povolotsky_, Apr 06 2015
%H Michael De Vlieger, <a href="/A204455/b204455.txt">Table of n, a(n) for n = 1..10000</a>
%F a(n) = A099985(n)/2 = A007947(2*n)/2.
%F a(n) = A000265(A007947(n)) = A007947(A000265(n)). - _Charles R Greathouse IV_, Jan 19 2012
%F Multiplicative with a(p^e)=p for p <> 2 and a(2^e)=1. - _R. J. Mathar_, Jul 02 2013
%F a(n) = Sum_{d|n} phi(d)*mu(2d)^2. - _Ridouane Oudra_, Sep 02 2019
%F From _Richard L. Ollerton_, May 09 2021: (Start)
%F a(n) = Sum_{k=1..n} mu(2*n/gcd(n,k))^2.
%F a(n) = Sum_{k=1..n} mu(2*gcd(n,k))^2*phi(gcd(n,k))/phi(n/gcd(n,k)). (End)
%F Sum_{k=1..n} a(k) ~ c * n^2, where c = (2/5) * Product_{p prime} (1 - 1/(p*(p+1))) = (2/5) * A065463 = 0.281776... . - _Amiram Eldar_, Nov 19 2022
%F a(n) = Sum_{d divides n, d odd} mu(d)^2 * phi(d). - _Peter Bala_, Feb 01 2024
%e a(5)=5 because 5 is a single odd prime.
%e a(9)=3 because 9=3*3 has as squarefree part 3.
%e a(1)=1 because 1 is a power of 2, having no odd primes as a factor.
%p A204455 := proc(n)
%p local p;
%p numtheory[factorset](n) minus {2} ;
%p mul(p,p=%) ;
%p end proc:
%p seq(A204455(n),n=1..40) ; # _R. J. Mathar_, Jan 25 2017
%t f[n_] := Select[First /@ FactorInteger@ n, PrimeQ@ # && OddQ@ # &]; Times @@@ (f /@ Range@ 120) (* _Michael De Vlieger_, Apr 08 2015 *)
%o (PARI) a(n) = {my(f = factor(n)); prod(k=1, #f~, if (f[k,1] % 2, f[k,1], 1));} \\ _Michel Marcus_, Apr 07 2015
%o (PARI) a(n) = factorback(setminus(factorint(n)[, 1]~, [2])) \\ _Jianing Song_, Aug 09 2022
%Y Cf. A099985, A099984, A007947, A000265, A065463, A106609.
%K nonn,mult
%O 1,3
%A _Wolfdieter Lang_, Jan 19 2012
|
