A006881 -id:A006881 - cp's OEIS frontend

A010553 a(n) = tau(tau(n)).

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

Offset: 1

N. J. A. Sloane

Ramanujan (1915) posed the problem of finding the extreme large values of a(n). Buttkewitz et al. determined the maximal order of log a(n).

Every number eventually appears. Sequence A193987 gives the least term where each number appears. - T. D. Noe, Aug 10 2011

S. Ramanujan, Highly composite numbers. Proc. London Math. Soc., series 2, 14 (1915), 347-409. Republished in Collected papers of Srinivasa Ramanujan, AMS Chelsea Publ., Providence, RI, 2000, pp. 78-128.

Indranil Ghosh, Table of n, a(n) for n = 1..10000 (first 2000 terms from Enrique Pérez Herrero)
M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards, Applied Math. Series 55, Tenth Printing, 1972.
R. Bellman and H. N. Shapiro, On a problem in additive number theory, Annals Math., 49 (1948), 333-340. See Eq. 1.5.
Yvonne Buttkewitz, Christian Elsholtz, Kevin Ford, Jan-Christoph Schlage-Puchta, A problem of Ramanujan, Erdos and Katai on the iterated divisor function, arXiv:1108.1815 [math.NT], Aug 08 2011.

Cf. A000005, A036450, A193987 (least number k such that tau(tau(k)) = n), A335831.

Maple
```
with(numtheory): f := n->tau(tau(n));
```

