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.

Previous Showing 11-20 of 32 results. Next

A183085 Numbers with 61 divisors.

Original entry on oeis.org

1152921504606846976, 42391158275216203514294433201, 867361737988403547205962240695953369140625, 508021860739623365322188197652216501772434524836001
Offset: 1

Views

Author

Omar E. Pol, Jul 31 2011

Keywords

Comments

Also, 60th powers of primes.
The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime.

Crossrefs

Programs

  • PARI
    a(n)=prime(n)^60

Formula

a(n) = A000040(n)^(61-1) = A000040(n)^60.
A000005(a(n)) = 61.

A030633 Numbers with 15 divisors.

Original entry on oeis.org

144, 324, 400, 784, 1936, 2025, 2500, 2704, 3969, 4624, 5625, 5776, 8464, 9604, 9801, 13456, 13689, 15376, 16384, 21609, 21904, 23409, 26896, 29241, 29584, 30625, 35344, 42849, 44944, 55696, 58564, 59536, 60025, 68121, 71824, 75625
Offset: 1

Views

Author

Keywords

Comments

Numbers of the form p^14 (subset of A010802) or p^2*q^4 (A189988) where p and q are distinct primes. - R. J. Mathar, Mar 01 2010

Crossrefs

Programs

  • Mathematica
    Select[Range[300000],DivisorSigma[0,#]==15&] (* Vladimir Joseph Stephan Orlovsky, May 05 2011 *)
  • PARI
    is(n)=numdiv(n)==15 \\ Charles R Greathouse IV, Jun 19 2016
    
  • Python
    from math import isqrt
    from sympy import primepi, primerange, integer_nthroot
    def A030633(n):
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            kmin = kmax >> 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        def f(x): return n+x-sum(primepi(isqrt(x//p**4)) for p in primerange(integer_nthroot(x,4)[0]+1))+primepi(integer_nthroot(x,6)[0])-primepi(integer_nthroot(x,14)[0])
        return bisection(f,n,n) # Chai Wah Wu, Feb 22 2025

Formula

From Amiram Eldar, Jul 03 2022: (Start)
A000005(a(n)) = 15.
Sum_{n>=1} 1/a(n) = P(2)*P(4) - P(6) + P(14) = 0.0178111..., where P is the prime zeta function. (End)

A261700 Numbers with 101 divisors.

Original entry on oeis.org

1267650600228229401496703205376, 515377520732011331036461129765621272702107522001, 7888609052210118054117285652827862296732064351090230047702789306640625, 3234476509624757991344647769100216810857203198904625400933895331391691459636928060001
Offset: 1

Views

Author

Omar E. Pol, Aug 28 2015

Keywords

Comments

Also, 100th powers of primes.
The n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime.

Examples

			a(1) = 2^100, a(2) = 3^100, a(3) = 5^100, a(4) = 7^100.
		

Crossrefs

Programs

Formula

a(n) = A000040(n)^(101-1) = A000040(n)^100.
A000005(a(n)) = 101.

A319075 Square array T(n,k) read by antidiagonal upwards in which row n lists the n-th powers of primes, hence column k lists the powers of the k-th prime, n >= 0, k >= 1.

Original entry on oeis.org

1, 2, 1, 4, 3, 1, 8, 9, 5, 1, 16, 27, 25, 7, 1, 32, 81, 125, 49, 11, 1, 64, 243, 625, 343, 121, 13, 1, 128, 729, 3125, 2401, 1331, 169, 17, 1, 256, 2187, 15625, 16807, 14641, 2197, 289, 19, 1, 512, 6561, 78125, 117649, 161051, 28561, 4913, 361, 23, 1, 1024, 19683, 390625, 823543, 1771561, 371293
Offset: 0

Views

Author

Omar E. Pol, Sep 09 2018

Keywords

Comments

If n = p - 1 where p is prime, then row n lists the numbers with p divisors.
The partial sums of column k give the column k of A319076.

Examples

			The corner of the square array is as follows:
         A000079 A000244 A000351  A000420    A001020    A001022     A001026
A000012        1,      1,      1,       1,         1,         1,          1, ...
A000040        2,      3,      5,       7,        11,        13,         17, ...
A001248        4,      9,     25,      49,       121,       169,        289, ...
A030078        8,     27,    125,     343,      1331,      2197,       4913, ...
A030514       16,     81,    625,    2401,     14641,     28561,      83521, ...
A050997       32,    243,   3125,   16807,    161051,    371293,    1419857, ...
A030516       64,    729,  15625,  117649,   1771561,   4826809,   24137569, ...
A092759      128,   2187,  78125,  823543,  19487171,  62748517,  410338673, ...
A179645      256,   6561, 390625, 5764801, 214358881, 815730721, 6975757441, ...
...
		

Crossrefs

Other rows n: A030635 (n=16), A030637 (n=18), A137486 (n=22), A137492 (n=28), A139571 (n=30), A139572 (n=36), A139573 (n=40), A139574 (n=42), A139575 (n=46), A173533 (n=52), A183062 (n=58), A183085 (n=60), A261700 (n=100).
Main diagonal gives A093360.
Second diagonal gives A062457.
Third diagonal gives A197987.
Removing the 1's we have A182944/ A182945.

Programs

  • PARI
    T(n, k) = prime(k)^n;

Formula

T(n,k) = A000040(k)^n, n >= 0, k >= 1.

A166546 Natural numbers n such that d(n) + 1 is prime.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 17, 18, 19, 20, 21, 22, 23, 26, 27, 28, 29, 31, 32, 33, 34, 35, 37, 38, 39, 41, 43, 44, 45, 46, 47, 48, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 67, 68, 69, 71, 72, 73, 74, 75, 76, 77, 79, 80, 82, 83, 84, 85, 86, 87, 89, 90
Offset: 1

Views

Author

Giovanni Teofilatto, Oct 16 2009

Keywords

Comments

Natural numbers n such that d(d(n)+1)= 2. - Giovanni Teofilatto, Oct 26 2009
The complement is the union of A001248, A030514, A030516, A030626, A030627, A030629, A030631, A030632, A030633 etc. - R. J. Mathar, Oct 26 2009

Crossrefs

Cf. A000005.
Cf. A073915. - R. J. Mathar, Oct 26 2009

Programs

  • Magma
    [n: n in [1..100] | IsPrime(NumberOfDivisors(n)+1)]; // Vincenzo Librandi, Jan 20 2019
  • Mathematica
    Select[Range@90, PrimeQ[DivisorSigma[0, #] + 1] &] (* Vincenzo Librandi, Jan 20 2019 *)
  • PARI
    isok(n) = isprime(numdiv(n)+1); \\ Michel Marcus, Jan 20 2019
    

Formula

{1} U A000040 U A030513 U A030515 U A030628 U A030630 U A030634 U A030636 U A137485 U A137491 U A137493 U ... . - R. J. Mathar, Oct 26 2009

A280298 Numbers with 67 divisors.

Original entry on oeis.org

73786976294838206464, 30903154382632612361920641803529, 13552527156068805425093160010874271392822265625, 59768263894155949306790119265585619217025149412430681649, 539407797827634189900210968137750826278309533633974732577186113975161
Offset: 1

Views

Author

Omar E. Pol, Dec 31 2016

Keywords

Comments

Also, 66th powers of primes.
More generally, the n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. In this case, p = 67.

Examples

			a(1) = 2^66, a(2) = 3^66, a(3) = 5^66, a(4) = 7^66, a(5) = 11^66.
		

Crossrefs

Programs

  • Mathematica
    Array[Prime[#]^66 &, {5}] (* Michael De Vlieger, Dec 31 2016 *)
  • PARI
    a(n)=prime(n)^66

Formula

a(n) = A000040(n)^(67-1) = A000040(n)^66.
A000005(a(n)) = 67.

A280299 Numbers with 71 divisors.

Original entry on oeis.org

1180591620717411303424, 2503155504993241601315571986085849, 8470329472543003390683225006796419620513916015625, 143503601609868434285603076356671071740077383739246066639249, 7897469567994392174328988784504809847540729881935024059662581894710332201
Offset: 1

Views

Author

Omar E. Pol, Dec 31 2016

Keywords

Comments

Also, 70th powers of primes.
More generally, the n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. In this case, p = 71.

Examples

			a(1) = 2^70, a(2) = 3^70, a(3) = 5^70, a(4) = 7^70, a(5) = 11^70.
		

Crossrefs

Programs

  • Mathematica
    Array[Prime[#]^70 &, {5}] (* Michael De Vlieger, Dec 31 2016 *)
  • PARI
    a(n)=prime(n)^70

Formula

a(n) = A000040(n)^(71-1) = A000040(n)^70.
A000005(a(n)) = 71.

A280301 Numbers with 73 divisors.

Original entry on oeis.org

4722366482869645213696, 22528399544939174411840147874772641, 211758236813575084767080625169910490512847900390625, 7031676478883553279994550741476882515263791803223057265323201, 955593817727321453093807642925081991552428315714137911219172409259950196321
Offset: 1

Views

Author

Omar E. Pol, Dec 31 2016

Keywords

Comments

Also, 72nd powers of primes.
More generally, the n-th number with p divisors is equal to the n-th prime raised to power p-1, where p is prime. In this case, p = 73.

Examples

			a(1) = 2^72, a(2) = 3^72, a(3) = 5^72, a(4) = 7^72, a(5) = 11^72.
		

Crossrefs

Programs

  • Mathematica
    Array[Prime[#]^72 &, {5}] (* Michael De Vlieger, Dec 31 2016 *)
  • PARI
    a(n)=prime(n)^72

Formula

a(n) = A000040(n)^(73-1) = A000040(n)^72.
A000005(a(n)) = 73.

A377654 Numbers m^2 for which the center part (containing the diagonal) of its symmetric representation of sigma, SRS(m^2), has width 1 and area m.

Original entry on oeis.org

1, 9, 25, 49, 81, 121, 169, 289, 361, 441, 529, 625, 729, 841, 961, 1089, 1369, 1521, 1681, 1849, 2209, 2401, 2601, 2809, 3025, 3249, 3481, 3721, 4225, 4489, 4761, 5041, 5329, 6241, 6561, 6889, 7225, 7569, 7921, 8649, 9025, 9409, 10201, 10609, 11449, 11881, 12321, 12769, 13225, 14161, 14641, 15129, 15625
Offset: 1

Views

Author

Hartmut F. W. Hoft, Nov 03 2024

Keywords

Comments

Since for numbers m^2 in the sequence the width at the diagonal of SRS(m^2) is 1, the area m of its center part is odd so that this sequence is a proper subsequence of A016754 and since SRS(m^2) has an odd number of parts it is a proper subsequence of A319529. The smallest odd square not in this sequence is 225 = 15^2. SRS(225) is {113, 177, 113}, its center part has maximum width 2, its width at the diagonal is 1.
The k+1 parts of SRS(p^(2k)), p an odd prime and k >= 0, through the diagonal including the center part have areas (p^(2k-i) + p^i)/2 for 0 <= i <= k. They form a strictly decreasing sequence. Since p^(2k) has 2k+1 divisors and SRS(p^(2k)) has 2k+1 parts, all of width 1 (A357581), the even powers of odd primes form a proper subsequence of A244579. For the subsequence of squares of odd primes p, SRS(p^2) consists of the 3 parts { (p^2 + 1)/2, p, (p^2 + 1)/2 } see A001248, A247687 and A357581.
The areas of the parts of SRS(m^2) need not be in descending order through the diagonal as a(112) = 275^2 = 75625 with SRS(75625) = (37813, 7565, 3443, 1525, 715, 738, 275, 738, 715, 1525, 3443, 7565, 37813) demonstrates.
An equivalent description of the sequence is: The center part of SRS(m^2) has width 1, m is odd, and A249223(m^2, m-1) = 0.
Conjectures (true for all a(n) <= 10^8):
(1) The central part of SRS(a(n)) is the minimum of all parts of SRS(a(n)), 1 <= n.
(2) The terms in this sequence are the squares of the terms in A244579.

Examples

			The center part of SRS(a(3)) = SRS(25) has area 5, all 3 parts have width 1, and 25 with 3 divisors also belongs to A244579.
The center part of SRS(a(7)) = SRS(169) has area 13, all 3 parts have width 1, and 169 with 3 divisors also belongs to A244579.
The center part of SRS(a(10)) = SRS(441) has area 21 and width 1, but the maximum width of SRS(441) is 2. Number 441 has 9 divisors and SRS(441) has 7 parts while 21 has 4 divisors and SRS(21) has 4 parts so that 21 is in A244579 while 441 is not.
		

Crossrefs

Programs

  • Mathematica
    (* t237591 and partsSRS compute rows in A237270 and A237591, respectively *)
    (* t249223 and widthPattern are also defined in A376829 *)
    row[n_] := Floor[(Sqrt[8 n+1]-1)/2]
    t237591[n_] := Map[Ceiling[(n+1)/#-(#+1)/2]-Ceiling[(n+1)/(#+1)-(#+2)/2]&, Range[row[n]]]
    partsSRS[n_] := Module[{widths=t249223[n], legs=t237591[n], parts, srs}, parts=widths legs; srs=Map[Apply[Plus, #]&, Select[SplitBy[Join[parts, Reverse[parts]], #!=0&], First[#]!=0&]]; srs[[Ceiling[Length[srs]/2]]]-=Last[widths]; srs]
    t249223[n_] := FoldList[#1+(-1)^(#2+1)KroneckerDelta[Mod[n-#2 (#2+1)/2, #2]]&, 1, Range[2, row[n]]]
    widthPattern[n_] := Map[First, Split[Join[t249223[n], Reverse[t249223[n]]]]]
    centerQ[n_] := Module[{pS=partsSRS[n]}, Sqrt[n]==pS[[(Length[pS]+1)/2]]]/;OddQ[n]
    widthQ[n_] := Module[{wP=SplitBy[widthPattern[n], #!=0&]}, wP[[(Length[wP]+1)/2]]]=={1}/;OddQ[n]
    a377654[m_, n_] := Select[Map[#^2&, Range[m, n, 2]], centerQ[#]&&widthQ[#]&]/;OddQ[m]
    a377654[1, 125]

A259417 Even powers of the odd primes listed in increasing order.

Original entry on oeis.org

1, 9, 25, 49, 81, 121, 169, 289, 361, 529, 625, 729, 841, 961, 1369, 1681, 1849, 2209, 2401, 2809, 3481, 3721, 4489, 5041, 5329, 6241, 6561, 6889, 7921, 9409, 10201, 10609, 11449, 11881, 12769, 14641, 15625, 16129, 17161, 18769, 19321, 22201, 22801, 24649
Offset: 1

Views

Author

Hartmut F. W. Hoft, Jun 26 2015

Keywords

Comments

Each of the following sequences, p^(q-1) with p >= 2 and q > 2 primes, except their respective first elements, powers of 2, is a subsequence:
A001248(p) = p^2, A030514(p) = p^4, A030516(p) = p^6,
A030629(p) = p^10, A030631(p) = p^12, A030635(p) = p^16,
A030637(p) = p^18, A137486(p) = p^22, A137492(p) = p^28,
A139571(p) = p^30, A139572(p) = p^36, A139573(p) = p^40,
A139574(p) = p^42, A139575(p) = p^46, A173533(p) = p^52,
A183062(p) = p^58, A183085(p) = p^60.
See also the link to the OEIS Wiki.
The sequences A053182(n)^2, A065509(n)^4, A163268(n)^6 and A240693(n)^10 are subsequences of this sequence.
The odd numbers in A023194 are a subsequence of this sequence.

Examples

			a(11) = 5^4 = 625 is followed by a(12) = 3^6 = 729 since no even power of an odd prime falls between them.
		

Programs

  • Mathematica
    a259417[bound_] := Module[{q, h, column = {}}, For[q = Prime[2], q^2 <= bound, q = NextPrime[q], For[h = 1, q^(2*h) <= bound, h++, AppendTo[column, q^(2*h)]]]; Prepend[Sort[column], 1]]
    a259417[25000] (* data *)
    With[{upto=25000},Select[Union[Flatten[Table[Prime[Range[2,Floor[ Sqrt[ upto]]]]^n,{n,0,Log[2,upto],2}]]],#<=upto&]] (* Harvey P. Dale, Nov 25 2017 *)

Formula

Sum_{n>=1} 1/a(n) = 1 + Sum_{k>=1} (P(2*k) - 1/2^(2*k)) = 1.21835996432366585110..., where P is the prime zeta function. - Amiram Eldar, Jul 10 2022
Previous Showing 11-20 of 32 results. Next