A085986 -id:A085986 - cp's OEIS frontend

A195086 Numbers k such that (number of prime factors of k counted with multiplicity) less (number of distinct prime factors of k) = 2.

8, 24, 27, 36, 40, 54, 56, 88, 100, 104, 120, 125, 135, 136, 152, 168, 180, 184, 189, 196, 225, 232, 248, 250, 252, 264, 270, 280, 296, 297, 300, 312, 328, 343, 344, 351, 375, 376, 378, 396, 408, 424, 440, 441, 450, 456, 459, 468, 472, 484, 488

Offset: 1

Harvey P. Dale, Sep 08 2011

nonn

From Amiram Eldar, Nov 07 2020: (Start)
Numbers whose powerful part (A057521) is either a cube of a prime (A030078) or a square of a squarefree semiprime (A085986).
The asymptotic density of this sequence is (6/Pi^2) * (Sum_{p prime} 1/(p^2*(p+1)) + Sum_{p=4} (-1)^(k+1)*(k-1)*P(k) + (Sum_{k>=2} (-1)^k*P(k))^2)/2 = 0.0963023158..., where P is the prime zeta function. (End)

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Alfred Rényi, On the density of certain sequences of integers, Publications de l'Institut Mathématique, Vol. 8 (1955), pp. 157-162.
Index entries for sequences computed from exponents in factorization of n

Subsequence of A048108.
Subsequences: A030078 and A085986.
Cf. A001221, A001222, A025487, A057521, A060687, A195069, A195087, A195088, A195089, A195090, A195091, A195092, A195093, A046660, A257851, A261256, A264959.

a195086 n = a195086_list !! (n-1)
a195086_list = filter ((== 2) . a046660) [1..]
-- Reinhard Zumkeller, Nov 29 2015

Select[Range[500],PrimeOmega[#]-PrimeNu[#]==2&]

is(n)=bigomega(n)-omega(n)==2 \\ Charles R Greathouse IV, Sep 14 2015

is(n)=my(f=factor(n)[,2]); vecsum(f)==#f+2 \\ Charles R Greathouse IV, Aug 01 2016

A001222(a(n)) - A001221(a(n)) = 2.

A046660(a(n)) = 2. - Reinhard Zumkeller, Nov 29 2015

A179668 Products of the 8th power of a prime and a distinct prime (p^8*q).

Original entry on oeis.org

768, 1280, 1792, 2816, 3328, 4352, 4864, 5888, 7424, 7936, 9472, 10496, 11008, 12032, 13122, 13568, 15104, 15616, 17152, 18176, 18688, 20224, 21248, 22784, 24832, 25856, 26368, 27392, 27904, 28928, 32512, 32805, 33536, 35072, 35584, 38144, 38656, 40192

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 23 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, Prime Signatures

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,8}; Select[Range[40000], f]
With[{nn=40},Take[Union[#[[1]]^8 #[[2]]&/@Flatten[Permutations/@Subsets[ Prime[Range[nn]],{2}],1]],nn]] (* Harvey P. Dale, Jan 20 2016 *)

list(lim)=my(v=List(),t);forprime(p=2,(lim\2)^(1/8),t=p^8;forprime(q=2,lim\t,if(p==q,next);listput(v,t*q)));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from sympy import primepi, primerange, integer_nthroot
def A179668(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**8) for p in primerange(integer_nthroot(x,8)[0]+1))+primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

A179692 Numbers of the form p^9*q where p and q are distinct primes.

Original entry on oeis.org

1536, 2560, 3584, 5632, 6656, 8704, 9728, 11776, 14848, 15872, 18944, 20992, 22016, 24064, 27136, 30208, 31232, 34304, 36352, 37376, 39366, 40448, 42496, 45568, 49664, 51712, 52736, 54784, 55808, 57856, 65024, 67072, 70144, 71168, 76288, 77312, 80384, 83456

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, Prime Signatures

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688, A179689, A179690, A179691.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,9}; Select[Range[90000], f]

list(lim)=my(v=List(),t);forprime(p=2, (lim\2)^(1/9), t=p^9;forprime(q=2, lim\t, if(p==q, next);listput(v,t*q))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 24 2011

from sympy import primepi, integer_nthroot, primerange
def A179692(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**9) for p in primerange(integer_nthroot(x,9)[0]+1))+primepi(integer_nthroot(x,10)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

A376936 Powerful numbers divisible by cubes of 2 distinct primes.

Original entry on oeis.org

216, 432, 648, 864, 1000, 1296, 1728, 1944, 2000, 2592, 2744, 3375, 3456, 3888, 4000, 5000, 5184, 5400, 5488, 5832, 6912, 7776, 8000, 9000, 9261, 10000, 10125, 10368, 10584, 10648, 10800, 10976, 11664, 13500, 13824, 15552, 16000, 16200, 16875, 17496, 17576, 18000

Offset: 1

Michael De Vlieger, Oct 16 2024

Numbers m with coreful divisors d, m/d such that neither d | m/d nor m/d | d, i.e., numbers m such that there exists a divisor pair (d, m/d) such that rad(d) = rad(m/d) but gcd(d, m/d) > 1 is neither d nor m/d, where rad = A007947. Divisors in each pair must be dissimilar and each in A126706.
Proper subset of A320966.
Contains A372695, A177493, and A162142. Does not contain A085986.

			216 is in the sequence since rad(12) | rad(18), but 12 does not divide 18 and 18 does not divide 12.
432 is a term since rad(18) | rad(24), but 18 does not divide 24 and 24 does not divide 18.
Table of coreful divisors d, a(n)/d such that neither d | a(n)/d nor a(n)/d | d for select a(n)
   n |   a(n)   divisor pairs d X a(n)/d
  ---------------------------------------------------------------------------
   1 |   216:   12 X 18;
   2 |   432:   18 X 24;
   3 |   648:   12 X 54;
   4 |   864:   24 X 36, 18 X 48;
   5 |  1000:   20 X 50;
   6 |  1296:   24 X 54;
   7 |  1728:   18 X 96, 36 X 48;
   8 |  1944:   12 X 162, 36 X 54;
   9 |  2000:   40 X 50;
  10 |  2592:   24 X 108, 48 X 54;
  11 |  2744:   28 X 98;
  12 |  3375:   45 X 75;
  13 |  3456:   18 X 192, 36 X 96, 48 X 72;
  22 |  7776:   24 X 324, 48 X 162, 54 X 144, 72 X 108;
  58 | 31104:   48 X 648, 54 X 576, 96 X 324, 108 X 288, 144 X 216, 162 X 192

Michael De Vlieger, Table of n, a(n) for n = 1..10000
Michael De Vlieger, Notes on this sequence

Cf. A007947, A030078, A085986, A126706, A162142, A177493, A286708, A307958, A320966, A361430, A370329, A372695, A376514.

Cf. A082020, A082695.

Union@ Select[
  Flatten@ Table[a^2*b^3, {b, Surd[#, 3]}, {a, Sqrt[#/b^3]}] &[20000],
  Length@ Select[FactorInteger[#][[All, -1]], # > 2 &] >= 2 &]

Sum_{n>=1} 1/a(n) = zeta(2)*zeta(3)/zeta(6) - (15/Pi^2) * (1 + Sum_{prime} 1/((p-1)*(p^2+1))) = 0.021194288968234037106579437374641326044... . - Amiram Eldar, Nov 08 2024

A074985 Squares of semiprimes (A001358).

Original entry on oeis.org

16, 36, 81, 100, 196, 225, 441, 484, 625, 676, 1089, 1156, 1225, 1444, 1521, 2116, 2401, 2601, 3025, 3249, 3364, 3844, 4225, 4761, 5476, 5929, 6724, 7225, 7396, 7569, 8281, 8649, 8836, 9025, 11236, 12321, 13225, 13924, 14161, 14641, 14884, 15129

Offset: 1

Jani Melik, Oct 07 2002

Disjoint union of 4th powers of primes, A030514, and squares of squarefree semiprimes, A085986. - M. F. Hasler, Nov 12 2021

			4 is divisible by 2 (twice) and 4*4 = 16.
6 is divisible by exactly 2 and 3 and 6*6 = 36.

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

Cf. A001358, A085548, A085964.

Cf. A030514 (4th powers of primes), A085986 (squares of squarefree semiprimes).

a074985 = a000290 . a001358  -- Reinhard Zumkeller, Aug 02 2012

readlib(issqr): ts_kv_sp := proc(n); if (numtheory[bigomega](n)=4 and issqr(n)='true') then RETURN(n); fi; end: seq(ts_kv_sp(i), i=1..50000);

Select[Range[200],PrimeOmega[#]==2&]^2 (* Harvey P. Dale, Oct 03 2011 *)

is(n)=if(issquare(n,&n), isprimepower(n)==2 || factor(n)[,2]==[1,1]~, 0) \\ Charles R Greathouse IV, Oct 16 2015

list(lim)=lim=sqrtint(lim\1); my(v=List()); forprime(p=2, sqrtint(lim), forprime(q=p, lim\p, listput(v, (p*q)^2))); Set(v) \\ Charles R Greathouse IV, Nov 13 2021

a(n) ~ (n log n/log log n)^2. - Charles R Greathouse IV, Oct 16 2015

Sum_{n>=1} 1/a(n) = (P(2)^2 + P(4))/2 = (A085548^2 + A085964)/2 = 0.1407604343..., where P is the prime zeta function. - Amiram Eldar, Oct 30 2020

A179642 Product of exactly 5 primes, 3 of which are distinct.

Original entry on oeis.org

120, 168, 180, 252, 264, 270, 280, 300, 312, 378, 396, 408, 440, 450, 456, 468, 520, 552, 588, 594, 612, 616, 680, 684, 696, 700, 702, 728, 744, 750, 760, 828, 882, 888, 918, 920, 945, 952, 980, 984, 1026, 1032, 1044, 1064, 1100, 1116, 1128, 1144, 1160

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 21 2010

nonn

			120=2^3*3*5, 168=2^3*3*7, 180=2^2*3^2*5, 252=2^2*3^2*7, 264=2^3*3*11, 270=2*3^3*5

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, Prime Signatures

Union of A189975 and A179643.

Cf. A006881, A007304, A085986, A085987, A178739.

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,1,3} || Sort[Last/@FactorInteger[n]]=={1,2,2}; Select[Range[2000], f]

list(lim)=my(v=List(),t);forprime(p=2,(lim\6)^(1/3),forprime(q=2,sqrt(lim\p^3),if(p==q,next);t=p^3*q;forprime(r=q+1,lim\t,if(p==r,next);listput(v,t*r))));forprime(p=2,sqrt(lim\12),forprime(q=p+1,sqrt(lim\p^2\2),t=(p*q)^2;forprime(r=2,lim\t,if(p==r||q==r,next);listput(v,t*r))));vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 19 2011

A255231 The number of factorizations n = Product_i b_i^e_i, where all bases b_i are distinct, and all exponents e_i are distinct >=1.

Original entry on oeis.org

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

Offset: 1

Saverio Picozzi, Feb 18 2015

nonn

Not multiplicative: a(48) = a(2^4*3) = 5 <> a(2^4)*a(3) = 4*1 = 4. - R. J. Mathar, Nov 05 2016

			From _R. J. Mathar_, Nov 05 2016: (Start)
a(4)=2: 4^1 = 2^2.
a(8)=2: 8^1 = 2^3.
a(9)=2: 9^1 = 3^2.
a(12)=2: 12^1 = 2^2*3^1.
a(16)=4: 16^1 = 4^2 = 2^2*4^1 = 2^4.
a(18)=2: 18^1 = 2*3^2.
a(20)=2: 20^1 = 2^2*5^1.
a(24)=3: 24^1 = 2^2*6^1 = 2^3*3^1.
a(32)=5: 32^1 = 2^1*4^2 = 2^2*8^1 = 2^3*4^1 = 2^5.
a(36)=4: 36^1 = 6^2 = 3^2*4^1 = 2^2*9^1.
a(48)=5: 48^1 = 3^1*4^2 = 2^2*12^1 = 2^3*6^1 = 2^4*3^1.
a(60)=2 : 60^1 = 2^2*15^1.
a(64)=7: 64^1 = 8^2 = 4^3 = 2^2*16^1 = 2^3*8^1 = 2^4*4^1 = 2^6.
a(72)=6 : 72^1 = 3^2*8^1 = 2^1*6^2 = 2^2*18^1 = 2^3*9^1 = 2^3*3^2.
(End)

R. J. Mathar, Table of n, a(n) for n = 1..419
R. J. Mathar, Factorizations of integers into factors with distinct bases and exponents
Index to sequences related to prime signature

Cf. A000688 (b_i not necessarily distinct).

Cf. A001248, A005117, A030078, A030514, A054753, A065036, A085986, A085987, A143610, A178739.

# Count solutions for products if n = dvs_i^exps(i) where i=1..pividx are fixed
Apiv := proc(n,dvs,exps,pividx)
    local dvscnt, expscopy,i,a,expsrt,e ;
    dvscnt := nops(dvs) ;
    a := 0 ;
    if pividx > dvscnt then
        # have exhausted the exponent list: leave of the recursion
        # check that dvs_i^exps(i) is a representation
        if n = mul( op(i,dvs)^op(i,exps),i=1..dvscnt) then
            # construct list of non-0 exponents
            expsrt := [];
            for i from 1 to dvscnt do
                if op(i,exps) > 0 then
                    expsrt := [op(expsrt),op(i,exps)] ;
                end if;
            end do;
            # check that list is duplicate-free
            if nops(expsrt) = nops( convert(expsrt,set)) then
                return 1;
            else
                return 0;
            end if;
        else
            return 0 ;
        end if;
    end if;
    # need a local copy of the list to modify it
    expscopy := [] ;
    for i from 1 to nops(exps) do
        expscopy := [op(expscopy),op(i,exps)] ;
    end do:
    # loop over all exponents assigned to the next base in the list.
    for e from 0 do
        candf := op(pividx,dvs)^e ;
        if modp(n,candf) <> 0 then
            break;
        end if;
        # assign e to the local copy of exponents
        expscopy := subsop(pividx=e,expscopy) ;
        a := a+procname(n,dvs,expscopy,pividx+1) ;
    end do:
    return a;
end proc:
A255231 := proc(n)
    local dvs,dvscnt,exps ;
    if n = 1 then
        return 1;
    end if;
    # candidates for the bases are all divisors except 1
    dvs := convert(numtheory[divisors](n) minus {1},list) ;
    dvscnt := nops(dvs) ;
    # list of exponents starts at all-0 and is
    # increased recursively
    exps := [seq(0,e=1..dvscnt)] ;
    # take any subset of dvs for the bases, i.e. exponents 0 upwards
    Apiv(n,dvs,exps,1) ;
end proc:
seq(A255231(n),n=1..120) ; # R. J. Mathar, Nov 05 2016

a(n)=1 for all n in A005117. a(n)=2 for all n in A001248 and for all n in A054753 and for all n in A085987 and for all n in A030078. a(n)=3 for all n in A065036. a(n)=4 for all n in A085986 and for all n in A030514. a(n)=5 for all n in A178739, all n in A179644 and for all n in A050997. a(n)=6 for all n in A143610, all n in A162142 and all n in A178740. a(n)=7 for all n in A030516. a(n)=9 for all n in A189988 and all n in A189987. a(n)=10 for all n in A092759. a(n) = 11 for all n in A179664. a(n)=12 for all n in A179646. - R. J. Mathar, Nov 05 2016, May 20 2017

Values corrected. Incorrect comments removed. - R. J. Mathar, Nov 05 2016

A343443 If n = Product (p_j^k_j) then a(n) = Product (k_j + 2), with a(1) = 1.

Original entry on oeis.org

1, 3, 3, 4, 3, 9, 3, 5, 4, 9, 3, 12, 3, 9, 9, 6, 3, 12, 3, 12, 9, 9, 3, 15, 4, 9, 5, 12, 3, 27, 3, 7, 9, 9, 9, 16, 3, 9, 9, 15, 3, 27, 3, 12, 12, 9, 3, 18, 4, 12, 9, 12, 3, 15, 9, 15, 9, 9, 3, 36, 3, 9, 12, 8, 9, 27, 3, 12, 9, 27, 3, 20, 3, 9, 12, 12, 9, 27, 3, 18

Offset: 1

Ilya Gutkovskiy, Apr 15 2021

Inverse Moebius transform of A056671.

a(n) depends only on the prime signature of n (see formulas). - Bernard Schott, May 03 2021

Cf. A000005, A001221, A005361, A007425, A025847, A034444, A056671, A064549, A074816, A107758, A107759, A348018, A363194.

a[1] = 1; a[n_] := Times @@ ((#[[2]] + 2) & /@ FactorInteger[n]); Table[a[n], {n, 80}]
a[n_] := Sum[If[GCD[d, n/d] == 1, DivisorSigma[0, d], 0], {d, Divisors[n]}]; Table[a[n], {n, 80}]

a(n) = sumdiv(n, d, if(gcd(d, n/d)==1, numdiv(d))) \\ Andrew Howroyd, Apr 15 2021

for(n=1, 100, print1(direuler(p=2, n, (1 + X - X^2)/(1-X)^2)[n], ", ")) \\ Vaclav Kotesovec, Feb 11 2023

from math import prod
from sympy import factorint
def A343443(n): return prod(e+2 for e in factorint(n).values()) # Chai Wah Wu, Feb 21 2025

a(n) = 2^omega(n) * tau_3(n) / tau(n), where omega = A001221, tau = A000005 and tau_3 = A007425.
a(n) = Sum_{d|n, gcd(d, n/d) = 1} tau(d).
From Bernard Schott, May 03 2021: (Start)
a(p^k) = k+2 for p prime, or signature [k].
a(A006881(n)) =  9 for signature [1, 1].
a(A054753(n)) = 12 for signature [2, 1].
a(A065036(n)) = 15 for signature [3, 1].
a(A085986(n)) = 16 for signature [2, 2].
a(A178739(n)) = 18 for signature [4, 1].
a(A143610(n)) = 20 for signature [3, 2].
a(A007304(n)) = 27 for signature [1, 1, 1]. (End)
Dirichlet g.f.: zeta(s)^2 * Product_{primes p} (1 + 1/p^s - 1/p^(2*s)). - Vaclav Kotesovec, Feb 11 2023
From Amiram Eldar, Sep 01 2023: (Start)
a(n) = A000005(A064549(n)).
a(n) = A363194(A348018(n)). (End)

A179689 Numbers with prime signature {7,2}, i.e., of form p^7*q^2 with p and q distinct primes.

Original entry on oeis.org

1152, 3200, 6272, 8748, 15488, 21632, 36992, 46208, 54675, 67712, 107163, 107648, 123008, 175232, 215168, 236672, 264627, 282752, 312500, 359552, 369603, 445568, 476288, 574592, 632043, 645248, 682112, 703125, 789507, 798848, 881792, 1013888

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
OEIS Wiki, Numbers with same prime signature.
Will Nicholes, Prime Signatures

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688.

Cf. A085548, A085967, A085969.

a:= proc(n) option remember; local k;
      for k from 1+ `if` (n=1, 1, a(n-1))
        while sort (map (x-> x[2], ifactors(k)[2]), `>`)<>[7, 2]
      do od; k
    end:
seq (a(n), n=1..32);  # Alois P. Heinz, Jan 23 2011

f[n_]:=Sort[Last/@FactorInteger[n]]=={2,7}; Select[Range[10^6], f]

list(lim)=my(v=List(),t);forprime(p=2, (lim\4)^(1/7), t=p^7;forprime(q=2, sqrt(lim\t), if(p==q, next);listput(v,t*q^2))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A179689(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(isqrt(x//p**7)) for p in primerange(integer_nthroot(x,7)[0]+1))+primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

Sum_{n>=1} 1/a(n) = P(2)*P(7) - P(9) = A085548 * A085967 - A085969 = 0.001741..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

Title edited by Daniel Forgues, Jan 22 2011

A179696 Numbers with prime signature {7,1,1}, i.e., of form p^7qr with p, q and r distinct primes.

Original entry on oeis.org

1920, 2688, 4224, 4480, 4992, 6528, 7040, 7296, 8320, 8832, 9856, 10880, 11136, 11648, 11904, 12160, 14208, 14720, 15232, 15744, 16512, 17024, 18048, 18304, 18560, 19840, 20352, 20608, 21870, 22656, 23424, 23680, 23936, 25728, 25984, 26240, 26752, 27264

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes List of Prime Signatures
OEIS Wiki, Numbers with same prime signature.

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688, A179689, A179690, A179691, A179692, A179693, A179694, A179695.

a:= proc(n) option remember; local k;
      for k from 1+ `if` (n=1, 1, a(n-1))
        while sort (map (x-> x[2], ifactors(k)[2]), `>`)<>[7, 1, 1]
      do od; k
    end:
seq (a(n), n=1..40); # Alois P. Heinz, Jan 23 2011

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,1,7}; Select[Range[30000], f]

list(lim)=my(v=List(),t1,t2);forprime(p=2, (lim\6)^(1/7), t1=p^7;forprime(q=2, lim\t1, if(p==q, next);t2=t1*q;forprime(r=q+1, lim\t2, if(p==r,next);listput(v,t2*r)))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from math import isqrt
from sympy import primerange, primepi, integer_nthroot
def A179696(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((t:=primepi(s:=isqrt(y:=x//r**7)))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(integer_nthroot(x,7)[0]+1))+sum(primepi(x//p**8) for p in primerange(integer_nthroot(x,8)[0]+1))-primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

Title edited by Daniel Forgues, Jan 22 2011

cp's OEIS Frontend

A195086 Numbers k such that (number of prime factors of k counted with multiplicity) less (number of distinct prime factors of k) = 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Formula

A179668 Products of the 8th power of a prime and a distinct prime (p^8*q).

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

A179692 Numbers of the form p^9*q where p and q are distinct primes.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

A376936 Powerful numbers divisible by cubes of 2 distinct primes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A074985 Squares of semiprimes (A001358).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

Formula

A179642 Product of exactly 5 primes, 3 of which are distinct.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

A255231 The number of factorizations n = Product_i b_i^e_i, where all bases b_i are distinct, and all exponents e_i are distinct >=1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

A179696 Numbers with prime signature {7,1,1}, i.e., of form p^7qr with p, q and r distinct primes.