Table[Nest[DivisorSigma[0, #] &, n, 2], {n, 81}] (* Michael De Vlieger, Dec 24 2015 *)

A010553(n)=numdiv(numdiv(n)); \\ Enrique Pérez Herrero, Jul 13 2010

a(n) = A000005(A000005(n)). a(1) = 1, a(p) = 2 for p = primes (A000040), a(pq) = 3 for pq = product of two distinct primes (A006881), a(pq...z) = k + 1 for pq...z = product of k (k > 2) distinct primes p,q,...,z (A120944), a(p^k) = A000005(k+1) for p^k = prime powers (A000961(n) for n > 1), k = natural numbers (A000027). - Jaroslav Krizek, Jul 17 2009
a(A007947(n)) = 1 + A001221(n); (n>1). - Enrique Pérez Herrero, May 30 2010
Asymptotically, Max_{i<=n} log(tau(tau(i))) = sqrt(log(n))/log_2(n) * (c + O(log_3(n)/log_2(n)) where c = 8*Sum_{j>=1} log^2 (1 + 1/j)) ~ 2.7959802335... [Buttkewitz et al.].

A007422 Multiplicatively perfect numbers j: product of divisors of j is j^2.

Original entry on oeis.org

1, 6, 8, 10, 14, 15, 21, 22, 26, 27, 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, 125, 129, 133, 134, 141, 142, 143, 145, 146, 155, 158, 159, 161, 166, 177, 178, 183, 185, 187

Offset: 1

N. J. A. Sloane

Or, numbers j such that product of proper divisors of j is j.
If M(j) denotes the product of the divisors of j, then j is said to be k-multiplicatively perfect if M(j) = j^k. All such numbers are of the form p q^(k-1) or p^(2k-1). This statement is in Sandor's paper. Therefore all 2-multiplicatively perfect numbers are semiprime p*q or cubes p^3. - Walter Kehowski, Sep 13 2005
All 2-multiplicatively perfect numbers except 1 have 4 divisors (as implied by Kehowski) and the converse is also true that all numbers with 4 divisors are 2-multiplicatively perfect. - Howard Berman (howard_berman(AT)hotmail.com), Oct 24 2008
Also 1 followed by numbers j such that A000005(j) = 4. - Nathaniel Johnston, May 03 2011
Fixed points of A007956. - Reinhard Zumkeller, Jan 26 2014

			The divisors of 10 are 1, 2, 5, 10 and 1 * 2 * 5 * 10 = 100 = 10^2.

Kenneth Ireland and Michael Ira Rosen, A Classical Introduction to Modern Number Theory. Springer-Verlag, NY, 1982, p. 19.
Edmund Landau, Elementary Number Theory, translation by Jacob E. Goodman of Elementare Zahlentheorie (Vol. I_1 (1927) of Vorlesungen ueber Zahlentheorie), by Edmund Landau, with added exercises by Paul T. Bateman and E. E. Kohlbecker, Chelsea Publishing Co., New York, 1958, pp. 31-32.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

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.
József Sándor, Multiplicatively perfect numbers, J. Ineq. Pure Appl. Math., Vol. 2, No. 1 (2001), Article 3, 6 pp.
Eric Weisstein's World of Mathematics, Divisor Product.
Eric Weisstein's World of Mathematics, Multiplicative Perfect Number.

Cf. A030513 (same as this sequence but without the 1), A027751, A006881 (subsequence), A030078 (subsequence), A084110, A084116, A236473.

a007422 n = a007422_list !! (n-1)
a007422_list = [x | x <- [1..], a007956 x == x]
-- Reinhard Zumkeller, Jan 26 2014

IsA007422:=func< n | &*Divisors(n) eq n^2 >; [ n: n in [1..200] | IsA007422(n) ]; // Klaus Brockhaus, May 04 2011

k:=2: MPL:=[]: for z from 1 to 1 do for n from 1 to 5000 do if convert(divisors(n),`*`) = n^k then MPL:=[op(MPL),n] fi od; od; MPL; # Walter Kehowski, Sep 13 2005
# second Maple program:
q:= n-> n=1 or numtheory[tau](n)=4:
select(q, [$1..200])[];  # Alois P. Heinz, Dec 17 2021

Select[Range[200], Times@@Divisors[#] == #^2 &]  (* Harvey P. Dale, Mar 27 2011 *)

is(n)=n==1 || numdiv(n) == 4 \\ Charles R Greathouse IV, Oct 15 2015

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A007422(n):
    def f(x): return int(n-1+x-primepi(integer_nthroot(x,3)[0])+(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 16 2024

A084110(a(n)) = 1, see also A084116. - Reinhard Zumkeller, May 12 2003

The number of terms not exceeding x is N(x) ~ x * log(log(x))/log(x) (Chau, 2004). - Amiram Eldar, Jun 29 2022

Some numbers were omitted - thanks to Erich Friedman for pointing this out.

A119899 Integers i such that bigomega(i) (A001222) and tau(i) (A000005) are both even.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 24, 26, 33, 34, 35, 38, 39, 40, 46, 51, 54, 55, 56, 57, 58, 60, 62, 65, 69, 74, 77, 82, 84, 85, 86, 87, 88, 90, 91, 93, 94, 95, 96, 104, 106, 111, 115, 118, 119, 122, 123, 126, 129, 132, 133, 134, 135, 136, 140, 141, 142, 143, 145, 146, 150

Offset: 1

Antti Karttunen, Jun 04 2006

nonn

Also numbers whose alternating sum of prime indices is < 0. Equivalently, numbers with even bigomega whose conjugate prime indices are not all even. This is the intersection of A028260 and A000037. - Gus Wiseman, Jun 20 2021

			From _Gus Wiseman_, Jun 20 2021: (Start)
The sequence of terms together with their prime indices begins:
       6: {1,2}          51: {2,7}          86: {1,14}
      10: {1,3}          54: {1,2,2,2}      87: {2,10}
      14: {1,4}          55: {3,5}          88: {1,1,1,5}
      15: {2,3}          56: {1,1,1,4}      90: {1,2,2,3}
      21: {2,4}          57: {2,8}          91: {4,6}
      22: {1,5}          58: {1,10}         93: {2,11}
      24: {1,1,1,2}      60: {1,1,2,3}      94: {1,15}
      26: {1,6}          62: {1,11}         95: {3,8}
      33: {2,5}          65: {3,6}          96: {1,1,1,1,1,2}
      34: {1,7}          69: {2,9}         104: {1,1,1,6}
      35: {3,4}          74: {1,12}        106: {1,16}
      38: {1,8}          77: {4,5}         111: {2,12}
      39: {2,6}          82: {1,13}        115: {3,9}
      40: {1,1,1,3}      84: {1,1,2,4}     118: {1,17}
      46: {1,9}          85: {3,7}         119: {4,7}
(End)

Superset: A119847. Subset: A006881. The intersection of A028260 and A000037.
Positions of negative terms in A316524.
The partitions with these Heinz numbers are counted by A344608.
Complement of A344609.
Cf. A000041, A000070, A000097, A027187, A103919, A116406, A239829, A239830, A343941, A344607, A344651.

Select[Range[200],And@@EvenQ[{PrimeOmega[#],DivisorSigma[0,#]}]&] (* Harvey P. Dale, Jan 24 2013 *)

A176504 a(n) = m + k where prime(m)*prime(k) = semiprime(n).

Original entry on oeis.org

2, 3, 4, 4, 5, 5, 6, 6, 6, 7, 7, 8, 7, 9, 8, 10, 8, 9, 8, 10, 11, 12, 9, 11, 13, 9, 14, 10, 15, 12, 10, 13, 16, 11, 17, 14, 12, 18, 11, 10, 19, 15, 16, 12, 20, 17, 21, 11, 13, 22, 14, 23, 18, 13, 24, 12, 19, 25, 20, 15, 12, 26, 21, 27, 14, 16, 28, 13, 22, 29, 17, 15, 30, 23, 13, 31

Offset: 1

Juri-Stepan Gerasimov, Apr 19 2010

nonn

			From _Gus Wiseman_, Dec 04 2020: (Start)
A semiprime (A001358) is a product of any two prime numbers. The sequence of all semiprimes together with their prime indices and weights begins:
   4: 1 + 1 = 2
   6: 1 + 2 = 3
   9: 2 + 2 = 4
  10: 1 + 3 = 4
  14: 1 + 4 = 5
  15: 2 + 3 = 5
  21: 2 + 4 = 6
  22: 1 + 5 = 6
  25: 3 + 3 = 6
  26: 1 + 6 = 7
(End)

A056239 is the version for not just semiprimes.
A087794 gives the product of the same two indices.
A176506 gives the difference of the same two indices.
A338904 puts the n-th semiprime in row a(n).
A001358 lists semiprimes.
A006881 lists squarefree semiprimes.
A338898/A338912/A338913 give the prime indices of semiprimes.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes, with product/sum/difference A339361/A339362/A338900.
Cf. A001222, A046315, A065516, A084126, A084127, A100484, A112798, A115392, A128301, A338906/A338907.

From R. J. Mathar, Apr 20 2010: (Start)
isA001358 := proc(n) numtheory[bigomega](n) = 2 ; end proc:
A001358 := proc(n) option remember ; if n = 1 then return 4 ; else for a from procname(n-1)+1 do if isA001358(a) then return a; end if; end do; end if; end proc:
A084126 := proc(n) min(op(numtheory[factorset](A001358(n)))) ; end proc:
A084127 := proc(n) max(op(numtheory[factorset](A001358(n)))) ; end proc:
A176504 := proc(n) numtheory[pi](A084126(n)) + numtheory[pi](A084127(n)) ; end proc: seq(A176504(n),n=1..80) ; (End)

Table[If[SquareFreeQ[n],Total[PrimePi/@First/@FactorInteger[n]],2*PrimePi[Sqrt[n]]],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

a(n) = A056239(A001358(n)) = A338912(n) + A338913(n). - Gus Wiseman, Dec 04 2020

sqrt(n/(log n log log n)) << a(n) << n/log log n. - Charles R Greathouse IV, Apr 17 2024

Entries checked by R. J. Mathar, Apr 20 2010

A087794 Products of prime-indices of factors of semiprimes.

Original entry on oeis.org

1, 2, 4, 3, 4, 6, 8, 5, 9, 6, 10, 7, 12, 8, 12, 9, 16, 14, 15, 16, 10, 11, 18, 18, 12, 20, 13, 21, 14, 20, 24, 22, 15, 24, 16, 24, 27, 17, 28, 25, 18, 26, 28, 32, 19, 30, 20, 30, 30, 21, 33, 22, 32, 36, 23, 36, 34, 24, 36, 36, 35, 25, 38, 26, 40, 39, 27, 40, 40, 28, 42, 44, 29

Offset: 1

Reinhard Zumkeller, Oct 09 2003

nonn

A semiprime (A001358) is a product of any two prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798. - Gus Wiseman, Dec 04 2020

			A001358(20)=57=3*19=A000040(2)*A000040(8), therefore a(20)=2*8=16.
From _Gus Wiseman_, Dec 04 2020: (Start)
The sequence of all semiprimes together with the products of their prime indices begins:
   4: 1 * 1 = 1
   6: 1 * 2 = 2
   9: 2 * 2 = 4
  10: 1 * 3 = 3
  14: 1 * 4 = 4
  15: 2 * 3 = 6
  21: 2 * 4 = 8
  22: 1 * 5 = 5
  25: 3 * 3 = 9
  26: 1 * 6 = 6
(End)

A003963 is the version for not just semiprimes.
A176504 gives the sum of the same two indices.
A176506 gives the difference of the same two indices.
A339361 is the squarefree case.
A001358 lists semiprimes.
A006881 lists squarefree semiprimes.
A289182/A115392 list the positions of odd/even terms of A001358.
A338898/A338912/A338913 give the prime indices of semiprimes.
A338899/A270650/A270652 give the prime indices of squarefree semiprimes.
A338904 groups semiprimes by weight.
Cf. A001222, A046315, A056239, A065516, A084126, A084127, A112798, A128301, A300912, A318990, A320655, A338909, A339362.

Table[If[SquareFreeQ[n],Times@@PrimePi/@First/@FactorInteger[n],PrimePi[Sqrt[n]]^2],{n,Select[Range[100],PrimeOmega[#]==2&]}] (* Gus Wiseman, Dec 04 2020 *)

a(n) = A049084(A084126(n))*A049084(A084127(n)) = A003963(A001358(n)).

a(n) = A003963(A001358(n)) = A338912(n) * A338913(n). - Gus Wiseman, Dec 04 2020

A320663 Number of non-isomorphic multiset partitions of weight n using singletons or pairs.

Original entry on oeis.org

1, 1, 4, 7, 21, 40, 106, 216, 534, 1139, 2715, 5962, 14012, 31420, 73484, 167617, 392714, 908600, 2140429, 5015655, 11905145, 28228533, 67590229, 162067916, 391695348, 949359190, 2316618809, 5673557284, 13979155798, 34583650498, 86034613145, 214948212879

Offset: 0

Gus Wiseman, Oct 18 2018

nonn

			Non-isomorphic representatives of the a(1) = 1 through a(4) = 21 multiset partitions:
  {{1}}  {{1,1}}    {{1},{1,1}}    {{1,1},{1,1}}
         {{1,2}}    {{1},{2,2}}    {{1,1},{2,2}}
         {{1},{1}}  {{1},{2,3}}    {{1,2},{1,2}}
         {{1},{2}}  {{2},{1,2}}    {{1,2},{2,2}}
                    {{1},{1},{1}}  {{1,2},{3,3}}
                    {{1},{2},{2}}  {{1,2},{3,4}}
                    {{1},{2},{3}}  {{1,3},{2,3}}
                                   {{1},{1},{1,1}}
                                   {{1},{1},{2,2}}
                                   {{1},{1},{2,3}}
                                   {{1},{2},{1,2}}
                                   {{1},{2},{2,2}}
                                   {{1},{2},{3,3}}
                                   {{1},{2},{3,4}}
                                   {{1},{3},{2,3}}
                                   {{2},{2},{1,2}}
                                   {{1},{1},{1},{1}}
                                   {{1},{1},{2},{2}}
                                   {{1},{2},{2},{2}}
                                   {{1},{2},{3},{3}}
                                   {{1},{2},{3},{4}}

Andrew Howroyd, Table of n, a(n) for n = 0..50

Cf. A001055, A001222, A001358, A005117, A006881, A007716, A007717, A037143, A320462, A320655, A320656, A320664, A320665.

EulerT(v)={Vec(exp(x*Ser(dirmul(v,vector(#v,n,1/n))))-1, -#v)}
permcount(v) = {my(m=1,s=0,k=0,t); for(i=1,#v,t=v[i]; k=if(i>1&&t==v[i-1],k+1,1); m*=t*k;s+=t); s!/m}
gs(v) = {sum(i=2, #v, sum(j=1, i-1, my(g=gcd(v[i],v[j])); g*x^(2*v[i]*v[j]/g))) + sum(i=1, #v, my(r=v[i]); (1 + (1+r)%2)*x^r + ((1+r)\2)*x^(2*r))}
a(n)={my(s=0); forpart(p=n, s+=permcount(p)*EulerT(Vec(gs(p) + O(x*x^n), -n))[n]); s/n!} \\ Andrew Howroyd, Oct 26 2018

Terms a(11) and beyond from Andrew Howroyd, Oct 26 2018

A320891 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into squarefree semiprimes.

Original entry on oeis.org

4, 9, 16, 24, 25, 40, 49, 54, 56, 64, 81, 88, 96, 104, 121, 135, 136, 144, 152, 160, 169, 184, 189, 224, 232, 240, 248, 250, 256, 289, 296, 297, 324, 328, 336, 344, 351, 352, 361, 375, 376, 384, 400, 416, 424, 459, 472, 486, 488, 513, 528, 529, 536, 544, 560

Offset: 1

Gus Wiseman, Oct 23 2018

nonn

A squarefree semiprime (A006881) is a product of any two distinct primes.

Also numbers with an even number x of prime factors, whose greatest prime multiplicity exceeds x/2.

			A complete list of all factorizations of 24 is:
  (2*2*2*3),
  (2*2*6), (2*3*4),
  (2*12), (3*8), (4*6),
  (24).
All of these contain at least one number that is not a squarefree semiprime, so 24 belongs to the sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001055, A001358, A005117, A006881, A007717, A028260, A318871, A318953, A320655, A320656, A320892, A320893, A320894.

semfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[semfacs[n/d],Min@@#>=d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
Select[Range[100],And[EvenQ[PrimeOmega[#]],semfacs[#]=={}]&]

A320894 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into distinct squarefree semiprimes.

Original entry on oeis.org

4, 9, 16, 24, 25, 36, 40, 49, 54, 56, 64, 81, 88, 96, 100, 104, 121, 135, 136, 144, 152, 160, 169, 184, 189, 196, 216, 224, 225, 232, 240, 248, 250, 256, 289, 296, 297, 324, 328, 336, 344, 351, 352, 360, 361, 375, 376, 384, 400, 416, 424, 441, 459, 472, 484

Offset: 1

Gus Wiseman, Oct 23 2018

nonn

A squarefree semiprime (A006881) is a product of any two distinct primes.

			A complete list of all strict factorizations of 24 is: (2*3*4), (2*12), (3*8), (4*6), (24). All of these contain at least one number that is not a squarefree semiprime, so 24 belongs to the sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001055, A001358, A005117, A006881, A007717, A028260, A318871, A318953, A320655, A320656, A320891, A320892, A320893.

strsqfsemfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[strsqfsemfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],And[SquareFreeQ[#],PrimeOmega[#]==2]&]}]];
Select[Range[100],And[EvenQ[PrimeOmega[#]],strsqfsemfacs[#]=={}]&]

A338900 Difference between the two prime indices of the n-th squarefree semiprime.

Original entry on oeis.org

1, 2, 3, 1, 2, 4, 5, 3, 6, 1, 7, 4, 8, 5, 2, 6, 9, 10, 3, 7, 11, 1, 12, 4, 13, 8, 2, 9, 14, 5, 15, 10, 6, 16, 3, 17, 11, 12, 4, 18, 13, 19, 1, 7, 20, 8, 21, 14, 5, 22, 15, 23, 16, 9, 2, 24, 17, 25, 6, 10, 26, 3, 18, 27, 11, 7, 28, 19, 1, 29, 12, 20, 2, 21, 4

Offset: 1

Gus Wiseman, Nov 16 2020

nonn

A squarefree semiprime is a product of any two distinct prime numbers. A prime index of n is a number m such that the m-th prime number divides n. The multiset of prime indices of n is row n of A112798.

Is this sequence an anti-run, i.e., are there no adjacent equal parts? I have verified this conjecture up to n = 10^6. - Gus Wiseman, Nov 18 2020

A176506 is the not necessarily squarefree version.
A338899 has row-differences equal to this sequence.
A338901 gives positions of first appearances.
A001221 counts distinct prime indices.
A001222 counts prime indices.
A001358 lists semiprimes.
A002100 and A338903 count partitions using squarefree semiprimes.
A004526 counts 2-part partitions, with strict case A140106 (shifted left).
A005117 lists squarefree numbers.
A006881 lists squarefree semiprimes, with odds A046388 and evens A100484.
A065516 gives first differences of semiprimes.
A166237 gives first differences of squarefree semiprimes.
A270650 and A270652 give the prime indices of squarefree semiprimes.
A338912 and A338913 give the prime indices of semiprimes.
Cf. A000040, A056239, A112798, A167171, A320656, A320891, A320894, A320911, A338898, A338905, A338908.

-Subtract@@PrimePi/@First/@FactorInteger[#]&/@Select[Range[100],SquareFreeQ[#]&&PrimeOmega[#]==2&]

If the n-th squarefree semiprime is prime(x) * prime(y) with x < y, then a(n) = y - x.

a(n) = A270652(n) - A270650(n).

A002100 a(n) = number of partitions of n into semiprimes (more precisely, number of ways of writing n as a sum of products of 2 distinct primes).

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 1, 0, 1, 0, 2, 2, 2, 0, 2, 1, 3, 2, 3, 1, 4, 2, 4, 3, 5, 4, 7, 3, 6, 5, 8, 6, 10, 6, 10, 9, 12, 9, 15, 11, 16, 14, 18, 14, 22, 19, 25, 22, 27, 23, 33, 29, 36, 33, 40, 38, 49, 43, 53, 51, 61, 57, 71, 64, 77, 76, 89, 86, 102, 96, 113, 111, 128, 125

Offset: 1

N. J. A. Sloane

nonn

			a(20) = 2: 20 = 2*3 + 2*7 = 2*5 + 2*5.

L. M. Chawla and S. A. Shad, On a restricted partition function t(n) and its table, J. Natural Sciences and Mathematics, 9 (1969), 217-221. Math. Rev. 41 #6761.
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, Table of n, a(n) for n = 1..1000

Cf. A006881, A073576, A101048.

a002100 = p a006881_list where
   p _          0 = 1
   p ks'@(k:ks) m = if m < k then 0 else p ks' (m - k) + p ks m
-- Reinhard Zumkeller, Mar 21 2014

a[n_] := SeriesCoefficient[1/Product[If[SquareFreeQ[k] && PrimeNu[k] == 2, 1 - z^k, 1], {k, 1, n}], {z, 0, n}];
Array[a, 100] (* Jean-François Alcover, Nov 26 2020, after PARI *)

a(n)=polcoeff(1/prod(k=1,n,if(issquarefree(k)*if(omega(k)-2,0,1),1-z^k,1))+O(z^(n+1)),n)

More terms from Benoit Cloitre, Jun 01 2003

cp's OEIS Frontend

A010553 a(n) = tau(tau(n)).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A007422 Multiplicatively perfect numbers j: product of divisors of j is j^2.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

Python

Formula

Extensions

A119899 Integers i such that bigomega(i) (A001222) and tau(i) (A000005) are both even.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A176504 a(n) = m + k where prime(m)*prime(k) = semiprime(n).

Views

Author

Keywords

Examples

Crossrefs

Programs

Maple

Mathematica

Formula

Extensions

A087794 Products of prime-indices of factors of semiprimes.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A320663 Number of non-isomorphic multiset partitions of weight n using singletons or pairs.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

PARI

Extensions

A320891 Numbers with an even number of prime factors (counted with multiplicity) that cannot be factored into squarefree semiprimes.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs