login
Product of odd divisors of n.
10

%I #31 Jun 26 2022 02:16:30

%S 1,1,3,1,5,3,7,1,27,5,11,3,13,7,225,1,17,27,19,5,441,11,23,3,125,13,

%T 729,7,29,225,31,1,1089,17,1225,27,37,19,1521,5,41,441,43,11,91125,23,

%U 47,3,343,125,2601,13,53,729,3025,7,3249,29,59,225,61,31,250047,1,4225,1089

%N Product of odd divisors of n.

%C Product of rows of triangle A182469. - _Reinhard Zumkeller_, May 01 2012

%H Reinhard Zumkeller, <a href="/A136655/b136655.txt">Table of n, a(n) for n = 1..10000</a>

%F a(p) = p if p noncomposite; a(2^n) = 1; a(pq) = p^2 * q^2 when p, q are odd primes.

%F a(n) = sqrt(n^od(n)/2^ed(n)), where od(n) = number of odd divisors of n = tau(2*n)-tau(n) and ed(n) = number of even divisors of n = 2*tau(n)-tau(2*n). - _Vladeta Jovovic_, Jun 25 2008

%F Also a(n) = A007955(A000265(n)). - David Wilson, Jun 26 2008

%F a(n) = Product_{h == 1 mod 4 and h | n}*Product_{i == 3 mod 4 and i | n}.

%F a(n) = Product_{j == 1 mod 6 and j | n}*Product_{k == 5 mod 6 and k | n}.

%F a(n) = A140210(n)*A140211(n). - _R. J. Mathar_, Jun 27 2008

%F a(n) = A007955(n) / A125911(n).

%p with(numtheory); f:=proc(n) local t1,i,k; t1:=divisors(n); k:=1; for i in t1 do if i mod 2 = 1 then k:=k*i; fi; od; k; end; # _N. J. A. Sloane_, Jul 14 2008

%t Array[Times @@ Select[Divisors@ #, OddQ] &, 66] (* _Michael De Vlieger_, Aug 03 2017 *)

%t a[n_] := (oddpart = n/2^IntegerExponent[n, 2])^(DivisorSigma[0, oddpart]/2); Array[a, 100] (* _Amiram Eldar_, Jun 26 2022 *)

%o (Haskell)

%o a136655 = product . a182469_row -- _Reinhard Zumkeller_, May 01 2012

%o (PARI) a(n) = my(d=divisors(n)); prod(k=1, #d, if (d[k]%2, d[k], 1)); \\ _Michel Marcus_, Aug 04 2017

%Y Cf. A000265, A000593, A007955, A007956, A078701.

%Y Cf. A140210, A140211, A140213, A140214, A140215.

%Y Cf. A125911, A126192.

%Y Cf. A001227, A183063.

%K nonn,easy

%O 1,3

%A _Jonathan Vos Post_, Jun 25 2008

%E More terms from _N. J. A. Sloane_, Jul 14 2008

%E Edited by _N. J. A. Sloane_, Aug 29 2008 at the suggestion of _R. J. Mathar_