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.

Showing 1-10 of 24 results. Next

A048691 a(n) = d(n^2), where d(k) = A000005(k) is the number of divisors of k.

Original entry on oeis.org

1, 3, 3, 5, 3, 9, 3, 7, 5, 9, 3, 15, 3, 9, 9, 9, 3, 15, 3, 15, 9, 9, 3, 21, 5, 9, 7, 15, 3, 27, 3, 11, 9, 9, 9, 25, 3, 9, 9, 21, 3, 27, 3, 15, 15, 9, 3, 27, 5, 15, 9, 15, 3, 21, 9, 21, 9, 9, 3, 45, 3, 9, 15, 13, 9, 27, 3, 15, 9, 27, 3, 35, 3, 9, 15, 15, 9, 27, 3, 27
Offset: 1

Views

Author

Keywords

Comments

Inverse Moebius transform of A034444: Sum_{d|n} 2^omega(d), where omega(n) = A001221(n) is the number of distinct primes dividing n.
Number of elements in the set {(x,y): x|n, y|n, gcd(x,y)=1}.
Number of elements in the set {(x,y): lcm(x,y)=n}.
Also gives total number of positive integral solutions (x,y), order being taken into account, to the optical or parallel resistor equation 1/x + 1/y = 1/n. Indeed, writing the latter as X*Y=N, with X=x-n, Y=y-n, N=n^2, the one-to-one correspondence between solutions (X, Y) and (x, y) is obvious, so that clearly, the solution pairs (x, y) are tau(N)=tau(n^2) in number. - Lekraj Beedassy, May 31 2002
Number of ordered pairs of positive integers (a,c) such that n^2 - ac = 0. Therefore number of quadratic equations of the form ax^2 + 2nx + c = 0 where a,n,c are positive integers and each equation has two equal (rational) roots, -n/a. (If a and c are positive integers, but, instead, the coefficient of x is odd, it is impossible for the equation to have equal roots.) - Rick L. Shepherd, Jun 19 2005
Problem A1 on the 21st Putnam competition in 1960 (see John Scholes link) asked for the number of pairs of positive integers (x,y) such that xy/(x+y) = n: the answer is a(n); for n = 4, the a(4) = 5 solutions (x,y) are (5,20), (6,12), (8,8), (12,6), (20,5). - Bernard Schott, Feb 12 2023
Numbers k such that a(k)/d(k) is an integer are in A217584 and the corresponding quotients are in A339055. - Bernard Schott, Feb 15 2023

References

  • A. M. Gleason et al., The William Lowell Putnam Mathematical Competitions, Problems & Solutions:1938-1960 Soln. to Prob. 1 1960, p. 516, MAA, 1980.
  • Ross Honsberger, More Mathematical Morsels, Morsel 43, pp. 232-3, DMA No. 10 MAA, 1991.
  • Loren C. Larson, Problem-Solving Through Problems, Prob. 3.3.7, p. 102, Springer 1983.
  • Alfred S. Posamentier and Charles T. Salkind, Challenging Problems in Algebra, Prob. 9-9 pp. 143 Dover NY, 1988.
  • D. O. Shklarsky et al., The USSR Olympiad Problem Book, Soln. to Prob. 123, pp. 28, 217-8, Dover NY.
  • Wacław Sierpiński, Elementary Theory of Numbers, pp. 71-2, Elsevier, North Holland, 1988.
  • James J. Tattersall, Elementary Number Theory in Nine Chapters, Cambridge University Press, 1999, page 91.
  • Charles W. Trigg, Mathematical Quickies, Question 194, pp. 53, 168, Dover, 1985.

Crossrefs

Partial sums give A061503.
For similar LCM sequences, see A070919, A070920, A070921.
For the earliest occurrence of 2n-1 see A016017.

Programs

Formula

