A048691 -id:A048691 - cp's OEIS frontend

A048785 a(0) = 0; a(n) = tau(n^3), where tau = number of divisors (A000005).

0, 1, 4, 4, 7, 4, 16, 4, 10, 7, 16, 4, 28, 4, 16, 16, 13, 4, 28, 4, 28, 16, 16, 4, 40, 7, 16, 10, 28, 4, 64, 4, 16, 16, 16, 16, 49, 4, 16, 16, 40, 4, 64, 4, 28, 28, 16, 4, 52, 7, 28, 16, 28, 4, 40, 16, 40, 16, 16, 4, 112, 4, 16, 28, 19, 16, 64, 4, 28, 16

Offset: 0

N. J. A. Sloane

The inverse Mobius transform of A074816. - R. J. Mathar, Feb 09 2011

a(n) is also the number of ordered triples (i,j,k) of positive integers such that i|n, j|n, k|n and i,j,k are pairwise relatively prime. - Geoffrey Critzer, Jan 11 2015

			a(6) = 16 because there are 16 divisors of 6^3 = 216: 1, 2, 3, 4, 6, 8, 9, 12, 18, 24, 27, 36, 54, 72, 108, 216.
Also there are 16 ordered triples of divisors of 6 that are pairwise relatively prime: (1,1,1), (1,1,2), (1,1,3), (1,1,6), (1,2,1), (1,2,3), (1,3,1), (1,3,2), (1,6,1), (2,1,1), (2,1,3), (2,3,1), (3,1,1), (3,1,2), (3,2,1), (6,1,1).

Antti Karttunen, Table of n, a(n) for n = 0..10000

Cf. A000005, A001221, A048691, A074816, A144943, A353551.

seq(numtheory:-tau(n^3), n=0..100); # Robert Israel, Jan 11 2015

Join[{0,1},Table[Product[3 k + 1, {k, FactorInteger[n][[All, 2]]}], {n, 2, 69}]] (* Geoffrey Critzer, Jan 11 2015 *)
Join[{0},DivisorSigma[0,Range[70]^3]] (* Harvey P. Dale, Jan 23 2016 *)

A048785(n) = if(!n,n,numdiv(n^3)); \\ Antti Karttunen, May 19 2017

print1("0, "); for(n=1, 100, print1(direuler(p=2, n, (1 + 2*X)/(1 - X)^2)[n], ", ")) \\ Vaclav Kotesovec, May 15 2021
print1("0, "); for(n=1, 100, print1(direuler(p=2, n, (1 - 3*X^2 + 2*X^3)/(1 - X)^4)[n], ", ")) \\ Vaclav Kotesovec, Aug 20 2021

from math import prod
from sympy import factorint
def A048785(n): return 0 if n == 0 else prod(3*e+1 for e in factorint(n).values()) # Chai Wah Wu, May 10 2022

from sympy import divisor_count
def A048785(n): return divisor_count(n**3) # Karl-Heinz Hofmann, May 10 2022

a(n) = Sum_{d|n} 3^omega(d), where omega(x) is the number of distinct prime factors in the factorization of x. - Benoit Cloitre, Apr 14 2002
Multiplicative with a(p^e) = 3e+1. - Mitch Harris, Jun 09 2005
L.g.f.: -log(Product_{k>=1} (1 - x^k)^(3^omega(k)/k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, May 26 2018
For n>0, a(n) = Sum_{d|n} mu(d)^2*tau(d)*tau(n/d). - Ridouane Oudra, Nov 18 2019
Dirichlet g.f.: zeta(s)^2 * Product_{primes p} (1 + 2/p^s). - Vaclav Kotesovec, May 15 2021
Dirichlet g.f.: zeta(s)^4 * Product_{primes p} (1 - 3/p^(2*s) + 2/p^(3*s)). - Vaclav Kotesovec, Aug 20 2021

A343656 Array read by antidiagonals where A(n,k) is the number of divisors of n^k.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 2, 1, 1, 4, 3, 3, 1, 1, 5, 4, 5, 2, 1, 1, 6, 5, 7, 3, 4, 1, 1, 7, 6, 9, 4, 9, 2, 1, 1, 8, 7, 11, 5, 16, 3, 4, 1, 1, 9, 8, 13, 6, 25, 4, 7, 3, 1, 1, 10, 9, 15, 7, 36, 5, 10, 5, 4, 1, 1, 11, 10, 17, 8, 49, 6, 13, 7, 9, 2, 1, 1, 12, 11, 19, 9, 64, 7, 16, 9, 16, 3, 6, 1

Offset: 1

Gus Wiseman, Apr 28 2021

First differs from A343658 at A(4,2) = 5, A343658(4,2) = 6.

As a triangle, T(n,k) = number of divisors of k^(n-k).

			Array begins:
       k=0 k=1 k=2 k=3 k=4 k=5 k=6 k=7
  n=1:  1   1   1   1   1   1   1   1
  n=2:  1   2   3   4   5   6   7   8
  n=3:  1   2   3   4   5   6   7   8
  n=4:  1   3   5   7   9  11  13  15
  n=5:  1   2   3   4   5   6   7   8
  n=6:  1   4   9  16  25  36  49  64
  n=7:  1   2   3   4   5   6   7   8
  n=8:  1   4   7  10  13  16  19  22
  n=9:  1   3   5   7   9  11  13  15
Triangle begins:
  1
  1  1
  1  2  1
  1  3  2  1
  1  4  3  3  1
  1  5  4  5  2  1
  1  6  5  7  3  4  1
  1  7  6  9  4  9  2  1
  1  8  7 11  5 16  3  4  1
  1  9  8 13  6 25  4  7  3  1
  1 10  9 15  7 36  5 10  5  4  1
  1 11 10 17  8 49  6 13  7  9  2  1
  1 12 11 19  9 64  7 16  9 16  3  6  1
  1 13 12 21 10 81  8 19 11 25  4 15  2  1
For example, row n = 8 counts the following divisors:
  1  64  243  256  125  36  7  1
     32  81   128  25   18  1
     16  27   64   5    12
     8   9    32   1    9
     4   3    16        6
     2   1    8         4
     1        4         3
              2         2
              1         1

Seiichi Manyama, Antidiagonals n = 1..140, flattened

Columns k=1..9 of the array give A000005, A048691, A048785, A344327, A344328, A344329, A343526, A344335, A344336.
Row n = 6 of the array is A000290.
Diagonal n = k of the array is A062319.
Array antidiagonal sums (row sums of the triangle) are A343657.
Dominated by A343658.
A000312 = n^n.
A007318 counts k-sets of elements of {1..n}.
A009998(n,k) = n^k (as an array, offset 1).
A059481 counts k-multisets of elements of {1..n}.
Cf. A000169, A000272, A002064, A002109, A066959, A143773, A146291, A176029, A251683, A282935, A326358, A327527, A334996, A343652.

Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]

A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

A(n,k) = A000005(A009998(n,k)), where A009998(n,k) = n^k is the interpretation as an array.

A(n,k) = Sum_{d|n} k^omega(d). - Seiichi Manyama, May 15 2021

A061283 Smallest number with exactly 2n-1 divisors.

Original entry on oeis.org

1, 4, 16, 64, 36, 1024, 4096, 144, 65536, 262144, 576, 4194304, 1296, 900, 268435456, 1073741824, 9216, 5184, 68719476736, 36864, 1099511627776, 4398046511104, 3600, 70368744177664, 46656, 589824, 4503599627370496, 82944

Offset: 1

Labos Elemer, May 22 2001

nonn

The terms are always squares (because the divisors of a nonsquare N come in pairs, d and N/d, and so their number is always even - N. J. A. Sloane, Dec 26 2018).

			For n=15, a(15)=144 with 15 divisors: 1,2,3,4,6,8,9,12,16,18,24,36,48,72 and 144.

Amiram Eldar, Table of n, a(n) for n = 1..1661 (terms 1..1000 from T. D. Noe)

Cf. A000005, A000290, A005408, A003680, A016017, A037992, A055079, A048691.

Second bisection of A005179.

mp[1, m_] := {{}}; mp[n_, 1] := {{}}; mp[n_?PrimeQ, m_] := If[m < n, {}, {{n}}]; mp[n_, m_] := Join @@ Table[Map[Prepend[#, d] &, mp[n/d, d]], {d, Select[Rest[Divisors[n]], # <= m &]}]; mp[n_] := mp[n, n]; Table[mulpar = mp[2*n-1] - 1; Min[Table[Product[Prime[s]^mulpar[[j, s]], {s, 1, Length[mulpar[[j]]]}], {j, 1, Length[mulpar]}]], {n, 1, 100}] (* Vaclav Kotesovec, Apr 04 2021 *)

a(n) = Min{k | A000005(k)=2n-1}.
a((p+1)/2) = 2^(p-1) for odd prime p. [Corrected by Jianing Song, Aug 30 2021]
From Jianing Song, Aug 30 2021: (Start)
a(n) = A016017(n)^2.
a(n) <= 2^(2n-2), where the equality holds if and only if n=1 or 2n-1 is prime. (End)

A126098 Where records occur in A018892.

Original entry on oeis.org

1, 2, 4, 6, 12, 24, 30, 60, 120, 180, 210, 360, 420, 840, 1260, 1680, 2520, 4620, 7560, 9240, 13860, 18480, 27720, 55440, 83160, 110880, 120120, 180180, 240240, 360360, 720720, 1081080, 1441440, 1801800, 2042040, 2882880, 3063060, 4084080, 5405400, 6126120, 12252240, 18378360, 24504480

Offset: 1

N. J. A. Sloane, Mar 05 2007

nonn

Remarkably similar to but ultimately different from A018894. - Jorg Brown and N. J. A. Sloane, Mar 06 2007
This sequence represents "where records occur" for a number of sequences in addition to A018892 including the following: A015995, A015996, A015999, A016001, A016002, A016003, A016005, A016006, A016007, A016008, A016009, A048691, A048785, A063647, A117677, A144943. - Ray Chandler, Dec 04 2008
Subsequence of A025487. - Ray Chandler, Sep 05 2008
Also record-setting elements of tau(n^2) (just as A002182 gives the record-setting elements of tau(n)). The point is that A018892 is (tau(n^2) + 1)/2. As tau(n^2) is odd, the record-setting elements of A018892 are also the record setting elements of tau(n^2). - Allen Tracht, Jan 20 2009

Ray Chandler, Table of n, a(n) for n = 1..582
Jorg Brown, Comparison of records in sigma(n)/phi(n) and A018892

Cf. A018892, A126097. Equals A117010(n) + 1.

More terms from Jorg Brown (jorg(AT)google.com) and T. D. Noe, Mar 05 2007

a(27) corrected by hupo001(AT)gmail.com, Jan 10 2008

A343652 Number of maximal pairwise coprime sets of divisors of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 4, 1, 4, 1, 4, 2, 2, 1, 6, 2, 2, 3, 4, 1, 5, 1, 5, 2, 2, 2, 8, 1, 2, 2, 6, 1, 5, 1, 4, 4, 2, 1, 8, 2, 4, 2, 4, 1, 6, 2, 6, 2, 2, 1, 10, 1, 2, 4, 6, 2, 5, 1, 4, 2, 5, 1, 12, 1, 2, 4, 4, 2, 5, 1, 8, 4, 2, 1, 10, 2, 2

Offset: 1

Gus Wiseman, Apr 25 2021

nonn

Also the number of maximal pairwise coprime sets of divisors > 1 of n. For example, the a(n) sets for n = 12, 30, 36, 60, 120 are:
  {6}    {30}     {6}    {30}     {30}
  {12}   {2,15}   {12}   {60}     {60}
  {2,3}  {3,10}   {18}   {2,15}   {120}
  {3,4}  {5,6}    {36}   {3,10}   {2,15}
         {2,3,5}  {2,3}  {3,20}   {3,10}
                  {2,9}  {4,15}   {3,20}
                  {3,4}  {5,6}    {3,40}
                  {4,9}  {5,12}   {4,15}
                         {2,3,5}  {5,6}
                         {3,4,5}  {5,12}
                                  {5,24}
                                  {8,15}
                                  {2,3,5}
                                  {3,4,5}
                                  {3,5,8}

			The a(n) sets for n = 12, 30, 36, 60, 120:
  {1,6}    {1,30}     {1,6}    {1,30}     {1,30}
  {1,12}   {1,2,15}   {1,12}   {1,60}     {1,60}
  {1,2,3}  {1,3,10}   {1,18}   {1,2,15}   {1,120}
  {1,3,4}  {1,5,6}    {1,36}   {1,3,10}   {1,2,15}
           {1,2,3,5}  {1,2,3}  {1,3,20}   {1,3,10}
                      {1,2,9}  {1,4,15}   {1,3,20}
                      {1,3,4}  {1,5,6}    {1,3,40}
                      {1,4,9}  {1,5,12}   {1,4,15}
                               {1,2,3,5}  {1,5,6}
                               {1,3,4,5}  {1,5,12}
                                          {1,5,24}
                                          {1,8,15}
                                          {1,2,3,5}
                                          {1,3,4,5}
                                          {1,3,5,8}

The case of pairs is A063647.
The case of triples is A066620.
The non-maximal version counting empty sets and singletons is A225520.
The non-maximal version with no 1's is A343653.
The non-maximal version is A343655.
The version for subsets of {1..n} is A343659.
The case without 1's or singletons is A343660.
A018892 counts pairwise coprime unordered pairs of divisors.
A048691 counts pairwise coprime ordered pairs of divisors.
A048785 counts pairwise coprime ordered triples of divisors.
A084422, A187106, A276187, and A320426 count pairwise coprime sets.
A100565 counts pairwise coprime unordered triples of divisors.
A305713 counts pairwise coprime non-singleton strict partitions.
A324837 counts minimal subsets of {1...n} with least common multiple n.
A325683 counts maximal Golomb rulers.
A326077 counts maximal pairwise indivisible sets.
Cf. A005361, A007359, A051026, A062319, A067824, A074206, A146291, A285572, A325859, A326359, A326496, A326675, A343654.

fasmax[y_]:=Complement[y,Union@@Most@*Subsets/@y];
Table[Length[fasmax[Select[Subsets[Divisors[n]],CoprimeQ@@#&]]],{n,100}]

a(n) = A343660(n) + A005361(n).

A016017 Smallest k such that 1/k can be written as a sum of exactly 2 unit fractions in n ways.

Original entry on oeis.org

1, 2, 4, 8, 6, 32, 64, 12, 256, 512, 24, 2048, 36, 30, 16384, 32768, 96, 72, 262144, 192, 1048576, 2097152, 60, 8388608, 216, 768, 67108864, 288, 1536, 536870912, 1073741824, 120, 576, 8589934592, 6144, 34359738368, 68719476736, 180, 864

Offset: 1

Robert G. Wilson v

nonn

From Jianing Song, Aug 30 2021: (Start)
a(n) is the smallest number whose square has exactly 2n-1 divisors.
a(n) is the earliest occurrence of 2n-1 in A048691. (End)

			a(1)=1 and a(2)=2 because 1/2 = 1/3 + 1/6 = 1/4 + 1/4.
a(3)=4 because 1/4 = 1/5 + 1/20 = 1/6 + 1/12 = 1/8 + 1/8.
a(4)=8 because 1/8 = 1/9 + 1/72 = 1/10 + 1/40 = 1/12 + 1/24 = 1/16 + 1/16.
a(5)=6 because 1/6 = 1/7 + 1/42 = 1/8 + 1/24 = 1/9 + 1/18 = 1/10 + 1/15 = 1/12 + 1/12.

David W. Wilson, Table of n, a(n) for n = 1..1000 (terms 122, 365 and 608 corrected by _Amiram Eldar_, Nov 08 2024)

Identical to A071571 shifted right.

Cf. A000005, A000290, A005408, A005179, A003680, A037992, A055079, A048691.

f[j_, n_] := (Times @@ (j(Last /@ FactorInteger[n]) + 1) + j - 1)/j; t = Table[0, {50}]; Do[a = f[2, n]; If[a < 51 && t[[a]] == 0, t[[a]] = n; Print[{a, n}]], {n, 2^30}] (* Robert G. Wilson v, Aug 03 2005 *)

a(n) = {k = 1; while (numdiv(k^2) != (2*n-1), k++); return (k); }; \\ Amiram Eldar, Jan 07 2019 after Michel Marcus at A071571

a(n+1) <= 2^n.
From Labos Elemer, May 22 2001: (Start)
a(n) = sqrt(A061283(n)).
a(n) = sqrt(Min{k| A000005(k)=2n-1}).
a((p+1)/2) = 2^((p-1)/2) = 2^A005097(i) if p is the i-th odd prime. [Corrected by Jianing Song, Aug 30 2021] (End)
a(n) is the least k such that (tau(k^2) + 1)/2 = n. - Vladeta Jovovic, Aug 01 2001

Entry revised by N. J. A. Sloane, Aug 14 2005

Offset corrected by David W. Wilson, Dec 27 2018

A055205 Number of nonsquare divisors of n^2.

Original entry on oeis.org

0, 1, 1, 2, 1, 5, 1, 3, 2, 5, 1, 9, 1, 5, 5, 4, 1, 9, 1, 9, 5, 5, 1, 13, 2, 5, 3, 9, 1, 19, 1, 5, 5, 5, 5, 16, 1, 5, 5, 13, 1, 19, 1, 9, 9, 5, 1, 17, 2, 9, 5, 9, 1, 13, 5, 13, 5, 5, 1, 33, 1, 5, 9, 6, 5, 19, 1, 9, 5, 19, 1, 23, 1, 5, 9, 9, 5, 19, 1, 17, 4, 5, 1, 33, 5, 5, 5, 13, 1, 33, 5, 9, 5, 5, 5

Offset: 1

Labos Elemer, Jun 19 2000

Seems to be equal to the number of unordered pairs of coprime divisors of n. (Checked up to 2*10^14.) - Charles R Greathouse IV, May 03 2013
Outline of a proof for this observation, R. J. Mathar, May 05 2013: (Start)
i) To construct the divisors of n, write n=product_i p_i^e_i as the standard prime power decomposition, take any subset of the primes p_i (including the empty set representing the 1) and run with the associated list exponents from 0 up to their individual e_i.
To construct the *nonsquare* divisors of n, ensure that one or more of the associated exponents is/are odd. (The empty set is interpreted as 1^0 with even exponent.) To construct the nonsquare divisors of n^2, the principle remains the same, although the exponents may individually range from 0 up to 2*e_i.
The nonsquare divisor is therefore a nonempty product of prime powers (at least one) with odd exponents times  a (potentially empty) product of prime powers (of different primes) with even exponents.
The nonsquare divisors of n^2 have exponents from 0 up to 2*e_i, but the subset of exponents in the "even/square" factor has e_i candidates (range 2, 4, .., 2*e_i) and in the "odd/nonsquare" factor also only e_i candidates (range 1,3,5,2*e_i-1).
ii) To construct the pairs of coprime divisors of n, take any two non-intersecting subsets of  the set of p_i (possibly the empty subset which represents the factor 1), and let the exponents run from 1 up to their individual e_i in each of the two products.
iii) The bijection between the sets constructed in i) and ii) is given by mapping the two non-intersection prime sets onto each other, and observing that the numbers of compositions of exponents have the same orders in both cases.
(End)

			n = 8, d(64) = 7 and from the 7 divisors {1,4,16,64} are square and the remaining 3 = a(8).
n = 12, d(144) = 15, from which 6 divisors are squares {1,4,9,16,36,144} so a(12) = d(144)-d(12) = 9
a(60) = (number of terms of finite A171425) = 33. [_Reinhard Zumkeller_, Dec 08 2009]

Alois P. Heinz, Table of n, a(n) for n = 1..10000 (first 1000 terms from T. D. Noe)

Cf. A000005, A000290, A048691, A056595.

Cf. A001620, A013661, A082633, A201994, A306016.

a055205 n = length [d | d <- [1..n^2], n^2 `mod` d == 0, a010052 d == 0]
-- Reinhard Zumkeller, Aug 15 2011

Table[Count[Divisors[n^2], d_ /;  ! IntegerQ[Sqrt[d]]], {n, 1, 95}] (* Jean-François Alcover, Mar 22 2011 *)
Table[DivisorSigma[0,n^2]-DivisorSigma[0,n],{n,100}] (* Harvey P. Dale, Sep 02 2017 *)

a(n)=my(f=factor(n)[,2]);prod(i=1,#f,2*f[i]+1)-prod(i=1,#f,f[i]+1) \\ Charles R Greathouse IV, May 02 2013

a(n) = A000005(n^2)-A000005(n) because the number of square divisors of n^2 equals the number of divisors of n.
a(n) = A056595(A000290(n)).
a(n) = A048691(n) - A000005(n). - Reinhard Zumkeller, Dec 08 2009
Sum_{k=1..n} a(k) ~ (n/zeta(2)) * (log(n)^2/2 + c_1 * log(n) + c_2), where c_1 =  3*gamma - 2*zeta'(2)/zeta(2) - zeta(2) - 1 = 0.226634..., c_2 = 3*gamma^2 - (2*gamma - 1)*zeta(2) - 3*gamma_1 + (1 - 3*gamma)*(2*zeta'(2)/zeta(2) + 1) + (2*zeta'(2)/zeta(2))^2 - 2*zeta''(2)/zeta(2) = -0.0529271..., gamma is Euler's constant (A001620), and gamma_1 is the first Stieltjes constant (A082633). - Amiram Eldar, Dec 01 2023

A225520 The number of subsets of the set of divisors of n in which elements are pairwise coprime.

Original entry on oeis.org

2, 4, 4, 6, 4, 10, 4, 8, 6, 10, 4, 16, 4, 10, 10, 10, 4, 16, 4, 16, 10, 10, 4, 22, 6, 10, 8, 16, 4, 30, 4, 12, 10, 10, 10, 26, 4, 10, 10, 22, 4, 30, 4, 16, 16, 10, 4, 28, 6, 16, 10, 16, 4, 22, 10, 22, 10, 10, 4, 50, 4, 10, 16, 14, 10, 30, 4, 16, 10, 30, 4, 36

Offset: 1

R. J. Mathar, May 09 2013

nonn

Note that this is not 1+A048691(n); n=30 is a counterexample.

The number of all subsets of the set of divisors (without the restriction) is 2^A000005(n), which therefore is an upper bound of the current sequence.

			For n=6, the set of divisors is {1,2,3,6} and the a(6)=10 subsets with pairwise coprime entries are {}, {1}, {2}, {3}, {6}, {1,2}, {1,3}, {1,6}, {2,3} and {1,2,3}.

T. D. Noe, Table of n, a(n) for n = 1..359

Cf. A076078 (subsets with lcm equal to n), A084422 (subsets of 1 through n).

paircoprime := proc(s)
    local L,i,j ;
    L := convert(s,list) ;
    for i from 1 to nops(L)-1 do
        for j from i+1 to nops(L) do
            if igcd(op(i,L),op(j,L)) <> 1 then
                return false;
            end if;
        end do:
    end do:
    return true;
end proc:
A225520 := proc(n)
    local dvs,a,p ;
    dvs := numtheory[divisors](n) ;
    a := 0 ;
    for p in combinat[powerset](dvs) do
        if paircoprime(p) then
            a := a+1 ;
        end if;
    end do:
    a ;
end proc:

Table[Length[Select[Subsets[Divisors[n]], If[Length[#] < 2, True, If[Length[#] == 2, CoprimeQ @@ #, And @@ CoprimeQ @@ #]] &]], {n, 100}] (* T. D. Noe, May 09 2013 *)

A100565 a(n) = Card{(x,y,z) : x <= y <= z, x|n, y|n, z|n, gcd(x,y)=1, gcd(x,z)=1, gcd(y,z)=1}.

Original entry on oeis.org

1, 2, 2, 3, 2, 5, 2, 4, 3, 5, 2, 8, 2, 5, 5, 5, 2, 8, 2, 8, 5, 5, 2, 11, 3, 5, 4, 8, 2, 15, 2, 6, 5, 5, 5, 13, 2, 5, 5, 11, 2, 15, 2, 8, 8, 5, 2, 14, 3, 8, 5, 8, 2, 11, 5, 11, 5, 5, 2, 25, 2, 5, 8, 7, 5, 15, 2, 8, 5, 15, 2, 18, 2, 5, 8, 8, 5, 15, 2, 14, 5, 5, 2, 25, 5, 5, 5, 11, 2, 25, 5, 8, 5, 5, 5, 17

Offset: 1

Vladeta Jovovic, Nov 28 2004

First differs from A018892 at a(30) = 15, A018892(30) = 14.
First differs from A343654 at a(210) = 51, A343654(210) = 52.
Also a(n) = Card{(x,y,z) : x <= y <= z and lcm(x,y)=n, lcm(x,z)=n, lcm(y,z)=n}.
In words, a(n) is the number of pairwise coprime unordered triples of divisors of n. - Gus Wiseman, May 01 2021

			From _Gus Wiseman_, May 01 2021: (Start)
The a(n) triples for n = 1, 2, 4, 6, 8, 12, 24:
  (1,1,1)  (1,1,1)  (1,1,1)  (1,1,1)  (1,1,1)  (1,1,1)   (1,1,1)
           (1,1,2)  (1,1,2)  (1,1,2)  (1,1,2)  (1,1,2)   (1,1,2)
                    (1,1,4)  (1,1,3)  (1,1,4)  (1,1,3)   (1,1,3)
                             (1,1,6)  (1,1,8)  (1,1,4)   (1,1,4)
                             (1,2,3)           (1,1,6)   (1,1,6)
                                               (1,2,3)   (1,1,8)
                                               (1,3,4)   (1,2,3)
                                               (1,1,12)  (1,3,4)
                                                         (1,3,8)
                                                         (1,1,12)
                                                         (1,1,24)
(End)

Antti Karttunen, Table of n, a(n) for n = 1..10000

Cf. A000005, A070919, A086222, A086165, A048691, A063647.
Positions of 2's through 5's are A000040, A001248, A030078, A068993.
The version for subsets of {1..n} instead of divisors is A015617.
The version for pairs of divisors is A018892.
The ordered version is A048785.
The strict case is A066620.
The version for strict partitions is A220377.
A version for sets of divisors of any size is A225520.
The version for partitions is A307719 (no 1's: A337563).
The case of distinct parts coprime is A337600 (ordered: A337602).
A001399(n-3) = A069905(n) = A211540(n+2) counts 3-part partitions.
A007304 ranks 3-part strict partitions.
A014311 ranks 3-part compositions.
A014612 ranks 3-part partitions.
A051026 counts pairwise indivisible subsets of {1..n}.
A302696 lists Heinz numbers of pairwise coprime partitions.
A337461 counts 3-part pairwise coprime compositions.
Cf. A000961, A000977, A007360, A023022, A087087, A276187, A282935, A337601, A337603, A338331, A343652, A343654.

pwcop[y_]:=And@@(GCD@@#==1&/@Subsets[y,{2}]);
Table[Length[Select[Tuples[Divisors[n],3],LessEqual@@#&&pwcop[#]&]],{n,30}] (* Gus Wiseman, May 01 2021 *)

A100565(n) = (numdiv(n^3)+3*numdiv(n)+2)/6; \\ Antti Karttunen, May 19 2017

a(n) = (tau(n^3) + 3*tau(n) + 2)/6.

A366438 The number of divisors of the exponentially odd numbers (A268335).

Original entry on oeis.org

1, 2, 2, 2, 4, 2, 4, 4, 2, 2, 4, 4, 2, 2, 4, 4, 2, 8, 4, 4, 2, 8, 2, 6, 4, 4, 4, 2, 4, 4, 8, 2, 8, 2, 4, 2, 4, 2, 8, 4, 8, 4, 4, 2, 2, 4, 4, 8, 2, 4, 8, 2, 2, 4, 4, 8, 2, 4, 2, 4, 4, 4, 8, 2, 4, 4, 4, 4, 12, 2, 2, 8, 2, 8, 8, 4, 2, 2, 8, 4, 2, 8, 4, 4, 4, 16, 4

Offset: 1

Amiram Eldar, Oct 10 2023

1 is the only odd term in this sequence.

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

Cf. A000005, A268335, A366439.

Similar sequences: A048691, A072048, A076400, A358040, A363194, A363195.

f[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, If[AllTrue[e, OddQ], Times @@ (e + 1), Nothing]]; f[1] = 1; Array[f, 150]

lista(max) = for(k = 1, max, my(e = factor(k)[, 2], isexpodd = 1); for(i = 1, #e, if(!(e[i] % 2), isexpodd = 0; break)); if(isexpodd, print1(vecprod(apply(x -> x+1, e)), ", ")));

from math import prod
from itertools import count, islice
from sympy import factorint
def A366438_gen(): # generator of terms
    for n in count(1):
        f = factorint(n).values()
        if all(e&1 for e in f):
            yield prod(e+1 for e in f)
A366438_list = list(islice(A366438_gen(),30)) # Chai Wah Wu, Oct 10 2023

a(n) = A000005(A268335(n)).

cp's OEIS Frontend

A048785 a(0) = 0; a(n) = tau(n^3), where tau = number of divisors (A000005).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Python

Formula

A343656 Array read by antidiagonals where A(n,k) is the number of divisors of n^k.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A061283 Smallest number with exactly 2n-1 divisors.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A126098 Where records occur in A018892.

Views

Author

Keywords

Comments

Links

Crossrefs

Extensions

A343652 Number of maximal pairwise coprime sets of divisors of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A016017 Smallest k such that 1/k can be written as a sum of exactly 2 unit fractions in n ways.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A055205 Number of nonsquare divisors of n^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs