A056595 -id:A056595 - cp's OEIS frontend

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

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

A030078 Cubes of primes.

Original entry on oeis.org

8, 27, 125, 343, 1331, 2197, 4913, 6859, 12167, 24389, 29791, 50653, 68921, 79507, 103823, 148877, 205379, 226981, 300763, 357911, 389017, 493039, 571787, 704969, 912673, 1030301, 1092727, 1225043, 1295029, 1442897, 2048383, 2248091, 2571353, 2685619, 3307949

Offset: 1

Patrick De Geest

Numbers with exactly three factorizations: A001055(a(n)) = 3 (e.g., a(4) = 1*343 = 7*49 = 7*7*7). - Reinhard Zumkeller, Dec 29 2001
Intersection of A014612 and A000578. Intersection of A014612 and A030513. - Wesley Ivan Hurt, Sep 10 2013
Let r(n) = (a(n)-1)/(a(n)+1) if a(n) mod 4 = 1, (a(n)+1)/(a(n)-1) otherwise; then Product_{n>=1} r(n) = (9/7) * (28/26) * (124/126) * (344/342) * (1332/1330) * ... = 48/35. - Dimitris Valianatos, Mar 06 2020
There exist 5 groups of order p^3, when p prime, so this is a subsequence of A054397. Three of them are abelian: C_p^3, C_p^2 X C_p and C_p X C_p X C_p = (C_p)^3. For 8 = 2^3, the 2 nonabelian groups are D_8 and Q_8; for odd prime p, the 2 nonabelian groups are (C_p x C_p) : C_p, and C_p^2 : C_p (remark, for p = 2, these two semi-direct products are isomorphic to D_8). Here C, D, Q mean Cyclic, Dihedral, Quaternion groups of the stated order; the symbols X and : mean direct and semidirect products respectively. - Bernard Schott, Dec 11 2021

			a(3) = 125; since the 3rd prime is 5, a(3) = 5^3 = 125.

Edmund Landau, Elementary Number Theory, translation by Jacob E. Goodman of Elementare Zahlentheorie (Vol. I_1 (1927) of Vorlesungen über Zahlentheorie), by Edmund Landau, with added exercises by Paul T. Bateman and E. E. Kohlbecker, Chelsea Publishing Co., New York, 1958, pp. 31-32.

T. D. Noe, Table of n, a(n) for n = 1..1000
Xavier Gourdon and Pascal Sebah, Some Constants from Number theory.
Eric Weisstein's World of Mathematics, Prime Power.
Wikipedia, p-group, Classification.
Index to sequences related to prime signature

Other sequences that are k-th powers of primes are: A000040 (k=1), A001248 (k=2), this sequence (k=3), A030514 (k=4), A050997 (k=5), A030516 (k=6), A092759 (k=7), A179645 (k=8), A179665 (k=9), A030629 (k=10), A079395 (k=11), A030631 (k=12), A138031 (k=13), A030635 (k=16), A138032 (k=17), A030637 (k=18).
Cf. A060800, A131991, A000578, subsequence of A046099.
Subsequence of A007422 and of A054397.
Cf. A036966, A085541, A088453, A157289, A258600.

a030078 = a000578 . a000040
a030078_list = map a000578 a000040_list -- Reinhard Zumkeller, May 26 2012

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