a(n) = A000005(A000290(n)).
tau(n^2) = Sum_{d|n} mu(n/d)*tau(d)^2, where mu(n) = A008683(n), cf. A061391.
Multiplicative with a(p^e) = 2e+1. - Vladeta Jovovic, Jul 23 2001
Also a(n) = Sum_{d|n} (tau(d)*moebius(n/d)^2), Dirichlet convolution of A000005 and A008966. - Benoit Cloitre, Sep 08 2002
a(n) = A055205(n) + A000005(n). - Reinhard Zumkeller, Dec 08 2009
Dirichlet g.f.: (zeta(s))^3/zeta(2s). - R. J. Mathar, Feb 11 2011
a(n) = Sum_{d|n} 2^omega(d). Inverse Mobius transform of A034444. - Enrique Pérez Herrero, Apr 14 2012
G.f.: Sum_{k>=1} 2^omega(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Mar 10 2018
Sum_{k=1..n} a(k) ~ n*(6/Pi^2)*(log(n)^2/2 + log(n)*(3*gamma - 1) + 1 - 3*gamma + 3*gamma^2 - 3*gamma_1 + (2 - 6*gamma - 2*log(n))*zeta'(2)/zeta(2) + (2*zeta'(2)/zeta(2))^2 - 2*zeta''(2)/zeta(2)), where gamma is Euler's constant (A001620) and gamma_1 is the first Stieltjes constant (A082633). - Amiram Eldar, Jan 26 2023

Extensions

Additional comments from Vladeta Jovovic, Apr 29 2001

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

Views

Author

Gus Wiseman, Apr 28 2021

Keywords

Comments

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).

Examples

			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
		

Crossrefs

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}.

Programs

  • Mathematica
    Table[DivisorSigma[0,k^(n-k)],{n,10},{k,n}]
  • PARI
    A(n, k) = numdiv(n^k); \\ Seiichi Manyama, May 15 2021

Formula

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

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

Views

Author

N. J. A. Sloane, Mar 05 2007

Keywords

Comments

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

Crossrefs

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

Extensions

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

Views

Author

Gus Wiseman, Apr 25 2021

Keywords

Comments

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}

Examples

			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}
		

Crossrefs

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.

Programs

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

Formula

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

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

Views

Author

Vladeta Jovovic, Nov 28 2004

Keywords

Comments

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

Examples

			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)
		

Crossrefs

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.

Programs

  • Mathematica
    pwcop[y_]:=And@@(GCD@@#==1&/@Subsets[y,{2}]);
    Table[Length[Select[Tuples[Divisors[n],3],LessEqual@@#&&pwcop[#]&]],{n,30}] (* Gus Wiseman, May 01 2021 *)
  • PARI
    A100565(n) = (numdiv(n^3)+3*numdiv(n)+2)/6; \\ Antti Karttunen, May 19 2017

Formula

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

A061391 a(n) = t(n,3) = Sum_{d|n} tau(d^3), where tau(n) = number of divisors of n, cf. A000005.

Original entry on oeis.org

1, 5, 5, 12, 5, 25, 5, 22, 12, 25, 5, 60, 5, 25, 25, 35, 5, 60, 5, 60, 25, 25, 5, 110, 12, 25, 22, 60, 5, 125, 5, 51, 25, 25, 25, 144, 5, 25, 25, 110, 5, 125, 5, 60, 60, 25, 5, 175, 12, 60, 25, 60, 5, 110, 25, 110, 25, 25, 5, 300, 5, 25, 60, 70, 25, 125, 5, 60, 25, 125, 5, 264
Offset: 1

Views

Author

Vladeta Jovovic, Apr 29 2001

Keywords

Comments

Inverse Mobius transform of A048785. - R. J. Mathar, Feb 09 2011

Crossrefs

Cf. t(n, 0) = A000005(n), t(n, 1) = A007425(n), t(n, 2) = A035116(n).
Cf. A048691.

Programs

  • Mathematica
    f[p_, e_] := (3*e^2 + 5*e + 2)/2; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 16 2020 *)
  • PARI
    A061391 = n -> sumdiv(n, d, numdiv(d^3));
    for(n=1, 10000, write("b061391.txt", n, " ", A061391(n)));
    \\ Antti Karttunen, Jan 17 2017
    
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, (1 + 2*X)/(1 - X)^3)[n], ", ")) \\ Vaclav Kotesovec, May 15 2021
    
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, (1 - 3*X^2 + 2*X^3)/(1 - X)^5)[n], ", ")) \\ Vaclav Kotesovec, Aug 20 2021

Formula

t(n, k) = Sum_{d|n} tau(d^k) is multiplicative: if the canonical factorization of n = Product p^e(p) over primes then t(n, k) = Product t(p^e(p), k), t(p^e(p), k) = (1/2) *(k*e(p)+2)*(e(p)+1).
For k=2 we get an interesting identity: Sum_{d|n} tau(d^2)=(tau(n))^2, cf. A048691, A035116.
a(n) = Sum_{d|n} tau(n*d). - Benoit Cloitre, Nov 30 2002
G.f.: Sum_{n>=1} tau(n^3)*x^n/(1-x^n). - Joerg Arndt, Jan 01 2011
Dirichlet g.f.: zeta(s)^3 * Product_{primes p} (1 + 2/p^s). - Vaclav Kotesovec, May 15 2021
Dirichlet g.f.: zeta(s)^5 * Product_{primes p} (1 - 3/p^(2*s) + 2/p^(3*s)). - Vaclav Kotesovec, Aug 20 2021

