login
Search: a059404 -id:a059404
     Sort: relevance | references | number | modified | created      Format: long | short | data
Number of distinct exponents when n is factorized as a product of primes.
+10
146
0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 1, 1, 1, 2, 1
OFFSET
1,12
COMMENTS
First term greater than 2 is a(360) = 3.
From Michel Marcus, Apr 24 2016: (Start)
A006939(n) gives the least m such that a(m) = n.
A062770 is the sequence of integers m such that a(m) = 1. (End)
We define the k-th omega of n to be Omega(red^{k-1}(n)) where Omega = A001222 and red^{k} is the k-th functional iteration of A181819. The first two omegas are A001222 and A001221, while this sequence is the third, and A323022 is the fourth. The zeroth omega is not uniquely determined from prime signature, but one possible choice is A056239 (sum of prime indices). - Gus Wiseman, Jan 02 2019
Sanna (2020) proved that for each k>=1, the sequence of numbers n with A071625(n) = k has an asymptotic density A_k = (6/Pi^2) * Sum_{n>=1, n squarefree} rho_k(n)/psi(n), where psi is the Dedekind psi function (A001615), and rho_k(n) is defined by rho_1(n) = 1 if n = 1 and 0 otherwise, rho_{k+1}(n) = 0 if n = 1 and (1/(n-1)) * Sum_{d|n, d<n} rho_k(d) otherwise. - Amiram Eldar, Oct 18 2020
LINKS
E. T. Bell, Functions of coprime divisors of integers, Bull. Amer. Math. Soc. 43 (1937), 818-822.
Carlo Sanna, On the number of distinct exponents in the prime factorization of an integer, Proceedings - Mathematical Sciences, Indian Academy of Sciences, Vol. 130, No. 1 (2020), Article 27, alternative link, arXiv preprint, arXiv:1902.09224 [math.NT], 2019.
EXAMPLE
n = 5040 = 2^4*(3*5)^2*7, three different exponents arise:4,2 and 1; so a(5040)=3.
MATHEMATICA
ffi[x_] := Flatten[FactorInteger[x]];
lf[x_] := Length[FactorInteger[x]];
ep[x_] := Table[Part[ffi[x], 2*w], {w, 1, lf[x]}];
Table[Length[Union[ep[w]]], {w, 1, 256}]
(* Second program: *)
{0}~Join~Array[Length@ Union@ FactorInteger[#][[All, -1]] &, 104, 2] (* Michael De Vlieger, Apr 10 2019 *)
PROG
(PARI) a(n) = #Set(factor(n)[, 2]); \\ Michel Marcus, Mar 12 2015
(Python)
from sympy import factorint
def a(n): return len(set(factorint(n).values()))
print([a(n) for n in range(1, 106)]) # Michael S. Branicky, Sep 01 2022
KEYWORD
nonn
AUTHOR
Labos Elemer, May 29 2002
STATUS
approved
Powers of squarefree numbers.
+10
124
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 13, 14, 15, 16, 17, 19, 21, 22, 23, 25, 26, 27, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 46, 47, 49, 51, 53, 55, 57, 58, 59, 61, 62, 64, 65, 66, 67, 69, 70, 71, 73, 74, 77, 78, 79, 81, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97
OFFSET
1,2
COMMENTS
a(n) = A072775(n)^A072776(n); complement of A059404.
Essentially the same as A062770. - R. J. Mathar, Sep 25 2008
Numbers m such that in canonical prime factorization all prime exponents are identical: A124010(m,k) = A124010(m,1) for k = 2..A000005(m). - Reinhard Zumkeller, Apr 06 2014
Heinz numbers of uniform partitions. An integer partition is uniform if all parts appear with the same multiplicity. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). - Gus Wiseman, Apr 16 2018
LINKS
MATHEMATICA
Select[Range[100], Length[Union[FactorInteger[#][[All, 2]]]] == 1 &] (* Geoffrey Critzer, Mar 30 2015 *)
PROG
(Haskell)
import Data.Map (empty, findMin, deleteMin, insert)
import qualified Data.Map.Lazy as Map (null)
a072774 n = a072774_list !! (n-1)
(a072774_list, a072775_list, a072776_list) = unzip3 $
(1, 1, 1) : f (tail a005117_list) empty where
f vs'@(v:vs) m
| Map.null m || xx > v = (v, v, 1) :
f vs (insert (v^2) (v, 2) m)
| otherwise = (xx, bx, ex) :
f vs' (insert (bx*xx) (bx, ex+1) $ deleteMin m)
where (xx, (bx, ex)) = findMin m
-- Reinhard Zumkeller, Apr 06 2014
(PARI) is(n)=ispower(n, , &n); issquarefree(n) \\ Charles R Greathouse IV, Oct 16 2015
(Python)
from math import isqrt
from sympy import mobius, integer_nthroot
def A072774(n):
def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))-1
def f(x): return n-2+x-sum(g(integer_nthroot(x, k)[0]) for k in range(1, 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 19 2024
CROSSREFS
Cf. A072777 (subsequence), A005117, A072778, A329332 (tabular arrangement).
A subsequence of A242414.
KEYWORD
nonn
AUTHOR
Reinhard Zumkeller, Jul 10 2002
STATUS
approved
Positive integers which are neither squarefree integers nor prime powers.
+10
112
12, 18, 20, 24, 28, 36, 40, 44, 45, 48, 50, 52, 54, 56, 60, 63, 68, 72, 75, 76, 80, 84, 88, 90, 92, 96, 98, 99, 100, 104, 108, 112, 116, 117, 120, 124, 126, 132, 135, 136, 140, 144, 147, 148, 150, 152, 153, 156, 160, 162, 164, 168, 171, 172, 175, 176, 180, 184, 188
OFFSET
1,1
LINKS
EXAMPLE
45 is in the sequence because 45=3^2*5, i.e., neither squarefree nor a prime power.
MAPLE
with(numtheory): a:=proc(n) if mobius(n)=0 and nops(factorset(n))>1 then n else fi end: seq(a(n), n=1..230); # Emeric Deutsch, Feb 17 2007
MATHEMATICA
Select[Range[200], Max @@ Last /@ FactorInteger[ # ] >1 && Length[FactorInteger[ # ]] > 1 &] (* Ray Chandler, Feb 17 2007 *)
Select[Range[200], !SquareFreeQ[#]&&!PrimePowerQ[#]&] (* Harvey P. Dale, Aug 05 2023 *)
PROG
(PARI) isok(k) = !issquarefree(k) && !isprimepower(k); \\ Michel Marcus, Nov 02 2022
(Python)
from math import isqrt
from sympy import primepi, integer_nthroot, mobius
def A126706(n):
def f(x): return int(n+sum(primepi(integer_nthroot(x, k)[0]) for k in range(2, x.bit_length()))+sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1)))
m, k = n, f(n)
while m != k:
m, k = k, f(k)
return m # Chai Wah Wu, Aug 15 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
Leroy Quet, Feb 11 2007
EXTENSIONS
Extended by Emeric Deutsch and Ray Chandler, Feb 17 2007
STATUS
approved
Tek's flip: Write n as the product of distinct factors of the form prime(i)^(2^(j-1)) with i and j integers, and replace each such factor with prime(j)^(2^(i-1)).
+10
94
1, 2, 4, 3, 16, 8, 256, 6, 9, 32, 65536, 12, 4294967296, 512, 64, 5, 18446744073709551616, 18, 340282366920938463463374607431768211456, 48, 1024, 131072, 115792089237316195423570985008687907853269984665640564039457584007913129639936, 24, 81, 8589934592, 36, 768
OFFSET
1,2
COMMENTS
This is a multiplicative self-inverse permutation of the integers.
A225547 gives the fixed points.
From Antti Karttunen and Peter Munn, Feb 02 2020: (Start)
This sequence operates on the Fermi-Dirac factors of a number. As arranged in array form, in A329050, this sequence reflects these factors about the main diagonal of the array, substituting A329050[j,i] for A329050[i,j], and this results in many relationships including significant homomorphisms.
This sequence provides a relationship between the operations of squaring and prime shift (A003961) because each successive column of the A329050 array is the square of the previous column, and each successive row is the prime shift of the previous row.
A329050 gives examples of how significant sets of numbers can be formed by choosing their factors in relation to rows and/or columns. This sequence therefore maps equivalent derived sets by exchanging rows and columns. Thus odd numbers are exchanged for squares, squarefree numbers for powers of 2 etc.
Alternative construction: For n > 1, form a vector v of length A299090(n), where each element v[i] for i=1..A299090(n) is a product of those distinct prime factors p(i) of n whose exponent e(i) has the bit (i-1) "on", or 1 (as an empty product) if no such exponents are present. a(n) is then Product_{i=1..A299090(n)} A000040(i)^A048675(v[i]). Note that because each element of vector v is squarefree, it means that each exponent A048675(v[i]) present in the product is a "submask" (not all necessarily proper) of the binary string A087207(n).
This permutation effects the following mappings:
A000035(a(n)) = A010052(n), A010052(a(n)) = A000035(n). [Odd numbers <-> Squares]
A008966(a(n)) = A209229(n), A209229(a(n)) = A008966(n). [Squarefree numbers <-> Powers of 2]
(End)
From Antti Karttunen, Jul 08 2020: (Start)
Moreover, we see also that this sequence maps between A016825 (Numbers of the form 4k+2) and A001105 (2*squares) as well as between A008586 (Multiples of 4) and A028983 (Numbers with even sum of the divisors).
(End)
FORMULA
Multiplicative, with a(prime(i)^j) = A019565(j)^A000079(i-1).
a(prime(i)) = 2^(2^(i-1)).
From Antti Karttunen and Peter Munn, Feb 06 2020: (Start)
a(A329050(n,k)) = A329050(k,n).
a(A329332(n,k)) = A329332(k,n).
Equivalently, a(A019565(n)^k) = A019565(k)^n. If n = 1, this gives a(2^k) = A019565(k).
a(A059897(n,k)) = A059897(a(n), a(k)).
The previous formula implies a(n*k) = a(n) * a(k) if A059895(n,k) = 1.
a(A000040(n)) = A001146(n-1); a(A001146(n)) = A000040(n+1).
a(A000290(a(n))) = A003961(n); a(A003961(a(n))) = A000290(n) = n^2.
a(A000265(a(n))) = A008833(n); a(A008833(a(n))) = A000265(n).
a(A006519(a(n))) = A007913(n); a(A007913(a(n))) = A006519(n).
A007814(a(n)) = A248663(n); A248663(a(n)) = A007814(n).
A048675(a(n)) = A048675(n) and A048675(a(2^k * n)) = A048675(2^k * a(n)) = k + A048675(a(n)).
(End)
From Antti Karttunen and Peter Munn, Jul 08 2020: (Start)
For all n >= 1, a(2n) = A334747(a(n)).
In particular, for n = A003159(m), m >= 1, a(2n) = 2*a(n). [Note that A003159 includes all odd numbers]
(End)
EXAMPLE
7744 = prime(1)^2^(2-1)*prime(1)^2^(3-1)*prime(5)^2^(2-1).
a(7744) = prime(2)^2^(1-1)*prime(3)^2^(1-1)*prime(2)^2^(5-1) = 645700815.
MATHEMATICA
Array[If[# == 1, 1, Times @@ Flatten@ Map[Function[{p, e}, Map[Prime[Log2@ # + 1]^(2^(PrimePi@ p - 1)) &, DeleteCases[NumberExpand[e, 2], 0]]] @@ # &, FactorInteger[#]]] &, 28] (* Michael De Vlieger, Jan 21 2020 *)
PROG
(PARI)
A019565(n) = factorback(vecextract(primes(logint(n+!n, 2)+1), n));
a(n) = {my(f=factor(n)); for (i=1, #f~, my(p=f[i, 1]); f[i, 1] = A019565(f[i, 2]); f[i, 2] = 2^(primepi(p)-1); ); factorback(f); } \\ Michel Marcus, Nov 29 2019
(PARI)
A048675(n) = { my(f = factor(n)); sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2; };
A225546(n) = if(1==n, 1, my(f=factor(n), u=#binary(vecmax(f[, 2])), prods=vector(u, x, 1), m=1, e); for(i=1, u, for(k=1, #f~, if(bitand(f[k, 2], m), prods[i] *= f[k, 1])); m<<=1); prod(i=1, u, prime(i)^A048675(prods[i]))); \\ Antti Karttunen, Feb 02 2020
(Python)
from math import prod
from sympy import prime, primepi, factorint
def A225546(n): return prod(prod(prime(i) for i, v in enumerate(bin(e)[:1:-1], 1) if v == '1')**(1<<primepi(p)-1) for p, e in factorint(n).items()) # Chai Wah Wu, Mar 17 2023
CROSSREFS
Cf. A225547 (fixed points) and the subsequences listed there.
Transposes A329050, A329332.
An automorphism of positive integers under the binary operations A059895, A059896, A059897, A306697, A329329.
An automorphism of A059897 subgroups: A000379, A003159, A016754, A122132.
Permutes lists where membership is determined by number of Fermi-Dirac factors: A000028, A050376, A176525, A268388.
Sequences f that satisfy f(a(n)) = f(n): A048675, A064179, A064547, A097248, A302777, A331592.
Pairs of sequences (f,g) that satisfy a(f(n)) = g(a(n)): (A000265,A008833), (A000290,A003961), (A005843,A334747), (A006519,A007913), (A008586,A334748).
Pairs of sequences (f,g) that satisfy a(f(n)) = g(n), possibly with offset change: (A000040,A001146), (A000079,A019565).
Pairs of sequences (f,g) that satisfy f(a(n)) = g(n), possibly with offset change: (A000035, A010052), (A008966, A209229), (A007814, A248663), (A061395, A299090), (A087207, A267116), (A225569, A227291).
Cf. A331287 [= gcd(a(n),n)].
Cf. A331288 [= min(a(n),n)], see also A331301.
Cf. A331309 [= A000005(a(n)), number of divisors].
Cf. A331590 [= a(a(n)*a(n))].
Cf. A331591 [= A001221(a(n)), number of distinct prime factors], see also A331593.
Cf. A331740 [= A001222(a(n)), number of prime factors with multiplicity].
Cf. A331733 [= A000203(a(n)), sum of divisors].
Cf. A331734 [= A033879(a(n)), deficiency].
Cf. A331735 [= A009194(a(n))].
Cf. A331736 [= A000265(a(n)) = a(A008833(n)), largest odd divisor].
Cf. A335914 [= A038040(a(n))].
A self-inverse isomorphism between pairs of A059897 subgroups: (A000079,A005117), (A000244,A062503), (A000290\{0},A005408), (A000302,A056911), (A000351,A113849 U {1}), (A000400,A062838), (A001651,A252895), (A003586,A046100), (A007310,A000583), (A011557,A113850 U {1}), (A028982,A042968), (A053165,A065331), (A262675,A268390).
A bijection between pairs of sets: (A001248,A011764), (A007283,A133466), (A016825, A001105), (A008586, A028983).
Cf. also A336321, A336322 (compositions with another involution, A122111).
KEYWORD
nonn,mult
AUTHOR
Paul Tek, May 10 2013
EXTENSIONS
Name edited by Peter Munn, Feb 14 2020
"Tek's flip" prepended to the name by Antti Karttunen, Jul 08 2020
STATUS
approved
Fourth omega of n. Number of distinct multiplicities in the prime signature of n.
+10
43
0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
OFFSET
1,60
COMMENTS
The indices of terms greater than 1 are {60, 84, 90, 120, 126, 132, 140, 150, ...}.
First term greater than 2 is a(1801800) = 3. In general, the first appearance of k is a(A182856(k)) = k.
The prime signature of n (row n of A118914) is the multiset of prime multiplicities in n.
We define the k-th omega of n to be Omega(red^{k-1}(n)) where Omega = A001222 and red^{k} is the k-th functional iteration of A181819. The first three omegas are A001222, A001221, A071625, and this sequence is the fourth. The zeroth omega is not uniquely determined from prime signature, but one possible choice is A056239 (sum of prime indices).
LINKS
EXAMPLE
The prime signature of 1286485200 is {1, 1, 1, 2, 2, 3, 4}, in which 1 appears three times, two appears twice, and 3 and 4 both appear once, so there are 3 distinct multiplicities {1, 2, 3} and hence a(1286485200) = 3.
MATHEMATICA
red[n_]:=Times@@Prime/@Last/@If[n==1, {}, FactorInteger[n]];
Table[PrimeNu[red[red[n]]], {n, 200}]
PROG
(PARI) a(n) = my(e=factor(n)[, 2], s = Set(e), m=Map(), v=vector(#s)); for(i=1, #s, mapput(m, s[i], i)); for(i=1, #e, v[mapget(m, e[i])]++); #Set(v) \\ David A. Corneth, Jan 02 2019
(PARI)
A071625(n) = #Set(factor(n)[, 2]); \\ From A071625
A181819(n) = factorback(apply(e->prime(e), (factor(n)[, 2])));
A323022(n) = A071625(A181819(n)); \\ Antti Karttunen, Jan 03 2019
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jan 02 2019
EXTENSIONS
More terms from Antti Karttunen, Jan 03 2019
STATUS
approved
The number obtained by reversing the sequence of nonzero exponents in the prime factorization of n with respect to distinct primes present, as ordered by their indices.
+10
18
1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 18, 13, 14, 15, 16, 17, 12, 19, 50, 21, 22, 23, 54, 25, 26, 27, 98, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 250, 41, 42, 43, 242, 75, 46, 47, 162, 49, 20, 51, 338, 53, 24, 55, 686, 57, 58, 59, 150, 61, 62, 147, 64, 65
OFFSET
1,2
COMMENTS
Equivalent description nearer to the old name: a(n) is a number obtained by reversing the indices of the primes present in the prime factorization of n, from the smallest to the largest, while keeping the nonzero exponents of those same primes at their old positions.
This self-inverse permutation of natural numbers fixes the numbers in whose prime factorization the sequence of nonzero exponents form a palindrome: A242414.
Integers which are changed are A242416.
Considered as a function on partitions encoded by the indices of primes in the prime factorization of n (as in table A112798), this implements an operation which reverses the order of vertical line segments of the "steps" in Young (or Ferrers) diagram of a partition, but keeps the order of horizontal line segments intact. Please see the last example in the example section.
FORMULA
If n = p_a^e_a * p_b^e_b * ... * p_j^e_j * p_k^e_k, where p_a < ... < p_k are distinct primes of the prime factorization of n (sorted into ascending order), and e_a, ..., e_k are their nonzero exponents, then a(n) = p_a^e_k * p_b^e_j * ... * p_j^e_b * p_k^e_a.
a(n) = product(A027748(o(n)+1-k)^A124010(k): k=1..o(n)) = product(A027748(k)^A124010(o(n)+1-k): k=1..o(n)), where o(n) = A001221(n). - Reinhard Zumkeller, Apr 27 2013
From Antti Karttunen, Jun 01 2014: (Start)
Can be obtained also by composing/conjugating related permutations:
a(n) = A242415(A242419(n)) = A242419(A242415(n)).
(End)
EXAMPLE
a(24) = 54 as 24 = p_1^3 * p_2^1 = 2^3 * 3^1 and 54 = p_1^1 * p_2^3 = 2 * 3^3.
For n = 2200, we see that it encodes the partition (1,1,1,3,3,5) in A112798 as 2200 = p_1 * p_1 * p_1 * p_3 * p_3 * p_5 = 2^3 * 5^2 * 11. This in turn corresponds to the following Young diagram in French notation:
_
| |
| |
| |_ _
| |
| |_ _
|_ _ _ _ _|
Reversing the order of vertical line segment lengths (3,2,1) to (1,2,3), but keeping the order of horizontal line segment lengths as (1,2,2), we get a new Young diagram
_
| |_ _
| |
| |_ _
| |
| |
|_ _ _ _ _|
which represents the partition (1,3,3,5,5,5), encoded in A112798 by p_1 * p_3^2 * p_5^3 = 2 * 5^2 * 11^3 = 66550, thus a(2200) = 66550.
MAPLE
A069799 := proc(n) local e, j; e := ifactors(n)[2]:
mul (e[j][1]^e[nops(e)-j+1][2], j=1..nops(e)) end:
seq (A069799(i), i=1..40);
# Peter Luschny, Jan 17 2011
MATHEMATICA
f[n_] := Block[{a = Transpose[ FactorInteger[n]], m = n}, If[ Length[a] == 2, Apply[ Times, a[[1]]^Reverse[a[[2]] ]], m]]; Table[ f[n], {n, 1, 65}]
PROG
(Haskell)
a069799 n = product $
zipWith (^) (a027748_row n) (reverse $ a124010_row n)
-- Reinhard Zumkeller, Apr 27 2013
(MIT/GNU Scheme, with Aubrey Jaffer's SLIB Scheme library)
(require 'factor)
(define (A069799 n) (let ((pf (ifactor n))) (apply * (map expt (uniq pf) (reverse (multiplicities pf))))))
(define (ifactor n) (cond ((< n 2) (list)) (else (sort (factor n) <))))
(define (uniq lista) (let loop ((lista lista) (z (list))) (cond ((null? lista) (reverse! z)) ((and (pair? z) (equal? (car z) (car lista))) (loop (cdr lista) z)) (else (loop (cdr lista) (cons (car lista) z))))))
(define (multiplicities lista) (let loop ((mults (list)) (lista lista) (prev #f)) (cond ((not (pair? lista)) (reverse! mults)) ((equal? (car lista) prev) (set-car! mults (+ 1 (car mults))) (loop mults (cdr lista) prev)) (else (loop (cons 1 mults) (cdr lista) (car lista))))))
;; Antti Karttunen, May 24 2014
(PARI) a(n) = {my(f = factor(n)); my(g = f); my(nbf = #f~); for (i=1, nbf, g[i, 1] = f[nbf-i+1, 1]; ); factorback(g); } \\ Michel Marcus, Jul 02 2015
CROSSREFS
A242414 gives the fixed points and A242416 is their complement.
{A000027, A069799, A242415, A242419} form a 4-group.
The set of permutations {A069799, A105119, A225891} generates an infinite dihedral group.
KEYWORD
nonn
AUTHOR
Amarnath Murthy, Apr 13 2002
EXTENSIONS
Edited, corrected and extended by Robert G. Wilson v and Vladeta Jovovic, Apr 15 2002
Definition corrected by Reinhard Zumkeller, Apr 27 2013
Definition again reworded, Comments section edited and Young diagram examples added by Antti Karttunen, May 30 2014
STATUS
approved
Nonsquarefree numbers that are not squareful.
+10
12
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,1
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)
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
EXAMPLE
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.
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 *)
PROG
(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
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).
KEYWORD
nonn,easy
AUTHOR
Bernard Schott, Feb 24 2020
STATUS
approved
Numbers that require exactly five iterations to reach a fixed point under the x -> A181819(x) map.
+10
9
60, 84, 90, 120, 126, 132, 140, 150, 156, 168, 180, 198, 204, 220, 228, 234, 240, 252, 260, 264, 270, 276, 280, 294, 300, 306, 308, 312, 315, 336, 340, 342, 348, 350, 364, 372, 378, 380, 396, 408, 414, 420, 440, 444, 450, 456, 460, 468, 476, 480, 490, 492, 495
OFFSET
1,1
COMMENTS
In each case, 2 is the fixed point that is reached (1 is the other fixed point of the x -> A181819(x) map).
Includes all integers whose prime signature a) contains two or more distinct numbers, and b) contains no number that occurs the same number of times as any other number. The first member of this sequence that does not fit that description is 75675600, whose prime signature is (4,3,2,2,1,1).
A full characterization is: Numbers whose prime signature (1) has not all equal multiplicities but (2) the numbers of distinct parts appearing with each distinct multiplicity are all equal. For example, the prime signature of 2520 is {1,1,2,3}, which satisfies (1) but fails (2), as the numbers of distinct parts appearing with each distinct multiplicity are 1 (with multiplicity 2, the part being 1) and 2 (with multiplicity 1, the parts being 2 and 3). Hence the sequence does not contain 2520. - Gus Wiseman, Jan 02 2019
LINKS
Eric Weisstein's World of Mathematics, Fixed Point
Eric Weisstein's World of Mathematics, Map
EXAMPLE
1. 180 requires exactly five iterations under the x -> A181819(x) map to reach a fixed point (namely, 2). A181819(180) = 18; A181819(18) = 6; A181819(6) = 4; A181819(4) = 3; A181819(3) = 2 (and A181819(2) = 2).
2. The prime signature of 180 (2^2*3^2*5) is (2,2,1).
a. Two distinct numbers appear in (2,2,1) (namely, 1 and 2).
b. Neither 1 nor 2 appears in (2,2,1) the same number of times as any other number that appears there.
MATHEMATICA
Select[Range[1000], With[{sig=Sort[Last/@FactorInteger[#]]}, And[!SameQ@@Length/@Split[sig], SameQ@@Length/@Union/@GatherBy[sig, Length[Position[sig, #]]&]]]&] (* Gus Wiseman, Jan 02 2019 *)
CROSSREFS
Numbers n such that A182850(n) = 5. See also A182853, A182854.
Subsequence of A059404 and A182851. Includes A085987 and A179642 as subsequences.
KEYWORD
nonn
AUTHOR
Matthew Vandermast, Jan 04 2011
STATUS
approved
Numbers with exactly three distinct exponents in their prime factorization, or three distinct parts in their prime signature.
+10
9
360, 504, 540, 600, 720, 756, 792, 936, 1008, 1176, 1188, 1200, 1224, 1350, 1368, 1400, 1404, 1440, 1500, 1584, 1620, 1656, 1836, 1872, 1960, 2016, 2052, 2088, 2160, 2200, 2232, 2250, 2268, 2352, 2400, 2448, 2484, 2520, 2600, 2646, 2664, 2736, 2800, 2880, 2904
OFFSET
1,1
COMMENTS
Positions of 3's in A071625.
Numbers k such that A001221(A181819(k)) = 3.
The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} r(n)/((n-1)*psi(n)) = 0.030575..., where psi is the Dedekind psi function (A001615), and r(n) = Sum_{d|n, 1<d<n} 1/(d-1) (Sanna, 2020). - Amiram Eldar, Oct 18 2020
LINKS
Carlo Sanna, On the number of distinct exponents in the prime factorization of an integer, Proceedings - Mathematical Sciences, Indian Academy of Sciences, Vol. 130, No. 1 (2020), Article 27, alternative link.
EXAMPLE
1500 = 2^2 * 3^1 * 5^3 has three distinct exponents {1, 2, 3}, so belongs to the sequence.
52500 = 2^2 * 3^1 * 5^4 * 7^1 has three distinct exponents {1, 2, 4}, so belongs to the sequence.
MATHEMATICA
tom[n_]:=Length[Union[Last/@If[n==1, {}, FactorInteger[n]]]];
Select[Range[1000], tom[#]==3&]
PROG
(PARI) is(n) = #Set(factor(n)[, 2]) == 3 \\ David A. Corneth, Jan 02 2019
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jan 02 2019
STATUS
approved
Numbers with exactly four distinct exponents in their prime factorization, or four distinct parts in their prime signature.
+10
8
75600, 105840, 113400, 118800, 126000, 140400, 151200, 158760, 178200, 183600, 198000, 205200, 210600, 211680, 232848, 234000, 237600, 246960, 248400, 252000, 261360, 275184, 275400, 280800, 283500, 294000, 302400, 306000, 307800, 313200, 315000, 334800
OFFSET
1,1
COMMENTS
Positions of 4's in A071625.
Numbers k such that A001221(A181819(k)) = 4.
Is a(n) ~ c * n for some c? - David A. Corneth, Jan 09 2019
The asymptotic density of this sequence is (6/Pi^2) * Sum_{n>=2, n squarefree} r(n)/((n-1)*psi(n)) = 0.00035750... (corresponding to c = 2797.1... in the question above, whose answer is affirmative), where psi is the Dedekind psi function (A001615), and r(n) = Sum_{d_1|n, 1<d_1<n} (1/(d_1-1)) * Sum_{d_2|d_1, 1<d_2<d_1} 1/(d_2-1) (Sanna, 2020). - Amiram Eldar, Oct 18 2020
LINKS
Carlo Sanna, On the number of distinct exponents in the prime factorization of an integer, Proceedings - Mathematical Sciences, Indian Academy of Sciences, Vol. 130, No. 1 (2020), Article 27, alternative link.
EXAMPLE
126000 = 2^4 * 3^2 * 5^3 * 7^1 has four distinct exponents {1, 2, 3, 4}, so belongs to the sequence.
831600 = 2^4 * 3^3 * 5^2 * 7^1 * 11^1 has four distinct exponents {1, 2, 3, 4}, so belongs to the sequence.
MATHEMATICA
tom[n_]:=Length[Union[Last/@If[n==1, {}, FactorInteger[n]]]];
Select[Range[100000], tom[#]==4&]
PROG
(PARI) is(n) = #Set(factor(n)[, 2]) == 4 \\ David A. Corneth, Jan 09 2019
KEYWORD
nonn
AUTHOR
Gus Wiseman, Jan 02 2019
STATUS
approved

Search completed in 0.017 seconds