A036785 -id:A036785 - cp's OEIS frontend

A085986 Squares of the squarefree semiprimes (p^2*q^2).

36, 100, 196, 225, 441, 484, 676, 1089, 1156, 1225, 1444, 1521, 2116, 2601, 3025, 3249, 3364, 3844, 4225, 4761, 5476, 5929, 6724, 7225, 7396, 7569, 8281, 8649, 8836, 9025, 11236, 12321, 13225, 13924, 14161, 14884, 15129, 16641, 17689, 17956, 19881

Offset: 1

Alford Arnold, Jul 06 2003

This sequence is a member of a family of sequences directly related to A025487. First terms and known sequences are listed below: 1, A000007; 2, A000040; 4, A001248; 6, A006881; 8, A030078; 12, A054753; 16, A030514; 24, A065036; 30, A007304; 32, A050997; 36, this sequence; 48, ?; 60, ?; 64, ?; ....
Subsequence of A077448. The numbers in A077448 but not in here are 1, the squares of A046386, the squares of A067885, etc. - R. J. Mathar, Sep 12 2008
a(4)-a(3)=29 and a(3)+a(4)=421 are both prime.  There are no other cases where the sum and difference of two members of this sequence are both prime. - Robert Israel and J. M. Bergot, Oct 25 2019

			A006881 begins 6 10 14 15 ... so this sequence begins 36 100 196 225 ...

Subsequence of A036785 and of A077448.
Subsequence of A062503.
Cf. A025487.
Cf. A085548, A085964.

