login
Squarefree product of all odd primes dividing n, and 1 if n is a power of 2: A099985/2.
23

%I #66 Jul 22 2024 06:32: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