cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 21-27 of 27 results.

A069280 19-almost primes (generalization of semiprimes).

Original entry on oeis.org

524288, 786432, 1179648, 1310720, 1769472, 1835008, 1966080, 2654208, 2752512, 2883584, 2949120, 3276800, 3407872, 3981312, 4128768, 4325376, 4423680, 4456448, 4587520, 4915200, 4980736, 5111808, 5971968, 6029312, 6193152
Offset: 1

Views

Author

Rick L. Shepherd, Mar 13 2002

Keywords

Comments

Product of 19 not necessarily distinct primes.
Divisible by exactly 19 prime powers (not including 1).

Crossrefs

Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275 (r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), this sequence (r = 19), A069281 (r = 20). - Jason Kimberley, Oct 02 2011

Programs

  • PARI
    k=19; start=2^k; finish=8000000; v=[]; for(n=start,finish, if(bigomega(n)==k,v=concat(v,n))); v
    
  • Python
    from math import prod, isqrt
    from sympy import primerange, integer_nthroot, primepi
    def A069280(n):
        def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
        def f(x): return int(n-1+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,19)))
        kmin, kmax = 1,2
        while f(kmax) >= kmax:
            kmax <<= 1
        while True:
            kmid = kmax+kmin>>1
            if f(kmid) < kmid:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmax # Chai Wah Wu, Aug 23 2024

Formula

Product p_i^e_i with Sum e_i = 19.

A101695 a(n) = n-th n-almost prime.

Original entry on oeis.org

2, 6, 18, 40, 108, 224, 480, 1296, 2688, 5632, 11520, 25600, 53248, 124416, 258048, 540672, 1105920, 2228224, 4587520, 9830400, 19922944, 40894464, 95551488, 192937984, 396361728, 822083584, 1660944384, 3397386240, 6845104128
Offset: 1

Views

Author

Jonathan Vos Post, Dec 12 2004

Keywords

Comments

A k-almost-prime is a positive integer that has exactly k prime factors, counted with multiplicity.
This is the diagonalization of the set of sequences {j-almost prime(k)}. The cumulative sums of this sequence are in A101696. This is the diagonal just below A078841.

Examples

			a(1) = first 1-almost prime = first prime = A000040(1) = 2.
a(2) = 2nd 2-almost prime = 2nd semiprime = A001358(2) = 6.
a(3) = 3rd 3-almost prime = A014612(3) = 18.
a(4) = 4th 4-almost prime = A014613(4) = 40.
a(5) = 5th 5-almost prime = A014614(5) = 108.
		

Crossrefs

Programs

  • Maple
    A101695 := proc(n)
        local s,a ;
        s := 0 ;
        for a from 2^n do
            if numtheory[bigomega](a) = n then
                s := s+1 ;
                if s = n then
                    return a;
                end if;
            end if;
        end do:
    end proc: # R. J. Mathar, Aug 09 2012
  • Mathematica
    AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
    AlmostPrime[k_, n_] := Block[{e = Floor[ Log[2, n] + k], a, b}, a = 2^e; Do[b = 2^p; While[ AlmostPrimePi[k, a] < n, a = a + b]; a = a - b/2, {p, e, 0, -1}]; a + b/2]; AlmostPrime[1, 1] = 2; lst = {}; Do[ AppendTo[lst, AlmostPrime[n, n]], {n, 30}]; lst (* Robert G. Wilson v, Oct 07 2007 *)
  • Python
    from math import prod, isqrt
    from sympy import primerange, primepi, integer_nthroot
    def A101695(n):
        if n == 1: return 2
        def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
        def f(x): return int(n-1+x-sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,n)))
        kmin, kmax = 1,2
        while f(kmax) >= kmax:
            kmax <<= 1
        while True:
            kmid = kmax+kmin>>1
            if f(kmid) < kmid:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmax # Chai Wah Wu, Aug 23 2024

Formula

