A028233 -id:A028233 - cp's OEIS frontend

A006530 Gpf(n): greatest prime dividing n, for n >= 2; a(1)=1.

1, 2, 3, 2, 5, 3, 7, 2, 3, 5, 11, 3, 13, 7, 5, 2, 17, 3, 19, 5, 7, 11, 23, 3, 5, 13, 3, 7, 29, 5, 31, 2, 11, 17, 7, 3, 37, 19, 13, 5, 41, 7, 43, 11, 5, 23, 47, 3, 7, 5, 17, 13, 53, 3, 11, 7, 19, 29, 59, 5, 61, 31, 7, 2, 13, 11, 67, 17, 23, 7, 71, 3, 73, 37, 5, 19, 11, 13, 79, 5, 3, 41, 83, 7, 17, 43

Offset: 1

N. J. A. Sloane

The initial term a(1)=1 is purely conventional: The unit 1 is not a prime number, although it has been considered so in the past. 1 is the empty product of prime numbers, thus 1 has no largest prime factor. - Daniel Forgues, Jul 05 2011
Greatest noncomposite number dividing n, (cf. A008578). - Omar E. Pol, Aug 31 2013
Conjecture: Let a, b be nonzero integers and f(n) denote the maximum prime factor of a*n + b if a*n + b <> 0 and f(n)=0 if a*n + b=0 for any integer n. Then the set {n, f(n), f(f(n)), ...} is finite of bounded size. - M. Farrokhi D. G., Jan 10 2021

