|
%I #8 Aug 26 2018 18:08:09
%S 1,2,3,1,5,6,7,1,1,10,11,3,13,14,15,1,17,2,19,5,21,22,23,3,1,26,1,7,
%T 29,30,31,1,33,34,35,1,37,38,39,5,41,42,43,11,5,46,47,3,1,2,51,13,53,
%U 2,55,7,57,58,59,15,61,62,7,1,65,66,67,17,69,70,71,2
%N a(n) = n/(largest perfect power divisor of n).
%C Not all terms are squarefree numbers; for example, a(500) = 4.
%H Andrew Howroyd, <a href="/A304328/b304328.txt">Table of n, a(n) for n = 1..1000</a>
%F a(n) * A203025(n) = n.
%t Table[n/Last[Select[Divisors[n],#===1||GCD@@FactorInteger[#][[All,2]]>1&]],{n,100}]
%o (PARI) a(n)={my(m=1); fordiv(n, d, if(ispower(d), m=max(m,d))); n/m} \\ _Andrew Howroyd_, Aug 26 2018
%Y Cf. A000961, A001055, A001222, A001597, A001694, A005117, A007716, A007916, A052485, A052486, A059404, A066636, A091050, A160400, A203025, A294068, A303707, A304327, A304339.
%K nonn
%O 1,2
%A _Gus Wiseman_, May 10 2018
|
