A000188 -id:A000188 - cp's OEIS frontend

A056169 Number of unitary prime divisors of n.

0, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 1, 0, 2, 0, 1, 1, 3, 1, 0, 2, 2, 2, 0, 1, 2, 2, 1, 1, 3, 1, 1, 1, 2, 1, 1, 0, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 0, 2, 3, 1, 1, 2, 3, 1, 0, 1, 2, 1, 1, 2, 3, 1, 1, 0, 2, 1, 2, 2, 2, 2, 1, 1, 2, 2, 1, 2, 2, 2, 1, 1, 1, 1, 0, 1, 3, 1, 1, 3

Offset: 1

Labos Elemer, Jul 27 2000

The zeros of this sequences are the powerful numbers (A001694). There are no arbitrarily long subsequences with a given upper bound; for example, every sequence of 4 values includes one divisible by 2 but not 4, so there are no more than 3 consecutive zeros. Similarly, there can be no more than 23 consecutive values with none divisible by both 2 and 3 but neither 4 nor 9 (so a(n) >= 2), etc. In general, this gives an upper bound that is a (relatively) small multiple of the k-th primorial number (prime(k)#). One suspects that the actual upper bounds for such subsequences are quite a bit lower; e.g., Erdős conjectured that there are no three consecutive powerful numbers. - Franklin T. Adams-Watters, Aug 08 2006
In particular, for every A048670(k)*A002110(k) consecutive terms, at least one is greater than or equal to k. - Charlie Neder, Jan 03 2019
Following Catalan's conjecture (which became Mihăilescu's theorem in 2002), the first case of two consecutive zeros in this sequence is for a(8) and a(9), because 8 = 2^3 and 9 = 3^2, and there are no other consecutive zeros for consecutive powers. However, there are other pairs of consecutive zeros at powerful numbers (A001694, A060355). The next example is a(288) = a(289) = 0, because 288 = 2^5 * 3^2 and 289 = 17^2, then also a(675) and a(676). - Bernard Schott, Jan 06 2019
a(2k-1) is the number of primes p such that p || x + y and p^2 || x^(2k-1) + y^(2k-1) for some positive integers x and y. For any positive integers x, y and k > 1, there is no prime p such that p || x + y and p^2 || x^(2k) + y^(2k). - Jinyuan Wang, Apr 08 2020

			9 = 3^2 so a(9) = 0; 10 = 2 * 5 so a(10) = 2; 11 = 11^1 so a(11) = 1.

T. D. Noe, Table of n, a(n) for n = 1..10000
Eric Weisstein's World of Mathematics, Catalan's Conjecture.
Index entries for sequences computed from exponents in factorization of n.

Cf. A001221, A001694, A002110, A034444, A056170, A055231, A076445, A162642, A275812, A295659, A295662, A295664, A001694.
See also A060355.
Cf. A077761, A085548.

a056169 = length . filter (== 1) . a124010_row
-- Reinhard Zumkeller, Sep 10 2013

a:= n-> nops(select(i-> i[2]=1, ifactors(n)[2])):
seq(a(n), n=1..120);  # Alois P. Heinz, Mar 27 2017

