A051709 -id:A051709 - cp's OEIS frontend

A001248 Squares of primes.

4, 9, 25, 49, 121, 169, 289, 361, 529, 841, 961, 1369, 1681, 1849, 2209, 2809, 3481, 3721, 4489, 5041, 5329, 6241, 6889, 7921, 9409, 10201, 10609, 11449, 11881, 12769, 16129, 17161, 18769, 19321, 22201, 22801, 24649, 26569, 27889, 29929, 32041, 32761, 36481

Offset: 1

N. J. A. Sloane

Also 4, together with numbers n such that Sum_{d|n}(-1)^d = -A048272(n) = -3. - Benoit Cloitre, Apr 14 2002
Also, all solutions to the equation sigma(x) + phi(x) = 2x + 1. - Farideh Firoozbakht, Feb 02 2005
Unique numbers having 3 divisors (1, their square root, themselves). - Alexandre Wajnberg, Jan 15 2006
Smallest (or first) new number deleted at the n-th step in an Eratosthenes sieve. - Lekraj Beedassy, Aug 17 2006
Subsequence of semiprimes A001358. - Lekraj Beedassy, Sep 06 2006
Integers having only 1 factor other than 1 and the number itself. Every number in the sequence is a multiple of 1 factor other than 1 and the number itself. 4 : 2 is the only factor other than 1 and 4; 9 : 3 is the only factor other than 1 and 9; and so on. - Rachit Agrawal (rachit_agrawal(AT)daiict.ac.in), Oct 23 2007
The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. - Omar E. Pol, May 06 2008
There are 2 Abelian groups of order p^2 (C_p^2 and C_p x C_p) and no non-Abelian group. - Franz Vrabec, Sep 11 2008
Also numbers n such that phi(n) = n - sqrt(n). - Michel Lagneau, May 25 2012
For n > 1, n is the sum of numbers from A006254(n-1) to A168565(n-1). - Vicente Izquierdo Gomez, Dec 01 2012
A078898(a(n)) = 2. - Reinhard Zumkeller, Apr 06 2015
Let r(n) = (a(n) - 1)/(a(n) + 1); then Product_{n>=1} r(n) = (3/5) * (4/5) * (12/13) * (24/25) * (60/61) * ... = 2/5. - Dimitris Valianatos, Feb 26 2019
Numbers k such that A051709(k) = 1. - Jianing Song, Jun 27 2021

Ray Chandler, Table of n, a(n) for n = 1..10000 (first 5000 terms from N. J. A. Sloane)
Nicholas John Bizzell-Browning, LIE scales: Composing with scales of linear intervallic expansion, Ph. D. Thesis, Brunel Univ. (UK, 2024). See p. 143.
R. P. Boas and N. J. A. Sloane, Correspondence, 1974
Marius Coman, On the special relation between the numbers of the form 505+ 1008k and the squares of primes, 2015.
Brady Haran and Matt Parker, Squaring Primes, Numberphile video (2018).
Eric Weisstein's World of Mathematics, Prime Power.
OEIS Wiki, Index entries for number of divisors
Index to sequences related to prime signature

Cf. A000005, A000040, A001358, A001694, A002033, A005408, A005722, A006254.
Cf. A008864, A033273, A049001, A048272, A051709, A059956, A060800, A062799.
Cf. A078898, A082020, A085548, A087112, A095874, A168565, A192134, A211110.
Cf. A258599.
Subsequence of A000430, A001358, and of A190641.
Cf. A024450 (partial sums), A069482 (first differences).

a001248 n = a001248_list !! (n-1)
a001248_list = map (^ 2) a000040_list -- Reinhard Zumkeller, Sep 23 2011

[p^2: p in PrimesUpTo(300)]; // Vincenzo Librandi, Mar 27 2014

A001248:=n->ithprime(n)^2; seq(A001248(k), k=1..50); # Wesley Ivan Hurt, Oct 11 2013

