|
%I #21 Sep 08 2022 08:46:20
%S 1,1,3,1,5,3,7,1,3,5,11,3,13,7,5,1,17,3,19,5,7,11,23,3,5,13,3,7,29,5,
%T 31,1,11,17,7,3,37,19,13,5,41,7,43,11,5,23,47,3,7,5,17,13,53,3,11,7,
%U 19,29,59,5,61,31,7,1,13,11,67,17,23,7,71,3,73,37,5,19,11,13,79,5,3,41,83,7,17,43
%N Greatest odd noncomposite divisor of n.
%C The sequence can be constructed replacing the 2's with 1's in A006530.
%F a(n) = A006530(A000265(n)). - _Michel Marcus_, Mar 15 2018
%F a(n) = 1 for n in A000079 (powers of 2). - _Michel Marcus_, Apr 14 2018
%p A299766 := proc(n)
%p local dvs ;
%p dvs := sort(convert(numtheory[factorset](n) minus {2},list)) ;
%p if nops(dvs) >= 1 then
%p op(-1,dvs) ;
%p else
%p 1;
%p end if;
%p end proc: # _R. J. Mathar_, May 05 2018
%t Array[FactorInteger[#][[-1, 1]] /. 2 -> 1 &, 86] (* _Michael De Vlieger_, Mar 16 2018 *)
%o (PARI) a(n) = my(m = n >> valuation(n, 2)); if (m==1, 1, vecmax(factor(m)[,1])); \\ _Michel Marcus_, Mar 15 2018
%o (Magma) [#f eq 1 select 1 else f[#f][1] where f is Factorization(2*n): n in [1..90]]; // _Vincenzo Librandi_, Mar 20 2018
%Y Cf. A000265 (odd part), A008578 (noncomposite numbers), A006005 (odd noncomposite numbers), A006530 (greatest noncomposite divisor).
%K nonn,easy
%O 1,3
%A _Omar E. Pol_, Mar 12 2018
|