A066620 Number of unordered triples of distinct pairwise coprime divisors of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 7, 0, 0, 1, 1, 1, 4, 0, 1, 1, 3, 0, 7, 0, 2, 2, 1, 0, 4, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 13, 0, 1, 2, 0, 1, 7, 0, 2, 1, 7, 0, 6, 0, 1, 2, 2, 1, 7, 0, 4, 0, 1, 0, 13, 1, 1, 1, 3, 0, 13, 1, 2, 1, 1, 1, 5, 0, 2, 2, 4, 0, 7, 0
Offset: 1

Views

Author

K. B. Subramaniam (kb_subramaniambalu(AT)yahoo.com) and Amarnath Murthy, Dec 24 2001

Keywords

Comments

a(m) = a(n) if m and n have same factorization structure.

Examples

			a(24) = 3: the divisors of 24 are 1, 2, 3, 4, 6, 8, 12 and 24. The triples are (1, 2, 3), (1, 2, 9), (1, 3, 4).
a(30) = 7: the triples are (1, 2, 3), (1, 2, 5), (1, 3, 5), (2, 3, 5), (1, 3, 10), (1, 5, 6), (1, 2, 15).
		

References

  • Amarnath Murthy, Decomposition of the divisors of a natural number into pairwise coprime sets, Smarandache Notions Journal, vol. 12, No. 1-2-3, Spring 2001.pp 303-306.

Crossrefs

Positions of zeros are A000961.
Positions of ones are A006881.
The version for subsets of {1..n} instead of divisors is A015617.
The non-strict ordered version is A048785.
The version for pairs of divisors is A063647.
The non-strict version (3-multisets) is A100565.
The version for partitions is A220377 (non-strict: A307719).
A version for sets of divisors of any size is A225520.
A000005 counts divisors.
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.
A018892 counts unordered pairs of coprime divisors (ordered: A048691).
A051026 counts pairwise indivisible subsets of {1..n}.
A337461 counts 3-part pairwise coprime compositions.
A338331 lists Heinz numbers of pairwise coprime partitions.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Divisors[n],{3}],CoprimeQ@@#&]],{n,100}] (* Gus Wiseman, Apr 28 2021 *)
  • PARI
    A066620(n) = (numdiv(n^3)-3*numdiv(n)+2)/6; \\ After Jovovic's formula. - Antti Karttunen, May 27 2017
    
  • Python
    from sympy import divisor_count as d
    def a(n): return (d(n**3) - 3*d(n) + 2)/6 # Indranil Ghosh, May 27 2017

Formula

In the reference it is shown that if k is a squarefree number with r prime factors and m with (r+1) prime factors then a(m) = 4*a(k) + 2^k - 1.
a(n) = (tau(n^3)-3*tau(n)+2)/6. - Vladeta Jovovic, Nov 27 2004

Extensions

More terms from Vladeta Jovovic, Apr 03 2003
Name corrected by Andrey Zabolotskiy, Dec 09 2020
Name corrected by Gus Wiseman, Apr 28 2021 (ordered version is 6*a(n))

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

Original entry on oeis.org

1, 2, 2, 3, 2, 6, 2, 4, 3, 6, 2, 10, 2, 6, 6, 5, 2, 10, 2, 10, 6, 6, 2, 14, 3, 6, 4, 10, 2, 22, 2, 6, 6, 6, 6, 17, 2, 6, 6, 14, 2, 22, 2, 10, 10, 6, 2, 18, 3, 10, 6, 10, 2, 14, 6, 14, 6, 6, 2, 38, 2, 6, 10, 7, 6, 22, 2, 10, 6, 22, 2, 24, 2, 6, 10, 10, 6, 22, 2, 18, 5, 6, 2, 38, 6, 6
Offset: 1

Views

Author

Keywords

Crossrefs

Programs

Formula

a(n) = (2+A048785(n))/3. - R. J. Mathar, May 07 2021

Extensions

Definition corrected by Vladeta Jovovic, Sep 03 2005

A059907 a(n) = |{m : multiplicative order of n mod m = 2}|.

Original entry on oeis.org

