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 11 results. Next

A068140 Smaller of two consecutive numbers each divisible by a cube greater than one.

Original entry on oeis.org

80, 135, 296, 343, 351, 375, 512, 567, 624, 728, 783, 944, 999, 1160, 1215, 1375, 1376, 1431, 1592, 1624, 1647, 1808, 1863, 2024, 2079, 2240, 2295, 2375, 2400, 2456, 2511, 2624, 2672, 2727, 2888, 2943, 3087, 3104, 3159, 3320, 3375, 3429, 3536, 3591
Offset: 1

Views

Author

Amarnath Murthy, Feb 22 2002

Keywords

Comments

Cubeful numbers with cubeful successors. This is to cubes as A068781 is to squares. 1375 is the smallest of three consecutive numbers divisible by a cube, since 1375 = 5^3 * 11 and 1376 = 2^5 * 43 and 1377 = 3^4 * 17. What is the smallest of four consecutive numbers divisible by a cube? Of n consecutive numbers divisible by a cube? - Jonathan Vos Post, Sep 18 2007
22624 is the smallest of four consecutive numbers each divisible by a cube, with factorizations 2^5 * 7 * 101, 5^3 * 181, 2 * 3^3 * 419, and 11^3 * 17. - D. S. McNeil, Dec 10 2010
18035622 is the smallest of five consecutive numbers each divisible by a cube. 4379776620 is the smallest of six consecutive numbers each divisible by a cube. 1204244328624 is the smallest of seven consecutive numbers each divisible by a cube. - Donovan Johnson, Dec 13 2010
The sequence is the union, over all pairs of distinct primes (p,q), of numbers == 0 mod p^3 and == -1 mod q^3 or vice versa. - Robert Israel, Aug 13 2018
The asymptotic density of this sequence is 1 - 2/zeta(3) + Product_{p prime} (1 - 2/p^3) = 1 - 2 * A088453 + A340153 = 0.013077991848467056243... - Amiram Eldar, Feb 16 2021

Examples

			343 is a term as 343 = 7^3 and 344= 2^3 * 43.
		

Crossrefs

