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

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

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;

A362841 Numbers with at least one 5 in their prime signature.

Original entry on oeis.org

32, 96, 160, 224, 243, 288, 352, 416, 480, 486, 544, 608, 672, 736, 800, 864, 928, 972, 992, 1056, 1120, 1184, 1215, 1248, 1312, 1376, 1440, 1504, 1568, 1632, 1696, 1701, 1760, 1824, 1888, 1944, 1952, 2016, 2080, 2144, 2208, 2272, 2336, 2400, 2430, 2464, 2528, 2592, 2656, 2673, 2720, 2784, 2848, 2912, 2976
Offset: 1

Views

Author

R. J. Mathar, May 05 2023

Keywords

Comments

Contains all odd multiples of 2^5: Each second term of A174312 is in this sequence.
The asymptotic density of this sequence is 1 - Product_{p prime} (1 - 1/p^5 + 1/p^6) = 0.01863624892... . - Amiram Eldar, May 05 2023

Examples

			Contains 2^5, 2^5*3, 2^5*5, 2^5*7, 3^5, 2^5*3^2, 2^5*11, 2^5*13, 2^5*3*5, 2*3^5, etc.
		

Crossrefs

Cf. A038109 (at least one 2), A176297 (at least one 3), A050997 (subsequence), A178740 (subsequence), A179646 (subsequence), A179667 (subsequence), A179671 (subsequence), A174312.

Programs

  • Mathematica
    Select[Range[3000], MemberQ[FactorInteger[#][[;;, 2]], 5] &] (* Amiram Eldar, May 05 2023 *)

A355462 Powerful numbers divisible by exactly 2 distinct primes.

Original entry on oeis.org

36, 72, 100, 108, 144, 196, 200, 216, 225, 288, 324, 392, 400, 432, 441, 484, 500, 576, 648, 675, 676, 784, 800, 864, 968, 972, 1000, 1089, 1125, 1152, 1156, 1225, 1296, 1323, 1352, 1372, 1444, 1521, 1568, 1600, 1728, 1936, 1944, 2000, 2025, 2116, 2304, 2312, 2500
Offset: 1

Views

Author

Amiram Eldar, Jul 03 2022

Keywords

Comments

First differs from A286708 at n = 25.
Number of the form p^i * q^j, where p != q are primes and i,j > 1.
Numbers k such that A001221(k) = 2 and A051904(k) >= 2.
The possible values of the number of the divisors (A000005) of terms in this sequence is any composite number that is not 8 or twice a prime (A264828 \ {1, 8}).
675 = 3^3*5^2 and 676 = 2^2*13^2 are 2 consecutive integers in this sequence. There are no other such pairs below 10^22 (the lesser members of such pairs are terms of A060355).

Examples

			36 is a term since 36 = 2^2 * 3^2.
		

Crossrefs

Programs

  • Mathematica
    Select[Range[2500], Length[(e = FactorInteger[#][[;; , 2]])] == 2 && Min[e] > 1 &]
  • PARI
    is(n) = {my(f=factor(n)); #f~ == 2 && vecmin(f[,2]) > 1};

Formula

Sum_{n>=1} 1/a(n) = ((Sum_{p prime} (1/(p*(p-1))))^2 - Sum_{p prime} (1/(p^2*(p-1)^2)))/2 = 0.1583860791... .

A381314 Powerful numbers that have a single exponent in their prime factorization that equals 2.

Original entry on oeis.org

4, 9, 25, 49, 72, 108, 121, 144, 169, 200, 288, 289, 324, 361, 392, 400, 500, 529, 576, 675, 784, 800, 841, 961, 968, 972, 1125, 1152, 1323, 1352, 1369, 1372, 1568, 1600, 1681, 1849, 1936, 2025, 2209, 2304, 2312, 2500, 2704, 2809, 2888, 2916, 3087, 3136, 3200
Offset: 1

Views

Author

Amiram Eldar, Feb 19 2025

Keywords

Comments

Number of the form A036966(m)/p, m >= 2, where p is a prime divisor of A036966(m).

Crossrefs

Programs

  • Mathematica
    With[{max = 3200}, Select[Union@ Flatten@ Table[i^2 * j^3, {j, 1, max^(1/3)}, {i, 1, Sqrt[max/j^3]}], Count[FactorInteger[#][[;; , 2]], 2] == 1 &]]
  • PARI
    isok(k) = if(k == 1, 0, my(e = factor(k)[, 2]); vecmin(e) > 1 && #select(x -> (x==2), e) == 1);

Formula

Sum_{n>=1} 1/a(n) = Sum_{p prime}((p-1)/(p^3-p^2+1)) * Product_{p prime} (1 + 1/(p^2*(p-1))) = 0.53045141423939736076... .
Previous Showing 11-15 of 15 results.