Conjecture: lim_{ n->inf.} a(n+1)/a(n) = 2. - Robert G. Wilson v, Oct 07 2007, Nov 13 2007
Stronger conjecture: a(n)/(n * 2^n) is polylogarithmic in n. That is, there exist real numbers b < c such that (log n)^b < a(n)/(n * 2^n) < (log n)^c for large enough n. Probably b and c can be chosen close to 0. - Charles R Greathouse IV, Aug 28 2012

Extensions

a(21)-a(30) from Robert G. Wilson v, Feb 11 2006
a(12) corrected by N. J. A. Sloane, Nov 23 2007

A125149 a(n) is the least k such that the n-almost prime count is positive and equal to the (n-1)-almost prime count. a(0) = 1.

Original entry on oeis.org

1, 2, 10, 15495, 151165506066
Offset: 0

Views

Author

Keywords

Comments

Unlike any of the prime number races in which any particular form may lead or trail, this sequence demonstrates that although the count of numbers having k prime factors begins by trailing the count for k-1 prime factors, eventually they exchange positions in the race. This can be seen by looking at A126279 or A126280.
The fundamental theorem of arithmetic, or unique factorization theorem, states that every natural number greater than 1 either is itself a prime number, or can be written as a unique product of prime numbers. It had a proof sketched by Euclid, then corrected and completed in "Disquisitiones Arithmeticae" [Carl Friedrich Gauss, 1801]. It fails in many rings of algebraic integers [Ernst Kummer, 1843], a discovery initiating algebraic number theory. Counting the elements in the unique product of prime numbers classifies natural numbers into primes, semiprimes, 3-almost primes and so on. This sequence quantifies a previously undescribed structure to that classification.
We took the first k where the two relevant counts are the same. If instead we took the least k such that the n-almost prime count from k onwards exceeds the (n-1)-almost prime count, the sequence would begin: 3, 34, 15530, ... [see A180126].
The prime count and the semiprime count are identical for 1, 10, 15, 16, 22, 25, 29, 30, 33.
The semiprime count and the 3-almost prime count are identical for 1, 2, 3, 15495, 15496, 15497, 15498, 15508, 15524, 15525, 15529.
The numbers of 3-almost primes and 4-almost primes are equal at 151165506066 and 731 larger numbers, the last one being 151165607041. See A180126. - T. D. Noe, Aug 11 2010
Landau's asymptotic formula suggests that a(n) is about exp(exp(n-1)). - Charles R Greathouse IV, Mar 14 2011

Examples

			a(1) = 2 since 1 has no prime factors and 2 has one prime factor, therefore prime factor counts of 0 and 1 occur equally often in the first 2 integers.
a(2) = 10 since there are 4 primes {2, 3, 5 & 7} and 4 semiprimes {4, 6, 9 & 10} less than or equal to 10.
a(4) = 151165506066 since there are 32437255807 4-almost primes and 3-almost primes <= a(4).
		

Crossrefs

Sequences listing r-almost primes, that is, k such that A001222(k) = r: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275 (r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20).
Cf. A180126.

Programs

  • Mathematica
    AlmostPrimePi[k_Integer, n_] := Module[{a, i}, a[0] = 1; If[k == 1, PrimePi[n], Sum[PrimePi[n/Times @@ Prime[Array[a, k - 1]]] - a[k - 1] + 1, Evaluate[ Sequence @@ Table[{a[i], a[i - 1], PrimePi[(n/Times @@ Prime[Array[a, i - 1]])^(1/(k - i + 1))]}, {i, k - 1}]]]]]; (* Eric W. Weisstein, Feb 07 2006 *)
    f[n_] := Block[{k = 2^n}, While[AlmostPrimePi[n, k] < AlmostPrimePi[n - 1, k], k++ ]; k];

Extensions

Changed 33 to 34 in a comment. - T. D. Noe, Aug 11 2010
Edited by Peter Munn, Dec 17 2022

A101696 a(n) = sum(i=1,n)(i-th i-almost prime). Cumulative sums of A101695.

Original entry on oeis.org

2, 8, 26, 66, 174, 398, 878, 2174, 4862, 10494, 22014, 45054, 98302, 222718, 480766, 1021438, 2127358, 4355582, 8943102, 18773502, 38696446, 79590910, 175142398, 368080382, 764442110, 1586525694, 3247470078, 6644856318, 13489960446
Offset: 1

Views

Author

Jonathan Vos Post, Dec 12 2004, Sep 28 2006

Keywords

Comments

It seems that this sum can never be a prime after a(1) = 2, since the n-th n-almost prime is always even. The number of prime factors (with multiplicity) of a(n) is 1, 3, 2, 3, 3, 2, 2, 2, 4, 5, 4, 4, 3, 3, 5, 4, 3, 4, 7, 4, 2, 5, 5, 2, 3, 7, 4, 3, 4.
This is the diagonalization of the set of sequences {j-almost prime(k)}. The cumulative sums of this sequence are in A101696. a(1)=2 is prime. a(2)=8 is a 3-almost prime. a(3)=26 is a semiprime. a(4)=66 is a 3-almost prime. a(5)= 174 is a 3-almost prime. a(6)=398 is a semiprime. a(7)=878 is a semiprime. a(8)=2174 is a semiprime. a(9)=4862 is a 4-almost prime. a(10)=10494 is a 5-almost prime. a(11)=22014 is a 4-almost prime. a(12)=45054 is a 3-almost prime. a(13)=98302 is a 3-almost prime. a(14)=222718 is a 3-almost prime. a(15)=480766 is a 5-almost prime. a(16)=1021438 is a 4-almost prime. a(17)=2127358 is a 3-almost prime. a(18)=4355582 is a 4-almost prime. a(19)=8943102 is a 7-almost prime. a(20)=18773502 is a 4-almost prime. 21-almost numbers are not yet listed in the OEIS.

Examples

			a(1) = first 1-almost prime = first prime = A000040(1) = 2.
a(2) = 2 + 2nd 2-almost prime = 2 + A001358(2) = 2+ 6 = 8.
a(3) = a(2) + 3rd 3-almost prime = 8+A014612(3) = 8+18 = 26.
a(4) = a(3) + 4th 4-almost prime = 26+A014613(4) = 26+40 = 66.
a(5) = a(4) + 5th 5-almost prime = 66+A014614(5) = 66+108=174.
...
a(12) = a(11) + 12th 12-almost prime = 22014 + 23040 = 45054 (the first nontrivial palindrome in the sequence).
		

Crossrefs

Formula

a(1) = first 1-almost prime = first prime = A000040(1). a(2) = a(1) + 2nd 2-almost prime = a(1) + 2nd semiprime = A000040(1)+A001358(2). a(3) = a(2) + 3rd 3-almost prime = a(2) + A014612(3). a(4) = a(3) + 4th 4-almost prime = a(3) + A014613(4)... a(n) = a(n-1) + n-th n-almost prime.

Extensions

Edited by N. J. A. Sloane, Jul 03 2008 at the suggestion of R. J. Mathar

A337112 Smallest term of A337081 that has exactly n prime factors, or 0 if no such term exists.

Original entry on oeis.org

0, 4, 0, 90, 675, 1134, 6318, 4374, 32805, 255879, 1003833, 531441, 327544803, 20751953125, 225830078125, 91552734375, 1068115234375, 23651123046875, 316619873046875, 1697540283203125, 13256072998046875, 85353851318359375, 541210174560546875, 4518032073974609375, 58233737945556640625
Offset: 1

Views

Author

Matej Veselovac, Aug 16 2020

Keywords

Comments

a(n) is the smallest product of n primes that has unordered factorizations whose sums of factors are the same (is a term of A337080) and all of whose proper divisors have the complementary property: that every unordered factorization has a distinct sum of factors (i.e., all proper divisors are terms of A337037).

Crossrefs