Join[{0},Table[Count[Transpose[FactorInteger[n]][[2]],1],{n,2,110}]] (* Harvey P. Dale, Mar 15 2012 *)
Table[DivisorSum[n, 1 &, And[PrimeQ@ #, CoprimeQ[#, n/#]] &], {n, 105}] (* Michael De Vlieger, Nov 28 2017 *)

a(n)=my(f=factor(n)[,2]); sum(i=1,#f,f[i]==1) \\ Charles R Greathouse IV, Apr 29 2015

from sympy import factorint
def a(n):
    f=factorint(n)
    return 0 if n==1 else sum(1 for i in f if f[i]==1)
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Jun 19 2017

;; With memoization-macro definec.
(definec (A056169 n) (if (= 1 n) 0 (+ (if (= 1 (A067029 n)) 1 0) (A056169 (A028234 n))))) ;; Antti Karttunen, Nov 28 2017

A prime factor of n is unitary iff its exponent is 1 in prime factorization of n. In general, gcd(p, n/p) = 1 or = p.
Additive with a(p^e) = 1 if e = 1, 0 otherwise.
a(n) = #{k: A124010(n,k) = 1, k = 1..A001221}. - Reinhard Zumkeller, Sep 10 2013
From Antti Karttunen, Nov 28 2017: (Start)
a(1) = 0; for n > 1, a(n) = A063524(A067029(n)) + a(A028234(n)).
a(n) = A001221(A055231(n)) = A001222(A055231(n)).
a(n) = A001221(n) - A056170(n) = A001221(n) - A001221(A000188(n)).
a(n) = A001222(n) - A275812(n).
a(n) = A162642(n) - A295662(n).
a(n) <= A162642(n) <= a(n) + A295659(n).
a(n) <= A295664(n).
(End)
Sum_{k=1..n} a(k) ~ n * (log(log(n)) + B - C), where B is Mertens's constant (A077761) and C = Sum_{p prime} (1/p^2) = 0.452247... (A085548). - Amiram Eldar, Sep 28 2023

A055231 Powerfree part of n: product of primes that divide n only once.

Original entry on oeis.org

1, 2, 3, 1, 5, 6, 7, 1, 1, 10, 11, 3, 13, 14, 15, 1, 17, 2, 19, 5, 21, 22, 23, 3, 1, 26, 1, 7, 29, 30, 31, 1, 33, 34, 35, 1, 37, 38, 39, 5, 41, 42, 43, 11, 5, 46, 47, 3, 1, 2, 51, 13, 53, 2, 55, 7, 57, 58, 59, 15, 61, 62, 7, 1, 65, 66, 67, 17, 69, 70, 71, 1, 73, 74, 3, 19, 77, 78, 79, 5

Offset: 1

Labos Elemer, Jun 21 2000

The previous name was: Write n = K^2*F where F is squarefree and F = g*f where g = gcd(K,F) and f = F/g; then a(n) = f(n) = F(n)/g(n). Thus gcd(K^2,f) = 1.
Differs from A007913; they coincide if and only if g(n) = 1.
a(n) is the powerfree part of n; i.e., if n=Product(pi^ei) over all i (prime factorization) then a(n)=Product(pi^ei) over those i with ei=1; if n=b*c^2*d^3 then a(n) is minimum possible value of b. - Henry Bottomley, Sep 01 2000
Also denominator of n/rad(n)^2, where rad is the squarefree kernel of n (A007947), numerator: A062378. - Reinhard Zumkeller, Dec 10 2002
Largest unitary squarefree number dividing n (the unitary squarefree kernel of n). - Steven Finch, Mar 01 2004
From Bernard Schott, Dec 19 2022: (Start)
a(n) = 1 iff n is a squareful number (A001694).
1 < a(n) < n iff n is a nonsquarefree number that is not squareful (A332785).
a(n) = n iff n is a squarefree number (A005117). (End)

Antti Karttunen, Table of n, a(n) for n = 1..20000 (first 1000 terms from T. D. Noe)
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
Vaclav Kotesovec, Plot of Sum_{k=1..n} a(k) / n^2 for n = 1..1000000

Positions of 1's: A001694.
Cf. A008833, A007913, A007947, A000188, A057521, A055773 (computed for n!), A056169 (number of prime divisors), A056671 (number of divisors), A092261 (sum of divisors of the n-th term), A197863, A332785.
Cf. A005117 (subsequence).

A055231 := proc(n)
    a := 1 ;
    if n > 1 then
        for f in ifactors(n)[2] do
            if op(2, f) = 1 then
                a := a*op(1, f) ;
            end if;
        end do:
    end if;
    a ;
end proc: # R. J. Mathar, Dec 23 2011

rad[n_] := Times @@ First /@ FactorInteger[n]; a[n_] := Denominator[n/rad[n]^2]; Table[a[n], {n, 1, 80}] (* Jean-François Alcover, Jun 20 2013, after Reinhard Zumkeller *)
f[p_, e_] := If[e==1, p, 1]; a[n_] := Times @@ (f @@@ FactorInteger[n]); Array[a, 100] (* Amiram Eldar, Sep 07 2020 *)

A055231(n)={
   local(a=1);
   f=factor(n) ;
   for(i=1,matsize(f)[1],
         if( f[i,2] ==1, a *=  f[i,1]
         )
   ) ;
   a ;
} /* R. J. Mathar, Mar 12 2012 */

a(n) = {my(f=factor(n)); for (k=1, #f~, if (f[k,2] > 1, f[k,2] = 0);); factorback(f);} \\ Michel Marcus, Aug 27 2017

from math import prod
from sympy import factorint
def A055231(n): return prod(p for p, e in factorint(n).items() if e == 1) # Chai Wah Wu, Nov 14 2022

;; With memoization-macro definec.
(definec (A055231 n) (if (= 1 n) 1 (* (if (= 1 (A067029 n)) (A020639 n) 1) (A055231 (A028234 n))))) ;; Antti Karttunen, Nov 28 2017

a(n) = A007913(n)/gcd(A008833(n), A007913(n)).
a(n) = n/A057521(n).
Multiplicative with a(p) = p and a(p^e) = 1 for e > 1. - Vladeta Jovovic, Nov 01 2001
Dirichlet g.f.: zeta(s)*Product_{primes p} (1 + p^(1-s) - p^(-s) - p^(1-2s) + p^(-2s)). - R. J. Mathar, Dec 21 2011
a(n) = A007947(n)/A071773(n). - observed by Velin Yanev, Aug 27 2017, confirmed by Antti Karttunen, Nov 28 2017
a(1) = 1; for n > 1, a(n) = A020639(n)^A063524(A067029(n)) * a(A028234(n)). - Antti Karttunen, Nov 28 2017
a(n*m) = a(n)*a(m)/(gcd(n,a(m))*gcd(m,a(n))) for all n and m > 0 (conjectured). - Velin Yanev, Feb 06 2019. [This follows easily from the comment of Vladeta Jovovic. - N. J. A. Sloane, Mar 14 2019]
From Vaclav Kotesovec, Dec 19 2019: (Start)
Dirichlet g.f.: zeta(s-1) * zeta(s) * Product_{primes p} (1 - p^(1-3*s) + p^(2-3*s) - p^(2-2*s) + p^(-2*s) - p^(-s)).
Sum_{k=1..n} a(k) ~ c * Pi^2 * n^2 / 12, where c = Product_{primes p} (1 - 2/p^2 + 2/p^4 - 1/p^5) = 0.394913518073109872954607634745304266741971541072... (End)
a(n) = A197863(n)/n. - Amiram Eldar, Sep 01 2023

Name replaced with a simpler description (based on Henry Bottomley's comment) by Antti Karttunen, Nov 28 2017

Incorrect comments and example deleted by Peter Munn, Nov 30 2022

A036537 Numbers whose number of divisors is a power of 2.

Original entry on oeis.org

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

Offset: 1

Labos Elemer

nonn

Primes and A030513(d(x)=4) are subsets; d(16k+4) and d(16k+12) have the form 3Q, so x=16k+4 or 16k-4 numbers are missing.
A number m is a term if and only if all its divisors are infinitary, or A000005(m) = A037445(m). - Vladimir Shevelev, Feb 23 2017
All exponents in the prime number factorization of a(n) have the form 2^k-1, k >= 1. So it is an S-exponential sequence (see Shevelev link) with S={2^k-1}. Using Theorem 1, we obtain that a(n) ~ C*n, where C = Product((1-1/p)*(1 + Sum_{i>=1} 1/p^(2^i-1))). - Vladimir Shevelev Feb 27 2017
This constant is C = 0.687827... . - Peter J. C. Moses, Feb 27 2017
From Peter Munn, Jun 18 2022: (Start)
1 and numbers j*m^2, j squarefree, m >= 1, such that all prime divisors of m divide j, and m is in the sequence.
Equivalently, the nonempty set of numbers whose squarefree part (A007913) and squarefree kernel (A007947) are equal, and whose square part's square root (A000188) is in the set.
(End)

			383, 384, 385, 386 have 2, 16, 8, 4 divisors, respectively, so they are consecutive terms of this sequence.

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Vladimir Shevelev, S-exponential numbers, Acta Arithm., 175(2016), 385-395.
Eric Weisstein's World of Mathematics, Square Part, Squarefree Part

A005117, A030513, A058891, A175496, A336591 are subsequences.
Complement of A162643; subsequence of A002035. - Reinhard Zumkeller, Jul 08 2009
Subsequence of A162644, A337533.
Cf. A000005, A000188, A007913, A007947, A036538, A352780.
The closure of the squarefree numbers under application of A355038(.) and lcm.

a036537 n = a036537_list !! (n-1)
a036537_list = filter ((== 1) . a209229 . a000005) [1..]
-- Reinhard Zumkeller, Nov 15 2012

bi[ x_ ] := 1-Sign[ N[ Log[ 2, x ], 5 ]-Floor[ N[ Log[ 2, x ], 5 ] ] ]; ld[ x_ ] := Length[ Divisors[ x ] ]; Flatten[ Position[ Table[ bi[ ld[ x ] ], {x, 1, m} ], 1 ] ]
Select[Range[110],IntegerQ[Log[2,DivisorSigma[0,#]]]&] (* Harvey P. Dale, Nov 20 2016 *)

is(n)=n=numdiv(n);n>>valuation(n,2)==1 \\ Charles R Greathouse IV, Mar 27 2013

isok(m) = issquarefree(m) || (omega(m) == omega(core(m)) && isok(core(m,1)[2])); \\ Peter Munn, Jun 18 2022

from itertools import count, islice
from sympy import factorint
def A036537_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:all(map(lambda m:not((k:=m+1)&-k)^k,factorint(n).values())),count(max(startvalue,1)))
A036537_list = list(islice(A036537_gen(),30)) # Chai Wah Wu, Jan 04 2023

A209229(A000005(a(n))) = 1. - Reinhard Zumkeller, Nov 15 2012
a(n) << n. - Charles R Greathouse IV, Feb 25 2017
m is in the sequence iff for k >= 0, A352780(m, k+1) | A352780(m, k)^2. - Peter Munn, Jun 18 2022

A019554 Smallest number whose square is divisible by n.

Original entry on oeis.org

1, 2, 3, 2, 5, 6, 7, 4, 3, 10, 11, 6, 13, 14, 15, 4, 17, 6, 19, 10, 21, 22, 23, 12, 5, 26, 9, 14, 29, 30, 31, 8, 33, 34, 35, 6, 37, 38, 39, 20, 41, 42, 43, 22, 15, 46, 47, 12, 7, 10, 51, 26, 53, 18, 55, 28, 57, 58, 59, 30, 61, 62, 21, 8, 65, 66, 67, 34, 69, 70, 71, 12, 73, 74, 15, 38, 77

Offset: 1

R. Muller

A note on square roots of numbers: we can write sqrt(n) = b*sqrt(c) where c is squarefree. Then b = A000188(n) is the "inner square root" of n, c = A007913(n), and b*c = A019554(n) = "outer square root" of n.
Instead of the terms "inner square root" and "outer square root", we may use the terms "lower square root" and "upper square root", respectively. Upper k-th roots have been studied by Broughan (2002, 2003, 2006). - Petros Hadjicostas, Sep 15 2019
The number of times each number k appears in this sequence is A034444(k). The first time k appears is at position A102631(k). - N. J. A. Sloane, Jul 28 2021

T. D. Noe, Table of n, a(n) for n = 1..10000
Henry Bottomley, Some Smarandache-type multiplicative sequences.
Kevin A. Broughan, Restricted divisor sums, Acta Arithmetica, 101(2) (2002), 105-114. See also here for another copy.
Kevin A. Broughan, Relationship between the integer conductor and k-th root functions, Int. J. Pure Appl. Math. 5(3) (2003), 253-275.
Kevin A. Broughan, Relaxations of the ABC Conjecture using integer k'th roots, New Zealand J. Math. 35(2) (2006), 121-136.
Florentin Smarandache, Collected Papers, Vol. II, Tempus Publ. Hse, Bucharest, 1996.
Eric Weisstein's World of Mathematics, Smarandache Ceil Function.

Cf. A000188 (inner square root), A053150 (inner 3rd root), A019555 (outer 3rd root), A053164 (inner 4th root), A053166 (outer 4th root), A015052 (outer 5th root), A015053 (outer 6th root).

Cf. also A007913, A007947, A008833, A015049, A034444, A053143, A102631, A346602, A005117, A133466, A195085.

a019554 n = product $ zipWith (^)
            (a027748_row n) (map ((`div` 2) . (+ 1)) $ a124010_row n)
-- Reinhard Zumkeller, Apr 13 2013
(Python 3.8+)
from math import prod
from sympy import factorint
def A019554(n): return n//prod(p**(q//2) for p, q in factorint(n).items()) # Chai Wah Wu, Aug 18 2021

with(numtheory):A019554 := proc(n) local i: RETURN(op(mul(i,i=map(x->x[1]^ceil(x[2]/2),ifactors(n)[2])))); end;

Flatten[Table[Select[Range[n],Divisible[#^2,n]&,1],{n,100}]] (* Harvey P. Dale, Oct 17 2011 *)
f[p_, e_] := p^Ceiling[e/2]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 18 2020 *)

a(n)=n/core(n,1)[2] \\ Charles R Greathouse IV, Feb 24 2011

Replace any square factors in n by their square roots.
Multiplicative with a(p^e) = p^ceiling(e/2).
Dirichlet series:
   Sum_{n>=1} a(n)/n^s = zeta(2*s-1)*zeta(s-1)/zeta(2*s-2), (Re(s) > 2);
   Sum_{n>=1} (1/a(n))/n^s = zeta(2*s+1)*zeta(s+1)/zeta(2*s+2), (Re(s) > 0).
a(n) = n/A000188(n).
a(n) = denominator of n/n^(3/2). - Arkadiusz Wesolowski, Dec 04 2011
a(n) = Product_{k=1..A001221(n)} A027748(n,k)^ceiling(A124010(n,k)/2). - Reinhard Zumkeller, Apr 13 2013
Sum_{k=1..n} a(k) ~ 3*zeta(3)*n^2 / Pi^2. - Vaclav Kotesovec, Sep 18 2020
Sum_{k=1..n} 1/a(k) ~ 3*log(n)^2/(2*Pi^2) + (9*gamma/Pi^2 - 36*zeta'(2)/Pi^4)*log(n) + 6*gamma^2/Pi^2 - 108*gamma*zeta'(2)/Pi^4 + 432*zeta'(2)^2/Pi^6 - 36*zeta''(2)/Pi^4 - 15*sg1/Pi^2, where gamma is the Euler-Mascheroni constant A001620 and sg1 is the first Stieltjes constant (see A082633). - Vaclav Kotesovec, Jul 27 2021
a(n) = sqrt(n*A007913(n)). - Jianing Song, May 08 2022
a(n) = sqrt(A053143(n)). - Amiram Eldar, Sep 02 2023
From Mia Boudreau, Jul 17 2025: (Start)
a(n^2) = n.
a(A005117(n)) = A005117(n).
a(A133466(n)) = A133466(n)/2.
a(A195085(n)) = A195085(n)/3. (End)

A055229 Greatest common divisor of largest square dividing n and squarefree part of n.

Original entry on oeis.org

1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 3, 1, 2, 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, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1

Offset: 1

Labos Elemer, Jun 21 2000

Record values occur at cubes of squarefree numbers: a(A062838(n)) = A005117(n) and a(m) < A005117(n) for m < A062838(n). - Reinhard Zumkeller, Apr 09 2010

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000

Cf. A000188, A007913, A007947, A008833.

Cf. A027748, A005117, A062838, A060476, A124010, A220218.

a055229 n = product $ zipWith (^) ps (map (flip mod 2) es) where
   (ps, es) = unzip $
              filter ((> 1) . snd) $ zip (a027748_row n) (a124010_row n)
-- Reinhard Zumkeller, Oct 27 2015

a[n_] := With[{sf = Times @@ Power @@@ ({#[[1]], Mod[#[[2]], 2]}& /@ FactorInteger[n])}, GCD[sf, n/sf]]; Table[a[n], {n, 1, 105}] (* Jean-François Alcover, Feb 05 2014 *)

a(n)=my(c=core(n));gcd(c,n/c) \\ Charles R Greathouse IV, Nov 20 2012

a(n) = gcd[A008833(n), A007913(n)].
Multiplicative with a(p^e)=1 for even e, a(p)=1, a(p^e)=p for odd e>1. - Vladeta Jovovic, Apr 30 2002
A220218(a(n)) = 1; A060476(a(n)) > 1 for n > 1. - Reinhard Zumkeller, Nov 30 2015
a(n) = core(n)*rad(n/core(n))/rad(n), where core = A007913 and rad = A007947. - Conjecture by Velin Yanev, proof by David J. Seal, Sep 19 2017
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = Product_{p prime} ((p^3 + p^2 + p - 1)/(p^2 * (p + 1))) = 1.2249749939341923764... . - Amiram Eldar, Oct 08 2022

A319057 Minimum sum of a strict factorization of n into factors > 1.

Original entry on oeis.org

0, 2, 3, 4, 5, 5, 7, 6, 9, 7, 11, 7, 13, 9, 8, 10, 17, 9, 19, 9, 10, 13, 23, 9, 25, 15, 12, 11, 29, 10, 31, 12, 14, 19, 12, 11, 37, 21, 16, 11, 41, 12, 43, 15, 14, 25, 47, 12, 49, 15, 20, 17, 53, 14, 16, 13, 22, 31, 59, 12, 61, 33, 16, 14, 18, 16, 67, 21, 26

Offset: 1

Gus Wiseman, Sep 09 2018

nonn

a(n) >= A001414(n), with equality iff n is squarefree or four times a squarefree number (i.e., A000188(n) <= 2). - Charlie Neder, Sep 10 2018

			The strict factorizations of 48 are (48), (2*24), (3*16), (4*12), (6*8), (2*3*8), (2*4*6), with sums 48, 26, 19, 16, 14, 13, 12 respectively, so a(48) = 12.

Charlie Neder, Table of n, a(n) for n = 1..1000

Cf. A001055, A007716, A045778, A056239, A162247, A215366, A246868, A318871, A318953, A318954.

strfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[strfacs[n/d],Min@@#>d&],{d,Rest[Divisors[n]]}]];
Table[Min[Total/@strfacs[n]],{n,100}]

A059896 The set of Fermi-Dirac factors of A(n,k) is the union of the Fermi-Dirac factors of n and k. Symmetric square array read by antidiagonals.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 4, 6, 6, 4, 5, 8, 3, 8, 5, 6, 10, 12, 12, 10, 6, 7, 6, 15, 4, 15, 6, 7, 8, 14, 6, 20, 20, 6, 14, 8, 9, 8, 21, 24, 5, 24, 21, 8, 9, 10, 18, 24, 28, 30, 30, 28, 24, 18, 10, 11, 10, 27, 8, 35, 6, 35, 8, 27, 10, 11, 12, 22, 30, 36, 40, 42, 42, 40, 36, 30, 22, 12, 13, 24

Offset: 1

Marc LeBrun, Feb 06 2001

Every positive integer, m, is the product of a unique subset, S(m), of the numbers listed in A050376 (primes, squares of primes etc.) The Fermi-Dirac factors of m are the members of S(m). So T(n,k) is the product of the members of (S(n) U S(k)).
Old name: Table a(i,j) = product prime(k)^(Ei(k) OR Ej(k)) where Ei and Ej are the vectors of exponents in the prime factorizations of i and j; OR is the bitwise operation on binary representation of the exponents.
Analogous to LCM, with OR replacing MAX.
A003418-analog seems to be A066616. - Antti Karttunen, Apr 12 2017
Considered as a binary operation, the result is the lowest common multiple of the squarefree parts of its operands multiplied by the square of the operation's result when applied to the square roots of the square parts of its operands. - Peter Munn, Mar 02 2022

			A(864,1944) = A(2^5*3^3,2^3*3^5) = 2^(5 OR 3) * 3^(3 OR 5) = 2^7*3^7 = 279936.
The top left 12 X 12 corner of the array:
   1,  2,  3,  4,  5,  6,  7,  8,   9,  10,  11,  12
   2,  2,  6,  8, 10,  6, 14,  8,  18,  10,  22,  24
   3,  6,  3, 12, 15,  6, 21, 24,  27,  30,  33,  12
   4,  8, 12,  4, 20, 24, 28,  8,  36,  40,  44,  12
   5, 10, 15, 20,  5, 30, 35, 40,  45,  10,  55,  60
   6,  6,  6, 24, 30,  6, 42, 24,  54,  30,  66,  24
   7, 14, 21, 28, 35, 42,  7, 56,  63,  70,  77,  84
   8,  8, 24,  8, 40, 24, 56,  8,  72,  40,  88,  24
   9, 18, 27, 36, 45, 54, 63, 72,   9,  90,  99, 108
  10, 10, 30, 40, 10, 30, 70, 40,  90,  10, 110, 120
  11, 22, 33, 44, 55, 66, 77, 88,  99, 110,  11, 132
  12, 24, 12, 12, 60, 24, 84, 24, 108, 120, 132,  12

Antti Karttunen, Table of n, a(n) for n = 1..10440; the first 144 antidiagonals of the array
Eric Weisstein's World of Mathematics, Square Part.
Eric Weisstein's World of Mathematics, Squarefree Part.

Cf. A003418, A003990, A007947, A028233, A028234, A066616, A284576.
Sequences used in a definition of this sequence: A003986, A000188/A007913/A008833, A052330/A052331.
Has simple/very significant relationships with A003961, A059895/A059897, A225546, A267116.

a[i_, i_] := i;
a[i_, j_] := Module[{f1 = FactorInteger[i], f2 = FactorInteger[j], e1, e2}, e1[] = 0; Scan[(e1[#[[1]]] = #[[2]])&, f1]; e2[] = 0; Scan[(e2[#[[1]]] = #[[2]])&, f2]; Times @@ (#^BitOr[e1[#], e2[#]]& /@ Union[f1[[All, 1]], f2[[All, 1]]])];
Table[a[i - j + 1, j], {i, 1, 15}, {j, 1, i}] // Flatten (* Jean-François Alcover, Jun 19 2018 *)

A059896(n,k) = if(n==k, n, lcm(core(n),core(k)) * A059896(core(n,1)[2],core(k,1)[2])^2) \\ Peter Munn, Mar 07 2022

(define (A059896 n) (A059896bi (A002260 n) (A004736 n)))
(define (A059896bi a b) (let loop ((a a) (b b) (m 1)) (cond ((= 1 a) (* m b)) ((= 1 b) (* m a)) ((equal? (A020639 a) (A020639 b)) (loop (A028234 a) (A028234 b) (* m (expt (A020639 a) (A003986bi (A067029 a) (A067029 b)))))) ((< (A020639 a) (A020639 b)) (loop (/ a (A028233 a)) b (* m (A028233 a)))) (else (loop a (/ b (A028233 b)) (* m (A028233 b)))))))
;; Antti Karttunen, Apr 11 2017

From Antti Karttunen, Apr 11 2017: (Start)
A(x,y) = A059895(x,y) * A059897(x,y).
A(x,y) * A059895(x,y) = x*y.
(End).
From Peter Munn, Mar 02 2022: (Start)
OR denotes the bitwise operation (A003986).
Limited multiplicative property: if gcd(n_1*k_1, n_2*k_2) = 1 then A(n_1*n_2, k_1*k_2) = A(n_1, k_1) * A(n_2, k_2).
For prime p, A(p^e_1, p^e_2) = p^(e_1 OR e_2).
A(n, A(m, k)) = A(A(n, m), k).
A(n, k) = A(k, n).
A(n, 1) = A(n, n) = n.
A(n^2, k^2) = A(n, k)^2.
A(n, k) = A(A007913(n), A007913(k)) * A(A008833(n), A008833(k)) = lcm(A007913(n), A007913(k)) * A(A000188(n), A000188(k))^2.
A007947(A(n, k)) = A007947(n*k).
Isomorphism: A(A052330(n), A052330(k)) = A052330(n OR k).
Equivalently, A(n, k) = A052330(A052331(n) OR A052331(k)).
A(A003961(n), A003961(k)) = A003961(A(n, k)).
A(A225546(n), A225546(k)) = A225546(A(n, k)).
(End)

New name from Peter Munn, Mar 02 2022

A323308 The number of exponential semiproper divisors of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 4, 1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1

Offset: 1

Amiram Eldar, Jan 10 2019

An exponential semiproper divisor of n is a divisor d such that rad(d) = rad(n) and gcd(d/rad(n), n/d) = 1, where rad(n) is the largest squarefree divisor of n (A007947).
a(n) is also the number of divisors of n that are squares of squarefree numbers (A062503). - Amiram Eldar, Oct 08 2022
a(n) is also the number of unitary divisors of n that are powerful (A001694). - Amiram Eldar, Feb 18 2023
The smallest integer that has exactly 2^n exponential semiproper divisors is A061742(n). - Bernard Schott, Feb 20 2023

Amiram Eldar, Table of n, a(n) for n = 1..10000
Nicusor Minculete, A new class of divisors: the exponential semiproper divisors, Bulletin of the Transilvania University of Brasov, Mathematics, Informatics, Physics, Series III, Vol. 7, No. 1 (2014), pp. 37-46.
Index entries for sequences computed from exponents in factorization of n.
Index entries for sequences related to divisors of numbers.

Cf. A000188, A001694, A007947, A008833, A034444, A057521, A061742, A062503, A082020, A323309, A323310.

f[p_, e_] := If[e==1, 1, 2]; a[1]=1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100]

a(n) = my(f=factor(n)); for (k=1, #f~, f[k,1] = min(f[k,2], 2); f[k,2] = 1); factorback(f); \\ Michel Marcus, Jan 11 2019

a(n) = A034444(n/A007947(n)).
Multiplicative with a(p^e) = 1 for e = 1 and 2 otherwise.
Asymptotic mean: Limit_{n->oo} (1/n) * Sum_{k=1..n} a(k) = 15/Pi^2 = 1.5198177546... (A082020). - Amiram Eldar, Nov 08 2020
a(n) = Sum_{d^2|n} mu(d)^2. - Wesley Ivan Hurt, Feb 13 2022
Dirichlet g.f.: zeta(s) * zeta(2*s) / zeta(4*s). - Werner Schulte, Dec 29 2022
a(n) = A034444(A000188(n)) = A034444(A008833(n)) (the number of unitary divisors of the largest square dividing n). - Amiram Eldar, Sep 03 2023
a(n) = A034444(A057521(n)) (the number of unitary divisors of the powerful part of n). - Amiram Eldar, Oct 03 2023

A056624 Number of unitary square divisors of n.

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 2, 2, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 1, 2, 2, 4, 1, 1, 1, 1, 1

Offset: 1

Labos Elemer, Aug 08 2000

Unitary analog of A046951.

The number of exponential divisors (A322791) of n that are cubefree (A004709). - Amiram Eldar, Jun 03 2025

			n=256, it has 5 square divisors of which only 2,{1,256} are unitary, 3 divisors are not.
n=124 has 2 (1 and 4) square divisors, both of them unitary a(124) = 2.
n=108 has 12 divisors, 4 square divisors: {1,4,9,36} of which 1 and 4 are unitary, 9 and 36 are not. So a(108)=2. The largest unitary square divisor of 108 is 4 with 1 prime divisor so a(108) = 2^1 = 2.

Antti Karttunen, Table of n, a(n) for n = 1..10000
Index entries for sequences computed from exponents in factorization of n.

Cf. A000188, A004709, A008833, A034444, A046952, A055229, A056626, A059841, A162641, A322791, A350388.

isA056624 := (n, d) -> igcd(n, d) = d and igcd(n/d, d) = d and igcd(n/d^2, d) = 1:
a := n -> nops(select(k -> isA056624(n, k), [seq(1..n)])):  # Peter Luschny, Jun 13 2025

Table[DivisorSum[n, 1 &, And[IntegerQ@ Sqrt@ #, CoprimeQ[#, n/#]] &], {n, 105}] (* Michael De Vlieger, Jul 28 2017 *)
f[p_, e_] := 2^(1 - Mod[e, 2]); a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Jan 03 2022 *)

a(n) = sumdiv(n, d, if(gcd(d, n/d)==1, issquare(d))); \\ Michel Marcus, Jul 28 2017

from sympy import factorint
def A056624(n): return 1<Chai Wah Wu, Aug 03 2024

def is_A056624(n, d): return gcd(n, d) == d and gcd(n//d, d) == d and gcd(n//(d*d), d) == 1
def a(n): return len([k for k in range(1, n+1) if is_A056624(n, k)])
print([a(n) for n in range(1, 106)])  # Peter Luschny, Jun 13 2025

(define (A056624 n) (if (= 1 n) n (* (A000079 (A059841 (A067029 n))) (A056624 (A028234 n))))) ;; Antti Karttunen, Jul 28 2017

a(n) = 2^r, where r is the number of prime factors of the largest unitary square divisor of n.
Multiplicative with a(p^e) = 2^(1-(e mod 2)). - Vladeta Jovovic, Dec 13 2002
Dirichlet g.f.: zeta(s)*zeta(2*s)/zeta(3*s). - Werner Schulte, Apr 03 2018
Sum_{k=1..n} a(k) ~ n*Pi^2/(6*zeta(3)) + sqrt(n)*zeta(1/2)/zeta(3/2). - Vaclav Kotesovec, Feb 07 2019
a(n) = 2^A162641(n). - Amiram Eldar, Sep 26 2022
a(n) = A034444(A350388(n)). - Amiram Eldar, Sep 09 2023

More terms from Vladeta Jovovic, Dec 13 2002

A293442 Multiplicative with a(p^e) = A019565(e).

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 6, 3, 4, 2, 6, 2, 4, 4, 5, 2, 6, 2, 6, 4, 4, 2, 12, 3, 4, 6, 6, 2, 8, 2, 10, 4, 4, 4, 9, 2, 4, 4, 12, 2, 8, 2, 6, 6, 4, 2, 10, 3, 6, 4, 6, 2, 12, 4, 12, 4, 4, 2, 12, 2, 4, 6, 15, 4, 8, 2, 6, 4, 8, 2, 18, 2, 4, 6, 6, 4, 8, 2, 10, 5, 4, 2, 12, 4, 4, 4, 12, 2, 12, 4, 6, 4, 4, 4, 20, 2, 6, 6, 9, 2, 8, 2, 12, 8

Offset: 1

Antti Karttunen, Oct 31 2017

From Peter Munn, Apr 06 2021: (Start)
a(n) is determined by the prime signature of n.
Compare with the multiplicative, self-inverse A225546, which also maps 2^e to the squarefree number A019565(e). However, this sequence maps p^e to the same squarefree number for every prime p, whereas A225546 maps the e-th power of progressively larger primes to progressively greater powers of A019565(e).
Both sequences map powers of squarefree numbers to powers of squarefree numbers.
(End)

Sequences used in a definition of this sequence: A000188, A003961, A019565, A028234, A059895, A067029, A162642.
Sequences with related definitions: A225546, A293443, A293444.
Cf. also A293214.
Sequences used to express relationship between terms of this sequence: A006519, A007913, A008833, A064989, A334747.
Sequences related via this sequence: (A001222, A048675, A064547), (A007814, A162642), (A087207, A267116), (A248663, A268387).

f[n_] := If[n == 1, 1, Apply[Times, Prime@ Flatten@ Position[Reverse@ IntegerDigits[Last@ #, 2], 1]] * f[n/Apply[Power, #]] &@ FactorInteger[n][[1]]]; Array[f, 105] (* Michael De Vlieger, Oct 31 2017 *)

a(1) = 1; for n > 1, a(n) = A019565(A067029(n)) * a(A028234(n)).
Other identities. For all n >= 1:
a(a(n)) = A293444(n).
A048675(a(n)) = A001222(n).
A001222(a(n)) = A064547(n) = A048675(A293444(n)).
A007814(a(n)) = A162642(n).
A087207(a(n)) = A267116(n).
A248663(a(n)) = A268387(n).
From Peter Munn, Mar 14 2021: (Start)
Alternative definition: a(1) = 1; a(2) = 2; a(n^2) = A003961(a(n)); a(A003961(n)) = a(n); if A059895(n, k) = 1, a(n*k) = a(n) * a(k).
For n >= 3, a(n) < n.
a(2n) = A334747(a(A006519(n))) * a(n/A006519(n)), where A006519(n) is the largest power of 2 dividing n.
a(2n+1) = a(A064989(2n+1)).
a(n) = a(A007913(n)) * a(A008833(n)) = 2^A162642(n) * A003961(a(A000188(n))).
(End)

cp's OEIS Frontend

A056169 Number of unitary prime divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Scheme

Formula

A055231 Powerfree part of n: product of primes that divide n only once.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Scheme

Formula

Extensions

A036537 Numbers whose number of divisors is a power of 2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

PARI

Python

Formula

A019554 Smallest number whose square is divisible by n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

A055229 Greatest common divisor of largest square dividing n and squarefree part of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Mathematica

PARI

Formula

A319057 Minimum sum of a strict factorization of n into factors > 1.

Views

Author

Keywords

Comments

Examples

Links