A082695 -id:A082695 - cp's OEIS frontend

A003557 n divided by largest squarefree divisor of n; if n = Product p(k)^e(k) then a(n) = Product p(k)^(e(k)-1), with a(1) = 1.

1, 1, 1, 2, 1, 1, 1, 4, 3, 1, 1, 2, 1, 1, 1, 8, 1, 3, 1, 2, 1, 1, 1, 4, 5, 1, 9, 2, 1, 1, 1, 16, 1, 1, 1, 6, 1, 1, 1, 4, 1, 1, 1, 2, 3, 1, 1, 8, 7, 5, 1, 2, 1, 9, 1, 4, 1, 1, 1, 2, 1, 1, 3, 32, 1, 1, 1, 2, 1, 1, 1, 12, 1, 1, 5, 2, 1, 1, 1, 8, 27, 1, 1, 2, 1, 1, 1, 4, 1, 3, 1, 2, 1, 1, 1, 16, 1, 7

Offset: 1

Marc LeBrun

a(n) is the size of the Frattini subgroup of the cyclic group C_n - Ahmed Fares (ahmedfares(AT)my-deja.com), Jun 07 2001.
Also of the Frattini subgroup of the dihedral group with 2*n elements. - Sharon Sela (sharonsela(AT)hotmail.com), Jan 01 2002
Number of solutions to x^m==0 (mod n) provided that n < 2^(m+1), i.e. the sequence of sequences A000188, A000189, A000190, etc. converges to this sequence. - Henry Bottomley, Sep 18 2001
a(n) is the number of nilpotent elements in the ring Z/nZ. - Laszlo Toth, May 22 2009
The sequence of partial products of a(n) is A085056(n). - Peter Luschny, Jun 29 2009
The first occurrence of n in this sequence is at A064549(n). - Franklin T. Adams-Watters, Jul 25 2014
From Hal M. Switkay, Jul 03 2025: (Start)
For n > 1, a(n) is a proper divisor of n. Thus the sequence n, a(n), a(a(n)), ... eventually becomes 1. This yields a minimal factorization of n as a product of squarefree numbers (A005117), each factor dividing all larger factors, in a factorization that is conjugate to the minimal factorization of n as a product of prime powers (A000961), as follows.
Let f(n,0) = n, and let f(n,k) = a(f(n,k-1)) for k > 0. A051903(n) is the minimal value of k such that f(n,k) = 1. A051903(n) <= log(n)/log(2). Since n/a(n) = A007947(n) is always squarefree by definition, n is a product of squarefree factors in the form Product_{i=1..A051903(n)} [f(n,i-1)/f(n,i)].
The two factorizations correspond to conjugate partitions of bigomega(n) = A001222(n). (End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Steven R. Finch, Idempotents and Nilpotents Modulo n, arXiv:math/0605019 [math.NT], 2006-2017.

Cf. A007947, A062378, A062379, A064549, A300717 (Möbius transform), A326306 (inv. Möbius transf.), A328572.
Sequences that are multiples of this sequence (the other factor of a pointwise product is given in parentheses): A000010 (A173557), A000027 (A007947), A001615 (A048250), A003415 (A342001), A007434 (A345052), A057521 (A071773).
Cf. A082695 (Dgf at s=2), A065487 (Dgf at s=3).
Cf. A000961, A001222, A005117, A051903.

a003557 n = product $ zipWith (^)
                      (a027748_row n) (map (subtract 1) $ a124010_row n)
-- Reinhard Zumkeller, Dec 20 2013

using Nemo
function A003557(n)
    n < 4 && return 1
    q = prod([p for (p, e) ∈ Nemo.factor(fmpz(n))])
    return n == q ? 1 : div(n, q)
end
[A003557(n) for n in 1:90] |> println  # Peter Luschny, Feb 07 2021

[(&+[(Floor(k^n/n)-Floor((k^n-1)/n)): k in [1..n]]): n in [1..100]]; // G. C. Greubel, Nov 02 2018

A003557 := n -> n/ilcm(op(numtheory[factorset](n))):
seq(A003557(n), n=1..98); # Peter Luschny, Mar 23 2011
seq(n / NumberTheory:-Radical(n), n = 1..98); # Peter Luschny, Jul 20 2021

Prepend[ Array[ #/Times@@(First[ Transpose[ FactorInteger[ # ] ] ])&, 100, 2 ], 1 ] (* Olivier Gérard, Apr 10 1997 *)

a(n)=n/factorback(factor(n)[,1]) \\ Charles R Greathouse IV, Nov 17 2014

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

from sympy.ntheory.factor_ import core
from sympy import divisors
def a(n): return n / max(i for i in divisors(n) if core(i) == i)
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Apr 16 2017

from math import prod
from sympy import primefactors
def A003557(n): return n//prod(primefactors(n)) # Chai Wah Wu, Nov 04 2022

def A003557(n) : return n*mul(1/p for p in prime_divisors(n))
[A003557(n) for n in (1..98)] # Peter Luschny, Jun 10 2012

Multiplicative with a(p^e) = p^(e-1). - Vladeta Jovovic, Jul 23 2001
a(n) = n/rad(n) = n/A007947(n) = sqrt(J_2(n)/J_2(rad(n))), where J_2(n) is A007434. - Enrique Pérez Herrero, Aug 31 2010
a(n) = (J_k(n)/J_k(rad(n)))^(1/k), where J_k is the k-th Jordan Totient Function: (J_2 is A007434 and J_3 A059376). - Enrique Pérez Herrero, Sep 03 2010
Dirichlet convolution of A000027 and A097945. - R. J. Mathar, Dec 20 2011
a(n) = A000010(n)/|A023900(n)|. - Eric Desbiaux, Nov 15 2013
a(n) = Product_{k = 1..A001221(n)} (A027748(n,k)^(A124010(n,k)-1)). - Reinhard Zumkeller, Dec 20 2013
a(n) = Sum_{k=1..n}(floor(k^n/n)-floor((k^n-1)/n)). - Anthony Browne, May 11 2016
a(n) = e^[Sum_{k=2..n} (floor(n/k)-floor((n-1)/k))*(1-A010051(k))*Mangoldt(k)] where Mangoldt is the Mangoldt function. - Anthony Browne, Jun 16 2016
a(n) = Sum_{d|n} mu(d) * phi(d) * (n/d), where mu(d) is the Moebius function and phi(d) is the Euler totient function (rephrases formula of Dec 2011). - Daniel Suteu, Jun 19 2018
G.f.: Sum_{k>=1} mu(k)*phi(k)*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Nov 02 2018
Dirichlet g.f.: Product_{primes p} (1 + 1/(p^s - p)). - Vaclav Kotesovec, Jun 24 2020
From Richard L. Ollerton, May 07 2021: (Start)
a(n) = Sum_{k=1..n} mu(n/gcd(n,k))*gcd(n,k).
a(n) = Sum_{k=1..n} mu(gcd(n,k))*(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). (End)
a(n) = A001615(n)/A048250(n) = A003415/A342001(n) = A057521(n)/A071773(n). - Antti Karttunen, Jun 08 2021

Secondary definition added to the name by Antti Karttunen, Jun 08 2021

A001615 Dedekind psi function: n * Product_{p|n, p prime} (1 + 1/p).

Original entry on oeis.org

1, 3, 4, 6, 6, 12, 8, 12, 12, 18, 12, 24, 14, 24, 24, 24, 18, 36, 20, 36, 32, 36, 24, 48, 30, 42, 36, 48, 30, 72, 32, 48, 48, 54, 48, 72, 38, 60, 56, 72, 42, 96, 44, 72, 72, 72, 48, 96, 56, 90, 72, 84, 54, 108, 72, 96, 80, 90, 60, 144, 62, 96, 96, 96, 84, 144, 68, 108, 96

Offset: 1

N. J. A. Sloane

Number of primitive sublattices of index n in generic 2-dimensional lattice; also index of Gamma_0(n) in SL_2(Z).
A generic 2-dimensional lattice L =  consists of all vectors of the form mV + nW, (m,n integers). A sublattice S =  has index |ad-bc| and is primitive if gcd(a,b,c,d) = 1. The generic lattice L has precisely a(2) = 3 sublattices of index 2, namely <2V,W>,  and  (which = ) and so on for other indices.
The sublattices of index n are in 1-to-1 correspondence with matrices [a b; 0 d] with a>0, ad=n, b in [0..d-1]. The number of these is Sum_{d|n} = sigma(n), which is A000203. A sublattice is primitive if gcd(a,b,d) = 1; the number of these is n * product_{p|n} (1+1/p), which is the present sequence.
SL_2(Z) = Gamma is the group of all 2 X 2 matrices [a b; c d] where a,b,c,d are integers with ad-bc = 1 and Gamma_0(N) is usually defined as the subgroup of this for which N|c. But conceptually Gamma is best thought of as the group of (positive) automorphisms of a lattice , its typical element taking V -> aV + bW, W -> cV + dW and then Gamma_0(N) can be defined as the subgroup consisting of the automorphisms that fix the sublattice  of index N. - J. H. Conway, May 05 2001
Dedekind proved that if n = k_i*j_i for i in I represents all the ways to write n as a product, and e_i=gcd(k_i,j_i), then a(n)= sum(k_i / (e_i * phi(e_i)), i in I ) [cf. Dickson, History of the Theory of Numbers, Vol. 1, p. 123].
Also a(n)= number of cyclic subgroups of order n in an Abelian group of order n^2 and type (1,1) (Fricke). - Len Smiley, Dec 04 2001
The polynomial degree of the classical modular equation of degree n relating j(z) and j(nz) is psi(n) (Fricke). - Michael Somos, Nov 10 2006; clarified by Katherine E. Stange, Mar 11 2022
The Mobius transform of this sequence is A063659. - Gary W. Adamson, May 23 2008
The inverse Mobius transform of this sequence is A060648. - Vladeta Jovovic, Apr 05 2009
The Dirichlet inverse of this sequence is A008836(n) * A048250(n). - Álvar Ibeas, Mar 18 2015
The Riemann Hypothesis is true if and only if a(n)/n - e^gamma*log(log(n)) < 0 for any n > 30. - Enrique Pérez Herrero, Jul 12 2011
The Riemann Hypothesis is also equivalent to another inequality, see the Sole and Planat link. - Thomas Ordowski, May 28 2017
An infinitary analog of this sequence is the sum of the infinitary divisors of n (see A049417). - Vladimir Shevelev, Apr 01 2014
Problem: are there composite numbers n such that n+1 divides psi(n)? - Thomas Ordowski, May 21 2017
The sum of divisors d of n such that n/d is squarefree. - Amiram Eldar, Jan 11 2019
Psi(n)/n is a new maximum for each primorial (A002110) [proof in link: Patrick Sole and Michel Planat, Proposition 1 page 2]. - Bernard Schott, May 21 2020
From Jianing Song, Nov 05 2022: (Start)
a(n) is the number of subgroups of C_n X C_n that are isomorphic to C_n, where C_n is the cyclic group of order n. Proof: the number of elements of order n in C_n X C_n is A007434(n) (they are the elements of the form (a,b) in C_n X C_n where gcd(a,b,n) = 1), and each subgroup isomorphic to C_n contains phi(n) generators, so the number of such subgroups is A007434(n)/phi(n) = a(n).
The total number of order-n subgroups of C_n X C_n is A000203(n). (End)

			Let L = <V,W> be a 2-dimensional lattice. The 6 primitive sublattices of index 4 are generated by <4V,W>, <V,4W>, <4V,W+-V>, <2V+W,2W>, <2V,2W+V>. Compare A000203.
G.f. = x + 3*x^2 + 4*x^3 + 6*x^4 + 6*x^5 + 12*x^6 + 8*x^7 + 12*x^8 + 12*x^9 + ...

Tom Apostol, Intro. to Analyt. Number Theory, page 71, Problem 11, where this is called phi_1(n).
David A. Cox, "Primes of the Form x^2 + n y^2", Wiley, 1989, p. 228.
R. Fricke, Die elliptischen Funktionen und ihre Anwendungen, Teubner, 1922, Vol. 2, see p. 220.
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004. See Section B41, p. 147.
B. Schoeneberg, Elliptic Modular Functions, Springer-Verlag, NY, 1974, p. 79.
G. Shimura, Introduction to the Arithmetic Theory of Automorphic Functions, Princeton, 1971, see p. 25, Eq. (1).
N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

T. D. Noe and N. J. A. Sloane, Table of n, a(n) for n = 1..10000
O. Bordelles and B. Cloitre, An Alternating Sum Involving the Reciprocal of Certain Multiplicative Functions, J. Int. Seq. 16 (2013) #13.6.3.
Harriet Fell, Morris Newman, and Edward Ordman, Tables of genera of groups of linear fractional transformations, J. Res. Nat. Bur. Standards Sect. B 67B 1963 61-68.
M. Hampejs, N. Holighaus, L. Toth and C. Wiesmeyr, Representing and counting the subgroups of the group Z_m X Z_n, arXiv:1211.1797 [math.GR], 2012.
W. Hürlimann, Dedekind's arithmetic function and primitive four squares counting functions, Journal of Algebra, Number Theory: Advances and Applications 14:2 (2015), 73-88.
F. A. Lewis et al., Problem 4002, Amer. Math. Monthly, Vol. 49, No. 9, Nov. 1942, pp. 618-619.
E. Pérez Herrero, Recycling Hardy & Wright, Average Order of Dedekind Psi Function, Psychedelic Geometry Blogspot.
Michel Planat, Riemann hypothesis from the Dedekind psi function, arXiv:1010.3239 [math.GM], 2010.
Patrick Sole and Michel Planat, Extreme values of the Dedekind Psi function, to appear in Journal of Combinatorics and Number Theory, arXiv:1011.1825 [math.NT], 2010-2011.
Eric Weisstein's World of Mathematics, Dedekind Function
Wikipedia, Dedekind psi function
Index entries for "core" sequences
Index entries for sequences related to sublattices

Other sequences that count lattices/sublattices: A000203 (with primitive condition removed), A003050 (hexagonal lattice instead), A003051, A054345, A160889, A160891.
Cf. A002110, A027748, A124010, A019269.
Cf. A301594.
Cf. A063659 (Möbius transform), A082020 (average order), A156303 (Euler transform), A173290 (partial sums), A175836 (partial products), A203444 (range).
Cf. A210523 (record values).
Algebraic combinations with other core sequences: A000082, A033196, A175732, A291784, A344695.
Sequences of the form n^k * Product_ {p|n, p prime} (1 + 1/p^k) for k=0..10: A034444 (k=0), this sequence (k=1), A065958 (k=2), A065959 (k=3), A065960 (k=4), A351300 (k=5), A351301 (k=6), A351302 (k=7), A351303 (k=8), A351304 (k=9), A351305 (k=10).
Cf. A082695 (Dgf at s=3), A339925 (Dgf at s=4).

import Data.Ratio (numerator)
a001615 n = numerator (fromIntegral n * (product $
            map ((+ 1) . recip . fromIntegral) $ a027748_row n))
-- Reinhard Zumkeller, Jun 03 2013, Apr 12 2012

m:=75; R:=PowerSeriesRing(Integers(), m); Coefficients(R!( (&+[MoebiusMu(k)^2*x^k/(1-x^k)^2: k in [1..2*m]]) )); // G. C. Greubel, Nov 23 2018

A001615 := proc(n) n*mul((1+1/i[1]),i=ifactors(n)[2]) end; # Mark van Hoeij, Apr 18 2012

Join[{1}, Table[n Times @@ (1 + 1/Transpose[FactorInteger[n]][[1]]), {n, 2, 100}]] (* T. D. Noe, Jun 11 2006 *)
Table[DirichletConvolve[j, MoebiusMu[j]^2, j, n], {n, 100}] (* Jan Mangaldan, Aug 22 2013 *)
a[n_] := n Sum[MoebiusMu[d]^2/d, {d, Divisors[n]}]; (* Michael Somos, Jan 10 2015 *)
Table[n Product[1 + 1/p, {p, Select[Divisors[n], PrimeQ]}], {n, 1, 100}] (* Vaclav Kotesovec, May 08 2021 *)
Table[n DivisorSum[n, MoebiusMu[#]^2/# &], {n, 20}] (* Eric W. Weisstein, Mar 09 2025 *)

{a(n) = if( n<1, 0, direuler( p=2, n, (1 + X) / (1 - p*X)) [n])};

{a(n) = if( n<1, 0, n * sumdiv( n, d, moebius(d)^2 / d))}; /* Michael Somos, Nov 10 2006 */

a(n)=my(f=factor(n)); prod(i=1,#f~, f[i,1]^f[i,2] + f[i,1]^(f[i,2]-1)) \\ Charles R Greathouse IV, Aug 22 2013

a(n) = n * sumdivmult(n, d, issquarefree(d)/d) \\ Charles R Greathouse IV, Sep 09 2014

from math import prod
from sympy import primefactors
def A001615(n):
    plist = primefactors(n)
    return n*prod(p+1 for p in plist)//prod(plist) # Chai Wah Wu, Jun 03 2021

def A001615(n) : return n*mul(1+1/p for p in prime_divisors(n))
[A001615(n) for n in (1..69)] # Peter Luschny, Jun 10 2012

Dirichlet g.f.: zeta(s) * zeta(s-1) / zeta(2*s). - Michael Somos, May 19 2000
Multiplicative with a(p^e) = (p+1)*p^(e-1). - David W. Wilson, Aug 01 2001
a(n) = A003557(n)*A048250(n) = n*A000203(A007947(n))/A007947(n). - Labos Elemer, Dec 04 2001
a(n) = n*Sum_{d|n} mu(d)^2/d, Dirichlet convolution of A008966 and A000027. - Benoit Cloitre, Apr 07 2002
a(n) = Sum_{d|n} mu(n/d)^2 * d. - Joerg Arndt, Jul 06 2011
From Enrique Pérez Herrero, Aug 22 2010: (Start)
a(n) = J_2(n)/J_1(n) = J_2(n)/phi(n) = A007434(n)/A000010(n), where J_k is the k-th Jordan Totient Function.
a(n) = (1/phi(n))*Sum_{d|n} mu(n/d)*d^(b-1), for b=3. (End)
a(n) = n / Sum_{d|n} mu(d)/a(d). - Enrique Pérez Herrero, Jun 06 2012
a(n^k)= n^(k-1) * a(n). - Enrique Pérez Herrero, Jan 05 2013
If n is squarefree, then a(n) = A049417(n) = A000203(n). - Vladimir Shevelev, Apr 01 2014
a(n) = Sum_{d^2 | n} mu(d) * A000203(n/d^2). - Álvar Ibeas, Dec 20 2014
The average order of a(n) is 15*n/Pi^2. - Enrique Pérez Herrero, Jan 14 2012. See Apostol. - N. J. A. Sloane, Sep 04 2017
G.f.: Sum_{k>=1} mu(k)^2*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Oct 25 2018
a(n) = Sum_{d|n} 2^omega(d) * phi(n/d), Dirichlet convolution of A034444 and A000010. - Daniel Suteu, Mar 09 2019
From Richard L. Ollerton, May 07 2021: (Start)
a(n) = Sum_{k=1..n} 2^omega(gcd(n,k)).
a(n) = Sum_{k=1..n} 2^omega(n/gcd(n,k))*phi(gcd(n,k))/phi(n/gcd(n,k)). (End)
a(n) = abs(A158523(n)) = A158523(n) * A008836(n). - Enrique Pérez Herrero, Nov 07 2022
a(n) = (1/n) * Sum_{d|n} mu(n/d)*sigma(d^2). - Ridouane Oudra, Mar 26 2025

More terms from Olivier Gérard, Aug 15 1997

A005361 Product of exponents of prime factorization of n.

Original entry on oeis.org

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

Offset: 1

Jeffrey Shallit and Olivier Gérard

a(n) depends only on prime signature of n (cf. A025487, A052306). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).
There was a comment here that said "a(n) is the number of nilpotents elements in the ring Z/nZ", but this is false, see A003557.
a(n) is the number of square-full divisors of n. a(n) is also the number of divisors d of n such that d and n have the same prime factors, i.e., A007947(d) = A007947(n). - Laszlo Toth, May 22 2009
Number of divisors u of n such that u|(u^n/n). Row lengths in triangle of A284318. - Juri-Stepan Gerasimov, Apr 05 2017

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)
Imanuel Chen and Michael Z. Spivey, Integral Generalized Binomial Coefficients of Multiplicative Functions, Preprint 2015; Summer Research Paper 238, Univ. Puget Sound.
P. Erdős and T. Motzkin, Problem 5735, Amer. Math. Monthly, 78 (1971), 680-681. (Incorrect solution!)
Vaclav Kotesovec, Graph - the asymptotic ratio (100000 terms)
J. Knopfmacher, A prime-divisor function, Proc. Amer. Math. Soc., 40 (1973), 373-377.
MathStackExchange, Dirichlet series for zeta(s)*zeta(2*s)*zeta(3*s)*zeta(6*s)^(-1).
H. N. Shapiro, Problem 5735, Amer. Math. Monthly, 97 (1990), 937.
D. Suryanarayana and R. Sitaramachandra Rao, The number of square-full divisors of an integer, Proc. Amer. Math. Soc., Vol. 34, No. 1 (1972), pp. 79-80.
Index entries for sequences computed from exponents in factorization of n

Cf. A000005, A002110, A005117 (indices of ones), A028234, A052306, A067029, A072411, A082695, A284318.
Cf. A360969, A360970, A361132.
Cf. A340065 (Dgf at s=2).

a005361 = product . a124010_row -- Reinhard Zumkeller, Jan 09 2012

A005361 := proc(n)
    local a, p ;
    a := 1 ;
    for p in ifactors(n)[2] do
       a := a*op(2, p) ;
    end do:
    a ;
end proc:
seq(A005361(n),n=1..30) ; # R. J. Mathar, Nov 20 2012
# second Maple program:
a:= n-> mul(i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Feb 18 2020

Prepend[ Array[ Times @@ Last[ Transpose[ FactorInteger[ # ] ] ]&, 100, 2 ], 1 ]
Array[Times@@Transpose[FactorInteger[#]][[2]]&,80] (* Harvey P. Dale, Aug 15 2012 *)

for(n=1,100, f=factor(n); print1(prod(i=1,omega(f), f[i,2]),",")) \\ edited by M. F. Hasler, Feb 18 2020

a(n)=factorback(factor(n)[,2]) \\ Charles R Greathouse IV, Nov 07 2014

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

from math import prod
from sympy import factorint
def a(n): return prod(factorint(n).values())
print([a(n) for n in range(1, 91)]) # Michael S. Branicky, Jul 04 2022

(define (A005361 n) (if (= 1 n) 1 (* (A067029 n) (A005361 (A028234 n))))) ;; Antti Karttunen, Mar 06 2017

n = Product (p_j^k_j) -> a(n) = Product (k_j).
Dirichlet g.f.: zeta(s)*zeta(2s)*zeta(3s)/zeta(6s).
Multiplicative with a(p^e) = e. - David W. Wilson, Aug 01 2001
a(n) = Sum_{d dividing n} floor(rad(d)/rad(n)) where rad(n) is A007947. - Enrique Pérez Herrero, Nov 06 2009
For n > 1: a(n) = Product_{k=1..A001221(n)} A124010(n,k). - Reinhard Zumkeller, Aug 27 2011
a(n) = tau(n/rad(n)), where tau is A000005 and rad is A007947. - Anthony Browne, May 11 2016
a(n) = Sum_{k=1..n}(floor(cos^2(Pi*k^n/n))*floor(cos^2(Pi*n/k))). - Anthony Browne, May 11 2016
From Antti Karttunen, Mar 06 2017: (Start)
For all n >= 1, a(prime^n) = n, a(A002110(n)) = a(A005117(n)) = 1. [From Crossrefs section.]
a(1) = 1; for n > 1, a(n) = A067029(n) * a(A028234(n)).
(End)
Let (b(n)) be multiplicative with b(p^e) = -1 + ( (floor((e-1)/3)+floor(e/3)) mod 4 ) for p prime and e > 0, then b(n) is the Dirichlet inverse of (a(n)). - Werner Schulte, Feb 23 2018
Sum_{i=1..k} a(i) ~ (zeta(2)*zeta(3)/zeta(6)) * k (Suryanarayana and Sitaramachandra Rao, 1972). - Amiram Eldar, Apr 13 2020
More precise asymptotics: Sum_{k=1..n} a(k) ~ 315*zeta(3)*n / (2*Pi^4) + zeta(1/2)*zeta(3/2)*sqrt(n) / zeta(3) + 6*zeta(1/3)*zeta(2/3)*n^(1/3) / Pi^2 [Knopfmacher, 1973]. - Vaclav Kotesovec, Jun 13 2020

A059956 Decimal expansion of 6/Pi^2.

Original entry on oeis.org

6, 0, 7, 9, 2, 7, 1, 0, 1, 8, 5, 4, 0, 2, 6, 6, 2, 8, 6, 6, 3, 2, 7, 6, 7, 7, 9, 2, 5, 8, 3, 6, 5, 8, 3, 3, 4, 2, 6, 1, 5, 2, 6, 4, 8, 0, 3, 3, 4, 7, 9, 2, 9, 3, 0, 7, 3, 6, 5, 4, 1, 9, 1, 3, 6, 5, 0, 3, 8, 7, 2, 5, 7, 7, 3, 4, 1, 2, 6, 4, 7, 1, 4, 7, 2, 5, 5, 6, 4, 3, 5, 5, 3, 7, 3, 1, 0, 2, 5, 6, 8, 1, 7, 3, 3

Offset: 0

Jason Earls, Mar 01 2001

"6/Pi^2 is the probability that two randomly selected numbers will be coprime and also the probability that a randomly selected integer is 'squarefree.'" [Hardy and Wright] - C. Pickover.
In fact, the probability that any k randomly selected numbers will be coprimes is 1/Sum_{n>=1} n^(-k) = 1/zeta(k). - Robert G. Wilson v [corrected by Ilya Gutkovskiy, Aug 18 2018]
6/Pi^2 is also the diameter of a circle whose circumference equals the ratio of volume of a cuboid to the inscribed ellipsoid. 6/Pi^2 is also the diameter of a circle whose circumference equals the ratio of surface area of a cube to the inscribed sphere. - Omar E. Pol, Oct 08 2011
6/(Pi^2 * n^2) is the probability that two randomly selected positive integers will have a greatest common divisor equal to n, n >= 1. - Geoffrey Critzer, May 28 2013
Equals lim_{n->oo} (Sum_{k=1..n} phi(k)/k)/n, i.e., the limit mean value of phi(k)/k, where phi(k) is Euler's totient function. Proof is trivial using the formula for Sum_{k=1..n} phi(k)/k listed at the Wikipedia link. For the limit mean value of k/phi(k), see A082695. - Stanislav Sykora, Nov 14 2014
This is the probability that a random point on a square lattice is visible from the origin, i.e., there is no other lattice point that lies on the line segment between this point and the origin. - Amiram Eldar, Jul 08 2020

			.6079271018540266286632767792583658334261526480...

Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, vol. 94, Cambridge University Press, Section 1.4, p. 18.
Hardy and Wright, 'An Introduction to the Theory of Numbers'. See Theorems 332 and 333.
C. Pickover, Wonders of Numbers, Oxford University Press, NY, 2001, p. 359.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See p. 184.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, pages 118-119.
David Wells, The Penguin Dictionary of Curious and Interesting Numbers. Penguin Books, NY, 1986, Revised edition 1987, p. 28.

Harry J. Smith, Table of n, a(n) for n = 0..20000
Persi Diaconis and Paul Erdős, On the distribution of the greatest common divisor, in A Festschrift for Herman Rubin, pp. 56-61, IMS Lecture Notes Monogr. Ser., 45, Inst. Math. Statist., Beachwood, OH, 2004.
C. A. Pickover, "Wonders of Numbers, Adventures in Mathematics, Mind and Meaning," Zentralblatt review.
Harry J. Smith, XPCalc. [Broken link]
Eric Weisstein's World of Mathematics, Hafner-Sarnak-McCurley Constant.
Eric Weisstein's World of Mathematics, Relatively Prime.
Eric Weisstein's World of Mathematics, Squarefree.
Wikipedia, Euler's totient function.
Index entries for transcendental numbers.

See A002117 for further references and links.

Cf. A005117 (squarefree numbers), A013661, A082695.

R:= RealField(100); 6/(Pi(R))^2; // G. C. Greubel, Mar 09 2018

evalf(1/Zeta(2)) ; # R. J. Mathar, Mar 27 2013

RealDigits[ 6/Pi^2, 10, 105][[1]]
RealDigits[1/Zeta[2], 10, 111][[1]] (* Robert G. Wilson v, Jan 20 2017 *)

default(realprecision, 20080); x=60/Pi^2; for (n=0, 20000, d=floor(x); x=(x-d)*10; write("b059956.txt", n, " ", d)); \\ Harry J. Smith, Jun 30 2009

Equals 1/A013661.

6/Pi^2 = Product_{k>=1} (1 - 1/prime(k)^2) = Sum_{k>=1} mu(k)/k^2. - Vladeta Jovovic, May 18 2001

A003958 If n = Product p(k)^e(k) then a(n) = Product (p(k)-1)^e(k).

Original entry on oeis.org

1, 1, 2, 1, 4, 2, 6, 1, 4, 4, 10, 2, 12, 6, 8, 1, 16, 4, 18, 4, 12, 10, 22, 2, 16, 12, 8, 6, 28, 8, 30, 1, 20, 16, 24, 4, 36, 18, 24, 4, 40, 12, 42, 10, 16, 22, 46, 2, 36, 16, 32, 12, 52, 8, 40, 6, 36, 28, 58, 8, 60, 30, 24, 1, 48, 20, 66, 16, 44, 24, 70, 4, 72, 36, 32, 18, 60, 24, 78, 4, 16

Offset: 1

Marc LeBrun

Completely multiplicative.

Dirichlet inverse of A097945. - R. J. Mathar, Aug 29 2011

Daniel Forgues, Table of n, a(n) for n = 1..100000 (first 1000 terms from T. D. Noe)
Vaclav Kotesovec, Graph - the asymptotic ratio (10^7 terms).
Index to divisibility sequences

Cf. A003959, A168065, A168066, A027746, A006093, A027748, A124010.

a003958 1 = 1
a003958 n = product $ map (subtract 1) $ a027746_row n
-- Reinhard Zumkeller, Apr 09 2012, Mar 02 2012

a:= n-> mul((i[1]-1)^i[2], i=ifactors(n)[2]):
seq(a(n), n=1..80);  # Alois P. Heinz, Sep 13 2017

DirichletInverse[f_][1] = 1/f[1]; DirichletInverse[f_][n_] := DirichletInverse[f][n] = -1/f[1]*Sum[ f[n/d]*DirichletInverse[f][d], {d, Most[ Divisors[n]]}]; muphi[n_] := MoebiusMu[n]*EulerPhi[n]; Table[ DirichletInverse[ muphi][n], {n, 1, 81}] (* Jean-François Alcover, Dec 12 2011, after R. J. Mathar *)
a[1] = 1; a[n_] := (fi = FactorInteger[n]; Times @@ ((fi[[All, 1]] - 1)^fi[[All, 2]])); Table[a[n], {n, 1, 50}] (* G. C. Greubel, Jun 10 2016 *)

a(n)=if(n<1,0,direuler(p=2,n,1/(1-p*X+X))[n]) /* Ralf Stephan */

from math import prod
from sympy import factorint
def a(n): return prod((p-1)**e for p, e in factorint(n).items())
print([a(n) for n in range(1, 82)]) # Michael S. Branicky, Feb 27 2022

Multiplicative with a(p^e) = (p-1)^e. - David W. Wilson, Aug 01 2001
a(n) = A000010(n) iff n is squarefree (see A005117). - Reinhard Zumkeller, Nov 05 2004
a(n) = abs(A125131(n)). - Tom Edgar, May 26 2014
Sum_{k=1..n} a(k) ~ c * n^2, where c = Pi^4 / (315 * zeta(3)) = 1/(2*A082695) = 0.25725505075419... - Vaclav Kotesovec, Jun 14 2020
Dirichlet g.f.: Product_{p prime} 1 / (1 - p^(1-s) + p^(-s)). - Ilya Gutkovskiy, Feb 27 2022
Dirichlet g.f.: zeta(s-1) * zeta(s) * Product_{primes p} (1 + (p^(1-s) - 2) / (1 - p + p^s)), (with a product that converges for s=2). - Vaclav Kotesovec, Feb 11 2023

Definition reedited (from formula) by Daniel Forgues, Nov 17 2009

A286708 Powerful numbers (A001694) that are not prime powers (A000961).

Original entry on oeis.org

36, 72, 100, 108, 144, 196, 200, 216, 225, 288, 324, 392, 400, 432, 441, 484, 500, 576, 648, 675, 676, 784, 800, 864, 900, 968, 972, 1000, 1089, 1125, 1152, 1156, 1225, 1296, 1323, 1352, 1372, 1444, 1521, 1568, 1600, 1728, 1764, 1800, 1936, 1944, 2000, 2025, 2116, 2304, 2312, 2500, 2592, 2601, 2700, 2704, 2744

Offset: 1

Ilya Gutkovskiy, May 13 2017

nonn

If a prime p divides a(n) then p^2 must also divide a(n) and number of distinct primes dividing a(n) > 1.

Intersection of A001694 and A024619.

			-------------------------------
| n | a(n) | prime            |
|   |      | factorization    |
|------------------------------
| 1 | 36   | {{2, 2}, {3, 2}} |
| 2 | 72   | {{2, 3}, {3, 2}} |
| 3 | 100  | {{2, 2}, {5, 2}} |
| 4 | 108  | {{2, 2}, {3, 3}} |
| 5 | 144  | {{2, 4}, {3, 2}} |
| 6 | 196  | {{2, 2}, {7, 2}} |
| 7 | 200  | {{2, 3}, {5, 2}} |
| 8 | 216  | {{2, 3}, {3, 3}} |
| 9 | 225  | {{3, 2}, {5, 2}} |
-------------------------------
a(n) = p_1^e_1*p_2^e_2*... : {{p_1, e_1}, {p_2, e_2}, ...}.

Michael De Vlieger, Table of n, a(n) for n = 1..10000 (first 5997 terms from Robert Israel)
Eric Weisstein's World of Mathematics, Prime Power.
Eric Weisstein's World of Mathematics, Powerful Number.
Index entries for sequences related to powerful numbers

Cf. A000961, A001221, A001694, A024619, A082695, A136141, A131605.

N:= 10000:
S:= {1}: P:= {1}:
p:= 1:
do
  p:= nextprime(p);
  if p^2 > N then break fi;
  S:= map(s -> (s, seq(s*p^k, k = 2 .. floor(log[p](N/s)))), S);
  P:= P union {seq(p^k, k=2..floor(log[p](N)))}:
od:
sort(convert(S minus P, list)); # Robert Israel, May 14 2017

Select[Range@2750, Min@FactorInteger[#][[All, 2]] > 1 && ! PrimePowerQ[#] &]
(* Second program *)
nn = 2^25; Select[Rest@ Union@ Flatten@ Table[a^2*b^3, {b, nn^(1/3)}, {a, Sqrt[nn/b^3]}], ! PrimePowerQ[#] &] (* Michael De Vlieger, Jun 22 2022 *)

from sympy import primefactors, factorint
print([n for n in range(4,2745) if len(primefactors(n)) > 1 and min(list(factorint(n).values())) > 1]) # Karl-Heinz Hofmann, Feb 07 2023

from math import isqrt
from sympy import integer_nthroot, primepi, mobius
def A286708(n):
    def squarefreepi(n): return int(sum(mobius(k)*(n//k**2) for k in range(1, isqrt(n)+1)))
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x):
        c, l = n+x, 0
        j = isqrt(x)
        while j>1:
            k2 = integer_nthroot(x//j**2,3)[0]+1
            w = squarefreepi(k2-1)
            c -= j*(w-l)
            l, j = w, isqrt(x//k2**3)
        c -= squarefreepi(integer_nthroot(x,3)[0])-l
        return c+1+sum(primepi(integer_nthroot(x, k)[0]) for k in range(2, x.bit_length()))
    return bisection(f,n,n) # Chai Wah Wu, Sep 10 2024

Sum_{n>=1} 1/a(n) = zeta(2)*zeta(3)/zeta(6) - Sum_{p prime} 1/(p*(p-1)) - 1 = A082695 - A136141 - 1 = 0.17043976777096407719... - Amiram Eldar, Feb 12 2021

A014197 Number of numbers m with Euler phi(m) = n.

Original entry on oeis.org

2, 3, 0, 4, 0, 4, 0, 5, 0, 2, 0, 6, 0, 0, 0, 6, 0, 4, 0, 5, 0, 2, 0, 10, 0, 0, 0, 2, 0, 2, 0, 7, 0, 0, 0, 8, 0, 0, 0, 9, 0, 4, 0, 3, 0, 2, 0, 11, 0, 0, 0, 2, 0, 2, 0, 3, 0, 2, 0, 9, 0, 0, 0, 8, 0, 2, 0, 0, 0, 2, 0, 17, 0, 0, 0, 0, 0, 2, 0, 10, 0, 2, 0, 6, 0, 0, 0, 6, 0, 0, 0, 3

Offset: 1

N. J. A. Sloane

Carmichael conjectured that there are no 1's in this sequence. - Jud McCranie, Oct 10 2000
Number of cyclotomic polynomials of degree n. - T. D. Noe, Aug 15 2003
Let v == 0 (mod 24), w = v + 24, and v < k < q < w, where k and q are integer. It seems that, for most values of v, there is no b such that b = a(k) + a(q) and b > a(v) + a(w). The first case where b > a(v) + a(w) occurs at v = 888: b = a(896) + a(900) = 15 + 4, b > a(888) + a(912), or 19 > 8 + 7. The first case where v < n < w and a(n) > a(v) + a(w) occurs at v = 2232: a(2240) > a(2232) + a(2256), or 27 > 7 + 8. - Sergey Pavlov, Feb 05 2017
One elementary result relating to phi(m) is that if m is odd, then phi(m)=phi(2m) because 1 and 2 both have phi value 1 and phi is multiplicative. - Roderick MacPhee, Jun 03 2017

Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section B39, pp. 144-146.
Joe Roberts, Lure of The Integers, The Mathematical Association of America, 1992, entry 32, page 182.

T. D. Noe, Table of n, a(n) for n = 1..10000
Max A. Alekseyev, Computing the Inverses, their Power Sums, and Extrema for Euler's Totient and Other Multiplicative Functions, Journal of Integer Sequences, Vol. 19 (2016), Article 16.5.2.
Max Alekseyev, PARI/GP Scripts for Miscellaneous Math Problems (invphi.gp).
R. D. Carmichael, On Euler’s phi-function, Bull. Amer. Math. Soc. 13 (1907), 241-243.
R. D. Carmichael, Note on Euler’s phi-function, Bull. Amer. Math. Soc. 28 (1922), 109-110.
Kevin Ford, The number of solutions of phi(x)=m, arXiv:math/9907204 [math.NT], 1999.
S. Sivasankaranarayana Pillai, On some functions connected with phi(n), Bull. Amer. Math. Soc. 35 (1929), 832-836.
Primefan, Totient Answers For The First 1000 Integers.
Eric Weisstein's World of Mathematics, Totient Function.
Eric Weisstein's World of Mathematics, Totient Valence Function.

Cf. A000010, A002202, A032446 (bisection), A049283, A051894, A055506, A057635, A057826, A058277 (nonzero terms), A058341, A063439, A066412, A070243 (partial sums), A070633, A071386 (positions of odd terms), A071387, A071388 (positions of primes), A071389 (where prime(n) occurs for the first time), A082695, A097942 (positions of records), A097946, A120963, A134269, A219930, A280611, A280709, A280712, A296655 (positions of positive even terms), A305353, A305656, A319048, A322019.
For records see A131934.
Column 1 of array A320000.

a := function(n)
local S, T, R, max, i, k, r;
S:=[];
for i in DivisorsInt(n)+1 do
    if IsPrime(i)=true then
        S:=Concatenation(S,[i]);
    fi;
od;
T:=[];
for k in [1..Size(S)] do
    T:=Concatenation(T,[S[k]/(S[k]-1)]);
od;
max := n*Product(T);
R:=[];
for r in [1..Int(max)] do
    if Phi(r)=n then
        R:=Concatenation(R,[r]);
    fi;
od;
return Size(R);
end; # Miles Englezou, Oct 22 2024

[#EulerPhiInverse(n): n in [1..100]]; // Marius A. Burtea, Sep 08 2019

with(numtheory): A014197:=n-> nops(invphi(n)): seq(A014197(n), n=1..200);

a[1] = 2; a[m_?OddQ] = 0; a[m_] := Module[{p, nmax, n, k}, p = Select[ Divisors[m]+1, PrimeQ]; nmax = m*Times @@ (p/(p - 1)); n = m; k = 0; While[n <= nmax, If[EulerPhi[n] == m, k++]; n++]; k]; Array[a, 92] (* Jean-François Alcover, Dec 09 2011, updated Apr 25 2016 *)
With[{nn = 116}, Function[s, Function[t, Take[#, nn] &@ ReplacePart[t, Map[# -> Length@ Lookup[s, #] &, Keys@ s]]]@ ConstantArray[0, Max@ Keys@ s]]@ KeySort@ PositionIndex@ Array[EulerPhi, Floor[nn^(3/2)] + 10]] (* Michael De Vlieger, Jul 19 2017 *)

A014197(n,m=1) = { n==1 && return(1+(m<2)); my(p,q); sumdiv(n, d, if( d>=m && isprime(d+1), sum( i=0,valuation(q=n\d,p=d+1), A014197(q\p^i,p))))} \\ M. F. Hasler, Oct 05 2009

a(n) = invphiNum(n); \\ Amiram Eldar, Nov 15 2024 using Max Alekseyev's invphi.gp

from sympy import totient, divisors, isprime, prod
def a(m):
    if m == 1: return 2
    if m % 2: return 0
    X = (x + 1 for x in divisors(m))
    nmax=m*prod(i/(i - 1) for i in X if isprime(i))
    n=m
    k=0
    while n<=nmax:
        if totient(n)==m:k+=1
        n+=1
    return k
print([a(n) for n in range(1, 51)]) # Indranil Ghosh, Jul 18 2017, after Mathematica code

Dirichlet g.f.: Sum_{n>=1} a(n)*n^-s = zeta(s)*Product_(1+1/(p-1)^s-1/p^s). - Benoit Cloitre, Apr 12 2003
Limit_{n->infinity} (1/n) * Sum_{k=1..n} a(k) = zeta(2)*zeta(3)/zeta(6) = 1.94359643682075920505707036... (see A082695). - Benoit Cloitre, Apr 12 2003
From Christopher J. Smyth, Jan 08 2017: (Start)
Euler transform = Product_{n>=1} (1-x^n)^(-a(n)) = g.f. of A120963.
Product_{n>=1} (1+x^n)^a(n)
= Product_{n>=1} ((1-x^(2n))/(1-x^n))^a(n)
= Product_{n>=1} (1-x^n)^(-A280712(n))
= Euler transform of A280712 = g.f. of A280611.
(End)
a(A000010(n)) = A066412(n). - Antti Karttunen, Jul 18 2017
From Antti Karttunen, Dec 04 2018: (Start)
a(A000079(n)) = A058321(n).
a(A000142(n)) = A055506(n).
a(A017545(n)) = A063667(n).
a(n) = Sum_{d|n} A008683(n/d)*A070633(d).
a(n) = A056239(A322310(n)).
(End)

A036689 Product of a prime and the previous number.

Original entry on oeis.org

2, 6, 20, 42, 110, 156, 272, 342, 506, 812, 930, 1332, 1640, 1806, 2162, 2756, 3422, 3660, 4422, 4970, 5256, 6162, 6806, 7832, 9312, 10100, 10506, 11342, 11772, 12656, 16002, 17030, 18632, 19182, 22052, 22650, 24492, 26406, 27722, 29756, 31862, 32580, 36290, 37056, 38612, 39402, 44310

Offset: 1

Felice Russo

Records in A002618. - Artur Jasinski, Jan 23 2008

Also records in A174857. - Vladimir Shevelev, Mar 31 2010

			2*1, 3*2, 5*4, 7*6, 11*10, 13*12, 17*16, ...

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Index to sequences related to prime powers

Cf. A000040, A001248, A002618, A005596, A053650, A053192, A053193, A053650, A082695, A117495, A136141.
Twice the terms of A008837.
Subsequence of A002378 (oblong numbers).
Column 1 of A257251. (Row 1 of A257252.)
Column 2 of A379010.

a036689 n = a036689_list !! (n-1)
a036689_list = zipWith (*) a000040_list $ map pred a000040_list
-- Reinhard Zumkeller, Sep 17 2011

[ n*(n-1): n in PrimesUpTo(220) ]; // Bruno Berselli, Apr 11 2011

A036689 := proc(n) local p ; p := ithprime(n) ; p*(p-1) ; end proc: # R. J. Mathar, Apr 11 2011

Table[Prime[n] EulerPhi[Prime[n]], {n, 100}] (* Artur Jasinski, Jan 23 2008 *)
Table[Prime[n] (Prime[n] - 1), {n, 1, 50}] (* Bruno Berselli, Apr 22 2014 *)
#(#-1)&/@Prime[Range[50]] (* Harvey P. Dale, Sep 08 2019 *)

forprime(p=2,1e3,print1(p^2-p", ")) \\ Charles R Greathouse IV, Jun 10 2011

A036689(n) = ((p->(p-1)*p)(prime(n))); \\ Antti Karttunen, Dec 14 2024

(define (A036689 n) (* (A000040 n) (- (A000040 n) 1))) ;; Antti Karttunen, May 01 2015

a(n) = prime(n) * (prime(n) - 1).
a(n) = phi(prime(n)^2) = A000010(A001248(n)).
a(n) = prime(n) * phi(prime(n)). - Artur Jasinski, Jan 23 2008
From Reinhard Zumkeller, Sep 17 2011: (Start)
a(n) = A000040(n) * A006093(n) = A001248(n) - A000040(n).
A006530(a(n)) = A000040(n). (End)
a(n) = A009262(prime(n)). - Enrique Pérez Herrero, May 12 2012
a(n) = prime(n)! mod (prime(n)^2). - J. M. Bergot, Apr 10 2014
a(n) = 2*A008837(n). - Antti Karttunen, May 01 2015
Sum_{n>=1} 1/a(n) = A136141. - Amiram Eldar, Nov 09 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(2)*zeta(3)/zeta(6) (A082695).
Product_{n>=1} (1 - 1/a(n)) = A005596. (End)

Deleted two incorrect comments. - N. J. A. Sloane, May 07 2020

A064549 a(n) = n * Product_{primes p|n} p.

Original entry on oeis.org

1, 4, 9, 8, 25, 36, 49, 16, 27, 100, 121, 72, 169, 196, 225, 32, 289, 108, 361, 200, 441, 484, 529, 144, 125, 676, 81, 392, 841, 900, 961, 64, 1089, 1156, 1225, 216, 1369, 1444, 1521, 400, 1681, 1764, 1849, 968, 675, 2116, 2209, 288, 343, 500, 2601, 1352

Offset: 1

Henry Bottomley, Oct 16 2001

Index of first occurrence of n in A003557. - Franklin T. Adams-Watters, Jul 25 2014

			a(12) = 72 since 12 = 2^2*3 and 12*2*3 = 72.

Seiichi Manyama, Table of n, a(n) for n = 1..10000 (terms 1..1000 from Harry J. Smith)
Nadia Heninger, E. M. Rains and N. J. A. Sloane, On the Integrality of n-th Roots of Generating Functions, J. Combinatorial Theory, Series A, Vol. 113, No. 8 (2006), pp. 1732-1745; arXiv preprint, arXiv:math/0509316 [math.NT], 2005-2006.

A permutation of the powerful numbers A001694.

Cf. A003557 (a left inverse), A007947, A057521, A078310, A082695, A202535.

a064549 n = a007947 n * n  -- Reinhard Zumkeller, Jul 23 2013

[n^2/( (&+[Floor(k^n/n)-Floor((k^n - 1)/n) : k in [1..n]]) ): n in [1..50]]; // G. C. Greubel, Nov 02 2018

a:= n -> n * convert(numtheory:-factorset(n), `*`):
seq(a(n),n=1..100); # Robert Israel, Jul 25 2014

a[n_] := n * Times @@ FactorInteger[n][[All, 1]]; Array[a, 100] (* Jean-François Alcover, Feb 17 2017 *)
Table[n*Product[If[PrimeQ[d], d, 1], {d, Divisors[n]}], {n, 1, 100}] (* Vaclav Kotesovec, Jun 15 2019 *)

popf(n)= { local(f,p=1); f=factor(n); for(i=1, matsize(f)[1], p*=f[i, 1]); return(p) } { for (n=1, 1000, write("b064549.txt", n, " ", n*popf(n)) ) } \\ Harry J. Smith, Sep 18 2009

A064549(n) = { my(f=factor(n)); for (i=1, #f~, f[i, 2]++); factorback(f); }; \\ Antti Karttunen, Aug 30 2018

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

Multiplicative with a(p^k)=p^(k+1) when k>0.
a(n) = n*A007947(n) = n^2/A003557(n).
Dirichlet convolution of A000027 and A202535. - R. J. Mathar, Dec 20 2011
a(n) = A078310(n) - 1. - Reinhard Zumkeller, Jul 23 2013
A003557(a(n)) = n; a(A003557(n)) = A057521(n). - Antti Karttunen, Aug 30 2018
G.f.: Sum_{k>=1} mu(k)^2*phi(k)*k*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Nov 02 2018
From Vaclav Kotesovec, Jun 24 2020: (Start)
Dirichlet g.f.: zeta(s-2) * zeta(s-1) * Product_{primes p} (1 + p^(3-2*s) - p^(4-2*s) - p^(1-s)).
Sum_{k=1..n} a(k) ~ c * n^3, where c = A065463/3 = A065464*Pi^2/18 = 0.234814...
(End)
Sum_{k>=1} 1/a(k) = zeta(2)*zeta(3)/zeta(6) = A082695. - Vaclav Kotesovec, Sep 19 2020
Sum_{k>=1} (-1)^(k+1)/a(k) = zeta(2)*zeta(3)/(3*zeta(6)) = (1/3) * A082695. - Amiram Eldar, Nov 18 2020

A005597 Decimal expansion of the twin prime constant C_2 = Product_{ p prime >= 3 } (1-1/(p-1)^2).

Original entry on oeis.org

6, 6, 0, 1, 6, 1, 8, 1, 5, 8, 4, 6, 8, 6, 9, 5, 7, 3, 9, 2, 7, 8, 1, 2, 1, 1, 0, 0, 1, 4, 5, 5, 5, 7, 7, 8, 4, 3, 2, 6, 2, 3, 3, 6, 0, 2, 8, 4, 7, 3, 3, 4, 1, 3, 3, 1, 9, 4, 4, 8, 4, 2, 3, 3, 3, 5, 4, 0, 5, 6, 4, 2, 3, 0, 4, 4, 9, 5, 2, 7, 7, 1, 4, 3, 7, 6, 0, 0, 3, 1, 4, 1, 3, 8, 3, 9, 8, 6, 7, 9, 1, 1, 7, 7, 9

Offset: 0

N. J. A. Sloane

C_2 = Product_{ p prime > 2} (p * (p-2) / (p-1)^2) is the 2-tuple case of the Hardy-Littlewood prime k-tuple constant (part of First H-L Conjecture): C_k = Product_{ p prime > k} (p^(k-1) * (p-k) / (p-1)^k).
Although C_2 is commonly called the twin prime constant, it is actually the prime 2-tuple constant (prime pair constant) which is relevant to prime pairs (p, p+2m), m >= 1.
The Hardy-Littlewood asymptotic conjecture for Pi_2m(n), the number of prime pairs (p, p+2m), m >= 1, with p <= n, claims that Pi_2m(n) ~ C_2(2m) * Li_2(n), where Li_2(n) = Integral_{2, n} (dx/log^2(x)) and C_2(2m) = 2 * C_2 * Product_{p prime > 2, p | m} (p-1)/(p-2), which gives: C_2(2) = 2 * C_2 as the prime pair (p, p+2) constant, C_2(4) = 2 * C_2 as the prime pair (p, p+4) constant, C_2(6) = 2* (2/1) * C_2 as the prime pair (p, p+6) constant, C_2(8) = 2 * C_2 as the prime pair (p, p+8) constant, C_2(10) = 2 * (4/3) * C_2 as the prime pair (p, p+10) constant, C_2(12) = 2 * (2/1) * C_2 as the prime pair (p, p+12) constant, C_2(14) = 2 * (6/5) * C_2 as the prime pair (p, p+14) constant, C_2(16) = 2 * C_2 as the prime pair (p, p+16) constant, ... and, for i >= 1, C_2(2^i) = 2 * C_2 as the prime pair (p, p+2^i) constant.
C_2 also occurs as part of other Hardy-Littlewood conjectures related to prime pairs, e.g., the Hardy-Littlewood conjecture concerning the distribution of the Sophie Germain primes (A156874) on primes p such that 2p+1 is also prime.
Another constant related to the twin primes is Viggo Brun's constant B (sometimes also called the twin primes Viggo Brun's constant B_2) A065421, where B_2 = Sum (1/p + 1/q) as (p,q) runs through the twin primes.
Reciprocal of the Selberg-Delange constant A167864. See A167864 for additional comments and references. - Jonathan Sondow, Nov 18 2009
C_2 = Product_{prime p>2} (p-2)p/(p-1)^2 is an analog for primes of Wallis' product 2/Pi = Product_{n=1 to oo} (2n-1)(2n+1)/(2n)^2. - Jonathan Sondow, Nov 18 2009
One can compute a cubic variant, product_{primes >2} (1-1/(p-1)^3) = 0.855392... = (2/3) * 0.6601618...* 1.943596... by multiplying this constant with 2/3 and A082695. - R. J. Mathar, Apr 03 2011
Cohen (1998, p. 7) referred to this number as the "twin prime and Goldbach constant" and noted that, conjecturally, the number of twin prime pairs (p,p+2) with p <= X tends to 2*C_2*X/log(X)^2 as X tends to infinity. - Artur Jasinski, Feb 01 2021

			0.6601618158468695739278121100145557784326233602847334133194484233354056423...

Henri Cohen, Number Theory, Volume II: Analytic and Modern Tools, GTM Vol. 240, Springer, 2007; see pp. 208-209.
Richard Crandall and Carl Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2001; see p. 11.
Steven R. Finch, Mathematical Constants, Encyclopedia of Mathematics and its Applications, Vol. 94, Cambridge University Press, 2003, pp. 84-93, 133.
R. K. Guy, Unsolved Problems in Number Theory, Section A8.
G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, ch. 22.20.
Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 194, 263-264.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Harry J. Smith, Table of n, a(n) for n = 0..1001
Folkmar Bornemann, PRIMES Is in P: Breakthrough for "Everyman", Notices Amer. Math. Soc., Vol. 50, No. 5 (May 2003), p. 549.
Paul S. Bruckman, Problem H-576, Advanced Problems and Solutions, The Fibonacci Quarterly, Vol. 39, No. 4 (2001), p. 379; General IZE, Solution to Problem H-576 by the proposer, ibid., Vol. 40, No. 4 (2002), pp. 383-384.
C. K. Caldwell, The Prime Glossary, twin prime constant.
Henri Cohen, High-precision computation of Hardy-Littlewood constants, (1998).
Henri Cohen, High-precision computation of Hardy-Littlewood constants. [pdf copy, with permission]
Steven R. Finch, Mathematical Constants, Errata and Addenda, arXiv:2001.00578 [math.HO], 2020, Sec. 2.1.
Philippe Flajolet and Ilan Vardi, Zeta Function Expansions of Classical Constants, Unpublished manuscript. 1996.
Daniel A. Goldston, Timothy Ngotiaoco and Julian Ziegler Hunts, The tail of the singular series for the prime pair and Goldbach problems, Functiones et Approximatio Commentarii Mathematici, Vol. 56, No. 1 (2017), pp. 117-141; arXiv preprint, arXiv:1409.2151 [math.NT], 2014.
R. J. Mathar, Hardy-Littlewood constants embedded into infinite products over all positive integers, arXiv:0903.2514 [math.NT], 2009-2011, constant T_1^(2).
B. H. Mayoh, The 2nd Goldbach conjecture revisited, BIT 8 (1968) 128-133 Table 5.
G. Niklasch, Some number theoretical constants: 1000-digit values. [Cached copy]
G. Niklasch, Twin primes constant.
Simon Plouffe, The twin primes constant.
Simon Plouffe, Plouffe's Inverter, The twin primes constant.
Pascal Sebah, Numbers, constants and computation (gives 5000 digits).
Eric Weisstein's World of Mathematics, Twin Primes Constant.
Eric Weisstein's World of Mathematics, Twin Prime Conjecture.
Eric Weisstein's World of Mathematics, k-Tuple Conjecture.
Eric Weisstein's World of Mathematics, Prime Constellation.
John W. Wrench, Jr., Evaluation of Artin's constant and the twin-prime constant, Math. Comp., Vol. 15, No. 76 (1961), pp. 396-398.

Cf. A065645 (continued fraction), A065646 (denominators of convergents to twin prime constant), A065647 (numerators of convergents to twin prime constant), A062270, A062271, A114907, A065418 (C_3), A167864, A000010, A008683.

s[n_] := (1/n)*N[ Sum[ MoebiusMu[d]*2^(n/d), {d, Divisors[n]}], 160]; C2 = (175/256)*Product[ (Zeta[n]*(1 - 2^(-n))*(1 - 3^(-n))*(1 - 5^(-n))*(1 - 7^(-n)))^(-s[n]), {n, 2, 160}]; RealDigits[C2][[1]][[1 ;; 105]] (* Jean-François Alcover, Oct 15 2012, after PARI *)
digits = 105; f[n_] := -2*(2^n-1)/(n+1); C2 = Exp[NSum[f[n]*(PrimeZetaP[n+1] - 1/2^(n+1)), {n, 1, Infinity}, NSumTerms -> 5 digits, WorkingPrecision -> 5 digits]]; RealDigits[C2, 10, digits][[1]] (* Jean-François Alcover, Apr 16 2016, updated Apr 24 2018 *)

\p1000; 175/256*prod(k=2,500,(zeta(k)*(1-1/2^k)*(1-1/3^k)*(1-1/5^k)*(1-1/7^k))^(-sumdiv(k,d,moebius(d)*2^(k/d))/k))

prodeulerrat(1-1/(p-1)^2, 1, 3) \\ Amiram Eldar, Mar 12 2021

Equals Product_{k>=2} (zeta(k)*(1-1/2^k))^(-Sum_{d|k} mu(d)*2^(k/d)/k). - Benoit Cloitre, Aug 06 2003
Equals 1/A167864. - Jonathan Sondow, Nov 18 2009
Equals Sum_{k>=1} mu(2*k-1)/phi(2*k-1)^2, where mu is the Möbius function (A008683) and phi is the Euler totient function (A000010) (Bruckman, 2001). - Amiram Eldar, Jan 14 2022

More terms from Vladeta Jovovic, Nov 08 2001
Commented and edited by Daniel Forgues, Jul 28 2009, Aug 04 2009, Aug 12 2009
PARI code removed by D. S. McNeil, Dec 26 2010

cp's OEIS Frontend

A003557 n divided by largest squarefree divisor of n; if n = Product p(k)^e(k) then a(n) = Product p(k)^(e(k)-1), with a(1) = 1.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Julia

Magma

Maple

Mathematica

PARI

PARI

Python

Python

Sage

Formula

Extensions

A001615 Dedekind psi function: n * Product_{p|n, p prime} (1 + 1/p).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Sage

Formula

Extensions

A005361 Product of exponents of prime factorization of n.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

Python

Scheme

Formula

A059956 Decimal expansion of 6/Pi^2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Magma

Maple

Mathematica

PARI

Formula

A003958 If n = Product p(k)^e(k) then a(n) = Product (p(k)-1)^e(k).

Views

Author

Keywords