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 21-30 of 39 results. Next

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]

A190464 Numbers with prime factorization p^4*q^6.

Original entry on oeis.org

5184, 11664, 40000, 153664, 250000, 455625, 937024, 1265625, 1750329, 1827904, 1882384, 5345344, 8340544, 9529569, 10673289, 17909824, 20820969, 28344976, 37515625, 45265984, 59105344, 60886809, 73530625, 77228944, 95004009, 119946304, 143496441, 180848704, 204004089, 218803264
Offset: 1

Views

Author

Keywords

Comments

A subsequence of A175745 (Numbers with 35 divisors).
First different term in A175745 is 17179869184(=2^34).

Crossrefs

Programs

  • Mathematica
    f[n_]:=Sort[Last/@FactorInteger[n]]=={4,6}; Select[Range[50000000],f] (*and*) lst={};Do[Do[If[n!=m,AppendTo[lst,Prime[n]^6*Prime[m]^4]], {n,50}],{m,50}]; Take[Union@lst,50]
  • PARI
    list(lim)=my(v=List(),t);forprime(p=2, (lim\16)^(1/6), t=p^6;forprime(q=2, (lim\t)^(1/4), if(p==q, next);listput(v,t*q^4))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 24 2011

Formula

Sum_{n>=1} 1/a(n) = P(4)*P(6) - P(10) = A085964 * A085966 - P(10) = 0.000320..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

A275387 Numbers of ordered pairs of divisors d < e of n such that gcd(d, e) > 1.

Original entry on oeis.org

0, 0, 0, 1, 0, 2, 0, 3, 1, 2, 0, 8, 0, 2, 2, 6, 0, 8, 0, 8, 2, 2, 0, 18, 1, 2, 3, 8, 0, 15, 0, 10, 2, 2, 2, 24, 0, 2, 2, 18, 0, 15, 0, 8, 8, 2, 0, 32, 1, 8, 2, 8, 0, 18, 2, 18, 2, 2, 0, 44, 0, 2, 8, 15, 2, 15, 0, 8, 2, 15, 0, 49, 0, 2, 8, 8, 2, 15, 0, 32, 6, 2
Offset: 1

Views

Author

Michel Lagneau, Aug 03 2016

Keywords

Comments

Number of elements in the set {(x, y): x|n, y|n, x < y, gcd(x, y) > 1}.
Every element of the sequence is repeated indefinitely, for instance:
a(n)=0 if n prime;
a(n)=1 if n = p^2 for p prime (A001248);
a(n)=2 if n is a squarefree semiprime (A006881);
a(n)=3 if n = p^3 for p prime (A030078);
a(n)=6 if n = p^4 for p prime (A030514);
a(n)=8 if n is a number which is the product of a prime and the square of a different prime (A054753);
a(n)=10 if n = p^5 for p prime (A050997);
a(n)=15 if n is in the set {A007304} union {64} = {30, 42, 64, 66, 70,...} = {Sphenic numbers} union {64};
a(n)=18 if n is the product of the cube of a prime (A030078) and a different prime (see A065036);
a(n)=21 if n = p^7 for p prime (A092759);
a(n)=24 if n is square of a squarefree semiprime (A085986);
a(n)=32 if n is the product of the 4th power of a prime (A030514) and a different prime (see A178739);
a(n)=36 if n = p^9 for p prime (A179665);
a(n)=44 if n is the product of exactly four primes, three of which are distinct (A085987);
a(n)=45 if n is a number with 11 divisors (A030629);
a(n)=49 if n is of the form p^2*q^3, where p,q are distinct primes (A143610);
a(n)=50 if n is the product of the 5th power of a prime (A050997) and a different prime (see A178740);
a(n)=55 if n if n = p^11 for p prime(A079395);
a(n)=72 if n is a number with 14 divisors (A030632);
a(n)=80 if n is the product of four distinct primes (A046386);
a(n)=83 if n is a number with 15 divisors (A030633);
a(n)=89 if n is a number with prime factorization pqr^3 (A189975);
a(n)=96 if n is a number that are the cube of a product of two distinct primes (A162142);
a(n)=98 if n is the product of the 7th power of a prime and a distinct prime (p^7*q) (A179664);
a(n)=116 if n is the product of exactly 2 distinct squares of primes and a different prime (p^2*q^2*r) (A179643);
a(n)=126 if n is the product of the 5th power of a prime and different distinct prime of the 2nd power (p^5*q^2) (A179646);
a(n)=128 if n is the product of the 8th power of a prime and a distinct prime (p^8*q) (A179668);
a(n)=150 if n is the product of the 4th power of a prime and 2 different distinct primes (p^4*q*r) (A179644);
a(n)=159 if n is the product of the 4th power of a prime and a distinct prime of power 3 (p^4*q^3) (A179666).
It is possible to continue with a(n) = 162, 178, 209, 224, 227, 238, 239, 260, 289, 309, 320, 333,...

Examples

			a(12) = 8 because the divisors of 12 are {1, 2, 3, 4, 6, 12} and GCD(d_i, d_j)>1 for the 8 following pairs of divisors: (2,4), (2,6), (2,12), (3,6), (3,12), (4,6), (4,12) and (6,12).
		

Crossrefs

Programs

  • Maple
    with(numtheory):nn:=100:
    for n from 1 to nn do:
    x:=divisors(n):n0:=nops(x):it:=0:
    for i from 1 to n0 do:
      for j from i+1 to n0 do:
       if gcd(x[i],x[j])>1
        then
        it:=it+1:
        else
       fi:
      od:
    od:
      printf(`%d, `,it):
    od:
  • Mathematica
    Table[Sum[Sum[(1 - KroneckerDelta[GCD[i, k], 1]) (1 - Ceiling[n/k] + Floor[n/k]) (1 - Ceiling[n/i] + Floor[n/i]), {i, k - 1}], {k, n}], {n, 100}] (* Wesley Ivan Hurt, Jan 01 2021 *)
  • PARI
    a(n)=my(d=divisors(n)); sum(i=2,#d, sum(j=1,i-1, gcd(d[i],d[j])>1)) \\ Charles R Greathouse IV, Aug 03 2016
    
  • PARI
    a(n)=my(f=factor(n)[,2],t=prod(i=1,#f,f[i]+1)); t*(t-1)/2 - (prod(i=1,#f,2*f[i]+1)+1)/2 \\ Charles R Greathouse IV, Aug 03 2016

Formula

a(n) = A066446(n) - A063647(n).
a(n) = Sum_{d1|n, d2|n, d1Wesley Ivan Hurt, Jan 01 2021

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

A280346 Numbers with 79 divisors.

Original entry on oeis.org

302231454903657293676544, 16423203268260658146231467800709255289, 3308722450212110699485634768279851414263248443603515625, 827269706064171159838078900184013751038269841857389464208009274449, 1692892739326831320764318961708001178036611459414853872137348292520966629744627081
Offset: 1

Views

Author

Omar E. Pol, Jan 01 2017

Keywords

Comments

Also, 78th 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 = 79.

Examples

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

Crossrefs

Programs

  • Mathematica
    With[{p = 22}, Table[Prime[n]^(Prime@ p - 1), {n, 5}]] (* Michael De Vlieger, Jan 01 2017 *)
  • PARI
    a(n)=prime(n)^78

Formula

a(n) = A000040(n)^(79-1) = A000040(n)^78.
A000005(a(n)) = 79.

A280347 Numbers with 83 divisors.

Original entry on oeis.org

4835703278458516698824704, 1330279464729113309844748891857449678409, 2067951531382569187178521730174907133914530277252197265625, 1986274564260074954771227439341817016242885890299592103563430267952049, 24785642596484137367310393918366845247634028377292875541962916350799472426091085092921
Offset: 1

Views

Author

Omar E. Pol, Jan 01 2017

Keywords

Comments

Also, 82nd 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 = 83.

Examples

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

Crossrefs

Programs

  • Mathematica
    With[{p = 23}, Table[Prime[n]^(Prime@ p - 1), {n, 5}]] (* Michael De Vlieger, Jan 01 2017 *)
  • PARI
    a(n)=prime(n)^82

Formula

a(n) = A000040(n)^(83-1) = A000040(n)^82.
A000005(a(n)) = 83.

A280349 Numbers with 89 divisors.

Original entry on oeis.org

309485009821345068724781056, 969773729787523602876821942164080815560161, 32311742677852643549664402033982923967414535582065582275390625, 233683216210633558353880137011125430143959282107856711392134007594290612801, 43909277783870034878569768760415886733743786946105343887995366053338664170638348798300219681
Offset: 1

Views

Author

Omar E. Pol, Jan 01 2017

Keywords

Comments

Also, 88th 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 = 89.

Examples

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

Crossrefs

Programs

  • Mathematica
    With[{p = 24}, Table[Prime[n]^(Prime@ p - 1), {n, 5}]] (* Michael De Vlieger, Jan 01 2017 *)
  • PARI
    a(n)=prime(n)^88

Formula

a(n) = A000040(n)^(89-1) = A000040(n)^88.
A000005(a(n)) = 89.

A382292 Numbers k such that A382290(k) = 1.

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, 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, 432, 440, 448, 456, 459, 472, 480, 486, 488, 500
Offset: 1

Views

Author

Amiram Eldar, Mar 21 2025

Keywords

Comments

First differs from A374590 and A375432 at n = 25: A374590(25) = A375432(25) = 216 is not a term of this sequence.
Numbers k such that A382291(k) = 2, i.e., numbers whose number of infinitary divisors is twice the number of their unitary divisors.
Numbers whose prime factorization has a single exponent that is a sum of two distinct powers of 2 (A018900) and all the other exponents, if they exist, are powers of 2. Equivalently, numbers of the form p^e * m, where p is a prime, e is a term in A018900, and m is a term in A138302 that is coprime to p.
If k is a term then k^2 is also a term. If m is a term in A138302 that is coprime to k then k * m is also a term. The primitive terms, i.e., the terms that cannot be generated from smaller terms using these rules, are the numbers of the form p^(2^i+1), where p is prime and i >= 1.
Analogous to A060687, which is the sequence of numbers k with prime excess A046660(k) = 2.
The asymptotic density of this sequence is A271727 * Sum_{p prime} (((1 - 1/p)/f(1/p)) * Sum_{k>=1} 1/p^A018900(k)) = 0.11919967112489084407..., where f(x) = 1 - x^3 + Sum_{k>=2} (x^(2^k)-x^(2^k+1)).

Crossrefs

Subsequences (numbers of the form): A030078 (p^3), A050997 (p^5), A030516 (p^6), A179665 (p^9), A030629 (p^10), A030631 (p^12), A065036 (p^3*q), A178740 (p^5*q), A189987 (p^6*q), A179692 (p^9*q), A143610 (p^2*q^3), A179646 (p^5*q^2), A189990 (p^2*q^6), A179702 (p^4*q^5), A179666 (p^4*q^3), A190464 (p^4*q^6), A163569 (p^3*q^2*r), A189975 (p*q*r^3), A190115 (p^2*q^3*r^4), A381315, A048109.

Programs

  • Mathematica
    f[p_, e_] := DigitCount[e, 2, 1] - 1; q[1] = False; q[n_] := Plus @@ f @@@ FactorInteger[n] == 1; Select[Range[500], q]
  • PARI
    isok(k) = vecsum(apply(x -> hammingweight(x) - 1, factor(k)[, 2])) == 1;
Previous Showing 21-30 of 39 results. Next