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

A369938 Numbers whose maximal exponent in their prime factorization is a power of 2.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77
Offset: 1

Views

Author

Amiram Eldar, Feb 06 2024

Keywords

Comments

First differs from its subsequence A138302 \ {1} at n = 378: a(378) = 432 = 2^4 * 3^3 is not a term of A138302.
First differs from A096432, A220218 \ {1}, A335275 \ {1} and A337052 \ {1} at n = 56, and from A270428 \ {1} at n = 113.
Numbers k such that A051903(k) is a power of 2.
The asymptotic density of this sequence is 1/zeta(3) + Sum_{k>=2} (1/zeta(2^k+1) - 1/zeta(2^k)) = 0.87442038669659566330... .

Crossrefs

Programs

  • Mathematica
    pow2Q[n_] := n == 2^IntegerExponent[n, 2];
    Select[Range[2, 100], pow2Q[Max[FactorInteger[#][[;; , 2]]]] &]
    Select[Range[2,80],IntegerQ[Log2[Max[FactorInteger[#][[;;,2]]]]]&] (* Harvey P. Dale, Nov 06 2024 *)
  • PARI
    ispow2(n) = n >> valuation(n, 2) == 1;
    is(n) = n > 1 && ispow2(vecmax(factor(n)[, 2]));

A369939 Numbers whose maximal exponent in their prime factorization is a Fibonacci number.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71
Offset: 1

Views

Author

Amiram Eldar, Feb 06 2024

Keywords

Comments

First differs from its subsequence A115063 at n = 2448. a(2448) = 2592 = 2^5 * 3^4 is not a term of A115063.
First differs from A209061 at n = 62.
Numbers k such that A051903(k) is a Fibonacci number.
The asymptotic density of this sequence is 1/zeta(4) + Sum_{k>=5} (1/zeta(Fibonacci(k)+1) - 1/zeta(Fibonacci(k))) = 0.94462177878047854647... .

Crossrefs

Similar sequences: A368714, A369937, A369938.

Programs

  • Mathematica
    fibQ[n_] := Or @@ IntegerQ /@ Sqrt[5*n^2 + {-4, 4}];
    Select[Range[100], fibQ[Max[FactorInteger[#][[;; , 2]]]] &]
  • PARI
    isfib(n) = issquare(5*n^2 - 4) || issquare(5*n^2 + 4);
    is(n) = n == 1 || isfib(vecmax(factor(n)[, 2]));

A374590 Numbers whose maximum exponent in their prime factorization is an evil number (A001969).

Original entry on oeis.org

8, 24, 27, 32, 40, 54, 56, 64, 72, 88, 96, 104, 108, 120, 125, 135, 136, 152, 160, 168, 184, 189, 192, 200, 216, 224, 232, 243, 248, 250, 264, 270, 280, 288, 296, 297, 312, 320, 328, 343, 344, 351, 352, 360, 375, 376, 378, 392, 408, 416, 424, 440, 448, 456, 459
Offset: 1

Views

Author

Amiram Eldar, Jul 12 2024

Keywords

Comments

The asymptotic density of this sequence is Sum_{k in A001969} (1/zeta(k+1) - 1/zeta(k)) = 0.12101890210392912747... .

Crossrefs

Subsequence of A013929 and A262675 \ {1}.
Similar sequences: A368714, A369937, A369938, A369939, A374588, A374589.

Programs

  • Mathematica
    evilQ[n_] := EvenQ[DigitCount[n, 2, 1]]; q[n_] := evilQ[Max[FactorInteger[n][[;; , 2]]]]; Select[Range[500], q]
  • PARI
    is(n) = n > 1 && !(hammingweight(vecmax(factor(n)[, 2])) % 2);

A374325 The maximal exponent in the prime factorization of the numbers whose maximal exponent in their prime factorization is a square.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Amiram Eldar, Jul 04 2024

Keywords

Comments

First differs from A369935 at n = 95.
The first occurrence of k^2, for k = 0, 1, ..., is at 1, 2, 12, 335, 42563, ..., which is the position of 2^(k^2) at A369937.

Crossrefs

Similar sequences: A374324, A374326, A374327, A374328.

Programs

  • Mathematica
    f[n_] := Module[{e = If[n == 1, 0, Max[FactorInteger[n][[;; , 2]]]]}, If[IntegerQ@ Sqrt[e], e, Nothing]]; Array[f, 200]
  • PARI
    lista(kmax) = {my(e); print1(0, ", "); for(k = 2, kmax, e = vecmax(factor(k)[, 2]); if(issquare(e), print1(e, ", ")));}

Formula

a(n) = A374326(n)^2.
a(n) = A051903(A369937(n)).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=1} k^2 * d(k) / Sum_{k>=1} d(k) = 1.19949058780036416105..., where d(k) = 1/zeta(k^2+1) - 1/zeta(k^2) for k>=2, and d(1) = 1/zeta(2).

A374326 The square root of the maximal exponent in the prime factorization of the numbers whose maximal exponent in their prime factorization is a square.

Original entry on oeis.org

0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
Offset: 1

Views

Author

Amiram Eldar, Jul 04 2024

Keywords

Comments

First differs from A369936 at n = 95.
The first occurrence of k = 0, 1, ... is at 1, 2, 12, 335, 42563, ..., which is the position of 2^(k^2) at A369937.

Crossrefs

Similar sequences: A374324, A374325, A374327, A374328.

Programs

  • Mathematica
    f[n_] := Module[{s = Sqrt[If[n == 1, 0, Max[FactorInteger[n][[;; , 2]]]]]}, If[IntegerQ[s], s, Nothing]]; Array[f, 200]
  • PARI
    lista(kmax) = {my(e); print1(0, ", "); for(k = 2, kmax, e = vecmax(factor(k)[, 2]); if(issquare(e), print1(sqrtint(e), ", ")));}

Formula

a(n) = sqrt(A374325(n)).
a(n) = sqrt(A051903(A369937(n))).
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Sum_{k>=1} k * d(k) / Sum_{k>=1} d(k) = 1.06543556163434367736..., where d(k) = 1/zeta(k^2+1) - 1/zeta(k^2) for k>=2, and d(1) = 1/zeta(2).

A374588 Numbers whose maximum exponent in their prime factorization is a composite number.

Original entry on oeis.org

16, 48, 64, 80, 81, 112, 144, 162, 176, 192, 208, 240, 256, 272, 304, 320, 324, 336, 368, 400, 405, 432, 448, 464, 496, 512, 528, 560, 567, 576, 592, 624, 625, 648, 656, 688, 704, 720, 729, 752, 768, 784, 810, 816, 832, 848, 880, 891, 912, 944, 960, 976, 1008
Offset: 1

Views

Author

Amiram Eldar, Jul 12 2024

Keywords

Comments

Subsequence of A322448 and first differs from it at n = 138: A322448(138) = 2592 = 2^5 * 3^4 is not a term of this sequence.
The asymptotic density of this sequence is d = Sum_{k composite} (1/zeta(k+1) - 1/zeta(k)) = 0.05296279266796920306... . The asymptotic density of this sequence within the nonsquarefree numbers (A013929) is d / (1 - 1/zeta(2)) = 0.13508404411123191108... .

Crossrefs

Complement of A074661 within A013929.
Subsequence of A322448 and A322449 \ {1}.
Similar sequences: A368714, A369937, A369938, A369939, A374589, A374590.

Programs

  • Maple
    filter:= proc(n) local m;
      m:= max(ifactors(n)[2][..,2]);
      m > 1 and not isprime(m)
    end proc:
    select(filter, [$1..10000]); # Robert Israel, Jul 14 2024
  • Mathematica
    Select[Range[1200], CompositeQ[Max[FactorInteger[#][[;; , 2]]]] &]
  • PARI
    iscomposite(n) = n > 1 && !isprime(n);
    is(n) = n > 1 && iscomposite(vecmax(factor(n)[, 2]));

A374589 Numbers whose maximum exponent in their prime factorization is a powerful number larger than 1.

Original entry on oeis.org

16, 48, 80, 81, 112, 144, 162, 176, 208, 240, 256, 272, 304, 324, 336, 368, 400, 405, 432, 464, 496, 512, 528, 560, 567, 592, 624, 625, 648, 656, 688, 720, 752, 768, 784, 810, 816, 848, 880, 891, 912, 944, 976, 1008, 1040, 1053, 1072, 1104, 1134, 1136, 1168, 1200
Offset: 1

Views

Author

Amiram Eldar, Jul 12 2024

Keywords

Comments

Subsequence of A130897 and first differs from it at n = 115: A130897(115) = 2592 = 2^5 * 3^4 is not a term of this sequence.
The asymptotic density of this sequence is d = Sum_{k > 1 and in A001694} (1/zeta(k+1) - 1/zeta(k)) = 0.043523813088759413253... . The asymptotic density of this sequence within A130897 is d/(1 - A262276) = 0.98744988886705430331... .

Crossrefs

Subsequence of A013929, A130897 and A372405.
Similar sequences: A368714, A369937, A369938, A369939, A374588, A374590.

Programs

  • Mathematica
    powQ[n_] := Min[FactorInteger[n][[;; , 2]]] > 1; q[n_] := powQ[Max[ FactorInteger[n][[;; , 2]] ]]; Select[Range[1200], q]
  • PARI
    ispow(n) = n > 1 && ispowerful(n);
    is(n) = n > 1 && ispow(vecmax(factor(n)[, 2]))
Showing 1-7 of 7 results.