Programs

  • Maple
    isA068140 := proc(n)
        isA046099(n) and isA046099(n+1) ;
    end proc:
    for n from 1 to 4000 do
        if isA068140(n) then
            printf("%d,",n) ;
        end if;
    end do: # R. J. Mathar, Dec 08 2015
  • Mathematica
    a = b = 0; Do[b = Max[ Transpose[ FactorInteger[n]] [[2]]]; If[a > 2 && b > 2, Print[n - 1]]; a = b, {n, 2, 5000}]
    Select[Range[2, 6000], Max[Transpose[FactorInteger[ # ]][[2]]] > 2 && Max[Transpose[FactorInteger[ # + 1]][[2]]] > 2 &] (* Jonathan Vos Post, Sep 18 2007 *)
    SequencePosition[Table[If[AnyTrue[Rest[Divisors[n]],IntegerQ[Surd[#,3]]&],1,0],{n,3600}],{1,1}][[All,1]] (* Requires Mathematica version 10 or later *) (* Harvey P. Dale, Apr 18 2020 *)

Formula

{k such that k is in A046099 and k+1 is in A046099}. - Jonathan Vos Post, Sep 18 2007

Extensions

Edited and extended by Robert G. Wilson v, Mar 02 2002
Title edited, cross-references added by Matthew Vandermast, Dec 09 2010
Definition clarified by Harvey P. Dale, Apr 18 2020

A071318 Lesser of 2 consecutive numbers which are cubefree and not squarefree, i.e., numbers k such that both k and k+1 are in A067259.

Original entry on oeis.org

44, 49, 75, 98, 99, 116, 147, 171, 244, 260, 275, 315, 332, 363, 387, 475, 476, 507, 524, 531, 548, 549, 603, 604, 636, 692, 724, 725, 747, 764, 774, 819, 844, 845, 846, 867, 908, 924, 931, 963, 980, 1035, 1075, 1083, 1179, 1196, 1251, 1274, 1275, 1324
Offset: 1

Views

Author

Labos Elemer, May 29 2002

Keywords

Comments

The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 0, 5, 41, 407, 4125, 41215, 412331, 4123625, 41236308, ... . Apparently, the asymptotic density of this sequence exists and equals 0.041236... . - Amiram Eldar, Jan 18 2023
The asymptotic density of this sequence is Product_{p prime} (1 - 2/p^3) - 2 * Product_{p prime} (1 - 1/p^2 - 1/p^3) + Product_{p prime} (1 - 2/p^2) = 0.041236147082334172926... . - Amiram Eldar, Jan 05 2024

Examples

			75 is a term since 75 = 3*5^2 and 76 = 2^2*19.
		

Crossrefs

Programs

  • Haskell
    a071318 n = a071318_list !! (n-1)
    a071318_list = [x | x <- [1..],  a212793 x == 1, a008966 x == 0,
                        let y = x+1, a212793 y == 1, a008966 y == 0]
    -- Reinhard Zumkeller, May 27 2012
    
  • Mathematica
    With[{s = Select[Range[1350], And[MemberQ[#, 2], FreeQ[#, k_ /; k > 2]] &@ FactorInteger[#][[All, -1]] &]}, Function[t, Part[s, #] &@ Position[t, 1][[All, 1]]]@ Differences@ s] (* Michael De Vlieger, Jul 30 2017 *)
  • PARI
    isok(n) = (n>1) && (vecmax(factor(n)[, 2])==2) && (vecmax(factor(n+1)[, 2])==2); \\ Michel Marcus, Aug 02 2017

Formula

A051903(k) = A051903(k+1) = 2 when k is a term.

A340153 Decimal expansion of Product_{p prime} (1 - 2/p^3).

Original entry on oeis.org

6, 7, 6, 8, 9, 2, 7, 3, 7, 0, 0, 9, 8, 8, 1, 9, 9, 3, 6, 1, 0, 2, 3, 7, 3, 2, 6, 7, 2, 4, 3, 8, 9, 2, 1, 2, 7, 9, 7, 6, 7, 8, 3, 9, 7, 4, 5, 9, 7, 8, 8, 8, 4, 5, 2, 7, 3, 2, 9, 7, 8, 2, 3, 0, 4, 4, 3, 2, 6, 3, 2, 0, 4, 6, 0, 3, 5, 7, 8, 6, 0, 5, 1, 2, 8, 3, 2, 6, 8, 4, 8, 1, 1, 1, 1, 0, 8, 4, 4, 9, 3, 1, 7, 0, 8, 4
Offset: 0

Views

Author

Amiram Eldar, Dec 29 2020

Keywords

Comments

The asymptotic density of the sequence of cubefree numbers k such that k+1 is also cubefree (A340152) (Carlitz, 1932).

Examples

			0.67689273700988199361023732672438921279767839745978...
		

Crossrefs

Programs

  • Mathematica
    $MaxExtraPrecision = 500; m = 500; c = LinearRecurrence[{0, 0, 2}, {0, 0, -6}, m]; RealDigits[Exp[NSum[Indexed[c, n]*(PrimeZetaP[n])/n, {n, 2, m}, NSumTerms -> m, WorkingPrecision -> m]], 10, 100][[1]]
  • PARI
    prodeulerrat(1 - 2/p^3)

Extensions

More digits from Vaclav Kotesovec, Jan 16 2021

A349307 Numbers k such that A072911(k) = A072911(k+1) > 1.

Original entry on oeis.org

80, 135, 296, 343, 375, 567, 624, 728, 783, 944, 1160, 1431, 1592, 1624, 1647, 1808, 1863, 2024, 2240, 2295, 2456, 2511, 2624, 2727, 2888, 3087, 3320, 3429, 3536, 3591, 3624, 3752, 3992, 4023, 4184, 4239, 4375, 4400, 4455, 4624, 4671, 4887, 4912, 4913, 5048, 5103
Offset: 1

Views

Author

Amiram Eldar, Nov 14 2021

Keywords

Comments

Without the restriction that A072911(k) > 1, all the terms of A340152 would be in this sequence.
In contrast to A001274, which has only one known pair of consecutive terms (5186 and 5187), this sequence seems to have many pairs of consecutive terms. The smaller members of these pairs are 4912, 5750, 6858, ...

Examples

			80 is a term since A072911(80) = A072911(81) = 2.
		

Crossrefs

Subsequence of A068140.
Similar sequences: A001274, A287055, A293184, A326403, A349308.

Programs

  • Mathematica
    f[p_, e_] := EulerPhi[e]; ephi[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[5000], ephi[#] == ephi[# + 1] > 1 &]

A349233 Numbers k such that k and k+2 are consecutive cubefree numbers.

Original entry on oeis.org

7, 15, 23, 26, 31, 39, 47, 53, 55, 63, 71, 87, 95, 103, 107, 111, 119, 124, 127, 143, 151, 159, 161, 167, 175, 183, 188, 191, 199, 207, 215, 223, 231, 239, 242, 247, 249, 255, 263, 269, 271, 279, 287, 303, 311, 319, 323, 327, 335, 359, 367, 377, 383, 391, 399
Offset: 1

Views

Author

Amiram Eldar, Nov 11 2021

Keywords

Comments

The asymptotic density of this sequence is 0.142325864924778... (Mossinghoff et al., 2021).

Examples

			7 is a term since 7 and 7 + 2 = 9 = 3^2 are cubefree, and 7 + 1 = 8 = 2^3 is not.
		

Crossrefs

Programs

  • Mathematica
    cubeFreeQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], # < 3 &]; Select[Range[400], Boole[cubeFreeQ /@ (# + {0, 1, 2})] == {1, 0, 1} &]

A349234 Numbers k such that k and k+3 are consecutive cubefree numbers.

Original entry on oeis.org

79, 134, 295, 342, 350, 374, 511, 566, 623, 727, 782, 943, 998, 1159, 1214, 1430, 1591, 1623, 1646, 1807, 1862, 2023, 2078, 2239, 2294, 2374, 2399, 2455, 2510, 2623, 2671, 2726, 2887, 2942, 3086, 3103, 3158, 3319, 3374, 3428, 3535, 3590, 3623, 3751, 3806, 3967
Offset: 1

Views

Author

Amiram Eldar, Nov 11 2021

Keywords

Comments

The asymptotic density of this sequence is 0.0123046264590258... (Mossinghoff et al., 2021).

Examples

			79 is a term since 79 and 79 + 3 = 82 = 2*41 are cubefree, and 79 + 1 = 80 = 2^4*5 and 79 + 2 = 81 = 3^4 are not.
		

Crossrefs

Programs

  • Mathematica
    cubeFreeQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], # < 3 &]; Select[Range[4000], Boole[cubeFreeQ /@ (# + Range[0, 3])] == {1, 0, 0, 1} &]
    SequencePosition[Table[If[Max[FactorInteger[n][[All,2]]]<3,1,0],{n,4000}],{1,0,0,1}][[All,1]] (* Harvey P. Dale, May 08 2022 *)

A349235 Numbers k such that k and k+4 are consecutive cubefree numbers.

Original entry on oeis.org

1374, 4373, 4911, 5749, 6857, 13309, 13374, 16118, 21247, 24351, 25622, 28374, 31373, 32749, 33613, 40471, 41741, 48247, 49623, 49733, 52622, 55374, 57966, 58373, 59749, 75247, 76623, 79622, 82374, 85373, 86749, 90206, 94470, 98439, 102247, 103623, 106622, 107701
Offset: 1

Views

Author

Amiram Eldar, Nov 11 2021

Keywords

Comments

The asymptotic density of this sequence is 0.000379098586237504... (Mossinghoff et al., 2021).

Examples

			1374 is a term since 1374 = 2*3*229 and 1374 + 4 = 1378 = 2*13*53 are cubefree, and 1374 + 1 = 1375 = 5^3*11, 1374 + 2 = 1376 = 2^5*43 and 1374 + 3 = 1377 = 3^4*17 are not.
		

Crossrefs

Programs

  • Mathematica
    cubeFreeQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], # < 3 &]; Select[Range[10^5], Boole[cubeFreeQ /@ (# + Range[0, 4])] == {1, 0, 0, 0, 1} &]

A349308 Numbers k such that A321167(k) = A321167(k+1) > 1.

Original entry on oeis.org

80, 135, 296, 343, 375, 624, 728, 1160, 1431, 1592, 1624, 2240, 2295, 2456, 2511, 2624, 2727, 2888, 3429, 3591, 3624, 3752, 3992, 4023, 4184, 4671, 4887, 4913, 5048, 5144, 5264, 5319, 5480, 5696, 6183, 6344, 6375, 6591, 6615, 6776, 6858, 6859, 7479, 7624, 7640
Offset: 1

Views

Author

Amiram Eldar, Nov 14 2021

Keywords

Comments

Without the restriction that A321167(k) > 1, all the terms of A340152 would be in this sequence.
In contrast to A001274, which has only one known pair of consecutive terms (5186 and 5187), this sequence seems to have many pairs of consecutive terms. The smaller members of these pairs are 6858, 13375, 22625, ...

Examples

			80 is a term since A321167(80) = A321167(81) = 3.
		

Crossrefs

Subsequence of A068140.
Similar sequences: A001274, A287055, A293184, A326403, A349307.

Programs

  • Mathematica
    f[p_, e_] := p^e - 1; uphi[1] = 1; uphi[n_] := Times @@ f @@@ FactorInteger[n]; fe[p_, e_] := uphi[e]; euphi[n_] := Times @@ fe @@@ FactorInteger[n]; Select[Range[8000], euphi[#] == euphi[# + 1] > 1 &]

A367695 Numbers k such that k and k+1 are both exponentially odd numbers (A268335).

Original entry on oeis.org

1, 2, 5, 6, 7, 10, 13, 14, 21, 22, 23, 26, 29, 30, 31, 32, 33, 34, 37, 38, 39, 40, 41, 42, 46, 53, 54, 55, 56, 57, 58, 61, 65, 66, 69, 70, 73, 77, 78, 82, 85, 86, 87, 88, 93, 94, 95, 96, 101, 102, 103, 104, 105, 106, 109, 110, 113, 114, 118, 119, 122, 127, 128
Offset: 1

Views

Author

Amiram Eldar, Nov 27 2023

Keywords

Comments

The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 6, 48, 478, 4734, 47195, 471707, 4716892, 47168363, 471681183, 4716806520, ... . Apparently, the asymptotic density of this sequence exists and equals Product_{p prime} (1 - 2/(p*(p+1))) = 0.47168... (A307868).

Crossrefs

Subsequence of A268335.
Cf. A307868.
Subsequences: A007674, A325058.
Similar sequences: A071318, A121495, A340152, A367696.

Programs

  • Mathematica
    expOddQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], OddQ]; Select[Range[128], And @@ expOddQ /@ {#, # + 1} &]
  • PARI
    isexpodd(n) = {my(f = factor(n)); for(i=1, #f~, if (!(f[i, 2] % 2), return (0))); 1;}
    is(n) = isexpodd(n) && isexpodd(n+1)

A367696 Numbers k such that k and k+1 are both exponentially odious numbers (A270428).

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 25, 28, 29, 30, 33, 34, 35, 36, 37, 38, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 57, 58, 59, 60, 61, 62, 65, 66, 67, 68, 69, 70, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85
Offset: 1

Views

Author

Amiram Eldar, Nov 27 2023

Keywords

Comments

The numbers of terms not exceeding 10^k, for k = 1, 2, ..., are 8, 78, 762, 7615, 76113, 761127, 7611222, 76111895, 761119135, 7611190807, ... . Apparently, the asymptotic density of this sequence exists and equals 0.761119... .

Crossrefs

Subsequence of A270428.
Subsequences: A007674, A367697.
Similar sequences: A071318, A121495, A340152, A367695.

Programs

  • Mathematica
    expOdQ[n_] := AllTrue[FactorInteger[n][[;; , 2]], OddQ[DigitCount[#, 2, 1]] &]; Select[Range[100], And @@ expOdQ /@ {#, # + 1} &]
  • PARI
    isexpod(n) = {my(f = factor(n)); for(i=1, #f~, if (!(hammingweight(f[i, 2]) % 2), return (0))); 1;}
    is(n) = isexpod(n) && isexpod(n+1)
Showing 1-10 of 11 results. Next