M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 844.
D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Section IV.1.
H. L. Montgomery, Ten Lectures on the Interface Between Analytic Number Theory and Harmonic Analysis, Amer. Math. Soc., 1996, p. 210.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from T. D. Noe)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972 [alternative scanned copy].
K. Alladi and P. Erdős, On an additive arithmetic function, Pacific J. Math., Volume 71, Number 2 (1977), 275-294. MR 0447086 (56 #5401).
G. Back and M. Caragiu, The greatest prime factor and recurrent sequences, Fib. Q., 48 (2010), 358-362.
A. E. Brouwer, Two number theoretic sums, Stichting Mathematisch Centrum. Zuivere Wiskunde, Report ZW 19/74 (1974): 3 pages. [Cached copy, included with the permission of the author]
Paul Erdős, Andrew Granville, Carl Pomerance and Claudia Spiro, On the normal behavior of the iterates of some arithmetic functions, Analytic number theory, Birkhäuser Boston, 1990, pp. 165-204.
Paul Erdős, Andrew Granville, Carl Pomerance and Claudia Spiro, On the normal behavior of the iterates of some arithmetic functions, Analytic number theory, Birkhäuser Boston, 1990, pp. 165-204. [Annotated copy with A-numbers]
Mats Granvik, Mathematica program to check the conjectured formula, Feb 28 2021.
Nathan McNew, The Most Frequent Values of the Largest Prime Divisor Function, Exper. Math., 2017, Vol. 26, No. 2, 210-224.
OEIS Wiki, Greatest prime factor of n
H. P. Robinson, Letter to N. J. A. Sloane, Oct 1981
David Singmaster, Letter to N. J. A. Sloane, Oct 03 1982.
Eric Weisstein's World of Mathematics, Greatest Prime Factor
Index entries for "core" sequences

Cf. A000040, A020639 (smallest prime divisor), A034684, A028233, A034699, A053585.
See also A032742, A052126, A070087, A070089, A061395, A175723.
Cf. A046670 (partial sums), A104350 (partial products).
See A385503 for "popular" primes.

[ #f eq 0 select 1 else f[ #f][1] where f is Factorization(n): n in [1..86] ]; // Klaus Brockhaus, Oct 23 2008

with(numtheory,divisors); A006530 := proc(n) local i,t1,t2,t3,t4,t5; t1 := divisors(n); t2 := convert(t1,list); t3 := sort(t2); t4 := nops(t3); t5 := 1; for i from 1 to t4 do if isprime(t3[t4+1-i]) then return t3[t4+1-i]; fi; od; 1; end;
# alternative
A006530 := n->max(1,op(numtheory[factorset](n))); # Peter Luschny, Nov 02 2010

Table[ FactorInteger[n][[ -1, 1]], {n, 100}] (* Ray Chandler, Nov 12 2005 and modified by Robert G. Wilson v, Jul 16 2014 *)

A006530(n)=if(n>1,vecmax(factor(n)[,1]),1) \\ Edited to cover n=1. - M. F. Hasler, Jul 30 2015

from sympy import factorint
def a(n): return 1 if n == 1 else max(factorint(n))
print([a(n) for n in range(1, 87)]) # Michael S. Branicky, Aug 08 2022

def A006530(n): return list(factor(n))[-1][0] if n > 1 else 1
print([A006530(n) for n in range(1, 87)])  # Peter Luschny, Jan 07 2024

;; The following uses macro definec for the memoization (caching) of the results. A naive implementation of A020639 can be found under that entry. It could be also defined with definec to make it faster on the later calls. See http://oeis.org/wiki/Memoization#Scheme
(definec (A006530 n) (let ((spf (A020639 n))) (if (= spf n) spf (A006530 (/ n spf)))))
;; Antti Karttunen, Mar 12 2017

a(n) = A027748(n, A001221(n)) = A027746(n, A001222(n)); a(n)^A071178(n) = A053585(n). - Reinhard Zumkeller, Aug 27 2011
a(n) = A000040(A061395(n)). - M. F. Hasler, Jan 16 2015
a(n) = n + 1 - Sum_{k=1..n} (floor((k!^n)/n) - floor(((k!^n)-1)/n)). - Anthony Browne, May 11 2016
n/a(n) = A052126(n). - R. J. Mathar, Oct 03 2016
If A020639(n) = n [when n is 1 or a prime] then a(n) = n, otherwise a(n) = a(A032742(n)). - Antti Karttunen, Mar 12 2017
a(n) has average order Pi^2*n/(12 log n) [Brouwer]. See also A046670. - N. J. A. Sloane, Jun 26 2017

Edited by M. F. Hasler, Jan 16 2015

A020639 Lpf(n): least prime dividing n (when n > 1); a(1) = 1. Or, smallest prime factor of n, or smallest prime divisor of n.

Original entry on oeis.org

1, 2, 3, 2, 5, 2, 7, 2, 3, 2, 11, 2, 13, 2, 3, 2, 17, 2, 19, 2, 3, 2, 23, 2, 5, 2, 3, 2, 29, 2, 31, 2, 3, 2, 5, 2, 37, 2, 3, 2, 41, 2, 43, 2, 3, 2, 47, 2, 7, 2, 3, 2, 53, 2, 5, 2, 3, 2, 59, 2, 61, 2, 3, 2, 5, 2, 67, 2, 3, 2, 71, 2, 73, 2, 3, 2, 7, 2, 79, 2, 3, 2, 83, 2, 5, 2, 3, 2, 89, 2, 7, 2, 3, 2, 5, 2, 97

Offset: 1

David W. Wilson

Also, the largest number of distinct integers such that all their pairwise differences are coprime to n. - Max Alekseyev, Mar 17 2006
The unit 1 is not a prime number (although it has been considered so in the past). 1 is the empty product of prime numbers, thus 1 has no least prime factor. - Daniel Forgues, Jul 05 2011
a(n) = least m > 0 for which n! + m and n - m are not relatively prime. - Clark Kimberling, Jul 21 2012
For n > 1, a(n) = the smallest k > 1 that divides n. - Antti Karttunen, Feb 01 2014
For n > 1, records are at prime indices. - Zak Seidov, Apr 29 2015
The initials "lpf" might be mistaken for "largest prime factor" (A009190), using "spf" for "smallest prime factor" would avoid this. - M. F. Hasler, Jul 29 2015
n = 89 is the first index > 1 for which a(n) differs from the smallest k > 1 such that (2^k + n - 2)/k is an integer. - M. F. Hasler, Aug 11 2015
From Stanislav Sykora, Jul 29 2017: (Start)
For n > 1, a(n) is also the smallest k, 1 < k <= n, for which the binomial(n,k) is not divisible by n.
Proof: (A) When k and n are relatively prime then binomial(n,k) is divisible by n because k*binomial(n,k) = n*binomial(n-1,k-1). (B) When gcd(n,k) > 1, one of its prime factors is the smallest; let us denote it p, p <= k, and consider the binomial(n,p) = (1/p!)*Product_{i=0..p-1} (n-i). Since p is a divisor of n, it cannot be a divisor of any of the remaining numerator factors. It follows that, denoting as e the largest e > 0 such that p^e|n, the numerator is divisible by p^e but not by p^(e+1). Hence, the binomial is divisible by p^(e-1) but not by p^e and therefore not divisible by n. Applying (A), (B) to all considered values of k completes the proof. (End)
From Bob Selcoe, Oct 11 2017, edited by M. F. Hasler, Nov 06 2017: (Start)
a(n) = prime(j) when n == J (mod A002110(j)), n, j >= 1, where J is the set of numbers <= A002110(j) with smallest prime factor = prime(j). The number of terms in J is A005867(j-1). So:
a(n) = 2 when n == 0 (mod 2);
a(n) = 3 when n == 3 (mod 6);
a(n) = 5 when n == 5 or 25 (mod 30);
a(n) = 7 when n == 7, 49, 77, 91, 119, 133, 161 or 203 (mod 210);
etc. (End)
For n > 1, a(n) is the leftmost term, other than 0 or 1, in the n-th row of A127093. - Davis Smith, Mar 05 2019

D. S. Mitrinovic et al., Handbook of Number Theory, Kluwer, Section IV.1.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (terms 1..10000 from T. D. Noe)
A. E. Brouwer, Two number theoretic sums, Stichting Mathematisch Centrum. Zuivere Wiskunde, Report ZW 19/74 (1974): 3 pages. [Copy included with the permission of the author.]
OEIS Wiki, Least prime factor of n
David Singmaster, Letter to N. J. A. Sloane, Oct 3 1982.
Eric Weisstein's World of Mathematics, Least Prime Factor
Index entries for "core" sequences

Cf. A090368 (bisection).
Cf. A000040, A009190, A006530, A034684, A028233, A034699, A053585.
See also A032742, A055396, A068319, A088377, A007978, A053669, A117818.
Cf. A046669 (partial sums), A072486 (partial products).
Cf. A002110, A005867.
Cf. A127093.

a020639 n = spf a000040_list where
  spf (p:ps) | n < p^2      = n
             | mod n p == 0 = p
             | otherwise    = spf ps
-- Reinhard Zumkeller, Jul 13 2011

A020639 := proc(n) if n = 1 then 1; else min(op(numtheory[factorset](n))) ; end if; end proc: seq(A020639(n),n=1..20) ; # R. J. Mathar, Oct 25 2010

f[n_]:=FactorInteger[n][[1,1]]; Join[{1}, Array[f,120,2]]  (* Robert G. Wilson v, Apr 06 2011 *)
Join[{1}, Table[If[EvenQ[n], 2, FactorInteger[n][[1,1]]], {n, 2, 120}]] (* Zak Seidov, Nov 17 2013 *)
Riffle[Join[{1},Table[FactorInteger[n][[1,1]],{n,3,101,2}]],2] (* Harvey P. Dale, Dec 16 2021 *)

A020639(n) = { vecmin(factor(n)[,1]) } \\ [Will yield an error for n = 1.] - R. J. Mathar, Mar 02 2012

A020639(n)=if(n>1, if(n>n=factor(n,0)[1,1], n, factor(n)[1,1]), 1) \\ Avoids complete factorization if possible. Often the smallest prime factor can be found quickly even if it is larger than primelimit. If factoring takes too long for large n, use debugging level >= 3 (\g3) to display the smallest factor as soon as it is found. - M. F. Hasler, Jul 29 2015

from sympy import factorint
def a(n): return 1 if n == 1 else min(factorint(n))
print([a(n) for n in range(1, 98)]) # Michael S. Branicky, Dec 09 2021

def A020639_list(n) : return [1] + [prime_divisors(n)[0] for n in (2..n)]
A020639_list(97) # Peter Luschny, Jul 16 2012

[trial_division(n) for n in (1..100)] # Giuseppe Coppoletta, May 25 2016

(define (A020639 n) (if (< n 2) n (let loop ((k 2)) (cond ((zero? (modulo n k)) k) (else (loop (+ 1 k))))))) ;; Antti Karttunen, Feb 01 2014

A014673(n) = a(A032742(n)); A115561(n) = a(A054576(n)). - Reinhard Zumkeller, Mar 10 2006
A028233(n) = a(n)^A067029(n). - Reinhard Zumkeller, May 13 2006
a(n) = A027746(n,1) = A027748(n,1). - Reinhard Zumkeller, Aug 27 2011
For n > 1: a(n) = A240694(n,2). - Reinhard Zumkeller, Apr 10 2014
a(n) = A000040(A055396(n)) = n / A032742(n). - Antti Karttunen, Mar 07 2017
a(n) has average order n/(2 log n) [Brouwer] - N. J. A. Sloane, Sep 03 2017

Deleted wrong comment from M. Lagneau in 2012, following an observation by Gionata Neri. - M. F. Hasler, Aug 11 2015
Edited by M. F. Hasler, Nov 06 2017
Expanded definition to make this easier to find. - N. J. A. Sloane, Sep 21 2020

A047994 Unitary totient (or unitary phi) function uphi(n).

Original entry on oeis.org

1, 1, 2, 3, 4, 2, 6, 7, 8, 4, 10, 6, 12, 6, 8, 15, 16, 8, 18, 12, 12, 10, 22, 14, 24, 12, 26, 18, 28, 8, 30, 31, 20, 16, 24, 24, 36, 18, 24, 28, 40, 12, 42, 30, 32, 22, 46, 30, 48, 24, 32, 36, 52, 26, 40, 42, 36, 28, 58, 24, 60, 30, 48, 63, 48, 20, 66, 48, 44, 24, 70

Offset: 1

N. J. A. Sloane

A divisor d of n is called a unitary divisor if gcd(d, n/d) = 1. Define gcd*(k,n) to be the largest divisor d of k that is also a unitary divisor of n (that is, such that gcd(d, n/d) = 1). The unitary totient function a(n) = number of k with 1 <= k <= n such that gcd*(k,n) = 1. - N. J. A. Sloane, Aug 08 2021
Unitary convolution of A076479 and A000027. - R. J. Mathar, Apr 13 2011
Multiplicative with a(p^e) = p^e - 1. - N. J. A. Sloane, Apr 30 2013

			a(12) = a(3)*a(4) = 2*3 = 6.

T. D. Noe, Table of n, a(n) for n = 1..10000
Eckford Cohen, Arithmetical functions associated with the unitary divisors of an integer, Math. Zeitschr. 74 (1960) 66-80
Steven R. Finch, Unitarism and infinitarism.
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
M. Lal, Iterates of the unitary totient function, Math. Comp., 28 (1974), 301-302.
R. J. Mathar, Survey of Dirichlet Series of Multiplicative Arithmetic Functions, arXiv:1106.4038 [math.NT], 2011, Remark 43.
László Tóth, On the Bi-Unitary Analogues of Euler's Arithmetical Function and the Gcd-Sum Function, JIS 12 (2009), Article 09.5.2.
László Tóth, A survey of gcd-sum functions, J. Int. Seq. 13 (2010), Article 10.8.1.

Cf. A000010, A000203, A001221, A003557, A049865, A003271, A028233, A076479.

a047994 n = f n 1 where
   f 1 uph = uph
   f x uph = f (x `div` sppf) (uph * (sppf - 1)) where sppf = a028233 x
-- Reinhard Zumkeller, Aug 17 2011

A047994 := proc(n)
    local a, f;
    a := 1 ;
    for f in ifactors(n)[2] do
        a := a*(op(1,f)^op(2,f)-1) ;
    end do:
    a ;
end proc:
seq(A047994(n),n=1..20) ; # R. J. Mathar, Dec 22 2011

uphi[n_] := (Times @@ (Table[ #[[1]]^ #[[2]] - 1, {1} ] & /@ FactorInteger[n]))[[1]]; Table[ uphi[n], {n, 2, 75}] (* Robert G. Wilson v, Sep 06 2004 *)
uphi[n_] := If[n==1, 1, Product[{p, e} = pe; p^e-1, {pe, FactorInteger[n]}] ]; Array[uphi, 80] (* Jean-François Alcover, Nov 17 2018 *)

A047994(n)=my(f=factor(n)~); prod(i=1, #f, f[1, i]^f[2, i]-1);

for(n=1, 100, print1(direuler(p=2, n, (1 - 2*X + p*X^2)/(1-X)/(1-p*X))[n], ", ")) \\ Vaclav Kotesovec, Jun 15 2020

from math import prod
from sympy import factorint
def A047994(n): return prod(p**e-1 for p, e in factorint(n).items()) # Chai Wah Wu, Sep 24 2021

If n = Product p_i^e_i, uphi(n) = Product (p_i^e_i - 1).
a(n) = A000010(n)*A000203(A003557(n))/A003557(n). - Velin Yanev and Charles R Greathouse IV, Aug 23 2017
From Amiram Eldar, May 29 2020: (Start)
a(n) = Sum_{d|n, gcd(d, n/d) = 1} (-1)^omega(d) * n/d.
Sum_{d|n, gcd(d, n/d) = 1} a(d) = n.
a(n) >= phi(n) = A000010(n), with equality if and only if n is squarefree (A005117). (End)
Sum_{k=1..n} a(k) ~ c * Pi^2 * n^2 / 12, where c = A065464 = Product_{primes p} (1 - 2/p^2 + 1/p^3). - Vaclav Kotesovec, Jun 15 2020
Dirichlet g.f.: zeta(s-1) * zeta(s) * Product_{p prime} (1 - 2/p^s + 1/p^(2*s-1)). - Amiram Eldar, May 22 2025

More terms from Jud McCranie

A067029 Exponent of least prime factor in prime factorization of n, a(1)=0.

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 3, 2, 1, 1, 1, 5, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 2, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 2, 6, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1, 3, 1, 1

Offset: 1

Reinhard Zumkeller, Feb 17 2002

Even bisection is A001511: a(2n) = A007814(n) + 1. - Ralf Stephan, Jan 31 2004
Number of occurrences of the smallest part in the partition with Heinz number n. The Heinz number of a partition p = [p_1, p_2, ..., p_r] is defined as Product_{j=1..r} (p_j-th prime) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). Example: a(24)=3 because the partition with Heinz number 24 = 3*2*2*2 is [2,1,1,1]. - Emeric Deutsch, Oct 02 2015
Together with A028234 is useful for defining sequences that are multiplicative with a(p^e) = f(e), as recurrences of the form: a(1) = 1 and for n > 1, a(n) = f(A067029(n)) * a(A028234(n)). - Antti Karttunen, May 29 2017

			a(18) = a(2^1 * 3^2) = 1.

T. D. Noe, Table of n, a(n) for n = 1..10000
Project Euler, Problem 779: Prime factor and exponent, (2022).

Cf. A051903, A020639, A028233, A034684, A071178, first column of A124010, A247180.

a067029 = head . a124010_row
-- Reinhard Zumkeller, Jul 05 2013, Jun 04 2012

A067029 := proc(n)
    local f,lp,a;
    a := 0 ;
    lp := n+1 ;
    for f in ifactors(n)[2] do
        p := op(1,f) ;
        if p < lp then
            a := op(2,f) ;
            lp := p;
        fi;
    end do:
    a ;
end proc: # R. J. Mathar, Jul 08 2015
seq(ifelse(n = 1, 0, ifactors(n)[2][1][2]), n = 1..90); # Peter Luschny, Jun 15 2025

Join[{0},Table[FactorInteger[n][[1,2]],{n,2,100}]] (* Harvey P. Dale, Oct 14 2011 *)

a(n) = if (n==1, 0, factor(n)[1,2]); \\ Michel Marcus, May 15 2017

from sympy import factorint
def a(n):
    f=factorint(n)
    return 0 if n==1 else f[min(f)] # Indranil Ghosh, May 15 2017

;; 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
(define (A067029 n) (if (< n 2) 0 (let ((mp (A020639 n))) (let loop ((e 0) (n (/ n mp))) (cond ((integer? n) (loop (+ e 1) (/ n mp))) (else e)))))) ;;  Antti Karttunen, May 29 2017

a(n) = A124010(n,1). - Reinhard Zumkeller, Aug 27 2011
A028233(n) = A020639(n)^a(n). - Reinhard Zumkeller, May 13 2006
a(A247180(n)) = 1. - Reinhard Zumkeller, Nov 23 2014
Asymptotic mean: lim_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=1} (Product_{i=1..k-1} (1 - 1/prime(i)))/(prime(k)-1) = 1/(prime(1)-1) + (1-1/prime(1))*(1/(prime(2)-1) + (1-1/prime(2))*(1/(prime(3)-1) + (1-1/prime(3))*( ... ))) = 1.6125177915... - Amiram Eldar, Oct 26 2021

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

Original entry on oeis.org

1, 1, 1, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 7, 5, 1, 1, 9, 1, 5, 7, 11, 1, 3, 1, 13, 1, 7, 1, 15, 1, 1, 11, 17, 7, 9, 1, 19, 13, 5, 1, 21, 1, 11, 5, 23, 1, 3, 1, 25, 17, 13, 1, 27, 11, 7, 19, 29, 1, 15, 1, 31, 7, 1, 13, 33, 1, 17, 23, 35, 1, 9, 1, 37, 25, 19, 11, 39, 1, 5, 1, 41, 1, 21

Offset: 1

Marc LeBrun

Together with A067029 is useful for defining sequences that are multiplicative with a(p^e) = f(e), as recurrences of the form: a(1) = 1 and for n > 1, a(n) = f(A067029(n)) * a(A028234(n)). - Antti Karttunen, May 29 2017

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A001221, A001222, A002110, A005867, A020639, A028233, A067029.

a := List(List(List(List([1..10^3], Factors), Collected), i -> i[1]), j -> j[1]^j[2]);;
A028234 := List([1..Length(a)],i->i/a[i]); # Muniru A Asiru, Jan 27 2018

a028234 n = n `div` a028233 n  -- Reinhard Zumkeller, Mar 27 2013

a[n_] := n / Power @@ First[FactorInteger[n]]; Table[a[n], {n, 1, 84}] (* Jean-François Alcover, Jun 12 2012 *)

a(n) = {my(f = factor(n)); if (#f~, f[1, 1] = 1); factorback(f);} \\ Michel Marcus, Feb 11 2016

from sympy import factorint
def a(n):
    f = factorint(n)
    return 1 if n==1 else n/(min(f)**f[min(f)]) # Indranil Ghosh, May 12 2017

(define (A028234 n) (/ n (A028233 n))) ;; Needs also code from A020639 and A028233. - Antti Karttunen, May 29 2017

a(n) = n / A028233(n).
A001221(a(n)) = A001221(n)-1; A001222(a(n)) = A001222(n)-A067029(n). - Reinhard Zumkeller, May 13 2006
Sum_{k=1..n} a(k) ~ c * n^2, where c = (1/2) * Sum_{k>=0} A005867(k)/(prime(k+1)*(prime(k+1)+1)*A002110(k)) = 0.114813... . - Amiram Eldar, Nov 19 2022

Edited name to include a(1) = 1 by Franklin T. Adams-Watters, Jan 27 2018

A059897 Symmetric square array read by antidiagonals: A(n,k) is the product of all factors that occur in one, but not both, of the Fermi-Dirac factorizations of n and k.

Original entry on oeis.org

1, 2, 2, 3, 1, 3, 4, 6, 6, 4, 5, 8, 1, 8, 5, 6, 10, 12, 12, 10, 6, 7, 3, 15, 1, 15, 3, 7, 8, 14, 2, 20, 20, 2, 14, 8, 9, 4, 21, 24, 1, 24, 21, 4, 9, 10, 18, 24, 28, 30, 30, 28, 24, 18, 10, 11, 5, 27, 2, 35, 1, 35, 2, 27, 5, 11, 12, 22, 30, 36, 40, 42, 42, 40, 36, 30, 22, 12, 13, 24, 33

Offset: 1

Marc LeBrun, Feb 06 2001

Old name: Square array read by antidiagonals: T(i,j) = product prime(k)^(Ei(k) XOR Ej(k)) where Ei and Ej are the vectors of exponents in the prime factorizations of i and j; XOR is the bitwise operation on binary representation of the exponents.
Analogous to multiplication, with XOR replacing +.
From Peter Munn, Apr 01 2019: (Start)
(1) Defines an abelian group whose underlying set is the positive integers. (2) Every element is self-inverse. (3) For all n and k, A(n,k) is a divisor of n*k. (4) The terms of A050376, sometimes called Fermi-Dirac primes, form a minimal set of generators. In ordered form, it is the lexicographically earliest such set.
The unique factorization of positive integers into products of distinct terms of the group's lexicographically earliest minimal set of generators seems to follow from (1) (2) and (3).
From (1) and (2), every row and every column of the table is a self-inverse permutation of the positive integers. Rows/columns numbered by nonmembers of A050376 are compositions of earlier rows/columns.
It is a subgroup of the equivalent group over the nonzero integers, which has -1 as an additional generator.
As generated by A050376, the subgroup of even length words is A000379. The complementary set of odd length words is A000028.
The subgroup generated by A000040 (the primes) is A005117 (the squarefree numbers).
(End)
Considered as a binary operation, the result is (the squarefree part of the product of its operands) times the square of (the operation's result when applied to the square roots of the square parts of its operands). - Peter Munn, Mar 21 2022

			A(864,1944) = A(2^5*3^3,2^3*3^5) = 2^(5 XOR 3) * 3^(3 XOR 5) = 2^6 * 3^6 = 46656.
The top left 12 X 12 corner of the array:
   1,  2,  3,  4,  5,  6,  7,  8,   9,  10,  11,  12
   2,  1,  6,  8, 10,  3, 14,  4,  18,   5,  22,  24
   3,  6,  1, 12, 15,  2, 21, 24,  27,  30,  33,   4
   4,  8, 12,  1, 20, 24, 28,  2,  36,  40,  44,   3
   5, 10, 15, 20,  1, 30, 35, 40,  45,   2,  55,  60
   6,  3,  2, 24, 30,  1, 42, 12,  54,  15,  66,   8
   7, 14, 21, 28, 35, 42,  1, 56,  63,  70,  77,  84
   8,  4, 24,  2, 40, 12, 56,  1,  72,  20,  88,   6
   9, 18, 27, 36, 45, 54, 63, 72,   1,  90,  99, 108
  10,  5, 30, 40,  2, 15, 70, 20,  90,   1, 110, 120
  11, 22, 33, 44, 55, 66, 77, 88,  99, 110,   1, 132
  12, 24,  4,  3, 60,  8, 84,  6, 108, 120, 132,   1
From _Peter Munn_, Apr 04 2019: (Start)
The subgroup generated by {6,8,10}, the first three integers > 1 not in A050376, has the following table:
    1     6     8    10    12    15    20   120
    6     1    12    15     8    10   120    20
    8    12     1    20     6   120    10    15
   10    15    20     1   120     6     8    12
   12     8     6   120     1    20    15    10
   15    10   120     6    20     1    12     8
   20   120    10     8    15    12     1     6
  120    20    15    12    10     8     6     1
(End)

Antti Karttunen, Table of n, a(n) for n = 1..10440; the first 144 antidiagonals of the array
Eric Weisstein's World of Mathematics, Group, Square Part, Squarefree Part.

Cf. A000040, A003987, A003991, A028233, A028234, A050376, A059896, A089913, A207901, A268387, A284577, A302033.
Cf. A284567 (A000142 or A003418-analog for this operation).
Rows/columns: A073675 (2), A120229 (3), A120230 (4), A307151 (5), A307150 (6), A307266 (8), A307267 (24).
Particularly significant subgroups or cosets: A000028, A000379, A003159, A005117, A030229, A252895. See also the lists in A329050, A352273.
Sequences that relate this sequence to multiplication: A000188, A007913, A059895.

a[i_, i_] = 1;
a[i_, j_] := Module[{f1 = FactorInteger[i], f2 = FactorInteger[j], e1, e2}, e1[] = 0; Scan[(e1[#[[1]]] = #[[2]])&, f1]; e2[] = 0; Scan[(e2[#[[1]]] = #[[2]])&, f2]; Times @@ (#^BitXor[e1[#], e2[#]]& /@ Union[f1[[All, 1]], f2[[All, 1]]])];
Table[a[i - j + 1, j], {i, 1, 15}, {j, 1, i}] // Flatten (* Jean-François Alcover, Jun 19 2018 *)

T(n,k) = {if (n==1, return (k)); if (k==1, return (n)); my(fn=factor(n), fk=factor(k)); vp = setunion(fn[,1]~, fk[,1]~); prod(i=1, #vp, vp[i]^(bitxor(valuation(n, vp[i]), valuation(k, vp[i]))));} \\ Michel Marcus, Apr 03 2019

T(i, j) = {if(gcd(i, j) == 1, return(i * j)); if(i == j, return(1)); my(f = vecsort(concat(factor(i)~, factor(j)~)), t = 1, res = 1); while(t + 1 <= #f, if(f[1, t] == f[1, t+1], res *= f[1, t] ^ bitxor(f[2, t] , f[2, t+1]); t+=2; , res*= f[1, t]^f[2, t]; t++; ) ); if(t == #f, res *= f[1, #f] ^ f[2, #f]); res } \\ David A. Corneth, Apr 03 2019

A059897(n,k) = if(n==k, 1, core(n*k) * A059897(core(n,1)[2],core(k,1)[2])^2) \\ Peter Munn, Mar 21 2022

(define (A059897 n) (A059897bi (A002260 n) (A004736 n)))
(define (A059897bi a b) (let loop ((a a) (b b) (m 1)) (cond ((= 1 a) (* m b)) ((= 1 b) (* m a)) ((equal? (A020639 a) (A020639 b)) (loop (A028234 a) (A028234 b) (* m (expt (A020639 a) (A003987bi (A067029 a) (A067029 b)))))) ((< (A020639 a) (A020639 b)) (loop (/ a (A028233 a)) b (* m (A028233 a)))) (else (loop a (/ b (A028233 b)) (* m (A028233 b)))))))
;; Antti Karttunen, Apr 11 2017

For all x, y >= 1, A(x,y) * A059895(x,y)^2 = x*y. - Antti Karttunen, Apr 11 2017
From Peter Munn, Apr 01 2019: (Start)
A(n,1) = A(1,n) = n
A(n, A(m,k)) = A(A(n,m), k)
A(n,n) = 1
A(n,k) = A(k,n)
if i_1 <> i_2 then A(A050376(i_1), A050376(i_2)) = A050376(i_1) * A050376(i_2)
if A(n,k_1) = n * k_1 and A(n,k_2) = n * k_2 then A(n, A(k_1,k_2)) = n * A(k_1,k_2)
(End)
T(k, m) = k*m for coprime k and m. - David A. Corneth, Apr 03 2019
if A(n*m,m) = n, A(n*m,k) = A(n,k) * A(m,k) / k. - Peter Munn, Apr 04 2019
A(n,k) = A007913(n*k) * A(A000188(n), A000188(k))^2. - Peter Munn, Mar 21 2022

New name from Peter Munn, Mar 21 2022

A034699 Largest prime power factor of n.

Original entry on oeis.org

1, 2, 3, 4, 5, 3, 7, 8, 9, 5, 11, 4, 13, 7, 5, 16, 17, 9, 19, 5, 7, 11, 23, 8, 25, 13, 27, 7, 29, 5, 31, 32, 11, 17, 7, 9, 37, 19, 13, 8, 41, 7, 43, 11, 9, 23, 47, 16, 49, 25, 17, 13, 53, 27, 11, 8, 19, 29, 59, 5, 61, 31, 9, 64, 13, 11, 67, 17, 23, 7, 71, 9, 73, 37, 25, 19, 11, 13, 79

Offset: 1

David W. Wilson

n divides lcm(1, 2, ..., a(n)).
a(n) = A210208(n,A073093(n)) = largest term of n-th row in A210208. - Reinhard Zumkeller, Mar 18 2012
a(n) = smallest m > 0 such that n divides A003418(m). - Thomas Ordowski, Nov 15 2013
a(n) = n when n is a prime power (A000961). - Michel Marcus, Dec 03 2013
Conjecture: For all n between two consecutive prime numbers, all a(n) are different. - I. V. Serov, Jun 19 2019
Disproved with between p=prime(574) = 4177 and prime(575) = 4201, a(4180) = a(4199) = 19. See A308752. - Michel Marcus, Jun 19 2019
Conjecture: For any N > 0, there exist numbers n and m, N < n < n+a(n) <= m, such that all n..m are composite and a(n) = a(m). - I. V. Serov, Jun 21 2019
Conjecture: For all n between two consecutive prime numbers, all (-1)^n*a(n) are different. Checked up to 5*10^7. - I. V. Serov, Jun 23 2019
Disproved: between p = prime(460269635) = 10120168277 and p = prime(460269636) = 10120168507 the numbers n = 10120168284 and m = 10120168498 form a pair such that (-1)^n*a(n) = (-1)^m*a(m) = 107. - L. Joris Perrenet, Jan 05 2020
a(n) = cardinality of smallest set on which idempotence of order n+1 (f^{n+1} = f) differs from idempotence of order e for 2 <= e <= n (see von Eitzen link for proof); derivable from A245501. - Mark Bowron, May 22 2025

Antti Karttunen, Table of n, a(n) for n = 1..65537 (first 1000 terms from T. D. Noe)
Hagen von Eitzen, Least k to distinguish f^n = f on {1,...,k} from f^e = f for every 1 < e < n is largest prime power factor of n-1?, Math Stack Exchange.

Cf. A006530, A010055, A020639, A027750, A034684, A028233, A051283, A052128, A053585, A057110, A060818, A038610, A081805, A088387, A088388, A100574, A210208, A245501, A284600, A308752.

Numbers n ordered by a(n): A305325.

a034699 = last . a210208_row
-- Reinhard Zumkeller, Mar 18 2012, Feb 14 2012

f[n_] := If[n == 1, 1, Max[ #[[1]]^#[[2]] & /@ FactorInteger@n]]; Array[f, 79] (* Robert G. Wilson v, Sep 02 2006 *)
Array[Max[Power @@@ FactorInteger@ #] &, 79] (* Michael De Vlieger, Jul 26 2018 *)

a(n) = if(1==n,n,my(f=factor(n)); vecmax(vector(#f[, 1], i, f[i, 1]^f[i, 2]))); \\ Charles R Greathouse IV, Nov 20 2012, check for a(1) added by Antti Karttunen, Aug 06 2018

A034699(n) = if(1==n,n,fordiv(n, d, if(isprimepower(n/d), return(n/d)))); \\ Antti Karttunen, Aug 06 2018

from sympy import factorint
def A034699(n): return max((p**e for p, e in factorint(n).items()), default=1) # Chai Wah Wu, Apr 17 2023

If n = p_1^e_1 *...* p_k^e_k, p_1 < ... < p_k primes, then a(n) = Max_i p_i^e_i.
a(n) = A088387(n)^A088388(n). - Antti Karttunen, Jul 22 2018
a(n) = n/A284600(n) = n - A081805(n) = A034684(n) + A100574(n). - Antti Karttunen, Aug 06 2018
a(n) = a(m) iff m = d*a(n), where d is a divisor of A038610(a(n)). - I. V. Serov, Jun 19 2019

A008475 If n = Product (p_j^k_j) then a(n) = Sum (p_j^k_j) (a(1) = 0 by convention).

Original entry on oeis.org

0, 2, 3, 4, 5, 5, 7, 8, 9, 7, 11, 7, 13, 9, 8, 16, 17, 11, 19, 9, 10, 13, 23, 11, 25, 15, 27, 11, 29, 10, 31, 32, 14, 19, 12, 13, 37, 21, 16, 13, 41, 12, 43, 15, 14, 25, 47, 19, 49, 27, 20, 17, 53, 29, 16, 15, 22, 31, 59, 12, 61, 33, 16, 64, 18, 16, 67, 21, 26, 14, 71, 17, 73

Offset: 1

Olivier Gérard

For n>1, a(n) is the minimal number m such that the symmetric group S_m has an element of order n. - Ahmed Fares (ahmedfares(AT)my-deja.com), Jun 26 2001

If gcd(u,w) = 1, then a(u*w) = a(u) + a(w); behaves like logarithm; compare A001414 or A056239. - Labos Elemer, Mar 31 2003

			a(180) = a(2^2 * 3^2 * 5) = 2^2 + 3^2 + 5 = 18.

F. J. Budden, The Fascination of Groups, Cambridge, 1972; pp. 322, 573.
József Sándor, Dragoslav S. Mitrinovic and Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Chapter IV, p. 147.
T. Z. Xuan, On some sums of large additive number theoretic functions (in Chinese), Journal of Beijing normal university, No. 2 (1984), pp. 11-18.

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 10000 terms from T. D. Noe)
John Bamberg, Grant Cairns and Devin Kilminster, The crystallographic restriction, permutations and Goldbach's conjecture, Amer. Math. Monthly, Vol. 110, No. 3 (March 2003), pp. 202-209.
Roger B. Eggleton and William P. Galvin, Upper Bounds on the Sum of Principal Divisors of an Integer, Mathematics Magazine, Vol. 77, No. 3 (Jun., 2004), pp. 190-200.

Cf. A001414, A000961, A005117, A051613, A072691, A081402, A081403, A081404, A027748, A124010, A001221, A028233, A034684, A053585, A159077, A023888, A078771, A092509, A286875.

See A222416 for the variant with a(1)=1.

a008475 1 = 0
a008475 n = sum $ a141809_row n
-- Reinhard Zumkeller, Jan 29 2013, Oct 10 2011

A008475 := proc(n) local e,j; e := ifactors(n)[2]:
add(e[j][1]^e[j][2], j=1..nops(e)) end:
seq(A008475(n), n=1..60); # Peter Luschny, Jan 17 2010

f[n_] := Plus @@ Power @@@ FactorInteger@ n; f[1] = 0; Array[f, 73]

for(n=1,100,print1(sum(i=1,omega(n), component(component(factor(n),1),i)^component(component(factor(n),2),i)),","))

a(n)=local(t);if(n<1,0,t=factor(n);sum(k=1,matsize(t)[1],t[k,1]^t[k,2])) /* Michael Somos, Oct 20 2004 */

A008475(n) = { my(f=factor(n)); vecsum(vector(#f~,i,f[i,1]^f[i,2])); }; \\ Antti Karttunen, Nov 17 2017

from sympy import factorint
def a(n):
    f=factorint(n)
    return 0 if n==1 else sum([i**f[i] for i in f]) # Indranil Ghosh, May 20 2017

Additive with a(p^e) = p^e.
a(A000961(n)) = A000961(n); a(A005117(n)) = A001414(A005117(n)).
a(n) = Sum_{k=1..A001221(n)} A027748(n,k) ^ A124010(n,k) for n>1. - Reinhard Zumkeller, Oct 10 2011
a(n) = Sum_{k=1..A001221(n)} A141809(n,k) for n > 1. - Reinhard Zumkeller, Jan 29 2013
Sum_{k=1..n} a(k) ~ (Pi^2/12)* n^2/log(n) + O(n^2/log(n)^2) (Xuan, 1984). - Amiram Eldar, Mar 04 2021

A100995 If n is a prime power p^m, m >= 1, then m, otherwise 0.

Original entry on oeis.org

0, 1, 1, 2, 1, 0, 1, 3, 2, 0, 1, 0, 1, 0, 0, 4, 1, 0, 1, 0, 0, 0, 1, 0, 2, 0, 3, 0, 1, 0, 1, 5, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 2, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 6, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 4, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0

Offset: 1

Reinhard Zumkeller, Nov 26 2004

nonn

Calculate matrix powers: (A175992^1)/1 - (A175992^2)/2 + (A175992^3)/3 - (A175992^4)/4 + ... Then the nonzero values of a(n) are found as reciprocals in the first column. Compare this to the Taylor series for log(1+x) = (x)/1 - (x^2)/2 + (x^3)/3 - (x^4)/4 + ... Therefore it is natural to write 0, 1/1, 1/1, 1/2, 1/1, 0, 1/1, 1/3, 1/2, 0, 1/1, ... Raising n to a such power gives A014963. - Mats Granvik, Gary W. Adamson, Apr 04 2011
The Dirichlet series that generates the reciprocals of this sequence is the logarithm of the Riemann zeta function. - Mats Granvik, Gary W. Adamson, Apr 04 2011
Number of automorphisms of the finite field with n elements, or 0 if the field does not exist. For n=p^k where p is a prime and k is an integer, the automorphism group of the finite field with n elements is a cyclic group of order k generated by the Frobenius endomorphism. - Yancheng Lu, Jan 11 2021

Daniel Forgues, Table of n, a(n) for n = 1..100000

Cf. A028233, A069513, A010055.

a100995 n = f 0 n where
   f e 1 = e
   f e x = if r > 0 then 0 else f (e + 1) x'
           where (x', r) = divMod x p
   p = a020639 n
-- Reinhard Zumkeller, Mar 19 2013

f:= proc(n) local F;
    F:= ifactors(n)[2];
    if nops(F) = 1 then F[1][2]
    else 0
    fi
end proc:
map(f, [$1..100]); # Robert Israel, Jun 09 2015

ppm[n_]:=If[PrimePowerQ[n],FactorInteger[n][[1,2]],0]; Array[ppm,110] (* Harvey P. Dale, Mar 03 2014 *)
a=Table[Limit[Sum[If[Mod[n, k] == 0, MoebiusMu[n/k]/(n/k)^(s - 1)/(1 - 1/n^(s - 1)), 0], {k, 1, n}], s -> 1], {n, 1, 105}];
Numerator[a]*Denominator[a] (* Mats Granvik, Jun 09 2015 *)
a = FullSimplify[Table[MangoldtLambda[n]/Log[n], {n, 1, 105}]]
Numerator[a]*Denominator[a] (* Mats Granvik, Jun 09 2015 *)

{a(n) = my(t); if( n<1, 0, t = factor(n); if( [1,2] == matsize(t), t[1,2], 0))} /* Michael Somos, Aug 15 2012 */

{a(n) = my(t); if( n<1, 0, if( t = isprimepower(n), t))} /* Michael Somos, Aug 15 2012 */

A100994(n) = A014963(n)^a(n);
a(A000961(n)) = A025474(n).
a(n) = Sum_{d|n, gcd(d, n/d) = 1} (-1)^omega(n/d) * bigomega(d). - Ilya Gutkovskiy, Apr 15 2021

Edited by Daniel Forgues and N. J. A. Sloane, Aug 18 2009

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

Original entry on oeis.org

1, 2, 3, 4, 5, 3, 7, 8, 9, 5, 11, 3, 13, 7, 5, 16, 17, 9, 19, 5, 7, 11, 23, 3, 25, 13, 27, 7, 29, 5, 31, 32, 11, 17, 7, 9, 37, 19, 13, 5, 41, 7, 43, 11, 5, 23, 47, 3, 49, 25, 17, 13, 53, 27, 11, 7, 19, 29, 59, 5, 61, 31, 7, 64, 13, 11, 67, 17, 23, 7, 71, 9, 73, 37, 25, 19, 11, 13, 79

Offset: 1

Frederick Magata (frederick.magata(AT)uni-muenster.de), Jan 19 2000

Let p be the largest prime dividing n, a(n) is the largest power of p dividing n.

			a(42)=7 because 42=2*3*7, a(144)=9 because 144=16*9=2^4*3^2.

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)

Cf. A020639, A006530, A034684, A028233, A051119, A008475.

Different from A034699.

a053585 = last . a141809_row  -- Reinhard Zumkeller, Jan 29 2013

a:= n-> `if`(n=1, 1, (i->i[1]^i[2])(sort(ifactors(n)[2])[-1])):
seq(a(n), n=1..100);  # Alois P. Heinz, Nov 03 2023

Table[Power @@ Last @ FactorInteger @ n, {n, 79}] (* Jean-François Alcover, Apr 01 2011 *)

a(n)=if(n>1,my(f=factor(n)); f[#f~,1]^f[#f~,2],1) \\ Charles R Greathouse IV, Nov 10 2015

from sympy import factorint, primefactors
def a(n):
    if n==1: return 1
    p = primefactors(n)[-1]
    return p**factorint(n)[p] # Indranil Ghosh, May 19 2017

a(n) = A006530(n)^A071178(n). - Reinhard Zumkeller, Aug 27 2011

a(n) = A141809(n,A001221(n)). - Reinhard Zumkeller, Jan 29 2013

More terms from Andrew Gacek (andrew(AT)dgi.net), Apr 20 2000

cp's OEIS Frontend

A006530 Gpf(n): greatest prime dividing n, for n >= 2; a(1)=1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Python

SageMath

Scheme

Formula

Extensions

A020639 Lpf(n): least prime dividing n (when n > 1); a(1) = 1. Or, smallest prime factor of n, or smallest prime divisor of n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Sage

Sage

Scheme

Formula

Extensions

A047994 Unitary totient (or unitary phi) function uphi(n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A067029 Exponent of least prime factor in prime factorization of n, a(1)=0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Scheme

Formula

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

Views

Author

Keywords

Comments

Links

Crossrefs

Programs