0, 1, 2, 2, 5, 2, 6, 4, 6, 3, 12, 2, 10, 6, 8, 4, 13, 2, 18, 6, 10, 4, 16, 4, 12, 9, 12, 4, 26, 2, 20, 6, 8, 12, 20, 4, 15, 6, 16, 4, 32, 2, 24, 10, 10, 6, 20, 4, 26, 9, 18, 4, 26, 6, 32, 12, 12, 4, 28, 2, 20, 10, 12, 18, 25, 4, 24, 6, 26, 4, 52, 2, 18, 10, 12, 18, 26, 4, 40, 8, 14, 5, 28
Offset: 1

Views

Author

Vladeta Jovovic, Feb 08 2001

Keywords

Comments

The multiplicative order of a mod m, GCD(a,m) = 1, is the smallest natural number d for which a^d = 1 (mod m).

Examples

			a(2) = |{3}| = 1, a(3) = |{4,8}| = 2, a(4) = |{5,15}| = 2, a(5) = |{3,6,8,12,24}| = 5, a(6) = |{7,35}| = 2, a(7) = |{4,8,12,16,24,48}| = 6,...
		

Crossrefs

Programs

  • Maple
    with(numtheory):f := n->tau(n^2-1)-tau(n-1):for n from 1 to 100 do printf(`%d,`,f(n)) od:
  • Mathematica
    a[n_] := Subtract @@ DivisorSigma[0, {n^2-1, n-1}]; a[1] = 0; Array[a, 100] (* Amiram Eldar, Jan 25 2025 *)
  • PARI
    a(n) = if(n == 1, 0, numdiv(n^2-1) - numdiv(n-1)); \\ Amiram Eldar, Jan 25 2025

Formula

a(n) = tau(n^2-1)-tau(n-1), where tau(n) = number of divisors of n A000005. Generally, if b(n, r) = |{m : multiplicative order of n mod m = r}| then b(n, r) = Sum_{d|r} mu(d)*tau(n^(r/d)-1), where mu(n) = Moebius function A008683.

A343653 Number of non-singleton pairwise coprime nonempty sets of divisors > 1 of n.

Original entry on oeis.org

0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 2, 0, 1, 1, 0, 0, 2, 0, 2, 1, 1, 0, 3, 0, 1, 0, 2, 0, 7, 0, 0, 1, 1, 1, 4, 0, 1, 1, 3, 0, 7, 0, 2, 2, 1, 0, 4, 0, 2, 1, 2, 0, 3, 1, 3, 1, 1, 0, 13, 0, 1, 2, 0, 1, 7, 0, 2, 1, 7, 0, 6, 0, 1, 2, 2, 1, 7, 0, 4, 0, 1, 0, 13, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Apr 25 2021

Keywords

Comments

First differs from A066620 at a(210) = 36, A066620(210) = 35.

Examples

			The a(n) sets for n = 6, 12, 24, 30, 36, 60, 72, 96:
  {2,3}  {2,3}  {2,3}  {2,3}    {2,3}  {2,3}    {2,3}  {2,3}
         {3,4}  {3,4}  {2,5}    {2,9}  {2,5}    {2,9}  {3,4}
                {3,8}  {3,5}    {3,4}  {3,4}    {3,4}  {3,8}
                       {5,6}    {4,9}  {3,5}    {3,8}  {3,16}
                       {2,15}          {4,5}    {4,9}  {3,32}
                       {3,10}          {5,6}    {8,9}
                       {2,3,5}         {2,15}
                                       {3,10}
                                       {3,20}
                                       {4,15}
                                       {5,12}
                                       {2,3,5}
                                       {3,4,5}
		

Crossrefs

The case of pairs is A089233.
The version with 1's, empty sets, and singletons is A225520.
The version for subsets of {1..n} is A320426.
The version for strict partitions is A337485.
The version for compositions is A337697.
The version for prime indices is A337984.
The maximal case with 1's is A343652.
The version with empty sets is a(n) + 1.
The version with singletons is A343654(n) - 1.
The version with empty sets and singletons is A343654.
The version with 1's is A343655.
The maximal case 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.
A051026 counts pairwise indivisible subsets of {1..n}.
A100565 counts pairwise coprime unordered triples of divisors.
A305713 counts pairwise coprime non-singleton strict partitions.
A343659 counts maximal pairwise coprime subsets of {1..n}.

Programs

  • Mathematica
    Table[Length[Select[Subsets[Rest[Divisors[n]]],CoprimeQ@@#&]],{n,100}]
Showing 1-10 of 24 results. Next