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.

A059269 Numbers m for which the number of divisors, tau(m), is divisible by 3.

Original entry on oeis.org

4, 9, 12, 18, 20, 25, 28, 32, 36, 44, 45, 49, 50, 52, 60, 63, 68, 72, 75, 76, 84, 90, 92, 96, 98, 99, 100, 108, 116, 117, 121, 124, 126, 132, 140, 144, 147, 148, 150, 153, 156, 160, 164, 169, 171, 172, 175, 180, 188, 196, 198, 200, 204, 207, 212, 220, 224, 225, 228
Offset: 1

Views

Author

Avi Peretz (njk(AT)netvision.net.il), Jan 24 2001

Keywords

Comments

tau(n) is divisible by 3 iff at least one prime in the prime factorization of n has exponent of the form 3*m + 2. This sequence is an extension of the sequence A038109 in which the numbers has at least one prime with exponent 2 (the case of m = 0 here ) in their prime factorization.
The union of A211337 and A211338 is the complementary sequence to this one. - Douglas Latimer, Apr 12 2012
Numbers whose cubefree part (A050985) is not squarefree (A005117). - Amiram Eldar, Mar 09 2021

Examples

			a(7) = 28 is a term because the number of divisors of 28, d(28) = 6, is divisible by 3.
		

Crossrefs

Characteristic function: A353470.

Programs

  • Maple
    with(numtheory): for n from 1 to 1000 do if tau(n) mod 3 = 0 then printf(`%d,`,n) fi: od:
  • Mathematica
    Select[Range[230], Divisible[DivisorSigma[0, #], 3] &] (* Amiram Eldar, Jul 26 2020 *)
  • PARI
    is(n)=vecmax(factor(n)[,2]%3)==2 \\ Charles R Greathouse IV, Apr 10 2012
    
  • PARI
    is(n)=numdiv(n)%3==0 \\ Charles R Greathouse IV, Sep 18 2015

Formula

Conjecture: a(n) ~ k*n where k = 1/(1 - Product(1 - (p-1)/(p^(3*i)))) = 3.743455... where p ranges over the primes and i ranges over the positive integers. - Charles R Greathouse IV, Apr 13 2012
The asymptotic density of this sequence is 1 - zeta(3)/zeta(2) = 1 - 6*zeta(3)/Pi^2 = 0.2692370305... (Sathe, 1945). Therefore, the above conjecture, a(n) ~ k*n, is true, but k = 1/(1-6*zeta(3)/Pi^2) = 3.7141993349... - Amiram Eldar, Jul 26 2020
A001248 UNION A030515 UNION A030627 UNION A030630 UNION A030633 UNION A030636 UNION ... - R. J. Mathar, May 05 2023

Extensions

More terms from James Sellers, Jan 24 2001