Cf. A337113 (factors of terms).
Cf. A056472 (all factorizations of n).
Cf. r-almost primes: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275 (r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20).

Extensions

a(14) onward from David A. Corneth, Aug 26 2020

A123118 Partial products of A101695.

Original entry on oeis.org

2, 12, 216, 8640, 933120, 209018880, 100329062400, 130026464870400, 349511137571635200, 1968446726803449446400, 22676506292775737622528000, 522466704985552994823045120000, 27820307107070725868337506549760000
Offset: 1

Views

Author

Jonathan Vos Post, Sep 28 2006

Keywords

Comments

The number of prime factors (with multiplicity) of a(n) is T(n) = A000217(n) = n*(n+1)/2.

Examples

			a(1) = 2 = prime(1).
a(2) = 12 = 2 * 6 = prime(1) * semiprime(2) = 2^2 * 3.
a(3) = 216 = 2 * 6 * 18 = prime(1) * semiprime(2) * 3-almostprime(3) = 2^3 * 3^3.
a(4) = 8640 = 2 * 6 * 18 * 40 = prime(1) * semiprime(2) * 3-almostprime(3) * 4-almostprime(4) = 2^6 * 3^3 * 5.
a(15) = 893179304874387947794472921245209518407680000 = 2 * 6 * 18 * 40 * 108 * 224 * 480 * 1296 * 2688 * 5632 * 11520 * 23040 * 53248 * 124416 * 258048 = 2^88 * 3^23 * 5^4 * 7^3 * 11 * 13.
		

Crossrefs

Formula

a(n) = Prod[i=1..n] i-th i-almost prime = Prod[i=1..n] A101695(i).

A321590 Smallest number m that is a product of exactly n primes and is such that m-1 and m+1 are products of exactly n-1 primes.

Original entry on oeis.org

4, 50, 189, 1863, 10449, 447849, 4449249, 5745249, 3606422049, 16554218751, 105265530369, 1957645712385
Offset: 2

Views

Author

Zak Seidov, Nov 13 2018

Keywords

Comments

From Jon E. Schoenfield, Nov 15 2018: (Start)
If a(11) is odd, it is 16554218751.
If a(12) is odd, it is 105265530369.
If a(13) is odd, it is 1957645712385. (End)
a(11), a(12), and a(13) are indeed odd. - Giovanni Resta, Jan 04 2019
10^13 < a(14) <= 240455334218751, a(15) <= 2992278212890624. - Giovanni Resta, Jan 06 2019

Examples

			For n = 3, 50 = 2*5*5, and the numbers before and after 50 are 49 = 7*7 and 51 = 3*17.
		

Crossrefs

Cf. A078840.
Sequences listing r-almost primes, that is, the n such that A001222(n) = r: A000040 (r = 1), A001358 (r = 2), A014612 (r = 3), A014613 (r = 4), A014614 (r = 5), A046306 (r = 6), A046308 (r = 7), A046310 (r = 8), A046312 (r = 9), A046314 (r = 10), A069272 (r = 11), A069273 (r = 12), A069274 (r = 13), A069275(r = 14), A069276 (r = 15), A069277 (r = 16), A069278 (r = 17), A069279 (r = 18), A069280 (r = 19), A069281 (r = 20).

Programs

  • Mathematica
    a[n_] := Module[{o={0,0,0}, k=1}, While[o!={n-1,n,n-1}, o=Rest[AppendTo[o,PrimeOmega[k]]]; k++]; k-2]; Array[a,7,2] (* Amiram Eldar, Nov 14 2018 *)
  • PARI
    {for(n=2,10,for(k=2^n,10^12,if(n==bigomega(k) &&
    n-1==bigomega(k-1) && n-1==bigomega(k+1),print1(k", ");break())))}

Extensions

a(10) from Jon E. Schoenfield, Nov 14 2018
a(11)-a(13) from Giovanni Resta, Jan 04 2019
Previous Showing 21-27 of 27 results.