cp's OEIS Frontend

A028233 If n = p_1^e_1 * ... * p_k^e_k, p_1 < ... < p_k primes, then a(n) = p_1^e_1, with a(1) = 1.

%I A028233 #55 Nov 03 2023 06:44:01
%S A028233 1,2,3,4,5,2,7,8,9,2,11,4,13,2,3,16,17,2,19,4,3,2,23,8,25,2,27,4,29,2,
%T A028233 31,32,3,2,5,4,37,2,3,8,41,2,43,4,9,2,47,16,49,2,3,4,53,2,5,8,3,2,59,
%U A028233 4,61,2,9,64,5,2,67,4,3,2,71,8,73,2,3,4,7,2,79,16,81,2,83,4,5,2
%N A028233 If n = p_1^e_1 * ... * p_k^e_k, p_1 < ... < p_k primes, then a(n) = p_1^e_1, with a(1) = 1.
%C A028233 Highest power of smallest prime dividing n. - _Reinhard Zumkeller_, Apr 09 2015
%H A028233 Reinhard Zumkeller, <a href="/A028233/b028233.txt">Table of n, a(n) for n = 1..10000</a> (first 1000 terms from T. D. Noe)
%F A028233 a(n) = A020639(n)^A067029(n). - _Reinhard Zumkeller_, May 13 2006
%F A028233 a(n) = A141809(n,1). - _Reinhard Zumkeller_, Jun 04 2012
%F A028233 a(n) = n / A028234(n). - _Antti Karttunen_, May 29 2017
%e A028233 From _Muniru A Asiru_, Jan 27 2018: (Start)
%e A028233 If n=10, then a(10) = 2 since 10 = 2^1*5^1.
%e A028233 If n=16, then a(16) = 16 since 16 = 2^4.
%e A028233 If n=29, then a(29) = 29 since 29 = 29^1.
%e A028233 (End)
%p A028233 A028233 := proc(n)
%p A028233     local spf,pf;
%p A028233     if n = 1 then
%p A028233         return 1 ;
%p A028233     end if;
%p A028233     spf := A020639(n) ;
%p A028233     for pf in ifactors(n)[2] do
%p A028233         if pf[1] = spf then
%p A028233             return pf[1]^pf[2] ;
%p A028233         end if;
%p A028233     end do:
%p A028233 end proc: # _R. J. Mathar_, Jul 09 2016
%p A028233 # second Maple program:
%p A028233 a:= n-> `if`(n=1, 1, (i->i[1]^i[2])(sort(ifactors(n)[2])[1])):
%p A028233 seq(a(n), n=1..100);  # _Alois P. Heinz_, Jan 29 2018
%t A028233 a[n_] := Power @@ First[ FactorInteger[n]]; Table[a[n], {n, 1, 86}] (* _Jean-François Alcover_, Dec 01 2011 *)
%o A028233 (Haskell)
%o A028233 a028233 = head . a141809_row
%o A028233 -- _Reinhard Zumkeller_, Jun 04 2012, Aug 17 2011
%o A028233 (PARI) a(n)=if(n>1,n=factor(n);n[1,1]^n[1,2],1) \\ _Charles R Greathouse IV_, Apr 26 2012
%o A028233 (Python)
%o A028233 from sympy import factorint
%o A028233 def a(n):
%o A028233     f = factorint(n)
%o A028233     return 1 if n==1 else min(f)**f[min(f)] # _Indranil Ghosh_, May 12 2017
%o A028233 (Scheme)
%o A028233 ;; Naive implementation of A020639 is given under that entry. All of these functions could be also defined with definec to make them faster on the later calls. See http://oeis.org/wiki/Memoization#Scheme
%o A028233 (define (A028233 n) (if (< n 2) n (let ((lpf (A020639 n))) (let loop ((m lpf) (n (/ n lpf))) (cond ((not (zero? (modulo n lpf))) m) (else (loop (* m lpf) (/ n lpf)))))))) ;; _Antti Karttunen_, May 29 2017
%o A028233 (GAP) List(List(List(List([1..10^3], Factors), Collected), i -> i[1]), j -> j[1]^j[2]); # _Muniru A Asiru_, Jan 27 2018
%Y A028233 Cf. A020639, A006530, A034684, A034699, A053585.
%Y A028233 See also A028234.
%Y A028233 Cf. A008475.
%Y A028233 Cf. A141809.
%K A028233 nonn,nice,easy
%O A028233 1,2
%A A028233 _Marc LeBrun_
%E A028233 Name edited to include a(1) = 1 by _Franklin T. Adams-Watters_, Jan 27 2018