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 166 results. Next

A380473 Numbers k neither squarefree nor prime power (i.e., in A126706) such that A119288(k) <= A003557(k) < A053669(k) < A006530(k).

Original entry on oeis.org

126, 168, 198, 234, 264, 306, 312, 342, 408, 414, 456, 522, 552, 558, 666, 696, 738, 744, 774, 846, 888, 954, 984, 990, 1032, 1062, 1098, 1128, 1170, 1206, 1272, 1278, 1314, 1320, 1386, 1416, 1422, 1464, 1494, 1530, 1560, 1602, 1608, 1638, 1650, 1704, 1710, 1746
Offset: 1

Views

Author

Michael De Vlieger, Jul 22 2025

Keywords

Comments

Let rad = A007947, p = A119288, q = A053669, g = A006530, and r = A003557.
Numbers k in A126706 such that p <= r < q < g.
Terms are products k of a number s in A033845 and a number t in A007310 with at least one prime power factor p^m | k such that m > 1.

Examples

			Table of n, a(n) for select n:
   n    a(n)                       r   q
  --------------------------------------
   1    126 = 2 * 3^2 * 7          3   5
   2    168 = 2^3 * 3 * 7          4   5
   3    198 = 2 * 3^2 * 11         3   5
   4    234 = 2 * 3^2 * 13         3   5
   5    264 = 2^3 * 3 * 11         4   5
   6    306 = 2 * 3^2 * 17         3   5
   7    312 = 2^3 * 3 * 13         4   5
  24    990 = 2 * 3^2 * 5 * 11     3   7
  29   1170 = 2 * 3^2 * 5 * 13     3   7
  45   1650 = 2 * 3 * 5^2 * 11     5   7
  57   1980 = 2^2 * 3^2 * 5 * 11   6   7
  68   2340 = 2^2 * 3^2 * 5 * 13   6   7
		

Crossrefs