Prime[Range[30]]^2 (* Zak Seidov, Dec 07 2011 *)
Select[Range[40000], DivisorSigma[0, #] == 3 &] (* Carlos Eduardo Olivieri, Jun 01 2015 *)

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

A001248(n)=prime(n)^2  \\ M. F. Hasler, Sep 16 2012

from sympy import prime
def A001248(n): return prime(n)**2 # Chai Wah Wu, Aug 09 2024

[n^2 for n in prime_range(1,301)] # G. C. Greubel, May 02 2024

n such that A062799(n) = 2. - Benoit Cloitre, Apr 06 2002
A000005(a(n)^(k-1)) = A005408(k) for all k>0. - Reinhard Zumkeller, Mar 04 2007
a(n) = A000040(n)^(3-1)=A000040(n)^2, where 3 is the number of divisors of a(n). - Omar E. Pol, May 06 2008
A000005(a(n)) = 3 or A002033(a(n)) = 2. - Juri-Stepan Gerasimov, Oct 10 2009
A033273(a(n)) = 3. - Juri-Stepan Gerasimov, Dec 07 2009
For n > 2: (a(n) + 17) mod 12 = 6. - Reinhard Zumkeller, May 12 2010
A192134(A095874(a(n))) = A005722(n) + 1. - Reinhard Zumkeller, Jun 26 2011
For n > 2: a(n) = 1 (mod 24). - Zak Seidov, Dec 07 2011
A211110(a(n)) = 2. - Reinhard Zumkeller, Apr 02 2012
a(n) = A087112(n,n). - Reinhard Zumkeller, Nov 25 2012
a(n) = prime(n)^2. - Jon E. Schoenfield, Mar 29 2015
Product_{n>=1} a(n)/(a(n)-1) = Pi^2/6. - Daniel Suteu, Feb 06 2017
Sum_{n>=1} 1/a(n) = P(2) = 0.4522474200... (A085548). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(2)/zeta(4) = 15/Pi^2 (A082020).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(2) = 6/Pi^2 (A059956). (End)

A006881 Squarefree semiprimes: Numbers that are the product of two distinct primes.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 26, 33, 34, 35, 38, 39, 46, 51, 55, 57, 58, 62, 65, 69, 74, 77, 82, 85, 86, 87, 91, 93, 94, 95, 106, 111, 115, 118, 119, 122, 123, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166, 177, 178, 183, 185, 187, 194, 201, 202, 203, 205

Offset: 1

N. J. A. Sloane, Robert Munafo, Simon Plouffe

Numbers k such that phi(k) + sigma(k) = 2*(k+1). - Benoit Cloitre, Mar 02 2002
Numbers k such that tau(k) = omega(k)^omega(k). - Benoit Cloitre, Sep 10 2002 [This comment is false. If k = 900 then tau(k) = omega(k)^omega(k) = 27 but 900 = (2*3*5)^2 is not the product of two distinct primes. - Peter Luschny, Jul 12 2023]
Could also be called 2-almost primes. - Rick L. Shepherd, May 11 2003
From the Goldston et al. reference's abstract: "lim inf [as n approaches infinity] [(a(n+1) - a(n))] <= 26. If an appropriate generalization of the Elliott-Halberstam Conjecture is true, then the above bound can be improved to 6." - Jonathan Vos Post, Jun 20 2005
The maximal number of consecutive integers in this sequence is 3 - there cannot be 4 consecutive integers because one of them would be divisible by 4 and therefore would not be product of distinct primes. There are several examples of 3 consecutive integers in this sequence. The first one is 33 = 3 * 11, 34 = 2 * 17, 35 = 5 * 7; (see A039833). - Matias Saucedo (solomatias(AT)yahoo.com.ar), Mar 15 2008
Number of terms less than or equal to 10^k for k >= 0 is A036351(k). - Robert G. Wilson v, Jun 26 2012
Are these the numbers k whose difference between the sum of proper divisors of k and the arithmetic derivative of k is equal to 1? - Omar E. Pol, Dec 19 2012
Intersection of A001358 and A030513. - Wesley Ivan Hurt, Sep 09 2013
A237114(n) (smallest semiprime k^prime(n)+1) is a term, for n != 2. - Jonathan Sondow, Feb 06 2014
a(n) are the reduced denominators of p_2/p_1 + p_4/p_3, where p_1 != p_2, p_3 != p_4, p_1 != p_3, and the p's are primes. In other words, (p_2*p_3 + p_1*p_4) never shares a common factor with p_1*p_3. - Richard R. Forberg, Mar 04 2015
Conjecture: The sums of two elements of a(n) forms a set that includes all primes greater than or equal to 29 and all integers greater than or equal to 83 (and many below 83). - Richard R. Forberg, Mar 04 2015
The (disjoint) union of this sequence and A001248 is A001358. - Jason Kimberley, Nov 12 2015
A263990 lists the subsequence of a(n) where a(n+1)=1+a(n). - R. J. Mathar, Aug 13 2019

N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
Zervos, Marie: Sur une classe de nombres composés. Actes du Congrès interbalkanique de mathématiciens 267-268 (1935)

T. D. Noe, Table of n, a(n) for n = 1..10000
D. A. Goldston, S. W. Graham, J. Pimtz and Y. Yildirim, "Small Gaps Between Primes or Almost Primes", arXiv:math/0506067 [math.NT], March 2005.
G. T. Leavens and M. Vermeulen, 3x+1 search programs, Computers and Mathematics with Applications, 24 (1992), 79-99. (Annotated scanned copy)
R. J. Mathar, Series of reciprocal powers of k-almost primes arXiv:0803.0900, table 6 k=2 shows sum 1/a(n)^s.
Eric Weisstein's World of Mathematics, Semiprime
Index to sequences related to prime signature

Products of exactly k distinct primes, for k = 1 to 6: A000040, A006881. A007304, A046386, A046387, A067885.
Cf. A030229, A051709, A001221 (omega(n)), A001222 (bigomega(n)), A001358 (semiprimes), A005117 (squarefree), A007304 (squarefree 3-almost primes), A213952, A039833, A016105 (subsequences), A237114 (subsequence, n != 2).
Subsequence of A007422.
Cf. A259758 (subsequence), A036351, A363923.

a006881 n = a006881_list !! (n-1)
a006881_list = filter chi [1..] where
   chi n = p /= q && a010051 q == 1 where
      p = a020639 n
      q = n `div` p
-- Reinhard Zumkeller, Aug 07 2011

[n: n in [1..210] | EulerPhi(n) + DivisorSigma(1,n) eq 2*(n+1)]; // Vincenzo Librandi, Sep 17 2015

N:= 1001: # to get all terms < N
Primes:= select(isprime, [2,seq(2*k+1,k=1..floor(N/2))]):
{seq(seq(p*q,q=Primes[1..ListTools:-BinaryPlace(Primes,N/p)]),p=Primes)} minus {seq(p^2,p=Primes)};
# Robert Israel, Jul 23 2014
# Alternative, using A001221:
isA006881 := proc(n)
     if numtheory[bigomega](n) =2 and A001221(n) = 2 then
        true ;
    else
        false ;
    end if;
end proc:
A006881 := proc(n) if n = 1 then 6; else for a from procname(n-1)+1 do if isA006881(a) then return a; end if; end do: end if;
end proc: # R. J. Mathar, May 02 2010
# Alternative:
with(NumberTheory): isA006881 := n -> is(NumberOfPrimeFactors(n, 'distinct') = 2 and NumberOfPrimeFactors(n) = 2):
select(isA006881, [seq(1..205)]); # Peter Luschny, Jul 12 2023

mx = 205; Sort@ Flatten@ Table[ Prime[n]*Prime[m], {n, Log[2, mx/3]}, {m, n + 1, PrimePi[ mx/Prime[n]]}] (* Robert G. Wilson v, Dec 28 2005, modified Jul 23 2014 *)
sqFrSemiPrimeQ[n_] := Last@# & /@ FactorInteger@ n == {1, 1}; Select[Range[210], sqFrSemiPrimeQ] (* Robert G. Wilson v, Feb 07 2012 *)
With[{upto=250},Select[Sort[Times@@@Subsets[Prime[Range[upto/2]],{2}]],#<=upto&]] (* Harvey P. Dale, Apr 30 2018 *)

for(n=1,214,if(bigomega(n)==2&&omega(n)==2,print1(n,",")))

for(n=1,214,if(bigomega(n)==2&&issquarefree(n),print1(n,",")))

list(lim)=my(v=List()); forprime(p=2,sqrt(lim), forprime(q=p+1, lim\p, listput(v,p*q))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from sympy import factorint
def ok(n): f=factorint(n); return len(f) == 2 and sum(f[p] for p in f) == 2
print(list(filter(ok, range(1, 206)))) # Michael S. Branicky, Jun 10 2021

from math import isqrt
from sympy import primepi, primerange
def A006881(n):
    def f(x): return int(n+x+(t:=primepi(s:=isqrt(x)))+(t*(t-1)>>1)-sum(primepi(x//k) for k in primerange(1, s+1)))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m # Chai Wah Wu, Aug 15 2024

def A006881_list(n) :
    R = []
    for i in (6..n) :
        d = prime_divisors(i)
        if len(d) == 2 :
            if d[0]*d[1] == i :
                R.append(i)
    return R
A006881_list(205)  # Peter Luschny, Feb 07 2012

A000005(a(n)^(k-1)) = A000290(k) for all k>0. - Reinhard Zumkeller, Mar 04 2007
A109810(a(n)) = 4; A178254(a(n)) = 6. - Reinhard Zumkeller, May 24 2010
A056595(a(n)) = 3. - Reinhard Zumkeller, Aug 15 2011
a(n) = A096916(n) * A070647(n). - Reinhard Zumkeller, Sep 23 2011
A211110(a(n)) = 3. - Reinhard Zumkeller, Apr 02 2012
Sum_{n >= 1} 1/a(n)^s = (1/2)*(P(s)^2 - P(2*s)), where P is Prime Zeta. - Enrique Pérez Herrero, Jun 24 2012
A050326(a(n)) = 2. - Reinhard Zumkeller, May 03 2013
sopf(a(n)) = a(n) - phi(a(n)) + 1 = sigma(a(n)) - a(n) - 1. - Wesley Ivan Hurt, May 18 2013
d(a(n)) = 4. Omega(a(n)) = 2. omega(a(n)) = 2. mu(a(n)) = 1. - Wesley Ivan Hurt, Jun 28 2013
a(n) ~ n log n/log log n. - Charles R Greathouse IV, Aug 22 2013
A089233(a(n)) = 1. - Reinhard Zumkeller, Sep 04 2013
From Peter Luschny, Jul 12 2023: (Start)
For k > 1: k is a term <=> k^A001221(k) = k*A007947(k).
For k > 1: k is a term <=> k^A001222(k) = k*A007947(k).
For k > 1: k is a term <=> A363923(k) = k. (End)
a(n) ~ n log n/log log n. - Charles R Greathouse IV, Jan 13 2025

Name expanded (based on a comment of Rick L. Shepherd) by Charles R Greathouse IV, Sep 16 2015

A065387 a(n) = sigma(n) + phi(n).

Original entry on oeis.org

2, 4, 6, 9, 10, 14, 14, 19, 19, 22, 22, 32, 26, 30, 32, 39, 34, 45, 38, 50, 44, 46, 46, 68, 51, 54, 58, 68, 58, 80, 62, 79, 68, 70, 72, 103, 74, 78, 80, 106, 82, 108, 86, 104, 102, 94, 94, 140, 99, 113, 104, 122, 106, 138, 112, 144, 116, 118, 118, 184, 122, 126, 140

Offset: 1

Labos Elemer, Nov 05 2001

a(n) = 2n for n listed in A008578, the prime numbers at the beginning of the 20th century. When a(n) = a(n + 1), n is probably listed in A066198, numbers n where phi changes as fast as sigma (the only exceptions below 10000 are 2 and 854). - Alonso del Arte, Nov 16 2005
A. Makowski proved that n is prime if and only if a(n) = n * d(n), where d is A000005. - Charles R Greathouse IV, Mar 19 2012
If n is semiprime, a(n) = 2n+1+ceiling(sqrt(n))-floor(sqrt(n)). - Wesley Ivan Hurt, May 05 2015
Atanassov proves that a(n) >= n + A001414(n). - Charles R Greathouse IV, Dec 06 2016
a(n) = 2*n+1 iff n is square of prime (A001248), a(n) = 2*(n+1) iff n is squarefree semiprime (A006881). - Bernard Schott, Feb 09 2020

			a(10) = 22 because there are 4 coprimes to 10 below 10, the divisors of 10 add up to 18, and 4 + 18 = 22.

K. Atanassov, New integer functions, related to ψ and σ functions. IV., Bull. Number Theory Related Topics 12 (1988), pp. 31-35.
Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004. See Section B41, p. 149.
James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 162.

N. J. A. Sloane, Table of n, a(n) for n = 1..10000 (First 1000 terms from T. D. Noe.)
A. Makowski, Aufgaben 339, Elemente der Mathematik 15 (1960), pp. 39-40.

Cf. A000010, A000203, A065388, A015702, A051709, A011774.

See A292768 for partial sums, A051612 for sigma - phi.

[DivisorSigma(1,k)+EulerPhi(k):k in [1..65]]; // Marius A. Burtea, Feb 09 2020

with(numtheory); A065387:=n->phi(n) + sigma(n); seq(A065387(n), n=1..100); # Wesley Ivan Hurt, Apr 08 2014

Table[EulerPhi[n] + DivisorSigma[1,n], {n, 65}] (* Alonso del Arte *)
a[n_] := SeriesCoefficient[Sum[(1+MoebiusMu[k])*x^k/(1-x^k)^2, {k, 1, n}], {x, 0, n}]; Array[a, 63] (* Jean-François Alcover, Sep 29 2017, after Ilya Gutkovskiy *)

a(n) = sigma(n) + eulerphi(n) \\ Harry J. Smith, Oct 17 2009

[sigma(n,1)+euler_phi(n) for n in range(1, 64)] # Stefano Spezia, Jul 20 2025

a(n) = A000203(n) + A000010(n).
a(n) = A051709(n) + 2n. - N. J. A. Sloane, Jun 12 2004
G.f.: Sum_{k>=1} (mu(k) + 1)*x^k/(1 - x^k)^2. - Ilya Gutkovskiy, Sep 29 2017

A344753 a(n) = sigma(n) + psi(n) - 2n = Sum_{d|n, d

Original entry on oeis.org

0, 2, 2, 5, 2, 12, 2, 11, 7, 16, 2, 28, 2, 20, 18, 23, 2, 39, 2, 38, 22, 28, 2, 60, 11, 32, 22, 48, 2, 84, 2, 47, 30, 40, 26, 91, 2, 44, 34, 82, 2, 108, 2, 68, 60, 52, 2, 124, 15, 83, 42, 78, 2, 120, 34, 104, 46, 64, 2, 192, 2, 68, 74, 95, 38, 156, 2, 98, 54, 148, 2, 195, 2, 80, 94, 108, 38, 180, 2, 170, 67, 88, 2

Offset: 1

Antti Karttunen, May 28 2021

Sigma is the sum of divisors (A000203), and psi is Dedekind psi-function (A001615). Coincides with the latter only on perfect numbers (A000396).

Antti Karttunen, Table of n, a(n) for n = 1..10080
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537

Cf. A001065, A001615, A008683, A008966, A033879, A173557, A306927, A344705, A342001, A344705, A344754, A344755, A344997, A344998.
Cf. A000396 [where a(n) = A001615(n)], A005100 [a(n) < A001615(n)], A005101 [a(n) > A001615(n)].
Cf. also A051709, A345001.
Cf. A013661, A082020.

a[n_] := Sum[d + If[SquareFreeQ[n/d], d, 0], {d, Most[Divisors[n]]}];
Array[a, 100] (* Jean-François Alcover, Jun 12 2021 *)

PARI
```
A344753(n) = sumdiv(n,d,(d
				
```

a(n) = Sum_{d|n, dA008966(n/d) * d).
a(n) = A001065(n) + A306927(n).
a(n) = A001615(n) - A033879(n).
a(n) = A344705(n) + 2*A001065(n) - n.
For squarefree n, a(n) = 2*A001065(n).
a(n) = A344997(n) / A173557(n) = A344998(n) / A342001(n). - Antti Karttunen, Jun 06 2021
Sum_{k=1..n} a(k) = c * n^2 / 2 + O(n*log(n)), where c = Pi^2/6 + 15/Pi^2 - 2 = 1.164751... . - Amiram Eldar, Dec 08 2023

New primary definition added by Antti Karttunen, Jun 06 2021

A345001 a(n) = sigma(n) + n' - 2n, where n' is the arithmetic derivative of n (A003415) and sigma is the sum of divisors (A000203).

Original entry on oeis.org

-1, 0, -1, 3, -3, 5, -5, 11, 1, 5, -9, 20, -11, 5, 2, 31, -15, 24, -17, 26, 0, 5, -21, 56, -9, 5, 13, 32, -27, 43, -29, 79, -4, 5, -10, 79, -35, 5, -6, 78, -39, 53, -41, 44, 27, 5, -45, 140, -27, 38, -10, 50, -51, 93, -22, 100, -12, 5, -57, 140, -59, 5, 29, 191, -28, 73, -65, 62, -16, 63, -69, 207, -71, 5, 29, 68

Offset: 1

Antti Karttunen, Jun 05 2021

Coincides with A003415 only on perfect numbers (A000396).

Cf. A000203, A000396, A001065, A003415, A033879, A168036, A344999, A345002, A345003, A345004, A345005, A345049.
Cf. also A051709, A344753 (analogous sequences).
Cf. A013661, A136141.

A003415[n_] := If[n < 2, 0, Module[{f = FactorInteger[n]}, If[PrimeQ[n], 1, Total[n*f[[All, 2]]/f[[All, 1]]]]]];
a[n_] := DivisorSigma[1, n] + A003415[n] - 2 n;
Array[a, 80] (* Jean-François Alcover, Jun 12 2021 *)

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A345001(n) = (sigma(n)+A003415(n)-(2*n));

a(n) = A003415(n) - A033879(n) = A000203(n) + A003415(n) - 2*n.
a(n) = A001065(n) + A168036(n).
a(n) = A344999(n) / A048250(n) = A345049(n) / A173557(n).
Sum_{k=1..n} a(k) ~ c * n^2 / 2, where c = A013661 + A136141 - 2 = 0.418090735898... . - Amiram Eldar, Dec 08 2023

cp's OEIS Frontend

A345048 a(n) = A342001(n) * A051709(n).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A344994 Numbers k such that A173557(k) divides nonzero A051709(k).

Views

Author

Keywords

Crossrefs

Programs

PARI

A344996 a(n) = A048250(n) * A051709(n).

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

Formula

A344995 Numbers k such that A051709(k)/A173557(k) is a positive natural number and a divisor of k.

Views

Author

Keywords

Links

Crossrefs

Programs

PARI

A345054 Odd numbers k such that A173557(k) divides nonzero A051709(k).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A001248 Squares of primes.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

Python

SageMath

Formula

A006881 Squarefree semiprimes: Numbers that are the product of two distinct primes.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Python

Sage

Formula