[k^2:k in [1..150]| IsSquarefree(k) and #PrimeDivisors(k) eq 2]; // Marius A. Burtea, Oct 24 2019

f[n_]:=Sort[Last/@FactorInteger[n]]=={2,2}; Select[Range[20000], f] (* Vladimir Joseph Stephan Orlovsky, Aug 14 2009 *)
Select[Range[200],PrimeOmega[#]==2&&SquareFreeQ[#]&]^2 (* Harvey P. Dale, Mar 07 2013 *)

list(lim)=my(v=List(), x=sqrtint(lim\=1), t); forprime(p=2, x\2, t=p; forprime(q=2, min(x\t,p-1), listput(v, (t*q)^2))); Set(v) \\ Charles R Greathouse IV, Sep 22 2015

is(n)=factor(n)[,2]==[2,2]~ \\ Charles R Greathouse IV, Oct 19 2015

from math import isqrt
from sympy import primepi, primerange
def A085986(n):
    def f(x): return int(n+x+(t:=primepi(s:=isqrt(x)))+(t*(t-1)>>1)-sum(primepi(x//k) for k in primerange(1, s+1)))
    m, k = n, f(n)
    while m != k:
        m, k = k, f(k)
    return m**2 # Chai Wah Wu, Aug 18 2024

a(n) = A006881(n)^2.

Sum_{n>=1} 1/a(n) = (P(2)^2 - P(4))/2 = (A085548^2 - A085964)/2 = 0.063767..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

A038838 Numbers that are divisible by the square of an odd prime.

Original entry on oeis.org

9, 18, 25, 27, 36, 45, 49, 50, 54, 63, 72, 75, 81, 90, 98, 99, 100, 108, 117, 121, 125, 126, 135, 144, 147, 150, 153, 162, 169, 171, 175, 180, 189, 196, 198, 200, 207, 216, 225, 234, 242, 243, 245, 250, 252, 261, 270, 275, 279, 288, 289, 294, 297, 300, 306

Offset: 1

David W. Wilson

nonn

Condition 1 of Theorem 7.5 (Robinson, 1979) includes: "k is a multiple of a square of an odd prime." - Jonathan Vos Post, Aug 06 2007
If m is a term, every k*m with k > 1 is another term and the primitive terms are the square of odd primes. The subsequence of odd terms is A053850 while the even terms 18, 36, 50, 54, 72, 90, 98, ... are exactly twice the terms of this sequence. - Bernard Schott, Nov 20 2020
The asymptotic density of this sequence is 1 - 8/Pi^2 = 0.189430... - Amiram Eldar, Nov 21 2020

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
M. Beeler, R. W. Gosper, and R. Schroeppel, HAKMEM ITEM 45.
Raphael M. Robinson, Multiple tiling of n-dimensional space by unit cubes, Math. Z., Vol. 166 (1979), pp. 225-264.
Chuanming Zong, What is known about unit cubes, Bull. Amer. Math. Soc., Vol. 42, No. 2 (2005), pp. 181-211; Robinson theorem cited on p. 204.

Cf. A000040, A065091, A122132 (complement).
Cf. A013929 (supersequence of nonsquarefrees).
Subsequences: A001248 \ {2} (primitives), A053850 (odds), A036785 (divisible by the squares of two distinct primes).
Subsequence of A167662. - Reinhard Zumkeller, Nov 08 2009

{a(n) = my(m, c); if( n<1, 0, c=0; m=0; while( cMichael Somos, Aug 22 2006 */

list(lim)=my(v=List(),n,e,t); forfactored(k=9,lim\1, e=k[2][,2]; t=#e; n=k[1]; if(if(n%2 && t, vecmax(e)>1, t>1, vecmax(e[2..t])>1, 0), listput(v, k[1]))); Vec(v) \\ Charles R Greathouse IV, Jan 08 2018

{a(n)} = {j such that for some k>1 A001248(k)|j} = {j such that for some k>0 (A065091(k)^2)|j}. - Jonathan Vos Post, Aug 06 2007

A008966(A000265(a(n))) = 0. - Reinhard Zumkeller, Nov 08 2009

A366890 Irregular triangle, wherein row n lists in ascending order all numbers k whose arithmetic derivative k' is equal to the n-th primorial, A002110(n), and that have more than two prime factors with multiplicity. Rows of length zero are simply omitted, i.e., when A369000(n) = 0.

Original entry on oeis.org

1547371, 79332523, 1102527599503, 25336943536819, 25962012375103, 25970380120783, 66702554987143, 526285951027003, 927949814519899, 7777707036642079, 9584173681667203, 13082430772438171, 22101822021783739, 4958985803436403, 32006922970429003, 32076018550175863, 49806227168831659, 84682266449971639, 97995266657958403

Offset: 1

Antti Karttunen, Jan 09 2024

For n > 0, numbers k such that A003415(k) = A002110(n) and A001222(k) > 2.
Sequence as a whole is not listed in ascending order, even though each batch of solutions for each n for which A369000(n) > 0 are. For example, we have a(14) < a(13) because A003415(22101822021783739) = A002110(12), while A003415(4958985803436403) = A002110(13). See the examples.
Question: Are there any common terms with A036785, that is, with A368697?

			For rows n=1..6, 9 & 10 nothing is listed, as those rows are empty.
Row for n=7 has just one term: 1547371 (= 7^2 * 23 * 1373). Note that A003415(1547371) = 510510 = A002110(7).
Row for n=8 has just one term: 79332523 (= 17^2 * 277 * 991).
Row for n=11 has two terms:
  1102527599503 (= 11^2 * 11071 * 823033),
  25336943536819 (= 157 * 743 * 5749 * 37781).
Row for n=12 has nine terms:
  25962012375103 (= 7^2 * 8597 * 61630451),
  25970380120783 (= 7^2 * 41387 * 12806141),
  66702554987143 (= 19^2 * 167 * 1106416889),
  526285951027003 (= 73 * 3919 * 7013 * 262313),
  927949814519899 (= 269 * 271 * 1697 * 7501033),
  7777707036642079 (= 2203 * 2791 * 7349 * 172127),
  9584173681667203 (= 2131 * 5953 * 7901 * 95621),
  13082430772438171 (= 3109 * 5861 * 24421 * 29399),
  22101822021783739 (= 8783 * 11777 * 13921 * 15349).
Row for n=13 has 18 terms, and begins with:
  4958985803436403 (= 37^2 * 137 * 26440450451),
and ends with:
  3206697143570677543 (= 36899 * 41983 * 45233 * 45763).
Note that A003415(3206697143570677543) = 304250263527210 = A002110(13).

Antti Karttunen, Table of n, a(n) for n = 1..31
Antti Karttunen, PARI program

When the whole sequence is sorted into ascending order, equal to A327978 without any semiprime solutions (solutions in A001358), and also a subsequence of following sequences: A004709, A327862, A328234.
Cf. A001222, A002110, A003415, A116979, A351029, A368703.
Cf. also A036785, A368697, A328243, A369240.

PARI
```
\\ See the attached PARI-program
```

A375398 Numbers k such that the minima of maximal anti-runs in the weakly increasing sequence of prime factors of k (with multiplicity) are distinct.

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 18, 19, 21, 22, 23, 26, 29, 30, 31, 33, 34, 35, 37, 38, 39, 41, 42, 43, 46, 47, 50, 51, 53, 55, 57, 58, 59, 61, 62, 65, 66, 67, 69, 70, 71, 73, 74, 75, 77, 78, 79, 82, 83, 85, 86, 87, 89, 90, 91, 93, 94, 95, 97, 98

Offset: 1

Gus Wiseman, Aug 16 2024

nonn

First differs from A375402 in lacking 20.
An anti-run is a sequence with no adjacent equal parts.
The minima of maximal anti-runs in a sequence are obtained by splitting it into maximal anti-run subsequences and taking the least term of each.
Note the prime factors can alternatively be taken in weakly decreasing order.

			The prime factors of 300 are {2,2,3,5,5}, with maximal anti-runs ((2),(2,3,5),(5)), with minima (2,2,5), so 300 is not in the sequence.
The prime factors of 450 are {2,3,3,5,5}, with maximal anti-runs ((2,3),(3,5),(5)), with minima (2,3,5), so 450 is in the sequence.

Gus Wiseman, Sequences counting and ranking compositions by their leaders (for six types of runs).

A version for compositions is A374638, counted by A374518.
These are positions of strict rows in A375128, sums A374706, ranks A375400.
Partitions (or reversed partitions) of this type are counted by A375134.
For identical instead of distinct we have A375396, counted by A115029.
The complement is A375399, counted by A375404.
For maxima instead of minima we have A375402, counted by A375133.
The complement for maxima is A375403, counted by A375401.
A000041 counts integer partitions, strict A000009.
A003242 counts anti-run compositions, ranks A333489.
A number's prime factors (A027746, reverse A238689) have sum A001414, min A020639, max A006530.
A number's prime indices (A112798, reverse A296150) have sum A056239, min A055396, max A061395.
Both have length A001222, distinct A001221.
Cf. A034296, A036785, A046660, A065200, A065201, A358836, A374632, A374761, A374767, A374768, A375136.

Select[Range[100],UnsameQ@@Min /@ Split[Flatten[ConstantArray@@@FactorInteger[#]],UnsameQ]&]

A346485 Möbius transform of A342001, where A342001(n) = A003415(n)/A003557(n).

Original entry on oeis.org

0, 1, 1, 1, 1, 3, 1, 1, 1, 5, 1, 2, 1, 7, 6, 1, 1, 1, 1, 4, 8, 11, 1, 2, 1, 13, 1, 6, 1, 14, 1, 1, 12, 17, 10, 0, 1, 19, 14, 4, 1, 20, 1, 10, 4, 23, 1, 2, 1, 1, 18, 12, 1, 1, 14, 6, 20, 29, 1, 8, 1, 31, 6, 1, 16, 32, 1, 16, 24, 34, 1, 0, 1, 37, 2, 18, 16, 38, 1, 4, 1, 41, 1, 12, 20, 43, 30, 10, 1, 4, 18, 22, 32, 47

Offset: 1

Antti Karttunen, Aug 26 2021

nonn

Conjecture 1: After the initial zero, the positions of other zeros is given by A036785.

Conjecture 2: No negative terms. Checked up to n = 2^24.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537

Cf. A003415, A003557, A008683, A036785, A342001, A347232 [= a(A276086(n))].

Cf. also A300251, A300717, A347234, A347235, A347395.

A003415(n) = if(n<=1, 0, my(f=factor(n)); n*sum(i=1, #f~, f[i, 2]/f[i, 1]));
A003557(n) = (n/factorback(factorint(n)[, 1]));
A342001(n) = (A003415(n) / A003557(n));
A346485(n) = sumdiv(n,d,moebius(n/d)*A342001(d));

a(n) = Sum_{d|n} A008683(n/d) * A342001(d).
Dirichlet g.f.: Product_{p prime} (1+p^(1-s)-p^(-s)) * Sum_{p prime} p^s/((p^s-1)*(p^s+p-1)). - Sebastian Karlsson, May 08 2022
Sum_{k=1..n} a(k) ~ c * A065464 * n^2 / 2, where c = Sum_{j>=2} (1/2 + (-1)^j * (Fibonacci(j) - 1/2))*PrimeZetaP(j) = 0.4526952873143153104685540856936425315834753528741817723313791528384... - Vaclav Kotesovec, Mar 04 2023

A375396 Numbers not divisible by the square of any prime factor except (possibly) the least. Hooklike numbers.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 51, 52, 53, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71

Offset: 1

Gus Wiseman, Aug 16 2024

nonn

Also numbers k such that the minima of the maximal anti-runs in the weakly increasing sequence of prime factors of k (with multiplicity) are identical. Here, an anti-run is a sequence with no adjacent equal parts, and the minima of the maximal anti-runs in a sequence are obtained by splitting it into maximal anti-run subsequences and taking the least term of each. Note the prime factors can alternatively be taken in weakly decreasing order.
The complement is a superset of A036785 = products of a squarefree number and a prime power.
The asymptotic density of this sequence is (1/zeta(2)) * (1 + Sum_{p prime} (1/(p^2-p)) / Product_{primes q <= p} (1 + 1/q)) = 0.884855661165... . - Amiram Eldar, Oct 26 2024

			The prime factors of 300 are {2,2,3,5,5}, with maximal anti-runs {{2},{2,3,5},{5}}, with minima (2,2,5), so 300 is not in the sequence.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Gus Wiseman, Sequences counting and ranking compositions by their leaders (for six types of runs).

The complement is a superset of A036785.
For maxima instead of minima we have A065200, counted by A034296.
The complement for maxima is A065201, counted by A239955.
Partitions of this type are counted by A115029.
A version for compositions is A374519, counted by A374517.
Also positions of identical rows in A375128, sums A374706, ranks A375400.
The complement is A375397, counted by A375405.
For distinct instead of identical minima we have A375398, counts A375134.
The complement for distinct minima is A375399, counted by A375404.
A000041 counts integer partitions, strict A000009.
A003242 counts anti-run compositions, ranks A333489.
A011782 comps counts compositions.
A number's prime factors (A027746, reverse A238689) have sum A001414, min A020639, max A006530.
A number's prime indices (A112798, reverse A296150) have sum A056239, min A055396, max A061395.
Both have length A001222, distinct A001221.
See the formula section for the relationships with A005117, A028234.
Cf. A000005, A013661, A046660, A272919, A319066, A358905, A374686, A374704, A374742, A375133, A375136, A375401.

Select[Range[100],SameQ@@Min /@ Split[Flatten[ConstantArray@@@FactorInteger[#]],UnsameQ]&]

is(k) = if(k > 1, my(e = factor(k)[, 2]); vecprod(e) == e[1], 1); \\ Amiram Eldar, Oct 26 2024

{a(n)} = {k >= 1 : A028234(k) is in A005117}. - Peter Munn, May 09 2025

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

Gus Wiseman, Sep 02 2018

nonn

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)

			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.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A001055, A001221, A001222, A007716, A045778, A051185, A078374, A281116, A303140, A303283, A305843, A305854, A317748, A318715, A318717, A318721.

Cf. A005117, A036785, A056170, A154945, A190641, A324833, A338539.

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 *)

A338539 Numbers having exactly two non-unitary prime factors.

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, 936, 968, 972, 980, 1000, 1008, 1044, 1080, 1089

Offset: 1

Amiram Eldar, Nov 01 2020

nonn

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

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

Amiram Eldar, Table of n, a(n) for n = 1..10000
Carl Pomerance and Andrzej Schinzel, Multiplicative Properties of Sets of Residues, Moscow Journal of Combinatorics and Number Theory, Vol. 1, No. 1 (2011), pp. 52-66. See pp. 61-62.

Subsequence of A013929 and A036785.
Cf. A001221, A056170, A057521, A190641, A338540, A338541, A338542, A347892 (subsequence).
Cf. A154945 (eta_1), A324833 (eta_2).

Select[Range[1000], Count[FactorInteger[#][[;;,2]], _?(#1 > 1 &)] == 2 &]

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

Antti Karttunen, Oct 19 2021

nonn

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.

Cf. A003557, A003959, A003968, A007947, A327564, A348036, A348037, A348038, A348039.

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 *)

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));

A375397 Numbers divisible by the square of some prime factor other than the least. Non-hooklike numbers.

Original entry on oeis.org

18, 36, 50, 54, 72, 75, 90, 98, 100, 108, 126, 144, 147, 150, 162, 180, 196, 198, 200, 216, 225, 234, 242, 245, 250, 252, 270, 288, 294, 300, 306, 324, 338, 342, 350, 360, 363, 375, 378, 392, 396, 400, 414, 432, 441, 450, 468, 484, 486, 490, 500, 504, 507, 522

Offset: 1

Gus Wiseman, Aug 16 2024

nonn

Contains no squarefree numbers A005117 or prime powers A000961, but some perfect powers A131605.
Also numbers k such that the minima of the maximal anti-runs in the weakly increasing sequence of prime factors of k (with multiplicity) are not identical. Here, an anti-run is a sequence with no adjacent equal parts, and the minima of the maximal anti-runs in a sequence are obtained by splitting it into maximal anti-run subsequences and taking the least term of each. Note the prime factors can alternatively be taken in weakly decreasing order.
Includes all terms of A036785 = non-products of a squarefree number and a prime power.
The asymptotic density of this sequence is 1 - (1/zeta(2)) * (1 + Sum_{p prime} (1/(p^2-p)) / Product_{primes q <= p} (1 + 1/q)) = 0.11514433883... . - Amiram Eldar, Oct 26 2024

			The prime factors of 300 are {2,2,3,5,5}, with maximal anti-runs ((2),(2,3,5),(5)), with minima (2,2,5), so 300 is in the sequence.
The terms together with their prime indices begin:
    18: {1,2,2}
    36: {1,1,2,2}
    50: {1,3,3}
    54: {1,2,2,2}
    72: {1,1,1,2,2}
    75: {2,3,3}
    90: {1,2,2,3}
    98: {1,4,4}
   100: {1,1,3,3}
   108: {1,1,2,2,2}
   126: {1,2,2,4}
   144: {1,1,1,1,2,2}

Amiram Eldar, Table of n, a(n) for n = 1..10000

A superset of A036785.
The complement for maxima is A065200, counted by A034296.
For maxima instead of minima we have A065201, counted by A239955.
A version for compositions is A374520, counted by A374640.
Also positions of non-constant rows in A375128, sums A374706, ranks A375400.
The complement is A375396, counted by A115029.
The complement for distinct minima is A375398, counted by A375134.
For distinct instead of identical minima we have A375399, counts A375404.
Partitions of this type are counted by A375405.
A000041 counts integer partitions, strict A000009.
A003242 counts anti-run compositions, ranks A333489.
A number's prime factors (A027746, reverse A238689) have sum A001414, min A020639, max A006530.
A number's prime indices (A112798, reverse A296150) have sum A056239, min A055396, max A061395.
Both have length A001222, distinct A001221.
Cf. A000005, A013661, A046660, A272919, A319066, A358905, A374686, A374704, A374742, A375133, A375136, A375401.

Select[Range[100],!SameQ@@Min /@ Split[Flatten[ConstantArray@@@FactorInteger[#]],UnsameQ]&]

is(k) = if(k > 1, my(e = factor(k)[, 2]); vecprod(e) > e[1], 0); \\ Amiram Eldar, Oct 26 2024

Name edited by Peter Munn, May 08 2025

cp's OEIS Frontend

A085986 Squares of the squarefree semiprimes (p^2*q^2).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

PARI

PARI

Python

Formula

A038838 Numbers that are divisible by the square of an odd prime.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

PARI

Formula

A366890 Irregular triangle, wherein row n lists in ascending order all numbers k whose arithmetic derivative k' is equal to the n-th primorial, A002110(n), and that have more than two prime factors with multiplicity. Rows of length zero are simply omitted, i.e., when A369000(n) = 0.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

A375398 Numbers k such that the minima of maximal anti-runs in the weakly increasing sequence of prime factors of k (with multiplicity) are distinct.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A346485 Möbius transform of A342001, where A342001(n) = A003415(n)/A003557(n).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

Formula

A375396 Numbers not divisible by the square of any prime factor except (possibly) the least. Hooklike numbers.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A338539 Numbers having exactly two non-unitary prime factors.