%I #19 Mar 31 2023 14:14:43
%S 1,1,1,1,1,3,1,1,1,5,1,3,1,7,5,1,1,3,1,5,7,11,1,3,1,13,1,7,1,3,1,1,11,
%T 17,7,3,1,19,13,5,1,3,1,11,5,23,1,3,1,5,17,13,1,3,11,7,19,29,1,3,1,31,
%U 7,1,13,3,1,17,23,5,1,3,1,37,5,19,11,3,1,5,1,41,1,3,17,43,29,11,1,3,13,23
%N a(n) is the second smallest prime factor of n, or 1 if n is a prime power.
%C Least prime factor of {n divided by the maximal power of the least prime factor of n}. - after the original name of the sequence.
%C a(n) = A020639(A028234(n)).
%C a(n) = 1 iff n is a prime power: a(A000961(n))=1 and a(A024619(n))>1.
%H Antti Karttunen, <a href="/A119288/b119288.txt">Table of n, a(n) for n = 1..10000</a>
%F A010055(n) = 0^(a(n)-1). - _Reinhard Zumkeller_, May 13 2006
%t Join[{1},Table[Which[PrimePowerQ[n],1,True,FactorInteger[n][[2,1]]],{n,2,100}]] (* _Harvey P. Dale_, Feb 08 2020 *)
%o (PARI) a(n) = if (isprimepower(n) || (n==1), 1, my(f=factor(n)[,1]); f[2]); \\ _Michel Marcus_, Mar 01 2023
%o (Python)
%o from sympy import primefactors
%o def A119288(n): return 1 if len(s:=primefactors(n)) <= 1 else sorted(s)[1] # _Chai Wah Wu_, Mar 31 2023
%Y Cf. A014673, A020639, A028234, A292269.
%K nonn
%O 1,6
%A _Reinhard Zumkeller_, May 13 2006
%E Name changed by _Antti Karttunen_, Oct 04 2017