Array[Prime[ # ]^3&, 5! ] (* Vladimir Joseph Stephan Orlovsky, Sep 01 2008 *)

a(n)=prime(n)^3 \\ Charles R Greathouse IV, Mar 20 2013

from sympy import prime, primerange
def aupton(terms): return [p**3 for p in primerange(1, prime(terms)+1)]
print(aupton(35)) # Michael S. Branicky, Aug 27 2021

[p**3 for p in prime_range(100)] # Zerinvary Lajos, May 15 2007

n such that A062799(n) = 3. - Benoit Cloitre, Apr 06 2002
a(n) = A000040(n)^3. - Omar E. Pol, Jul 27 2009
A064380(a(n)) = A000010(a(n)). - Vladimir Shevelev, Apr 19 2010
A003415(a(n)) = A079705(n). - Reinhard Zumkeller, Jun 26 2011
A056595(a(n)) = 2. - Reinhard Zumkeller, Aug 15 2011
A000005(a(n)) = 4. - Wesley Ivan Hurt, Sep 10 2013
a(n) = A119959(n) * A008864(n) -1.- R. J. Mathar, Aug 13 2019
Sum_{n>=1} 1/a(n) = P(3) = 0.1747626392... (A085541). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(3)/zeta(6) (A157289).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(3) (A088453). (End)

A046951 a(n) is the number of squares dividing n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 3, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 3, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 3, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 4, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 2, 2, 1, 1, 1, 3, 3, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 3, 1, 2, 2, 4, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2

Offset: 1

Simon Colton (simonco(AT)cs.york.ac.uk)

Rediscovered by the HR automatic theory formation program.
a(n) depends only on prime signature of n (cf. A025487, A046523). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3, 1).
First differences of A013936. Average value tends towards Pi^2/6 = 1.644934... (A013661, A013679). - Henry Bottomley, Aug 16 2001
We have a(n) = A159631(n) for all n < 125, but a(125) = 2 < 3 = A159631(125). - Steven Finch, Apr 22 2009
Number of 2-generated Abelian groups of order n, if n > 1. - Álvar Ibeas, Dec 22 2014 [In other words, number of order-n abelian groups with rank <= 2. Proof: let b(n) be such number. A finite abelian group is the inner direct product of all Sylow-p subgroups, so {b(n)} is multiplicative. Obviously b(p^e) = floor(e/2)+1 (corresponding to the groups C_(p^r) X C_(p^(e-r)) for 0 <= r <= floor(e/2)), hence b(n) = a(n) for all n. - Jianing Song, Nov 05 2022]
Number of ways of writing n = r*s such that r|s. - Eric M. Schmidt, Jan 08 2015
The number of divisors of the square root of the largest square dividing n. - Amiram Eldar, Jul 07 2020
The number of unordered factorizations of n into cubefree powers of primes (1, primes and squares of primes, A166684). - Amiram Eldar, Jun 12 2025

			a(16) = 3 because the squares 1, 4, and 16 divide 16.
G.f. = x + x^2 + x^3 + 2*x^4 + x^5 + x^6 + x^7 + 2*x^8 + 2*x^9 + x^10 + ...

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Antonio Amariti, Claudius Klare, Domenico Orlando and Susanne Reffert, The M-theory origin of global properties of gauge theories, Nuclear Physics B, Vol. 901 (2015), pp. 318-337, arXiv preprint, arXiv:1507.04743 [hep-th], 2015 (see (A.13)).
Simon Colton, Refactorable Numbers - A Machine Invention, J. Integer Sequences, Vol. 2 (1999), Article 99.1.2.
Simon Colton, HR - Automatic Theory Formation in Pure Mathematics
Ian G. Connell, A number theory problem concerning finite groups and rings, Canad. Math. Bull, 7 (1964), 23-34. See delta(n).
Andrew V. Lelechenko, Average number of squares dividing mn, arXiv preprint arXiv:1407.1222 [math.NT], 2014.
Werner Georg Nowak and László Tóth, On the average number of subgroups of the group Z_m X Z_n, International Journal of Number Theory, Vol. 10, No. 2 (2014), pp. 363-374, arXiv preprint, arXiv:1307.1414 [math.NT], 2013.
N. J. A. Sloane, Transforms.
Index entries for sequences computed from exponents in factorization of n.

Cf. A000005, A000188, A004101, A005117 (positions of 1's), A008619, A008833, A013936 (partial sums), A038538, A046952, A052304, A056595, A159631, A007814, A010052, A027750, A239930, A007862, A046523, A064989, A065704, A130279, A156552, A166684, A278161, A307445 (Dirichlet inverse).
One more than A071325.
Differs from A096309 for the first time at n=32, where a(32) = 3, while A096309(32) = 2 (and also A185102(32) = 2).
Sum of the k-th powers of the square divisors of n for k=0..10: this sequence (k=0), A035316 (k=1), A351307 (k=2), A351308 (k=3), A351309 (k=4), A351310 (k=5), A351311 (k=6), A351313 (k=7), A351314 (k=8), A351315 (k=9), A351315 (k=10).
Sequences of the form n^k * Sum_{d^2|n} 1/d^k for k = 0..10: this sequence (k=0), A340774 (k=1), A351600 (k=2), A351601 (k=3), A351602 (k=4), A351603 (k=5), A351604 (k=6), A351605 (k=7), A351606 (k=8), A351607 (k=9), A351608 (k=10).
Cf. A082293 (a(n)==2), A082294 (a(n)==3).

a046951 = sum . map a010052 . a027750_row
-- Reinhard Zumkeller, Dec 16 2013

[#[d: d in Divisors(n)|IsSquare(d)]:n in [1..120]]; // Marius A. Burtea, Jan 21 2020

A046951 := proc(n)
    local a,s;
    a := 1 ;
    for p in ifactors(n)[2] do
        a := a*(1+floor(op(2,p)/2)) ;
    end do:
    a ;
end proc: # R. J. Mathar, Sep 17 2012
# Alternatively:
isbidivisible := (n, d) -> igcd(n, d) = d and igcd(n/d, d) = d:
a := n -> nops(select(k -> isbidivisible(n, k), [seq(1..n)])): # Peter Luschny, Jun 13 2025

a[n_] := Length[ Select[ Divisors[n], IntegerQ[Sqrt[#]]& ] ]; Table[a[n], {n, 1, 105}] (* Jean-François Alcover, Jun 26 2012 *)
Table[Length[Intersection[Divisors[n], Range[10]^2]], {n, 100}] (* Alonso del Arte, Dec 10 2012 *)
a[ n_] := If[ n < 1, 0, Sum[ Mod[ DivisorSigma[ 0, d], 2], {d, Divisors @ n}]]; (* Michael Somos, Jun 13 2014 *)
a[ n_] := If[ n < 2, Boole[ n == 1], Times @@ (Quotient[ #[[2]], 2] + 1 & /@ FactorInteger @ n)]; (* Michael Somos, Jun 13 2014 *)
a[ n_] := If[ n < 0, 0, SeriesCoefficient[ Sum[ x^k^2 / (1 - x^k^2), {k, Sqrt @ n}], {x, 0, n}]]; (* Michael Somos, Jun 13 2014 *)
f[p_, e_] := 1 + Floor[e/2]; a[1] = 1; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Sep 15 2020 *)

a(n)=my(f=factor(n));for(i=1,#f[,1],f[i,2]\=2);numdiv(factorback(f)) \\ Charles R Greathouse IV, Dec 11 2012

a(n) = direuler(p=2, n, 1/((1-X^2)*(1-X)))[n]; \\ Michel Marcus, Mar 08 2015

a(n)=factorback(apply(e->e\2+1, factor(n)[,2])) \\ Charles R Greathouse IV, Sep 17 2015

from math import prod
from sympy import factorint
def A046951(n): return prod((e>>1)+1 for e in factorint(n).values()) # Chai Wah Wu, Aug 04 2024

def is_bidivisible(n, d) -> bool: return gcd(n, d) == d and gcd(n//d, d) == d
def aList(n) -> list[int]: return [k for k in range(1, n+1) if is_bidivisible(n, k)]
print([len(aList(n)) for n in range(1, 126)])  # Peter Luschny, Jun 13 2025

(definec (A046951 n) (if (= 1 n) 1 (* (A008619 (A007814 n)) (A046951 (A064989 n)))))
(define (A008619 n) (+ 1 (/ (- n (modulo n 2)) 2)))
;; Antti Karttunen, Nov 14 2016

a(p^k) = A008619(k) = [k/2] + 1. a(A002110(n)) = 1 for all n. (This is true for any squarefree number, A005117). - Original notes clarified by Antti Karttunen, Nov 14 2016
a(n) = |{(i, j) : i*j = n AND i|j}| = |{(i, j) : i*j^2 = n}|. Also tau(A000188(n)), where tau = A000005.
Multiplicative with p^e --> floor(e/2) + 1, p prime. - Reinhard Zumkeller, May 20 2007
a(A130279(n)) = n and a(m) <> n for m < A130279(n); A008966(n)=0^(a(n) - 1). - Reinhard Zumkeller, May 20 2007
Inverse Moebius transform of characteristic function of squares (A010052). Dirichlet g.f.: zeta(s)*zeta(2s).
G.f.: Sum_{k > 0} x^(k^2)/(1 - x^(k^2)). - Vladeta Jovovic, Dec 13 2002
a(n) = Sum_{k=1..A000005(n)} A010052(A027750(n,k)). - Reinhard Zumkeller, Dec 16 2013
a(n) = Sum_{k = 1..n} ( floor(n/k^2) - floor((n-1)/k^2) ). - Peter Bala, Feb 17 2014
From Antti Karttunen, Nov 14 2016: (Start)
a(1) = 1; for n > 1, a(n) = A008619(A007814(n)) * a(A064989(n)).
a(n) = A278161(A156552(n)). (End)
G.f.: Sum_{k>0}(theta(q^k)-1)/2, where theta(q)=1+2q+2q^4+2q^9+2q^16+... - Mamuka Jibladze, Dec 04 2016
From  Antti Karttunen, Nov 12 2017: (Start)
a(n) = A000005(n) - A056595(n).
a(n) = 1 + A071325(n).
a(n) = 1 + A001222(A293515(n)). (End)
L.g.f.: -log(Product_{k>=1} (1 - x^(k^2))^(1/k^2)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Jul 30 2018
a(n) = Sum_{d|n} A000005(d) * A008836(n/d). - Torlach Rush, Jan 21 2020
a(n) = A000005(sqrt(A008833(n))). - Amiram Eldar, Jul 07 2020
a(n) = Sum_{d divides n} mu(core(d)^2), where core(n) = A007913(n). - Peter Bala, Jan 24 2024

Data section filled up to 125 terms and wrong claim deleted from Crossrefs section by Antti Karttunen, Nov 14 2016

A030514 a(n) = prime(n)^4.

Original entry on oeis.org

16, 81, 625, 2401, 14641, 28561, 83521, 130321, 279841, 707281, 923521, 1874161, 2825761, 3418801, 4879681, 7890481, 12117361, 13845841, 20151121, 25411681, 28398241, 38950081, 47458321, 62742241, 88529281, 104060401, 112550881, 131079601, 141158161

Offset: 1

Jeff Burch

Numbers with 5 divisors (1, p, p^2, p^3, p^4, where p is the n-th prime). - Alexandre Wajnberg, Jan 15 2006
Subsequence of A036967. - Reinhard Zumkeller, Feb 05 2008
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
The general product formula for even s is: product_{p = A000040} (p^s-1)/(p^s+1) = 2*Bernoulli(2s)/( binomial(2s, s)*Bernoulli^2(s)), where the infinite product is over all primes. Here, with s = 4, product_{n = 1, 2, ...} (a(n)-1)/(a(n)+1) = 6/7. In A030516, where s = 6, the product of the ratios is 691/715. For s = 8, the 8th row in A120458, the corresponding product of ratios is 7234/7293. - R. J. Mathar, Feb 01 2009
Except for the first three terms, all others are congruent to 1 mod 240. - Robert Israel, Aug 29 2014

R. J. Mathar, Table of n, a(n) for n = 1..1000
Eric Weisstein's World of Mathematics, Prime Power.
OEIS Wiki, Index entries for number of divisors
Index to sequences related to prime signature

Cf. A030078, A085964, A131991, A131992, A000005, A000040, A001248, A157290, A215267.

Cf. A258601.

a030514 = (^ 4) . a000040
a030514_list = map (^ 4) a000040_list
-- Reinhard Zumkeller, Jun 03 2015

[NthPrime(n)^4: n in [1..100] ]; // Vincenzo Librandi, Apr 22 2011

map(p -> p^4, select(isprime,[2,seq(2*i+1,i=1..100)])); # Robert Israel, Aug 29 2014

Array[Prime[#]^4 &, 5!] (* Vladimir Joseph Stephan Orlovsky, Sep 01 2008 *)

a(n)=prime(n)^4 \\ Charles R Greathouse IV, Mar 21 2013

[p**4 for p in prime_range(100)]
# Zerinvary Lajos, May 15 2007

a(n) = A000040(n)^(5-1) = A000040(n)^4, where 5 is the number of divisors of a(n). - Omar E. Pol, May 06 2008
A000005(a(n)) = 5. - Alexandre Wajnberg, Jan 15 2006
A056595(a(n)) = 2. - Reinhard Zumkeller, Aug 15 2011
Sum_{n>=1} 1/a(n) = P(4) = 0.0769931397... (A085964). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(4)/zeta(8) = 105/Pi^4 (A157290).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(4) = 90/Pi^4 (A215267). (End)

Description corrected by Eric W. Weisstein

A054753 Numbers which are the product of a prime and the square of a different prime (p^2 * q).

Original entry on oeis.org

12, 18, 20, 28, 44, 45, 50, 52, 63, 68, 75, 76, 92, 98, 99, 116, 117, 124, 147, 148, 153, 164, 171, 172, 175, 188, 207, 212, 236, 242, 244, 245, 261, 268, 275, 279, 284, 292, 316, 325, 332, 333, 338, 356, 363, 369, 387, 388, 404, 412, 423, 425, 428, 436, 452

Offset: 1

Henry Bottomley, Apr 25 2000

nonn

A178254(a(n)) = 4; union of A095990 and A096156. - Reinhard Zumkeller, May 24 2010
Numbers with prime signature (2,1) = union of numbers with ordered prime signature (1,2) and numbers with ordered prime signature (2,1) (restating second part of above comment). - Daniel Forgues, Feb 05 2011
A056595(a(n)) = 4. - Reinhard Zumkeller, Aug 15 2011
For k>1, Sum_{n>=1} 1/a(n)^k = P(k) * P(2*k) - P(3*k), where P is the prime zeta function. - Enrique Pérez Herrero, Jun 27 2012
Also numbers n with A001222(n)=3 and A001221(n)=2. - Enrique Pérez Herrero, Jun 27 2012
A089233(a(n)) = 2. - Reinhard Zumkeller, Sep 04 2013
Subsequence of the triprimes (A014612). If a(n) is even, then a(n)/2 is semiprime (A001358). - Wesley Ivan Hurt, Sep 08 2013
From Bernard Schott, Sep 16 2017: (Start)
These numbers are called "Nombres d'Einstein" on the French site "Diophante" (see link) because a(n) = m * c^2 where m and c are two different primes.
The numbers 44 = 2^2 * 11 and 45 = 3^2 * 5 are the two smallest consecutive "Einstein numbers"; 603, 604, 605 are the three smallest consecutive integers in this sequence. It's not possible to get more than five such consecutive numbers (proof in the link); the first set of five such consecutive numbers begins at the 17-digit number 10093613546512321. Where does the first sequence of four consecutive "Einstein numbers" begin? (End) [corrected by Jon E. Schoenfield, Sep 20 2017]
The first set of four consecutive integers in this sequence begins at the 11-digit number 17042641441. (Each such set must include two even numbers, one of which is of the form 2^2*q, the other of the form p^2*2; a quick search, taking the factorizations of consecutive integers before and after numbers of the latter form, shows that the number of sets of four consecutive k-digit integers in this sequence is 1, 7, 12, 18 for k = 11, 12, 13, 14, respectively.) - Jon E. Schoenfield, Sep 16 2017
The first 13 sets of 5 consecutive integers in this sequence have as their first terms 10093613546512321, 14414905793929921, 266667848769941521, 562672865058083521, 1579571757660876721, 1841337567664174321, 2737837351207392721, 4456162869973433521, 4683238426747860721, 4993613853242910721, 5037980611623036721, 5174116847290255921, 5344962129269790721. Each of these numbers except for the last is 7^2 times a prime; the last is 23^2 times a prime. - Jon E. Schoenfield, Sep 17 2017

			a(1) = 12 because 12 = 2^2*3 is the smallest number of the form p^2*q.

Reinhard Zumkeller, Table of n, a(n) for n = 1..1000
Guilhem Castagnos, Antoine Joux, Fabien Laguillaumie, and Phong Q. Nguyen, Factoring pq^2 with quadratic forms: nice cryptanalyses, Advances in Cryptology - ASIACRYPT 2009. Lecture Notes in Computer Science Volume 5912 (2009), pp. 469-486.
Diophante, A 350, Les Nombres d'Einstein (in French).
Mathematics Stack Exchange, Sequence of numbers with prime factorization pq^2
René Peralta and Eiji Okamoto, Faster factoring of integers of a special form (1996).
Index to sequences related to prime signature

Cf. A001221, A001222, A001358, A014612, A056595, A089233, A095990, A096156, A178254.
Numbers with 6 divisors (A030515) which are not 5th powers of primes (A050997).
Subsequence of A325241.  Supersequence of A096156.
Table giving for each subsequence the corresponding number of groups of order p^2*q, from Bernard Schott, Jan 23 2022
  -------------------------------------------------------------------------------
  | Subsequence  | A350638 | A143928 | A350115 | A349495 | A350245 | A350422 (*)|
  -------------------------------------------------------------------------------
  |A000001(p^2*q)| (q+9)/2 |    5    |    5    |    4    |    3    |      2     |
  -------------------------------------------------------------------------------
(*) A350422 equals disjoint union of A350332 (pA350421 (p>q).

Select[Range[12,452], {1,2}==Sort[Last/@FactorInteger[ # ]]&] (* Zak Seidov, Jul 19 2009 *)
With[{nn=60},Take[Union[Flatten[{#[[1]]#[[2]]^2,#[[1]]^2 #[[2]]}&/@ Subsets[ Prime[Range[nn]],{2}]]],nn]] (* Harvey P. Dale, Dec 15 2014 *)

is(n)=vecsort(factor(n)[,2])==[1,2]~ \\ Charles R Greathouse IV, Dec 30 2014

for(n=1, 1e3, if(numdiv(n) - bigomega(n) == 3, print1(n, ", "))) \\ Altug Alkan, Nov 24 2015

from sympy import factorint
def ok(n): return sorted(factorint(n).values()) == [1, 2]
print([k for k in range(453) if ok(k)]) # Michael S. Branicky, Dec 18 2021

from math import isqrt
from sympy import primepi, primerange, integer_nthroot
def A054753(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(x//p**2) for p in primerange(isqrt(x)+1))+primepi(integer_nthroot(x,3)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

Link added and incorrect Mathematica code removed by David Bevan, Sep 17 2011

A030516 Numbers with 7 divisors. 6th powers of primes.

Original entry on oeis.org

64, 729, 15625, 117649, 1771561, 4826809, 24137569, 47045881, 148035889, 594823321, 887503681, 2565726409, 4750104241, 6321363049, 10779215329, 22164361129, 42180533641, 51520374361, 90458382169, 128100283921

Offset: 1

Jeff Burch

These are the numbers p^6 with p prime. - Jorge Coveiro, Apr 13 2004

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

R. J. Mathar, Table of n, a(n) for n = 1..1000
OEIS Wiki, Index entries for number of divisors
Index to sequences related to prime signature

Subsequence of A046306.

Cf. A013664, A050997, A085966, A131993, A000005, A000040, A001248, A030514, A258603, A269404.

a030516 = (^ 6) . a000040
a030516_list = map (^ 6) a000040_list
-- Reinhard Zumkeller, Jun 03 2015

[NthPrime(n)^6: n in [1..400]];  // Bruno Berselli, Jul 30 2011

Array[Prime[ # ]^6&,60] (* Vladimir Joseph Stephan Orlovsky, Dec 14 2009 *)

for(n=1,1e3,print1(prime(n)^6", "));  \\ Bruno Berselli, Jul 30 2011

[p**6 for p in prime_range(100)]
# Zerinvary Lajos, May 16 2007

a(n) = A000040(n)^(7-1) = A000040(n)^6. - Omar E. Pol, May 06 2008
A056595(a(n)) = 3. - Reinhard Zumkeller, Aug 15 2011
Sum_{n>=1} 1/a(n) = P(6) = 0.0170700868... (A085966). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(6)/zeta(12) = 675675/(691*Pi^6) (A269404).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(6) = 945/Pi^6 = 1/A013664. (End)

A050997 Fifth powers of primes.

Original entry on oeis.org

32, 243, 3125, 16807, 161051, 371293, 1419857, 2476099, 6436343, 20511149, 28629151, 69343957, 115856201, 147008443, 229345007, 418195493, 714924299, 844596301, 1350125107, 1804229351, 2073071593, 3077056399, 3939040643, 5584059449, 8587340257, 10510100501

Offset: 1

Eric W. Weisstein

Numbers k such that A062799(k) = 5.

Let r(n) = (a(n)+1)/(a(n)-1) if a(n) mod 4 = 3, (a(n)-1)/(a(n)+1) otherwise; then Product_{n>=1} r(n) = (31/33) * (244/242) * (3124/3126) * (16808/16806) * ... = 246016/259875. - Dimitris Valianatos, Mar 09 2020

T. D. Noe, Table of n, a(n) for n = 1..1000
Xavier Gourdon and Pascal Sebah, Some Constants from Number theory.
Eric Weisstein's World of Mathematics, MathWorld: Prime Power.
Index to sequences related to prime signature

Cf. A000040, A001248, A030078, A030514, A085965, A131992, A131993, A013663, A157291.

Cf. A258602.

a050997 = (^ 5) . a000040
a050997_list = map (^ 5) a000040_list
-- Reinhard Zumkeller, Jun 03 2015

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

Array[Prime[ # ]^5 &, 30] (* Vladimir Joseph Stephan Orlovsky, May 01 2008 *)

PARI
```
vector(66,n,prime(n)^5)
```

A056595(a(n)) = 3. - Reinhard Zumkeller, Aug 15 2011
Sum_{n>=1} 1/a(n) = P(5) = 0.0357550174... (A085965). - Amiram Eldar, Jul 27 2020
From Amiram Eldar, Jan 23 2021: (Start)
Product_{n>=1} (1 + 1/a(n)) = zeta(5)/zeta(10) (A157291).
Product_{n>=1} (1 - 1/a(n)) = 1/zeta(5) = 1/A013663. (End)

A000430 Primes and squares of primes.

Original entry on oeis.org

2, 3, 4, 5, 7, 9, 11, 13, 17, 19, 23, 25, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223

Offset: 1

R. Muller

Also numbers n such that the product of proper divisors is < n.
See A050216 for lengths of blocks of consecutive primes. - Reinhard Zumkeller, Sep 23 2011
Numbers q > 1 such that d(q) < 4. Numbers k such that the number of ways of writing k = m + t is equal to the number of ways of writing k = r*s, where m|t and r|s. - Juri-Stepan Gerasimov, Oct 14 2017
Called multiplicatively deficient numbers by Chau (2004). - Amiram Eldar, Jun 29 2022

F. Smarandache, Definitions solved and unsolved problems, conjectures and theorems in number theory and geometry, edited by M. Perez, Xiquan Publishing House 2000
F. Smarandache, Sequences of Numbers Involved in Unsolved Problems, Hexis, Phoenix, 2006.

T. D. Noe, Table of n, a(n) for n = 1..1000
William Chau, The tau, sigma, rho functions, and some related numbers, Pi Mu Epsilon Journal, Vol. 11, No. 10 (Spring 2004), pp. 519-534; entire issue.
F. Smarandache, Sequences of Numbers Involved in Unsolved Problems.

Union of A000040 and A001248.

Cf. A007422, A010051, A010055, A032741, A046951, A050216, A056595, A058080, A064911, A084110, A084114, A293575.

a000430 n = a000430_list !! (n-1)
a000430_list = m a000040_list a001248_list where
   m (x:xs) (y:ys) | x < y = x : m xs (y:ys)
                   | x > y = y : m (x:xs) ys
-- Reinhard Zumkeller, Sep 23 2011

nn = 223; t = Union[Prime[Range[PrimePi[nn]]], Prime[Range[PrimePi[Sqrt[nn]]]]^2] (* T. D. Noe, Apr 11 2011 *)
Module[{upto=250,prs},prs=Prime[Range[PrimePi[upto]]];Select[Join[ prs,prs^2], #<=upto&]]//Sort (* Harvey P. Dale, Oct 08 2016 *)

is(n)=isprime(n) || (issquare(n,&n) && isprime(n)) \\ Charles R Greathouse IV, Sep 04 2013

from math import isqrt
from sympy import primepi
def A000430(n):
    def f(x): return n+x-primepi(x)-primepi(isqrt(x))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return int(m) # Chai Wah Wu, Aug 09 2024

A084114(a(n)) = 0, see also A084110. - Reinhard Zumkeller, May 12 2003
A109810(a(n)) = 2. - Reinhard Zumkeller, May 24 2010
A010051(a(n)) + A010055(a(n))*A064911(a(n)) = 1;
A056595(a(n)) = 1. - Reinhard Zumkeller, Aug 15 2011
A032741(a(n)) = A046951(a(n)); A293575(a(n)) = 0. - Juri-Stepan Gerasimov, Oct 14 2017
The number of terms not exceeding x is N(x) ~ (x + 2*sqrt(x))/log(x) (Chau, 2004). - Amiram Eldar, Jun 29 2022

cp's OEIS Frontend

A194095 Record values in A056595.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

A194096 Where record values occur in A056595.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

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

Extensions

A030078 Cubes of primes.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Mathematica

PARI

Python

Sage

Formula

A046951 a(n) is the number of squares dividing n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

Python

Python

Scheme

Formula

Extensions

A030514 a(n) = prime(n)^4.

Views