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

A036785 Numbers divisible by the squares of two distinct primes.

Original entry on oeis.org

36, 72, 100, 108, 144, 180, 196, 200, 216, 225, 252, 288, 300, 324, 360, 392, 396, 400, 432, 441, 450, 468, 484, 500, 504, 540, 576, 588, 600, 612, 648, 675, 676, 684, 700, 720, 756, 784, 792, 800, 828, 864, 882, 900, 936, 968, 972, 980, 1000, 1008, 1044
Offset: 1

Views

Author

Keywords

Comments

Not squarefree, not a nontrivial prime power and not in {squarefree} times {nontrivial prime powers}.
Numbers k such that A056170(k) > 1. The asymptotic density of this sequence is 1 - (6/Pi^2) * (1 + A154945) = 0.05668359058... - Amiram Eldar, Nov 01 2020

References

  • CRC Standard Mathematical Tables and Formulae, 30th ed., (1996) page 102-105.

Crossrefs

Equivalent sequence for 3 distinct primes: A318720.
Cf. A085986, A338539, A339245 (subsequences).
Subsequence of A038838.

Programs

  • Mathematica
    Select[Range@ 1050, And[Length@ # > 1, Total@ Boole@ Map[# > 1 &, #[[All, -1]]] > 1] &@ FactorInteger@ # &] (* Michael De Vlieger, Apr 25 2017 *)
    dstdpQ[n_]:=Length[Select[Sqrt[#]&/@Divisors[n],PrimeQ]]>1; Select[ Range[ 1100],dstdpQ] (* Harvey P. Dale, Jan 15 2020 *)
  • PARI
    is(n)=my(f=vecsort(factor(n)[,2],,4));#f>1&&f[2]>1 \\ Charles R Greathouse IV, Nov 15 2012

Extensions

More terms from Larry Reeves (larryr(AT)acm.org), Apr 03 2000
New name from Charles R Greathouse IV, Nov 15 2012

A318720 Numbers k such that there exists a strict relatively prime factorization of k in which no pair of factors is relatively prime.

Original entry on oeis.org

900, 1764, 1800, 2700, 3528, 3600, 4356, 4500, 4900, 5292, 5400, 6084, 6300, 7056, 7200, 8100, 8712, 8820, 9000, 9800, 9900, 10404, 10584, 10800, 11025, 11700, 12100, 12168, 12348, 12600, 12996, 13068, 13500, 14112, 14400, 14700, 15300, 15876, 16200, 16900
Offset: 1

Views

Author

Gus Wiseman, Sep 02 2018

Keywords

Comments

From Amiram Eldar, Nov 01 2020: (Start)
Also, numbers with more than two non-unitary prime divisors, i.e., numbers k such that A056170(k) > 2, or equivalently, numbers divisible by the squares of three distinct primes.
The complement of the union of A005117, A190641 and A338539.
The asymptotic density of this sequence is 1 - 6/Pi^2 - (6/Pi^2)*A154945 - (3/Pi^2)*(A154945^2 - A324833) = 0.0033907041... (End)

Examples

			900 is in the sequence because the factorization 900 = (6*10*15) is relatively prime (since the GCD of (6,10,15) is 1) but each of the pairs (6,10), (6,15), (10,15) has a common divisor > 1. Larger examples are:
1800 = (6*15*20) = (10*12*15).
9900 = (6*10*165) = (6*15*110) = (10*15*66).
5400 = (6*20*45) = (10*12*45) = (10*15*36) = (15*18*20).
60 is not in the sequence because all its possible factorizations (4 * 15, 3 * 4 * 5, etc.) contain at least one pair that is coprime, if not more than one prime.
		

Crossrefs

Programs

  • Mathematica
    strfacs[n_] := If[n <= 1, {{}}, Join@@Table[(Prepend[#1, d] &)/@Select[strfacs[n/d], Min@@#1 > d &], {d, Rest[Divisors[n]]}]]; Select[Range[10000], Function[n, Select[strfacs[n], And[GCD@@# == 1, And@@(GCD[##] > 1 &)@@@Select[Tuples[#, 2], Less@@# &]] &] != {}]]
    Select[Range[20000], Count[FactorInteger[#][[;;,2]], ?(#1 > 1 &)] > 2 &] (* _Amiram Eldar, Nov 01 2020 *)

A347960 Numbers k for which A348036(k) > A007947(k).

Original entry on oeis.org

36, 72, 100, 108, 144, 180, 196, 200, 216, 225, 252, 288, 300, 324, 360, 392, 396, 400, 432, 450, 468, 484, 500, 504, 540, 576, 588, 600, 612, 648, 675, 676, 684, 700, 720, 756, 784, 792, 800, 828, 864, 900, 936, 968, 972, 980, 1000, 1008, 1044, 1080, 1089, 1100, 1116, 1125, 1152, 1156, 1176, 1188, 1200, 1224, 1260, 1296
Offset: 1

Views

Author

Antti Karttunen, Oct 19 2021

Keywords

Comments

Numbers k such that A348039(k) > 1.
Numbers k such that A348037(k) < A003557(k).
Numbers k such that A327564(k) > A348038(k).
Differs from A036785 and A338539 for the first time at n=20, where a(n) = 450, as A036785(20) = A338539(20) = 441 is not included in this sequence.

Crossrefs

Programs

  • Mathematica
    f[p_, e_] := p*(p + 1)^(e - 1); q[n_] := GCD[n, Times @@ f @@@ (fct = FactorInteger[n])] > Times @@ First /@ fct; Select[Range[1300], q] (* Amiram Eldar, Oct 20 2021 *)
  • PARI
    A003968(n) = {my(f=factor(n)); for (i=1, #f~, p= f[i, 1]; f[i, 1] = p*(p+1)^(f[i, 2]-1); f[i, 2] = 1); factorback(f); }
    A007947(n) = factorback(factorint(n)[, 1]);
    A348036(n) = gcd(n, A003968(n));
    isA347960(n) = (A348036(n)>A007947(n));

A338540 Numbers having exactly three non-unitary prime factors.

Original entry on oeis.org

900, 1764, 1800, 2700, 3528, 3600, 4356, 4500, 4900, 5292, 5400, 6084, 6300, 7056, 7200, 8100, 8712, 8820, 9000, 9800, 9900, 10404, 10584, 10800, 11025, 11700, 12100, 12168, 12348, 12600, 12996, 13068, 13500, 14112, 14400, 14700, 15300, 15876, 16200, 16900, 17100
Offset: 1

Views

Author

Amiram Eldar, Nov 01 2020

Keywords

Comments

Numbers k such that A056170(k) = A001221(A057521(k)) = 3.
Numbers divisible by the squares of exactly three distinct primes.
Subsequence of A318720 and first differs from it at n = 123.
The asymptotic density of this sequence is (eta_1^3 - 3*eta_1*eta_2 + 2*eta_3)/Pi^2 = 0.0032920755..., where eta_j = Sum_{p prime} 1/(p^2-1)^j (Pomerance and Schinzel, 2011).

Examples

			900 = 2^2 * 3^2 * 5^2 is a term since it has exactly 3 prime factors, 2, 3 and 5, that are non-unitary.
		

Crossrefs

Subsequence of A013929, A318720 and A327877.
Cf. A154945 (eta_1), A324833 (eta_2), A324834 (eta_3).

Programs

  • Mathematica
    Select[Range[17000], Count[FactorInteger[#][[;;,2]], _?(#1 > 1 &)] == 3 &]

A338541 Numbers having exactly four non-unitary prime factors.

Original entry on oeis.org

44100, 88200, 108900, 132300, 152100, 176400, 213444, 217800, 220500, 260100, 264600, 298116, 304200, 308700, 324900, 326700, 352800, 396900, 426888, 435600, 441000, 456300, 476100, 485100, 509796, 520200, 529200, 544500, 573300, 592900, 596232, 608400, 617400
Offset: 1

Views

Author

Amiram Eldar, Nov 01 2020

Keywords

Comments

Numbers k such that A056170(k) = A001221(A057521(k)) = 4.
Numbers divisible by the squares of exactly four distinct primes.
The asymptotic density of this sequence is (eta_1^4 - 6*eta_1^2*eta_2 + 3*eta_2^2 + 8*eta_1*eta_3 - 6*eta_4)/(4*Pi^2) = 0.0000970457..., where eta_j = Sum_{p prime} 1/(p^2-1)^j (Pomerance and Schinzel, 2011).

Examples

			44100 = 2^2 * 3^2 * 5^2 * 7^2 is a term since it has exactly 4 prime factors, 2, 3, 5 and 7, that are non-unitary.
		

Crossrefs

Subsequence of A013929 and A318720.
Cf. A154945 (eta_1), A324833 (eta_2), A324834 (eta_3), A324835 (eta_4).

Programs

  • Mathematica
    Select[Range[620000], Count[FactorInteger[#][[;;,2]], _?(#1 > 1 &)] == 4 &]

A338542 Numbers having exactly five non-unitary prime factors.

Original entry on oeis.org

5336100, 7452900, 10672200, 12744900, 14905800, 15920100, 16008300, 18404100, 21344400, 22358700, 23328900, 25489800, 26680500, 29811600, 31472100, 31840200, 32016600, 36072036, 36808200, 37088100, 37264500, 37352700, 38234700, 39312900, 42380100, 42688800, 43956900
Offset: 1

Views

Author

Amiram Eldar, Nov 01 2020

Keywords

Comments

Numbers k such that A056170(k) = A001221(A057521(k)) = 5.
Numbers divisible by the squares of exactly five distinct primes.
The asymptotic density of this sequence is (eta_1^5 - 10*eta_1^3*eta_2 + 15*eta_1*eta_2^2 + 20*eta_1^2*eta_3 - 20*eta_2*eta_3 - 30*eta_1*eta_4 + 24*eta_5)/(20*Pi^2) = 0.0000015673..., where eta_j = Sum_{p prime} 1/(p^2-1)^j (Pomerance and Schinzel, 2011).

Examples

			5336100 = 2^2 * 3^2 * 5^2 * 7^2 * 11^2 is a term since it has exactly 5 prime factors, 2, 3, 5, 7 and 11, that are non-unitary.
		

Crossrefs

Subsequence of A013929, A318720 and A327877.
Cf. A154945 (eta_1), A324833 (eta_2), A324834 (eta_3), A324835 (eta_4), A324836 (eta_5).

Programs

  • Mathematica
    Select[Range[2*10^7], Count[FactorInteger[#][[;;,2]], _?(#1 > 1 &)] == 5 &]

A347892 Nonsquarefree numbers k such that A003968(k) is a multiple of k, where A003968 is multiplicative with a(p^e) = p*(p+1)^(e-1).

Original entry on oeis.org

36, 72, 180, 216, 252, 360, 396, 432, 468, 504, 612, 684, 792, 828, 864, 936, 1044, 1080, 1116, 1224, 1260, 1296, 1332, 1368, 1476, 1512, 1548, 1656, 1692, 1908, 1980, 2088, 2124, 2160, 2196, 2232, 2340, 2376, 2412, 2520, 2556, 2592, 2628, 2664, 2772, 2808, 2844, 2952, 2988, 3024, 3060, 3096, 3204, 3276, 3384, 3420
Offset: 1

Views

Author

Antti Karttunen, Oct 29 2021

Keywords

Comments

All terms k are multiples of 36 and A056170(k) = 2. Proof: Imagine k had a non-unitary prime divisor p^e, with p > 3 and e > 1. Then p^e divides A003968(k) only if k has also another non-unitary prime divisor q^h (with h > 1), such that p divides (q+1), which implies that q > p. But then q^h divides A003968(k) only if there is yet another non-unitary prime divisor r^i, such that r > q (and i > 1), and so on, which is clearly impossible by reductio ad infinitum. Therefore we should consider only the cases p=2 and p=3, because they are only primes that can occur as non-unitary prime factors in k, and at least either of them must occur with exponent larger than one, because every k is nonsquarefree. Let e = A007814(k) and h = A007949(k), so that 2^e and 3^h are the highest powers of 2 and 3 that divide k. Because A003968 changes "extra" 2's to 3's and extra 3's to 4's, it must follow that e >= h > e/2. Therefore, if e >= 2 (k is a multiple of 4), h must be at least 2. On the other hand, if h >= 2, then e also must be at least 2. In other words, if k is a multiple of 4, it must then also be a multiple of 9, and vice versa, thus k is a multiple of 36 and k has exactly two non-unitary prime divisors (2^e and 3^h, with e, h > 1), therefore this is a subsequence of A338539.

Crossrefs

Intersection of A013929 and A348499. Subsequence of A338539.

Programs

  • Mathematica
    f[p_, e_] := p*(p + 1)^(e - 1); s[n_] := Times @@ f @@@ FactorInteger[n]; Select[Range[3500], ! SquareFreeQ[#] && Divisible[s[#], #] &] (* Amiram Eldar, Oct 29 2021 *)
  • PARI
    A003968(n) = {my(f=factor(n)); for (i=1, #f~, p= f[i, 1]; f[i, 1] = p*(p+1)^(f[i, 2]-1); f[i, 2] = 1); factorback(f); }
    isA347892(n) = (!issquarefree(n) && !(A003968(n)%n));
Showing 1-7 of 7 results.