Programs

  • Mathematica
    a053669[x_] := Block[{q = 2}, While[Divisible[x, q], q = NextPrime[q] ]; q];
    s = Select[Range[2^12], Nor[SquareFreeQ[#], PrimePowerQ[#]] &];
    Select[s, And[#3 < #4 < #2[[-1, 1]], #2[[2, 1]] <= #3] & @@
      {#1, #2, #1/Apply[Times, #2[[All, 1]]], a053669[#1]} & @@
      {#, FactorInteger[#]} &]

Formula

Intersection of A364998 and A080259 = A364998 \ A055932 = A364998 \ A369540.

A380691 Number of divisors d | k, d < k/d, such that (d, k/d) are neither unitary nor both coreful, where k is neither squarefree nor prime power (in A126706).

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 2, 1, 1, 3, 1, 1, 2, 2, 2, 1, 1, 3, 1, 1, 3, 2, 2, 2, 1, 4, 1, 1, 2, 2, 3, 3, 1, 1, 4, 1, 2, 2, 2, 2, 2, 4, 1, 1, 2, 2, 1, 2, 4, 3, 1, 4, 1, 1, 1, 3, 5, 2, 1, 2, 5, 2, 2, 3, 2, 1, 3, 1, 4, 2, 4, 2, 2, 2, 2, 1, 6, 1, 1, 1, 2, 2, 5, 2, 1, 4, 1
Offset: 1

Views

Author

Michael De Vlieger, Feb 09 2025

Keywords

Comments

A divisor d | k is said to be coreful if rad(d) = rad(k), where rad = A007947.
In other words, half the number of divisors d | k such that both gcd(d, k/d) > 1 and rad(d) != rad(k/d).
Divisors d and k/d have at least 1 prime factor in common and at least one prime factor q divides one but not the other divisor. Thus, the reference domain S is the intersection of nonsquarefree numbers (k in A013929) and numbers that are not prime powers (k in A024619).
Let S = { prime p : p | d } and let T = { prime p : p | k/d }. Then this sequence counts divisor pairs (d, k/d), d < k/d, such that the symmetric difference of S and T is not empty. For instance, for k = 24 = 2*12 = 4*6, where, in both cases, the product P of the symmetric difference is 3. For k = 180 = 2*90 = 3*60 = 6*30 = 10*18 = 12*15, the products of symmetric differences are 15, 10, 5, 15, and 10, respectively. In the case of 10*18, it is evident that neither rad(10) = rad(180) nor rad(18) = rad(30).

Examples

			Table of n, a(n) listing divisors d and S(n)/d for select values of n:
    n  S(n) a(n)  d*S(n)/d
  ---------------------------------------------------------------------
    1    12   1   2*6
    2    18   1   3*6
    3    20   1   2*10
    4    24   2   2*12, 4*6
    5    28   1   2*14
    6    36   2   2*18, 3*12
    7    40   2   2*20, 4*10
   10    48   3   2*24, 4*12, 6*8
   26    96   4   2*48, 4*24, 6*16, 8*12
   57   180   5   2*90, 3*60, 6*30, 10*18, 12*15
   77   240   6   2*120, 4*60, 6*40, 8*30, 10*24, 12*20
  123   360   8   2*180, 3*120, 4*90, 6*60, 10*36, 12*30, 15*24, 18*20
		

Crossrefs

Programs

  • Mathematica
    Table[1/2*(DivisorSigma[0, k] - 2^PrimeNu[k] - Apply[Times, FactorInteger[k][[All, -1]] - 1]), {k, Select[Range[12, 240], Nor[PrimePowerQ[#], SquareFreeQ[#]] &] }]

Formula

Let tau = A000005, omega = A001221, rad = A007947, and S = A126706.
a(n) = card({ d | k : d < k/d, gcd(d, k/d) > 1, rad(d) != rad(k/d) }), k = S(n).
For k in S(n), a(n) = 1/2 * tau(k) - 2^omega(k) - Product_{p|k} m-1, where p^m | k but p^(m-1) does not divide k.
For k = S(n), a(n) = 1/2 * (A000005(k) - A034444(k) - A361430(k)).

A378108 Primes p such that neither p-1 nor p+1 are in A126706.

Original entry on oeis.org

2, 3, 5, 7, 31, 257, 131071, 618970019642690137449562111, 162259276829213363391578010288127
Offset: 1

Views

Author

Michael De Vlieger, Nov 26 2024

Keywords

Comments

Primes p < 11 are in the sequence since the smallest number in A126706 is 12.
Consider p > 11, odd primes; then both p-1 and p+1 are even. Let j and k be neighbors of p. One neighbor, j, is also divisible by 4, while the other neighbor k is not divisible by 2^m, m > 1. The latter statement implies k cannot be a perfect power q^m, p != q, m > 0, but q^m may divide k.
This sequence is that of primes where j = 2^m and k is squarefree.
Proper subset of A141453.
The neighbor k is also divisible by 3, since abs(p-k) = 1 and neither are divisible by 3. Therefore, 6 | k.

Examples

			17 = 2^4+1 is not in the sequence since 18 = 2 * 3^2.
31 = 2^5-1 is in the sequence since 30 = 2*3*5 is squarefree.
127 = 2^7-1 is not in the sequence because 126 = 2 * 3^2 * 7.
257 = 2^8+1 is in the sequence since 258 = 2*3*43 is squarefree.
8191 = 2^13-1 is not in the sequence because 8190 = 2 * 3^2 * 5 * 7 * 13.
65537 = 2^16+1 is not in the sequence since 65538 = 2 * 3^2 * ll * 331.
131071 = 2^17-1 is in the sequence since 131070 = 2 * 3 * 5 * 17 * 257, etc
		

Crossrefs

Programs

  • Mathematica
    Reap[Do[Which[
      And[PrimeQ[# + 1], SquareFreeQ[(# + 2)/6]], Sow[# + 1],
      And[PrimeQ[# - 1], SquareFreeQ[(# - 2)/6]], Sow[# - 1] ] &[2^i],
    {i, 0, 650}] ][[-1, 1]]

A361102 1 together with numbers having at least two distinct prime factors.

Original entry on oeis.org

1, 6, 10, 12, 14, 15, 18, 20, 21, 22, 24, 26, 28, 30, 33, 34, 35, 36, 38, 39, 40, 42, 44, 45, 46, 48, 50, 51, 52, 54, 55, 56, 57, 58, 60, 62, 63, 65, 66, 68, 69, 70, 72, 74, 75, 76, 77, 78, 80, 82, 84, 85, 86, 87, 88, 90, 91, 92, 93, 94, 95, 96, 98, 99, 100, 102, 104, 105, 106, 108, 110, 111, 112
Offset: 1

Views

Author

Keywords

Comments

This is the union of 1 and A024619. It is the sequence C used in the definition of A360519. Since C is central to the analysis of A360519 it deserves its own entry.
This has the same relationship to A024619 as A000469 does to A120944 for squarefree numbers.

Crossrefs

Programs

  • Maple
    isa := n -> is(irem(ilcm(seq(1..n-1)), n) = 0):
    aList := upto -> select(isa, [seq(1..upto)]):
    aList(112); # Peter Luschny, May 17 2023
  • Mathematica
    Select[Range[120], Not@*PrimePowerQ] (* Michael De Vlieger, May 17 2023 *)
  • Python
    from sympy import primepi, integer_nthroot
    def A361102(n):
        def f(x): return int(n+sum(primepi(integer_nthroot(x,k)[0]) for k in range(1,x.bit_length())))
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        return bisection(f) # Chai Wah Wu, Aug 31 2024
  • SageMath
    def A361102List(upto: int) -> list[int]:
        return sorted(Set(1..upto).difference(prime_powers(upto)))
    print(A361102List(112))  # Peter Luschny, May 17 2023
    

Formula

From Peter Luschny and Michael De Vlieger, May 17 2023: (Start)
The sequence is the complement of the prime powers in the positive integers, a = A000027 \ A246655.
k is in this sequence <=> k divides lcm(1, 2, ..., k-1). (End)
This sequence is {1} U { A120944 U A126706 } = {1} U A024619. - Michael De Vlieger, May 17 2023

Extensions

Offset set to 1 by Peter Luschny, May 17 2023

A131605 Perfect powers of nonprimes (m^k where m is a nonprime positive integer and k >= 2).

Original entry on oeis.org

1, 36, 100, 144, 196, 216, 225, 324, 400, 441, 484, 576, 676, 784, 900, 1000, 1089, 1156, 1225, 1296, 1444, 1521, 1600, 1728, 1764, 1936, 2025, 2116, 2304, 2500, 2601, 2704, 2744, 2916, 3025, 3136, 3249, 3364, 3375, 3600, 3844, 3969, 4225, 4356, 4624
Offset: 1

Views

Author

Daniel Forgues, May 27 2008

Keywords

Comments

Although 1 is a square, is a cube, and so on..., 1 is sometimes excluded from perfect powers since it is not a well-defined power of 1 (1 = 1^k for any k in [2, 3, 4, 5, ...])
From Michael De Vlieger, Aug 11 2025: (Start)
This sequence is A001597 \ A246547, i.e., perfect powers without proper prime powers.
Union of {1} with the intersection of A001597 and A126706, where A126706 is the sequence of numbers that are neither prime powers nor squarefree.
Union of {1} and A286708 \ A052486, i.e., powerful numbers that are not prime powers, without Achilles numbers, but including the empty product. (End)

Crossrefs

Programs

  • Mathematica
    With[{nn = 2^20}, {1}~Join~Select[Union@ Flatten@ Table[a^2*b^3, {b, Surd[nn, 3]}, {a, Sqrt[nn/b^3]}], And[Length[#2] > 1, GCD @@ #2 > 1] & @@ {#, FactorInteger[#][[;; , -1]]} &] ] (* Michael De Vlieger, Aug 11 2025 *)
  • PARI
    isok(n) = if (n == 1, return (1), return (ispower(n, ,&np) && (! isprime(np)))); \\ Michel Marcus, Jun 12 2013
    
  • Python
    from sympy import mobius, integer_nthroot, primepi
    def A131605(n):
        def f(x): return int(n-2+x+sum(mobius(k)*((a:=integer_nthroot(x,k)[0])-1)+primepi(a) for k in range(2,x.bit_length())))
        kmin, kmax = 1,2
        while f(kmax) >= kmax:
            kmax <<= 1
        while True:
            kmid = kmax+kmin>>1
            if f(kmid) < kmid:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmax # Chai Wah Wu, Aug 14 2024

Formula

Sum_{n>=1} 1/a(n) = 1 + A072102 - A136141 = 1.10130769935514973882... . - Amiram Eldar, Aug 15 2025

A360519 Let C consist of 1 together with all numbers with at least two distinct prime factors; this is the lexicographically earliest infinite sequence {a(n)} of distinct elements of C such that, for n>2, a(n) has a common factor with a(n-1) but not with a(n-2).

Original entry on oeis.org

1, 6, 10, 35, 21, 12, 20, 55, 33, 18, 14, 77, 99, 15, 40, 22, 143, 39, 24, 28, 91, 65, 30, 34, 119, 63, 36, 26, 221, 51, 42, 38, 95, 45, 48, 44, 187, 85, 50, 46, 69, 57, 76, 52, 117, 75, 70, 58, 87, 93, 62, 56, 105, 111, 74, 68, 153, 123, 82, 80, 115, 161, 84, 60, 145, 203, 98, 54, 129, 215, 100, 66, 141
Offset: 1

Views

Author

Keywords

Comments

In other words, C contains all positive numbers except powers of primes p^k, k>=1.
This is a modified version of the Enots Wolley sequence A336957. The modification ensures that the sequence does not contain the prime 2.
Let Ker(k), the kernel of k, denote the set of primes dividing k. Thus Ker(36) = {2,3}, Ker(1) = {}.
Theorem: a(1)=1, a(2)=6; thereafter, a(n) is the smallest number m not yet in the sequence such that
(i) Ker(m) intersect Ker(a(n-1)) is nonempty,
(ii) Ker(m) intersect Ker(a(n-2)) is empty, and
(iii) The set Ker(m) \ Ker(a(n-1)) is nonempty.
Conjecture: The sequence is a permutation of C.

Crossrefs

For a number of sequences related to this, see A361102 (the sequence C) and the following entries.

Programs

  • Maple
    with(numtheory);
    N:= 10^4: # to get a(1) to a(n) where a(n+1) is the first term > N
    B:= Vector(N, datatype=integer[4]):
    A[1]:=1; A[2]:=6;
    for n from 3 do
      for k from 10 to N do
        if B[k] = 0 and igcd(k, A[n-1]) > 1 and igcd(k, A[n-2]) = 1 then
              if nops(factorset(k) minus factorset(A[n-1])) > 0 then
           A[n]:= k;
           B[k]:= 1;
           break;
              fi;
        fi
      od:
      if k > N then break; fi;
    od:
    s1:=[seq(A[i], i=1..n-1)];
  • Mathematica
    nn = 2^12; c[_] = False;
    f[x_] := f[x] = Times @@ FactorInteger[x][[All, 1]];
    MapIndexed[
     Set[{a[First[#2]], c[#1]}, {#1, True}] &, {1, 6}];
     u = 10; i = a[1]; j = a[2];
    Do[k = u;
      While[Nand[! PrimePowerQ[k], ! c[k],
        CoprimeQ[i, k], ! CoprimeQ[j, k], ! Divisible[j, f[k]]], k++];
      Set[{a[n], c[k], i, j}, {k, True, j, f[k]}];
      If[k == u, While[Or[c[u], PrimePowerQ[u]], u++]]
      , {n, 3, nn}];
    Array[a, nn] (* Michael De Vlieger, Mar 03 2023 *)

A340596 Number of co-balanced factorizations of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 1, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 4, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 5, 1, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 4, 1, 1, 1, 3, 1, 4, 1, 2, 1, 1, 1, 5, 1, 2, 2, 4, 1, 1, 1, 3, 1, 1, 1, 5, 1, 1, 1, 4, 1, 1, 1, 2, 2, 1, 1, 8
Offset: 1

Views

Author

Gus Wiseman, Jan 15 2021

Keywords

Comments

We define a factorization of n into factors > 1 to be co-balanced if it has exactly A001221(n) factors.

Examples

			The a(n) co-balanced factorizations for n = 12, 24, 36, 72, 120, 144, 180:
  2*6    3*8     4*9     8*9     3*5*8     2*72     4*5*9
  3*4    4*6     6*6     2*36    4*5*6     3*48     5*6*6
         2*12    2*18    3*24    2*2*30    4*36     2*2*45
                 3*12    4*18    2*3*20    6*24     2*3*30
                         6*12    2*4*15    8*18     2*5*18
                                 2*5*12    9*16     2*6*15
                                 2*6*10    12*12    2*9*10
                                 3*4*10             3*3*20
                                                    3*4*15
                                                    3*5*12
                                                    3*6*10
		

Crossrefs

Positions of terms > 1 are A126706.
Positions of 1's are A303554.
The version for unlabeled multiset partitions is A319616.
The alt-balanced version is A340599.
The balanced version is A340653.
The cross-balanced version is A340654.
The twice-balanced version is A340655.
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
Other balance-related sequences:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A098124 counts balanced compositions.
- A106529 lists Heinz numbers of balanced partitions.
- A340597 lists numbers with an alt-balanced factorization.
- A340598 counts balanced set partitions.
- A340600 counts unlabeled balanced multiset partitions.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],Length[#]==PrimeNu[n]&]],{n,100}]
  • PARI
    A340596(n, m=n, om=omega(n)) = if(1==n,(0==om), sumdiv(n, d, if((d>1)&&(d<=m), A340596(n/d, d, om-1)))); \\ Antti Karttunen, Jun 10 2024

Extensions

Data section extended up to a(120) by Antti Karttunen, Jun 10 2024

A386482 a(1)=1, a(2)=2; thereafter a(n) is either the greatest number k < a(n-1) not already used such that gcd(k, a(n-1)) > 1, or if no such k exists then a(n) is the smallest number k > a(n-1) not already used such that gcd(k, a(n-1)) > 1.

Original entry on oeis.org

1, 2, 4, 6, 3, 9, 12, 10, 8, 14, 7, 21, 18, 16, 20, 15, 5, 25, 30, 28, 26, 24, 22, 11, 33, 27, 36, 34, 32, 38, 19, 57, 54, 52, 50, 48, 46, 44, 42, 40, 35, 45, 39, 13, 65, 60, 58, 56, 49, 63, 51, 17, 68, 66, 64, 62, 31, 93, 90, 88, 86, 84, 82, 80, 78, 76, 74, 72, 70, 55, 75, 69, 23, 92, 94, 47, 141, 138, 136, 134, 132, 130, 128
Offset: 1

Views

Author

N. J. A. Sloane, Aug 15 2025, based on email messages from Geoffrey Caveney

Keywords

Comments

Similar to the EKG sequence A064413, but whereas in that sequence a(n) is chosen to be as small as possible, here the primary goal is to choose a(n) to be less than a(n-1) and as close to it as possible. This sequence first differs from the EKG sequence at n = 8, where a(8) = k = 10 is closer to a(7) = 12 than A064413(8) = 8 is.
A significant difference from the EKG sequence is that the primes do not appear in their natural order. Also, it is not always true that a prime p is preceded by 2*p when it first appears. 4k+3 primes appear to be preceded by smaller multiples than 4k+1 primes.
It is conjectured that every positive number appears.
It is interesting to study what happens if the first two terms are taken to be 1,s, with s >= 2, or if the first s terms are taken to be 1,2,3,...,s, with s >= 2. Call two such sequences equivalent if they eventually merge. The 1,3 and 1,2,3 sequences merge with each other after half-a-dozen terms. But at present we do not know if they merge with the 1,2 sequence.
It appears that many sequences that start 1,s and 1,2,3,...,s with small s merge with one of the sequences 1,2 or 1,2,3 or 1,2,3,...,11.
[The preceding comments are from Geoffrey Caveney's emails.]
From Michael De Vlieger, Aug 15 2025: (Start)
There are long runs of terms with the same parity in this sequence. For example, beginning at a(481) = 948, there are 100 consecutive even terms. Starting with a(730076) = 1026330, there are 100869 consecutive even terms, followed by 36709 consecutive odd terms. Runs of even terms tend to be longer than those of odd.
There are long runs of first differences of -2 and -6 in this sequence, and that there appear to be three phases. The predominant (A) phase has a(n) = a(n-1)-2, the second (B) phase has a(n) = a(n-1)-6, and then there is a turbulent (C) phase [C] with varied differences.
Generally the even runs correspond to differences a(n)-a(n-1) = 2 and feature square-free terms separated by an odd number of terms in A126706. Phase [C] tends to be largely odd squarefree semiprimes and includes prime powers. (End)

References

  • Geoffrey Caveney, Emails to N. J. A. Sloane, Aug 13 2025 - Aug 15 2025.

Crossrefs

Cf. A064413 (EKG), A387072 (inverse), A387073 (record high points), A387074 (indices of record high points), A387075 (first differences), A387076 (primes in order of appearance), A387077 (indices of primes), A387078 (run lengths of consecutive odd and even terms), A387080 (variant that begins with 1,3).

Programs

  • Mathematica
    aList[n_] := Module[{an = 2, aset = <|2 -> True|>, m}, Reap[Sow[1]; Sow[an];
    Do[m = SelectFirst[Range[an - 1, 2, -1], ! KeyExistsQ[aset, #] && GCD[#, an] > 1 & ];
    If[MissingQ[m], m = NestWhile[# + 1 &, an + 1, !(! KeyExistsQ[aset, #] && GCD[#, an] > 1) & ]];
    aset[m] = True; an = m; Sow[an], {n - 2}]][[2, 1]]]; aList[83]  (* Peter Luschny, Aug 15 2025 *)
  • PARI
    \\ See Links section.
    
  • Python
    from math import gcd
    from itertools import count, islice
    def A386482_gen(): # generator of terms
        yield 1
        an, aset = 2, {2}
        while True:
            yield an
            m = next((k for k in range(an-1, 1, -1) if k not in aset and gcd(k, an) > 1), False)
            if not m: m = next(k for k in count(an+1) if k not in aset and gcd(k, an) > 1)
            an = m
            aset.add(an)
    print(list(islice(A386482_gen(), 83))) # Michael S. Branicky, Aug 15 2025

A332785 Nonsquarefree numbers that are not squareful.

Original entry on oeis.org

12, 18, 20, 24, 28, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 68, 75, 76, 80, 84, 88, 90, 92, 96, 98, 99, 104, 112, 116, 117, 120, 124, 126, 132, 135, 136, 140, 147, 148, 150, 152, 153, 156, 160, 162, 164, 168, 171, 172, 175, 176, 180, 184, 188, 189, 192, 198, 204, 207, 208, 212, 220, 224
Offset: 1

Views

Author

Bernard Schott, Feb 24 2020

Keywords

Comments

Sometimes nonsquarefree numbers are misnamed squareful numbers (see 1st comment of A013929). Indeed, every squareful number > 1 is nonsquarefree, but the converse is false. This sequence = A013929 \ A001694 and consists of these counterexamples.
This sequence is not a duplicate: the first 16 terms (<= 68) are the same first 16 terms of A059404, A323055, A242416 and A303946, then 72 is the 17th term of these 4 sequences. Also, the first 37 terms (<= 140) are the same first 37 terms of A317616 then 144 is the 38th term of this last sequence.
From Amiram Eldar, Sep 17 2023: (Start)
Called "hybrid numbers" by Jakimczuk (2019).
These numbers have a unique representation as a product of two numbers > 1, one is squarefree (A005117) and the other is powerful (A001694).
Equivalently, numbers k such that A055231(k) > 1 and A057521(k) > 1.
Equivalently, numbers that have in their prime factorization at least one exponent that is equal to 1 and at least one exponent that is larger than 1.
The asymptotic density of this sequence is 1 - 1/zeta(2) (A229099). (End)

Examples

			18 = 2 * 3^2 is nonsquarefree as it is divisible by the square 3^2, but it is not squareful because 2 divides 18 but 2^2 does not divide 18, hence 18 is a term.
72 = 2^3 * 3^2 is nonsquarefree as it is divisible by the square 3^2, but it is also squareful because primes 2 and 3 divide 72, and 2^2 and 3^2 divide also 72, so 72 is not a term.
		

Crossrefs

Cf. A005117 (squarefree), A013929 (nonsquarefree), A001694 (squareful), A052485 (not squareful).
Cf. A059404, A126706, A229099, A242416, A286708, A303946, A317616, A323055 (first terms are the same).

Programs

  • Maple
    filter:= proc(n) local F;
     F:= ifactors(n)[2][..,2];
     max(F) > 1 and min(F) = 1
    end proc:
    select(filter, [$1..1000]); # Robert Israel, Sep 15 2024
  • Mathematica
    Select[Range[225], Max[(e = FactorInteger[#][[;;,2]])] > 1 && Min[e] == 1 &] (* Amiram Eldar, Feb 24 2020 *)
  • PARI
    isok(m) = !issquarefree(m) && !ispowerful(m); \\ Michel Marcus, Feb 24 2020
    
  • Python
    from math import isqrt
    from sympy import mobius, integer_nthroot
    def A332785(n):
        def squarefreepi(n): return int(sum(mobius(k)*(n//k**2) for k in range(1, isqrt(n)+1)))
        def bisection(f,kmin=0,kmax=1):
            while f(kmax) > kmax: kmax <<= 1
            while kmax-kmin > 1:
                kmid = kmax+kmin>>1
                if f(kmid) <= kmid:
                    kmax = kmid
                else:
                    kmin = kmid
            return kmax
        def f(x):
            c, l, j = n-1+squarefreepi(integer_nthroot(x,3)[0])+squarefreepi(x), 0, isqrt(x)
            while j>1:
                k2 = integer_nthroot(x//j**2,3)[0]+1
                w = squarefreepi(k2-1)
                c += j*(w-l)
                l, j = w, isqrt(x//k2**3)
            return c-l
        return bisection(f,n,n) # Chai Wah Wu, Sep 14 2024

Formula

This sequence is A126706 \ A286708.
Sum_{n>=1} 1/a(n)^s = 1 + zeta(s) - zeta(s)/zeta(2*s) - zeta(2*s)*zeta(3*s)/zeta(6*s), s > 1. - Amiram Eldar, Sep 17 2023

A340654 Number of cross-balanced factorizations of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 3, 1, 1, 1, 2, 2, 1, 1, 4, 1, 2, 1, 2, 1, 3, 1, 3, 1, 1, 1, 3, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 6, 1, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 3, 1, 1, 1, 3, 1, 3, 1, 2, 1, 1, 1, 5, 1, 2, 2, 5, 1, 1, 1, 3, 1
Offset: 1

Views

Author

Gus Wiseman, Jan 15 2021

Keywords

Comments

We define a factorization of n into factors > 1 to be cross-balanced if either (1) it is empty or (2) the maximum image of A001222 over the factors is A001221(n).

Examples

			The cross-balanced factorizations for n = 12, 24, 36, 72, 144, 240:
  2*6   4*6     4*9     2*4*9     4*4*9       8*30
  3*4   2*2*6   6*6     2*6*6     4*6*6       12*20
        2*3*4   2*2*9   3*4*6     2*2*4*9     5*6*8
                2*3*6   2*2*2*9   2*2*6*6     2*4*30
                3*3*4   2*2*3*6   2*3*4*6     2*6*20
                        2*3*3*4   3*3*4*4     2*8*15
                                  2*2*2*2*9   3*4*20
                                  2*2*2*3*6   3*8*10
                                  2*2*3*3*4   4*5*12
                                              2*10*12
                                              2*3*5*8
                                              2*2*2*30
                                              2*2*3*20
                                              2*2*5*12
		

Crossrefs

Positions of terms > 1 are A126706.
Positions of 1's are A303554.
The co-balanced version is A340596.
The version for unlabeled multiset partitions is A340651.
The balanced version is A340653.
The twice-balanced version is A340655.
A001055 counts factorizations.
A045778 counts strict factorizations.
A316439 counts factorizations by product and length.
A320655 counts factorizations into semiprimes.
Other balance-related sequences:
- A010054 counts balanced strict partitions.
- A047993 counts balanced partitions.
- A098124 counts balanced compositions.
- A106529 lists Heinz numbers of balanced partitions.
- A340597 have an alt-balanced factorization.
- A340598 counts balanced set partitions.
- A340599 counts alt-balanced factorizations.
- A340652 counts unlabeled twice-balanced multiset partitions.
- A340656 have no twice-balanced factorizations.
- A340657 have a twice-balanced factorization.

Programs

  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[n],#=={}||PrimeNu[n]==Max[PrimeOmega/@#]&]],{n,100}]
  • PARI
    A340654(n, m=n, om=omega(n),mbo=0) = if(1==n,(mbo==om), sumdiv(n, d, if((d>1)&&(d<=m), A340654(n/d, d, om, max(mbo,bigomega(d)))))); \\ Antti Karttunen, Jun 19 2024

Extensions

Data section extended up to a(105) by Antti Karttunen, Jun 19 2024
Previous Showing 21-30 of 166 results. Next