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.

A374459 Nonsquarefree exponentially odd numbers.

Original entry on oeis.org

8, 24, 27, 32, 40, 54, 56, 88, 96, 104, 120, 125, 128, 135, 136, 152, 160, 168, 184, 189, 216, 224, 232, 243, 248, 250, 264, 270, 280, 296, 297, 312, 328, 343, 344, 351, 352, 375, 376, 378, 384, 408, 416, 424, 440, 456, 459, 472, 480, 486, 488, 512, 513, 520, 536
Offset: 1

Views

Author

Amiram Eldar, Jul 09 2024

Keywords

Comments

First differs from A301517 at n = 1213. A301517(1213) = 12500 = 2^2 * 5^5 is not an exponentially odd number.
Numbers whose exponents in their prime factorization are all odd and at least one of them is larger than 1.
All the squarefree numbers (A005117) are exponentially odd. Therefore, the sequence of exponentially odd numbers (A268335) is a disjoint union of the squarefree numbers and this sequence.
The asymptotic density of this sequence is A065463 - A059956 = 0.096515099145... .

Examples

			8 = 2^3 is a term since 3 is odd and larger than 1.
		

Crossrefs

Intersection of A013929 (or A046099) and A268335.
Subsequence of A301517.
Subsequences: A062838 \ {1}, A065036, A102838, A113850, A113852, A179671, A190011, A335988 \ {1}.

Programs

  • Mathematica
    q[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, AllTrue[e, OddQ] && ! AllTrue[e, # == 1 &]]; Select[Range[1000], q]
  • PARI
    is(k) = {my(e = factor(k)[, 2]); for(i = 1, #e, if(!(e[i] %2), return(0))); for(i = 1, #e, if(e[i] >1, return(1))); 0;}

Formula

a(n) = A268335(A374460(n)).
Sum_{n>=1} 1/a(n)^s = zeta(2*s) * (Product_{p prime} (1 + 1/p^s - 1/p^(2*s))) - zeta(s)/zeta(2*s) for s > 1.

A374536 a(n) is the least exponentially odd number that is nonsquarefree and is followed by exactly n successive exponentially odd numbers that are squarefree, or -1 if no such number exists.

Original entry on oeis.org

135, 24, 120, 27, 96, 88, 32, 40, 328, 168, 136, 104, 1288, 1161, 352, 488, 8, 783, 189, 952, 4520, 56, 11576, 67384, 5088, 1336, 35768, 16173, 53768, 80328, 128169, 28576, 247375, 208552, 2556192, 1486568, 3099368, 1653032, 910568, 7864008, 34242976, 14484152
Offset: 0

Views

Author

Amiram Eldar, Jul 11 2024

Keywords

Examples

			a(0) = 135 because 135 and 136 are successive nonsquarefree exponentially odd numbers with no squarefree number between them.
a(1) = 24 because 24 and 27 are successive nonsquarefree exponentially odd numbers with one squarefree number between them, 26.
a(2) = 120 because 120 and 125 are successive nonsquarefree exponentially odd numbers with two squarefree number between them, 122 and 123.
		

Crossrefs

Programs

  • Mathematica
    sq[k_] := Module[{e = FactorInteger[k][[;;, 2]]}, If[AnyTrue[e, EvenQ], 0, If[k == 1 || Max[e] == 1, 2, 1]]]; seq[len_, kmax_ : Infinity] := Module[{v = Table[0, {len}], c = 0, k = 1, k0 = 0, m, i = 1}, While[c < len && k < kmax, m = sq[k]; If[m > 0, If[m == 2, i++, If[k0 > 0, If[i <= len && v[[i]] == 0, c++; v[[i]] = k0]; i = 1];	k0 = k]]; k++]; v]; seq[10]
  • PARI
    issq(k) = {my(e = factor(k)[, 2]); for(i = 1, #e, if(!(e[i] % 2), return(0))); if(k == 1 || vecmax(e) == 1, 2, 1);}
    lista(len, kmax = oo) = {my(v = vector(len), c = 0, k = 1, k0 = 0, m, i = 1); while(c < len && k < kmax, m = issq(k); if(m > 0, if(m == 2, i++, if(k0 > 0, if(i <= len && v[i] == 0, c++; v[i] = k0); i = 1); k0 = k)); k++); v; }

Formula

a(n) = A268335(A374460(k)), where k is the least number such that A374459(k+1) - A374459(k) = n + 1.
Showing 1-2 of 2 results.