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

A067128 Ramanujan's largely composite numbers, defined to be numbers m such that d(m) >= d(k) for k = 1 to m-1.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 10, 12, 18, 20, 24, 30, 36, 48, 60, 72, 84, 90, 96, 108, 120, 168, 180, 240, 336, 360, 420, 480, 504, 540, 600, 630, 660, 672, 720, 840, 1080, 1260, 1440, 1680, 2160, 2520, 3360, 3780, 3960, 4200, 4320, 4620, 4680, 5040, 7560, 9240
Offset: 1

Views

Author

Amarnath Murthy, Jan 09 2002

Keywords

Comments

This sequence is a subsequence of A034287; are they identical? They match for m up to 1500000.
Identical to A034287 for the 105834 terms less than 10^150.
Every subsequence of terms, having the fixed greatest prime divisor prime(k), k=1,2,..., is finite. For a proof see A273015. The list of these subsequences begins {2,4,8}, {3,6,12,18,24,36,48,72,96,108}, ... - Vladimir Shevelev, May 13 2016
By a result of Erdős (1944), a(n+1) <= 2*a(n): see Erdős link. - David A. Corneth, May 20 2016
It appears that if n > 13, then a(n) = A363658(n). - Simon Jensen, Aug 31 2023
Out of the first 10000 terms of this sequence, 1766 are adjacent to a prime. - Dmitry Kamenetsky, Jul 02 2024

Examples

			8 is a term as d(8) = 4 and d(k) <= 4 for k = 1,...,7.
		

Crossrefs

For n with strictly increasing number of divisors, see A002182; A272314, A273011 (infinitary analog), subsequences A273015, A273016, A273018.
Number of divisors of a(n): A273353.

Programs

  • Maple
    isA067128 := proc(n)
        local nd,k ;
        nd := numtheory[tau](n) ;
        for k from 1 to n-1 do
            if numtheory[tau](k) > nd then
                return false ;
            end if;
        end do:
        true ;
    end proc:
    A067128 := proc(n)
        option remember;
        if n = 1 then
            1;
        else
            for a from procname(n-1)+1 do
                if isA067128(a) then
                    return a;
                end if;
            end do:
        end if;
    end proc:
    seq(A067128(n),n=1..60) ; # R. J. Mathar, Apr 15 2024
  • Mathematica
    For[n=1; max=0, True, n++, If[(d=DivisorSigma[0, n])>=max, Print[n]; max=d]]
    NestList[Function[last,
      NestWhile[# + 1 &, last + 1,
       DivisorSigma[0, #] < DivisorSigma[0, last] &]], 1, 70] (* Steven Lu, Nov 28 2022 *)
  • PARI
    is(n) = my(nd=numdiv(n)); for(k=1, n-1, if(numdiv(k) > nd, return(0))); return(1) \\ Felix Fröhlich, May 22 2016

Extensions

Edited by Dean Hickerson, Jan 15 2002 and by T. D. Noe, Nov 07 2002

A217854 Product of all divisors of n, positive or negative.

Original entry on oeis.org

-1, 4, 9, -64, 25, 1296, 49, 4096, -729, 10000, 121, 2985984, 169, 38416, 50625, -1048576, 289, 34012224, 361, 64000000, 194481, 234256, 529, 110075314176, -15625, 456976, 531441, 481890304, 841, 656100000000, 961, 1073741824, 1185921
Offset: 1

Views

Author

Keywords

Comments

a(n) is negative iff n is a square.

Crossrefs

Programs

  • Mathematica
    a[n_] := (-n)^DivisorSigma[0, n]; Array[a, 33] (* Amiram Eldar, Aug 31 2023 *)
  • PARI
    a(n) = (-n)^numdiv(n); \\ Michel Marcus, Aug 31 2023

Formula

a(n) = (-n)^tau(n) = (-n)^A000005(n).
a(n) = A158387(n) * A062758(n). - Andrew Howroyd, Aug 31 2023

A363657 Numbers m where A217854(m) is a record minimum.

Original entry on oeis.org

1, 4, 9, 16, 36, 100, 144, 324, 400, 576, 900, 1764, 2304, 3600, 7056, 8100, 14400, 28224, 32400, 44100, 57600, 108900, 112896, 129600, 176400, 396900, 435600, 518400, 608400, 705600, 1587600, 2822400, 5336100, 6350400, 14288400, 15681600, 17640000, 21344400
Offset: 1

Views

Author

Simon Jensen, Jun 13 2023

Keywords

Comments

(-m)^tau(m) < 0 and (-m)^tau(m) < (-k)^tau(k) for all positive k < m, where tau is the number of divisors function.
All terms are squares.
It is conjectured that if m is a term, then abs((-m)^tau(m)) <= abs((-k)^tau(k)) for some k < m. See the link.

Examples

			9 is a term since (-9)^tau(9) = (-9)^3 = -729 and -729 < (-k)^tau(k) for k = 1,...,8.
25 is not a term since (-25)^tau(5) = (-25)^3 = -15625 > (-16)^tau(16) = (-16)^5 = -1048576 and 16 < 25.
		

Crossrefs

Programs

  • PARI
    isok(m) = my(x=(-m)^numdiv(m)); for (k=1, m-1, if (x >= (-k)^numdiv(k), return(0))); return(1); \\ Michel Marcus, Jun 18 2023
Showing 1-3 of 3 results.