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-2 of 2 results.

A365866 Integers that are divisible by the cube of their least prime factor.

Original entry on oeis.org

8, 16, 24, 27, 32, 40, 48, 56, 64, 72, 80, 81, 88, 96, 104, 112, 120, 125, 128, 135, 136, 144, 152, 160, 168, 176, 184, 189, 192, 200, 208, 216, 224, 232, 240, 243, 248, 256, 264, 272, 280, 288, 296, 297, 304, 312, 320, 328, 336, 343, 344, 351, 352, 360, 368, 376
Offset: 1

Views

Author

Amiram Eldar, Sep 21 2023

Keywords

Comments

Numbers k such that A067029(k) >= 3.
The asymptotic density of terms with least prime factor prime(n) (within all the positive integers) is d(n) = (1/prime(n)^3) * Product_{k=1..(n-1)} (1-1/prime(k)). For example, for n = 1, 2, 3, 4, 5 and 6, d(n) = 1/8, 1/54, 1/375, 4/5145, 8/46585 and 16/169169.
The asymptotic density of this sequence is Sum_{n>=1} d(n) = 0.147333958520714364623977...

Examples

			8 is a term since 2 is the least prime factor of 8 and 8 is divisible by 2^3 = 8.
		

Crossrefs

Subsequence of A046099 and A283050.
Subsequences: A365867, A365868.

Programs

  • Mathematica
    Select[Range[400], FactorInteger[#][[1, -1]] >= 3 &]
  • PARI
    is(n) = n > 1 && factor(n)[1,2] >= 3;

A365868 Starts of runs of 3 consecutive integers that are divisible by the cube of their least prime factor.

Original entry on oeis.org

1375, 16119, 25623, 28375, 52623, 55375, 57967, 79623, 82375, 90207, 94471, 106623, 109375, 129623, 133623, 134431, 136375, 160623, 163375, 164295, 187623, 190375, 206143, 214623, 217375, 241623, 244375, 268623, 271375, 280231, 295623, 298375, 312471, 322623, 325375
Offset: 1

Views

Author

Amiram Eldar, Sep 21 2023

Keywords

Comments

Numbers k such that k, k+1 and k+2 are all terms of A365866.
Numbers of the form 4*k+2 are not terms of A365866. Therefore, there are no runs of 4 or more consecutive integers. Since the middle integer in each triple is divisible by 8, all the terms of this sequence are of the form 8*k+7.
The numbers of terms not exceeding 10^k, for k = 4, 5, ..., are , 1, 11, 109, 1092, 10899, 109125, 1091335, 10912897, ... . Apparently, the asymptotic density of this sequence exists and equals 0.0001091... .

Examples

			1375 is a term since 5 is the least prime factor of 1375 and 1375 is divisible by 5^3 = 125, 2 is the least prime factor of 1376 and 1376 is divisible by 2^3 = 8, and 3 is the least prime factor of 1377 and 1377 is divisible by 3^3 = 27.
		

Crossrefs

Cf. A067029.
Subsequence of A004771, A068140, A365866 and A365867.

Programs

  • Mathematica
    Select[8 * Range[41000] + 7, AllTrue[# + {0, 1, 2}, FactorInteger[#1][[1, -1]] >= 3 &] &]
  • PARI
    is(n) = factor(n)[1,2] >= 3;
    lista(kmax) = forstep(k = 7, kmax, 8, if(is(k) && is(k+1) && is(k+2), print1(k, ", ")));
Showing 1-2 of 2 results.