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

A321147 Odd exponential abundant numbers: odd numbers k whose sum of exponential divisors A051377(k) > 2*k.

Original entry on oeis.org

225450225, 385533225, 481583025, 538472025, 672624225, 705699225, 985646025, 1121915025, 1150227225, 1281998025, 1566972225, 1685513025, 1790559225, 1826280225, 2105433225, 2242496025, 2466612225, 2550755025, 2679615225, 2930852925, 2946861225, 3132081225
Offset: 1

Views

Author

Amiram Eldar, Oct 28 2018

Keywords

Comments

From Amiram Eldar, Jun 08 2020: (Start)
Exponential abundant numbers that are odd are relatively rare: there are 235290 even exponential abundant number smaller than the first odd term, i.e., a(1) = A129575(235291).
Odd exponential abundant numbers k such that k-1 or k+1 is also exponential abundant number exist (e.g. (73#/5#)^2-1 and (73#/5#)^2 are both exponential abundant numbers, where prime(k)# = A002110(k)). Which pair is the least?
The least exponential abundant number that is coprime to 6 is (31#/3#)^2 = 1117347505588495206025. In general, the least exponential abundant number that is coprime to A002110(k) is (A007708(k+1)#/A002110(k))^2. (End)
The asymptotic density of this sequence is Sum_{n>=1} f(A328136(n)) = 5.29...*10^(-9), where f(n) = (4/(Pi^2*n))*Product_{prime p|n}(p/(p+1)) if n is odd and 0 otherwise. - Amiram Eldar, Sep 02 2022

Examples

			225450225 is in the sequence since it is odd and A051377(225450225) = 484323840 > 2 * 225450225.
		

Crossrefs

The exponential version of A005231.
The odd subsequence of A129575.

Programs

  • Mathematica
    esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; s={};Do[If[esigma[n]>2n,AppendTo[s,n]],{n,1,10^10,2}]; s

A318167 Numbers k such that both k and k+1 are bi-unitary abundant numbers.

Original entry on oeis.org

21735, 21944, 43064, 49664, 58695, 76544, 106784, 135135, 144584, 160544, 188055, 209055, 227744, 256095, 262184, 300104, 345344, 348704, 382304, 387584, 407295, 409184, 414855, 437535, 498015, 520695, 560384, 567944, 611415, 679455, 687015, 705375, 709695
Offset: 1

Views

Author

Amiram Eldar, Aug 20 2018

Keywords

Comments

The bi-unitary version of A096399.

Examples

			21735 is in the sequence since both 21735 and 21736 are bi-unitary abundant numbers.
		

Crossrefs

Cf. A096399 (analog for sigma), A188999 (bi-unitary sigma).
Cf. A292982 (bi-unitary abundant), A293186 (odd bi-unitary abundant).

Programs

  • Mathematica
    f[n_] := Select[Divisors[n], Function[d, CoprimeQ[d, n/d]]]; bsigma[m_] := DivisorSum[m, # &, Last@Intersection[f@#, f[m/#]] == 1 &]; bAbundantQ[n_] := bsigma[n] > 2 n; seq={}; n=1; While[Length[seq]<32,If[bAbundantQ[n] && bAbundantQ [n+1],AppendTo[seq,n]];n++];seq
  • PARI
    a188999(n) = {f = factor(n); for (i=1, #f~, p = f[i, 1]; e = f[i, 2]; f[i, 1] = if (e % 2, (p^(e+1)-1)/(p-1), (p^(e+1)-1)/(p-1) -p^(e/2)); f[i, 2] = 1; ); factorback(f); }
    isok(n) = (a188999(n) > 2*n) && (a188999(n+1) > 2*(n+1)); \\ Michel Marcus, Aug 21 2018

A348275 Odd noninfinitary abundant numbers: the odd terms of A348274.

Original entry on oeis.org

99225, 1091475, 1289925, 1334025, 1576575, 1686825, 1715175, 1863225, 1885275, 2027025, 2061675, 2282175, 2304225, 2395575, 2401245, 2436525, 2480625, 2650725, 2723175, 2789325, 2877525, 2962575, 3031875, 3075975, 3132675, 3185325, 3186225, 3296475, 3353805, 3501225
Offset: 1

Views

Author

Amiram Eldar, Oct 09 2021

Keywords

Comments

The number of terms below 10^k, for k = 5, 6, ..., are 1, 113, 630, 7771, 73685, ... Apparently this sequence has an asymptotic density 0.000007...

Examples

			99225 is a term since A348271(99225) = 107207 > 99225.
		

Crossrefs

Cf. A348271.
Subsequence of A005231 and A348274.
Similar sequences: A094889, A127666, A129485, A293186, A321147.

Programs

  • Mathematica
    f[p_, e_] := Module[{b = IntegerDigits[e, 2], m}, m = Length[b]; Product[If[b[[j]] > 0, 1 + p^(2^(m - j)), 1], {j, 1, m}]]; isigma[1] = 1; isigma[n_] := Times @@ f @@@ FactorInteger[n]; s[n_] := DivisorSigma[1,n] - isigma[n]; Select[Range[1, 2*10^6, 2], s[#] > # &]

A333950 Odd recursive abundant numbers: odd numbers k such that A333926(k) > 2*k.

Original entry on oeis.org

1575, 2205, 3465, 4095, 5355, 5775, 5985, 6435, 6825, 7245, 8085, 8415, 8925, 9135, 9555, 9765, 11025, 11655, 12705, 12915, 13545, 14805, 15015, 16695, 17325, 18585, 19215, 19635, 20475, 21105, 21945, 22365, 22995, 23205, 24255, 24885, 25935, 26145, 26565, 26775
Offset: 1

Views

Author

Amiram Eldar, Apr 11 2020

Keywords

Examples

			1575 is a term since it is odd and A333926(1575) = 3224 > 2 * 1575.
		

Crossrefs

Intersection of A005408 and A333928.
Cf. A333926.
Analogous sequences: A005231, A094889 (nonunitary), A129485 (unitary), A127666 (infinitary), A293186 (bi-unitary), A321147 (exponential).

Programs

  • Mathematica
    recDivQ[n_, 1] = True; recDivQ[n_, d_] := recDivQ[n, d] = Divisible[n, d] && AllTrue[FactorInteger[d], recDivQ[IntegerExponent[n, First[#]], Last[#]] &]; recDivs[n_] := Select[Divisors[n], recDivQ[n, #] &]; f[p_, e_] := 1 + Total[p^recDivs[e]]; recDivSum[1] = 1; recDivSum[n_] := Times @@ (f @@@ FactorInteger[n]); Select[2*Range[15000] + 1, recDivSum[#] > 2*# &]

A360526 Odd numbers k such that A360522(k) > 2*k.

Original entry on oeis.org

15015, 19635, 21945, 23205, 25935, 26565, 31395, 33495, 33915, 35805, 39585, 41055, 42315, 42735, 45885, 47355, 49665, 50505, 51765, 54285, 55965, 58695, 61215, 64155, 68145, 70455, 72345, 77385, 80535, 82005, 83265, 84315, 91245, 95865, 102795, 112035, 116655
Offset: 1

Views

Author

Amiram Eldar, Feb 10 2023

Keywords

Comments

First differs from A112643, A129485, A249263 at n = 46: a(46) = 165165 is not a term of these sequences.

Examples

			15015 is a term since A360522(15015) = 32256 > 2*15015.
		

Crossrefs

Cf. A360522.
Subsequence of A005101, A005231 and A360525.

Programs

  • Mathematica
    f[p_, e_] := p^e + e; s[1] = 1; s[n_] := Times @@ f @@@ FactorInteger[n]; q[n_] := s[n] > 2*n; Select[Range[1, 10^5, 2], q]
  • PARI
    isab(n) = { my(f = factor(n)); prod(i = 1, #f~, f[i,1]^f[i,2] + f[i,2]) > 2*n;}
    is(n) = n%2 && isab(n);

A335053 Odd bi-unitary abundant numbers whose bi-unitary abundancy is closer to 2 than that of any smaller odd bi-unitary abundant number.

Original entry on oeis.org

945, 25515, 46035, 49875, 83265, 354585, 359205, 361515, 366135, 382305, 389235, 396165, 400785, 403095, 407715, 414645, 416955, 423885, 430815, 437745, 442365, 13351635, 132335385, 159030135, 1756753845, 6561644355, 10394173335, 13455037365, 37456183215
Offset: 1

Views

Author

Amiram Eldar, May 21 2020

Keywords

Comments

The bi-unitary abundancy of a number k is bsigma(k)/k, where bsigma(k) is the sum of bi-unitary divisors of k (A188999).

Examples

			The bi-unitary abundancies of the first terms are 2.031..., 2.005..., 2.0019..., 2.0018..., 2.0015..., ...
		

Crossrefs

Programs

  • Mathematica
    fun[p_, e_] := If[OddQ[e], (p^(e + 1) - 1)/(p - 1), (p^(e + 1) - 1)/(p - 1) - p^(e/2)]; bsigma[1] = 1; bsigma[n_] := Times @@ (fun @@@ FactorInteger[n]); seq = {}; r = 3; Do[s = bsigma[n]/n; If[s > 2 && s < r, AppendTo[seq, n]; r = s], {n, 1, 10^6, 2}]; seq

A339936 Odd coreful abundant numbers: the odd terms of A308053.

Original entry on oeis.org

99225, 165375, 231525, 297675, 496125, 694575, 826875, 893025, 1091475, 1157625, 1225125, 1289925, 1488375, 1620675, 1686825, 1819125, 1885275, 2083725, 2149875, 2282175, 2480625, 2546775, 2679075, 2811375, 2877525, 3009825, 3075975, 3142125, 3274425, 3472875
Offset: 1

Views

Author

Amiram Eldar, Dec 23 2020

Keywords

Comments

The asymptotic density of this sequence is Sum_{n>=1} f(A356871(n)) = 9.1348...*10^(-6), where f(n) = (4/(Pi^2*n))*Product_{prime p|n}(p/(p+1)) if n is odd and 0 otherwise. - Amiram Eldar, Sep 02 2022

Examples

			99225 is a term since it is odd and the sum of its coreful divisors is A057723(99225) = 201600 > 2 * 99225.
		

Crossrefs

Intersection of A005408 and A308053.
Subsequence of A321147.

Programs

  • Mathematica
    f[p_, e_] := (p^(e + 1) - 1)/(p - 1) - 1; s[1] = 1; s[n_] := Times @@ (f @@@ FactorInteger[n]); Select[Range[1, 10^6, 2], s[#] > 2*# &]

A339938 Odd non-coreful abundant numbers: the odd terms of A308127.

Original entry on oeis.org

15015, 19635, 21945, 23205, 25935, 26565, 31395, 33495, 33915, 35805, 39585, 41055, 42315, 42735, 45885, 47355, 49665, 50505, 51765, 54285, 55965, 58695, 61215, 64155, 68145, 70455, 72345, 75075, 77385, 80535, 82005, 83265, 84315, 91245, 95865, 102795, 105105
Offset: 1

Views

Author

Amiram Eldar, Dec 23 2020

Keywords

Comments

First differs from A112643, A129485 and A249263 at n = 28.

Examples

			15015 is a term since it is odd and the sum of its non-coreful divisors is A308135(15015) = 17241 > 15015.
		

Crossrefs

Intersection of A005408 and A308127.
Cf. A308135.

Programs

  • Mathematica
    f[p_, e_] := (p^(e + 1) - 1)/(p - 1); fc[p_, e_] := f[p, e] - 1; s[1] = 0; s[n_] := Times @@ (f @@@ FactorInteger[n]) - Times @@ (fc @@@ FactorInteger[n]); Select[Range[1, 10^5, 2], s[#] > # &]

A348605 Odd nonexponential abundant numbers: odd numbers k such that A160135(k) > k.

Original entry on oeis.org

8505, 10395, 12285, 15015, 16065, 17955, 19635, 21735, 21945, 23205, 25515, 25935, 26565, 28875, 31185, 31395, 33495, 33915, 34125, 35805, 36855, 39585, 41055, 42315, 42735, 45885, 47355, 48195, 49665, 50505, 51765, 53865, 54285, 55965, 56595, 58695, 61215, 64155
Offset: 1

Views

Author

Amiram Eldar, Oct 25 2021

Keywords

Comments

The odd terms of A348604.
The numbers of terms not exceeding 10^k, for k = 4, 5, ..., are 1, 51, 360, 4117, 39803, 418663, 4099004, ... Apparently this sequence has an asymptotic density 0.0004...

Examples

			8505 is a term since A160135(8505) = 8862 > 8505.
		

Crossrefs

Cf. A160135.
Subsequence of A005231 and A348604.
Similar sequences: A094889, A127666, A129485, A293186, A321147, A348275.

Programs

  • Mathematica
    esigma[n_] := Times @@ (Sum[First[#]^d, {d, Divisors[Last[#]]}] &) /@ FactorInteger[n]; q[n_] := DivisorSigma[1, n] - esigma[n] > n; Select[Range[1, 65000, 2], q]

A379031 Odd modified exponential abundant numbers: odd numbers k such that A241405(k) > 2*k.

Original entry on oeis.org

15015, 19635, 21945, 23205, 25935, 26565, 31395, 33495, 33915, 35805, 39585, 41055, 42315, 42735, 45885, 47355, 49665, 50505, 51765, 54285, 55965, 58695, 61215, 64155, 68145, 70455, 72345, 77385, 80535, 82005, 83265, 84315, 91245, 95865, 102795, 112035, 116655
Offset: 1

Views

Author

Amiram Eldar, Dec 14 2024

Keywords

Comments

First differs from its subsequences A112643 and A249263 at n = 51: a(51) = 195195 is not a term of these two sequences.
First differs from its subsequence A129485 at n = 363: a(363) = 2537535 is not a term of A129485.
First differs from A339938 at n = 28: A339938(28) = 75075 is not a term of this sequence.
First differs from A360526 at n = 46: A360526(46) = 165165 is not a term of this sequence.

Crossrefs

Intersection of A005408 and A379029.
Subsequence of A005231.
Subsequences: A112643, A129485, A249263.
Cf. A241405.

Programs

  • Mathematica
    f[p_, e_] := DivisorSum[e + 1, p^(# - 1) &]; mesigma[1] = 1; mesigma[n_] := Times @@ f @@@ FactorInteger[n]; meAbQ[n_] := mesigma[n] > 2*n; Select[Range[1, 10^5, 2], meAbQ]
  • PARI
    is(k) = if(!(k%2), 0, my(f=factor(k)); prod(i=1, #f~, sumdiv(f[i, 2]+1, d, f[i, 1]^(d-1))) > 2*k);
Showing 1-10 of 14 results. Next