A206778 -id:A206778 - cp's OEIS frontend

A005117 Squarefree numbers: numbers that are not divisible by a square greater than 1.

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

Offset: 1

N. J. A. Sloane

1 together with the numbers that are products of distinct primes.
Also smallest sequence with the property that a(m)*a(k) is never a square for k != m. - Ulrich Schimke (ulrschimke(AT)aol.com), Dec 12 2001
Numbers k such that there is only one Abelian group with k elements, the cyclic group of order k (the numbers such that A000688(k) = 1). - Ahmed Fares (ahmedfares(AT)my-deja.com), Apr 25 2001
Numbers k such that A007913(k) > phi(k). - Benoit Cloitre, Apr 10 2002
a(n) is the smallest m with exactly n squarefree numbers <= m. - Amarnath Murthy, May 21 2002
k is squarefree <=> k divides prime(k)# where prime(k)# = product of first k prime numbers. - Mohammed Bouayoun (bouyao(AT)wanadoo.fr), Mar 30 2004
Numbers k such that omega(k) = Omega(k) = A072047(k). - Lekraj Beedassy, Jul 11 2006
The LCM of any finite subset is in this sequence. - Lekraj Beedassy, Jul 11 2006
This sequence and the Beatty Pi^2/6 sequence (A059535) are "incestuous": the first 20000 terms are bounded within (-9, 14). - Ed Pegg Jr, Jul 22 2008
Let us introduce a function D(n) = sigma_0(n)/2^(alpha(1) + ... + alpha(r)), sigma_0(n) number of divisors of n (A000005), prime factorization of n = p(1)^alpha(1) * ... * p(r)^alpha(r), alpha(1) + ... + alpha(r) is sequence (A001222). Function D(n) splits the set of positive integers into subsets, according to the value of D(n). Squarefree numbers (A005117) has D(n)=1, other numbers are "deviated" from the squarefree ideal and have 0 < D(n) < 1. For D(n)=1/2 we have A048109, for D(n)=3/4 we have A060687. - Ctibor O. Zizka, Sep 21 2008
Numbers k such that gcd(k,k')=1 where k' is the arithmetic derivative (A003415) of k. - Giorgio Balzarotti, Apr 23 2011
Numbers k such that A007913(k) = core(k) = k. - Franz Vrabec, Aug 27 2011
Numbers k such that sqrt(k) cannot be simplified. - Sean Loughran, Sep 04 2011
Indices m where A057918(m)=0, i.e., positive integers m for which there are no integers k in {1,2,...,m-1} such that k*m is a square. - John W. Layman, Sep 08 2011
It appears that these are numbers j such that Product_{k=1..j} (prime(k) mod j) = 0 (see Maple code). - Gary Detlefs, Dec 07 2011. - This is the same claim as Mohammed Bouayoun's Mar 30 2004 comment above. To see why it holds: Primorial numbers, A002110, a subsequence of this sequence, are never divisible by any nonsquarefree number, A013929, and on the other hand, the index of the greatest prime dividing any n is less than n. Cf. A243291. - Antti Karttunen, Jun 03 2014
Conjecture: For each n=2,3,... there are infinitely many integers b > a(n) such that Sum_{k=1..n} a(k)*b^(k-1) is prime, and the smallest such an integer b does not exceed (n+3)*(n+4). - Zhi-Wei Sun, Mar 26 2013
The probability that a random natural number belongs to the sequence is 6/Pi^2, A059956 (see Cesàro reference). - Giorgio Balzarotti, Nov 21 2013
Booker, Hiary, & Keating give a subexponential algorithm for testing membership in this sequence without factoring. - Charles R Greathouse IV, Jan 29 2014
Because in the factorizations into prime numbers these a(n) (n >= 2) have exponents which are either 0 or 1 one could call the a(n) 'numbers with a fermionic prime number decomposition'. The levels are the prime numbers prime(j), j >= 1, and the occupation numbers (exponents) e(j) are 0 or 1 (like in Pauli's exclusion principle). A 'fermionic state' is then denoted by a sequence with entries 0 or 1, where, except for the zero sequence, trailing zeros are omitted. The zero sequence stands for a(1) = 1. For example a(5) = 6 = 2^1*3^1 is denoted by the 'fermionic state' [1, 1], a(7) = 10 by [1, 0, 1]. Compare with 'fermionic partitions' counted in A000009. - Wolfdieter Lang, May 14 2014
From Vladimir Shevelev, Nov 20 2014: (Start)
The following is an Eratosthenes-type sieve for squarefree numbers. For integers > 1:
1) Remove even numbers, except for 2; the minimal non-removed number is 3.
2) Replace multiples of 3 removed in step 1, and remove multiples of 3 except for 3 itself; the minimal non-removed number is 5.
3) Replace multiples of 5 removed as a result of steps 1 and 2, and remove multiples of 5 except for 5 itself; the minimal non-removed number is 6.
4) Replace multiples of 6 removed as a result of steps 1, 2 and 3 and remove multiples of 6 except for 6 itself; the minimal non-removed number is 7.
5) Repeat using the last minimal non-removed number to sieve from the recovered multiples of previous steps.
Proof. We use induction. Suppose that as a result of the algorithm, we have found all squarefree numbers less than n and no other numbers. If n is squarefree, then the number of its proper divisors d > 1 is even (it is 2^k - 2, where k is the number of its prime divisors), and, by the algorithm, it remains in the sequence. Otherwise, n is removed, since the number of its squarefree divisors > 1 is odd (it is 2^k-1).
(End)
The lexicographically least sequence of integers > 1 such that each entry has an even number of proper divisors occurring in the sequence (that's the sieve restated). - Glen Whitney, Aug 30 2015
0 is nonsquarefree because it is divisible by any square. - Jon Perry, Nov 22 2014, edited by M. F. Hasler, Aug 13 2015
The Heinz numbers of partitions with distinct parts. We define the Heinz number of a partition p = [p_1, p_2, ..., p_r] as Product_{j=1..r} prime(j) (concept used by Alois P. Heinz in A215366 as an "encoding" of a partition). For example, for the partition [1, 1, 2, 4, 10] the Heinz number is 2*2*3*7*29 = 2436. The number 30 (= 2*3*5) is in the sequence because it is the Heinz number of the partition [1,2,3]. - Emeric Deutsch, May 21 2015
It is possible for 2 consecutive terms to be even; for example a(258)=422 and a(259)=426. - Thomas Ordowski, Jul 21 2015. [These form a subsequence of A077395 since their product is divisible by 4. - M. F. Hasler, Aug 13 2015]
There are never more than 3 consecutive terms. Runs of 3 terms start at 1, 5, 13, 21, 29, 33, ... (A007675). - Ivan Neretin, Nov 07 2015
a(n) = product of row n in A265668. - Reinhard Zumkeller, Dec 13 2015
Numbers without excess, i.e., numbers k such that A001221(k) = A001222(k). - Juri-Stepan Gerasimov, Sep 05 2016
Numbers k such that b^(phi(k)+1) == b (mod k) for every integer b. - Thomas Ordowski, Oct 09 2016
Boreico shows that the set of square roots of the terms of this sequence is linearly independent over the rationals. - Jason Kimberley, Nov 25 2016 (reference found by Michael Coons).
Numbers k such that A008836(k) = A008683(k). - Enrique Pérez Herrero, Apr 04 2018
The prime zeta function P(s) "has singular points along the real axis for s=1/k where k runs through all positive integers without a square factor". See Wolfram link. - Maleval Francis, Jun 23 2018
Numbers k such that A007947(k) = k. - Kyle Wyonch, Jan 15 2021
The Schnirelmann density of the squarefree numbers is 53/88 (Rogers, 1964). - Amiram Eldar, Mar 12 2021
Comment from Isaac Saffold, Dec 21 2021: (Start)
Numbers k such that all groups of order k have a trivial Frattini subgroup [Dummit and Foote].
Let the group G have order n. If n is squarefree and n > 1, then G is solvable, and thus by Hall's Theorem contains a subgroup H_p of index p for all p | n. Each H_p is maximal in G by order considerations, and the intersection of all the H_p's is trivial. Thus G's Frattini subgroup Phi(G), being the intersection of G's maximal subgroups, must be trivial. If n is not squarefree, the cyclic group of order n has a nontrivial Frattini subgroup. (End)
Numbers for which the squarefree divisors (A206778) and the unitary divisors (A077610) are the same; moreover they are also the set of divisors (A027750). - Bernard Schott, Nov 04 2022
0 = A008683(a(n)) - A008836(a(n)) = A001615(a(n)) - A000203(a(n)). - Torlach Rush, Feb 08 2023
From Robert D. Rosales, May 20 2024: (Start)
Numbers n such that mu(n) != 0, where mu(n) is the Möbius function (A008683).
Solutions to the equation Sum_{d|n} mu(d)*sigma(d) = mu(n)*n, where sigma(n) is the sum of divisors function (A000203). (End)
a(n) is the smallest root of x = 1 + Sum_{k=1..n-1} floor(sqrt(x/a(k))) greater than a(n-1). - Yifan Xie, Jul 10 2024
Number k such that A001414(k) = A008472(k). - Torlach Rush, Apr 14 2025
To elaborate on the formula from Greathouse (2018), the maximum of a(n) - floor(n*Pi^2/6 + sqrt(n)/17) is 10 at indices n = 48715, 48716, 48721, and 48760. The maximum is 11, at the same indices, if floor is taken individually for the two addends and the square root. If the value is rounded instead, the maximum is 9 at 10 indices between 48714 and 48765. - M. F. Hasler, Aug 08 2025

Jean-Marie De Koninck, Ces nombres qui nous fascinent, Entry 165, p. 53, Ellipses, Paris, 2008.
David S. Dummit and Richard M. Foote, Abstract algebra. Vol. 1999. Englewood Cliffs, NJ: David S.Prentice Hall, 1991.
Ivan M. Niven and Herbert S. Zuckerman, An Introduction to the Theory of Numbers. 2nd ed., Wiley, NY, 1966, p. 251.
Michael Pohst and Hans J. Zassenhaus, Algorithmic Algebraic Number Theory, Cambridge Univ. Press, page 432.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Daniel Forgues, Table of n, a(n) for n = 1..60794 (first 10000 terms from T. D. Noe)
Zenon B. Batang, Squarefree integers and the abc conjecture, arXiv:2109.10226 [math.GM], 2021.
Andrew R. Booker, Ghaith A. Hiary and Jon P. Keating, Detecting squarefree numbers, Duke Mathematical Journal, Vol. 164, No. 2 (2015), pp. 235-275; arXiv preprint, arXiv:1304.6937 [math.NT], 2013-2015.
Iurie Boreico, Linear independence of radicals, The Harvard College Mathematics Review 2(1), 87-92, Spring 2008.
Ernesto Cesàro, La serie di Lambert in aritmetica assintotica, Rendiconto della Reale Accademia delle Scienze di Napoli, Serie 2, Vol. 7 (1893), pp. 197-204.
Henri Cohen, Francois Dress, and Mohamed El Marraki, Explicit estimates for summatory functions linked to the Möbius μ-function, Functiones et Approximatio Commentarii Mathematici 37 (2007), part 1, pp. 51-63.
H. Gent, Letter to N. J. A. Sloane, Nov 27 1975.
Andrew Granville, ABC means we can count squarefrees, International Mathematical Research Notices 19 (1998), 991-1009.
Pentti Haukkanen, Mika Mattila, Jorma K. Merikoski and Timo Tossavainen, Can the Arithmetic Derivative be Defined on a Non-Unique Factorization Domain?, Journal of Integer Sequences, Vol. 16 (2013), Article 13.1.2.
Aaron Krowne, squarefree number, PlanetMath.org.
Louis Marmet, First occurrences of squarefree gaps and an algorithm for their computation.
Louis Marmet, First occurrences of square-free gaps and an algorithm for their computation, arXiv preprint arXiv:1210.3829 [math.NT], 2012.
Srinivasa Ramanujan, Irregular numbers, J. Indian Math. Soc., Vol. 5 (1913), pp. 105-106.
Kenneth Rogers, The Schnirelmann density of the squarefree integers, Proceedings of the American Mathematical Society, Vol. 15, No. 4 (1964), pp. 515-516.
J. A. Scott, Square-freedom revisited, The Mathematical Gazette, Vol. 90, No. 517 (2006), pp. 112-113.
Vladimir Shevelev, Set of all densities of exponentially S-numbers, arXiv preprint arXiv:1511.03860 [math.NT], 2015.
O. Trifonov, On the Squarefree Problem II, Math. Balkanica, Vol. 3 (1989), Fasc. 3-4.
Eric Weisstein's World of Mathematics, Squarefree.
Eric Weisstein's World of Mathematics, Prime Zeta Function.
Wikipedia, Squarefree integer.
Index entries for "core" sequences.

Complement of A013929. Subsequence of A072774 and A209061.
Characteristic function: A008966 (mu(n)^2, where mu = A008683).
Subsequences: A000040, A002110, A235488.
Cf. A001414, A008472, A076259 (first differences), A173143 (partial sums), A000688, A003277, A013928, A020753, A020754, A020755, A030059, A030229, A033197, A034444, A039956, A048672, A053797, A057918, A059956, A071403, A072284, A120992, A133466, A136742, A136743, A160764, A243289, A243347, A243348, A243351, A215366, A046660, A265668, A265675.
Cf. A027750, A077610, A206778.
Subsequences: numbers j such that j*a(k) is squarefree where k > 1: A056911 (k = 2), A261034 (k = 3), A274546 (k = 5), A276378 (k = 6).

a005117 n = a005117_list !! (n-1)
a005117_list = filter ((== 1) . a008966) [1..]
-- Reinhard Zumkeller, Aug 15 2011, May 10 2011

[ n : n in [1..1000] | IsSquarefree(n) ];

with(numtheory); a := [ ]; for n from 1 to 200 do if issqrfree(n) then a := [ op(a), n ]; fi; od:
t:= n-> product(ithprime(k),k=1..n): for n from 1 to 113 do if(t(n) mod n = 0) then print(n) fi od; # Gary Detlefs, Dec 07 2011
A005117 := proc(n) option remember; if n = 1 then 1; else for a from procname(n-1)+1 do if numtheory[issqrfree](a) then return a; end if; end do: end if; end proc:  # R. J. Mathar, Jan 09 2013

Select[ Range[ 113], SquareFreeQ] (* Robert G. Wilson v, Jan 31 2005 *)
Select[Range[150], Max[Last /@ FactorInteger[ # ]] < 2 &] (* Joseph Biberstine (jrbibers(AT)indiana.edu), Dec 26 2006 *)
NextSquareFree[n_, k_: 1] := Block[{c = 0, sgn = Sign[k]}, sf = n + sgn; While[c < Abs[k], While[ ! SquareFreeQ@ sf, If[sgn < 0, sf--, sf++]]; If[ sgn < 0, sf--, sf++]; c++]; sf + If[ sgn < 0, 1, -1]]; NestList[ NextSquareFree, 1, 70] (* Robert G. Wilson v, Apr 18 2014 *)
Select[Range[250], MoebiusMu[#] != 0 &] (* Robert D. Rosales, May 20 2024 *)

bnd = 1000; L = vector(bnd); j = 1; for (i=1,bnd, if(issquarefree(i),L[j]=i; j=j+1)); L

{a(n)= local(m,c); if(n<=1,n==1, c=1; m=1; while( cMichael Somos, Apr 29 2005 */

list(n)=my(v=vectorsmall(n,i,1),u,j); forprime(p=2,sqrtint(n), forstep(i=p^2, n, p^2, v[i]=0)); u=vector(sum(i=1,n,v[i])); for(i=1,n,if(v[i],u[j++]=i)); u \\ Charles R Greathouse IV, Jun 08 2012

for(n=1, 113, if(core(n)==n, print1(n, ", "))); \\ Arkadiusz Wesolowski, Aug 02 2016

S(n) = my(s); forsquarefree(k=1,sqrtint(n),s+=n\k[1]^2*moebius(k)); s;
a(n) = my(min=1, max=231, k=0, sc=0); if(n >= 144, min=floor(zeta(2)*n - 5*sqrt(n)); max=ceil(zeta(2)*n + 5*sqrt(n))); while(min <= max, k=(min+max)\2; sc=S(k); if(abs(sc-n) <= sqrtint(n), break); if(sc > n, max=k-1, if(sc < n, min=k+1, break))); while(!issquarefree(k), k-=1); while(sc != n, my(j=1); if(sc > n, j = -1); k += j; sc += j; while(!issquarefree(k), k += j)); k; \\ Daniel Suteu, Jul 07 2022

first(n)=my(v=vector(n),i); forsquarefree(k=1,if(n<268293,(33*n+30)\20,(n*Pi^2/6+0.058377*sqrt(n))\1), if(i++>n, return(v)); v[i]=k[1]); v \\ Charles R Greathouse IV, Jan 10 2023

A5117=[1..3]; A005117(n)={if(n>#A5117, my(N=#A5117); A5117=Vec(A5117, max(n+999, N*5\4)); iferr(forsquarefree(k=A5117[N]+1, #A5117*Pi^2\6+sqrtint(#A5117)\17+11, A5117[N++]=k[1]),E,)); A5117[n]} \\ M. F. Hasler, Aug 08 2025

from sympy.ntheory.factor_ import core
def ok(n): return core(n, 2) == n
print(list(filter(ok, range(1, 114)))) # Michael S. Branicky, Jul 31 2021

from itertools import count, islice
from sympy import factorint
def A005117_gen(startvalue=1): # generator of terms >= startvalue
    return filter(lambda n:all(x == 1 for x in factorint(n).values()),count(max(startvalue,1)))
A005117_list = list(islice(A005117_gen(),20)) # Chai Wah Wu, May 09 2022

from math import isqrt
from sympy import mobius
def A005117(n):
    def f(x): return n+x-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, Jul 22 2024

Limit_{n->oo} a(n)/n = Pi^2/6 (see A013661). - Benoit Cloitre, May 23 2002
Equals A039956 UNION A056911. - R. J. Mathar, May 16 2008
A122840(a(n)) <= 1; A010888(a(n)) < 9. - Reinhard Zumkeller, Mar 30 2010
a(n) = A055229(A062838(n)) and a(n) > A055229(m) for m < A062838(n). - Reinhard Zumkeller, Apr 09 2010
A008477(a(n)) = 1. - Reinhard Zumkeller, Feb 17 2012
A055653(a(n)) = a(n); A055654(a(n)) = 0. - Reinhard Zumkeller, Mar 11 2012
A008966(a(n)) = 1. - Reinhard Zumkeller, May 26 2012
Sum_{n>=1} 1/a(n)^s = zeta(s)/zeta(2*s). - Enrique Pérez Herrero, Jul 07 2012
A056170(a(n)) = 0. - Reinhard Zumkeller, Dec 29 2012
A013928(a(n)+1) = n. - Antti Karttunen, Jun 03 2014
A046660(a(n)) = 0. - Reinhard Zumkeller, Nov 29 2015
Equals {1} UNION A000040 UNION A006881 UNION A007304 UNION A046386 UNION A046387 UNION A067885 UNION A123321 UNION A123322 UNION A115343 ... - R. J. Mathar, Nov 05 2016
|a(n) - n*Pi^2/6| < 0.058377*sqrt(n) for n >= 268293; this result can be derived from Cohen, Dress, & El Marraki, see links. - Charles R Greathouse IV, Jan 18 2018
From Amiram Eldar, Jul 07 2021: (Start)
Sum_{n>=1} (-1)^(a(n)+1)/a(n)^2 = 9/Pi^2.
Sum_{k=1..n} 1/a(k) ~ (6/Pi^2) * log(n).
Sum_{k=1..n} (-1)^(a(k)+1)/a(k) ~ (2/Pi^2) * log(n).
(all from Scott, 2006) (End)

A034444 a(n) is the number of unitary divisors of n (d such that d divides n, gcd(d, n/d) = 1).

Original entry on oeis.org

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

Offset: 1

N. J. A. Sloane

If n = Product p_i^a_i, d = Product p_i^c_i is a unitary divisor of n if each c_i is 0 or a_i.
Also the number of squarefree divisors of n. - Labos Elemer
Also number of divisors of the squarefree kernel of n: a(n) = A000005(A007947(n)). - Reinhard Zumkeller, Jul 19 2002
Also shadow transform of pronic numbers A002378.
For n >= 1 define an n X n (0,1) matrix A by A[i,j] = 1 if lcm(i,j) = n, A[i,j] = 0 if lcm(i,j) <> n for 1 <= i,j <= n. a(n) is the rank of A. - Yuval Dekel (dekelyuval(AT)hotmail.com), Aug 11 2003
a(n) is also the number of solutions to x^2 - x == 0 (mod n). - Yuval Dekel (dekelyuval(AT)hotmail.com), Sep 21 2003
a(n) is the number of squarefree divisors of n, but in general the set of unitary divisors of n is not the set of squarefree divisors (compare the rows of A077610 and A206778). - Jaroslav Krizek, May 04 2009
Row lengths of the triangles in A077610 and in A206778. - Reinhard Zumkeller, Feb 12 2012
a(n) is also the number of distinct residues of k^phi(n) (mod n), k=0..n-1. - Michel Lagneau, Nov 15 2012
a(n) is the number of irreducible fractions y/x that satisfy x*y=n (and gcd(x,y)=1), x and y positive integers. - Luc Rousseau, Jul 09 2017
a(n) is the number of (x,y) lattice points satisfying both x*y=n and (x,y) is visible from (0,0), x and y positive integers. - Luc Rousseau, Jul 10 2017
Conjecture: For any nonnegative integer k and positive integer n, the sum of the k-th powers of the unitary divisors of n is divisible by the sum of the k-th powers of the odd unitary divisors of n (note that this sequence lists the sum of the 0th powers of the unitary divisors of n). - Ivan N. Ianakiev, Feb 18 2018
a(n) is the number of one-digit numbers, k, when written in base n such that k and k^2 end in the same digit. - Matthew Scroggs, Jun 01 2018
Dirichlet convolution of A271102 and A000005. - Vaclav Kotesovec, Apr 08 2019
Conjecture: Let b(i; n), n > 0, be multiplicative sequences for some fixed integer i >= 0 with b(i; p^e) = (Sum_{k=1..i+1} A164652(i, k) * e^(k-1)) * (i+2) / (i!) for prime p and e > 0. Then we have Dirichlet generating functions: Sum_{n > 0} b(i; n) / n^s = (zeta(s))^(i+2) / zeta((i+2) * s). Examples for i=0 this sequence, for i=1 A226602, and for i=2 A286779. - Werner Schulte, Feb 17 2022
The smallest integer with 2^m unitary divisors, or equivalently, the smallest integer with 2^m squarefree divisors, is A002110(m). - Bernard Schott, Oct 04 2022

			a(12) = 4 because the four unitary divisors of 12 are 1, 3, 4, 12, and also because the four squarefree divisors of 12 are 1, 2, 3, 6.

R. K. Guy, Unsolved Problems in Number Theory, Sect. B3.

T. D. Noe, Table of n, a(n) for n = 1..10000
O. Bagdasar, On some functions involving the lcm and gcd of integer tuples, Scientific Publications of the State University of Novi Pazar, Appl. Maths. Inform. and Mech., Vol. 6, 2 (2014), 91--100.
Masum Billal, Number of Ways To Write as Product of Co-prime Numbers, arXiv:1909.07823 [math.GM], 2019.
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
Steven R. Finch, Mathematical Constants II, Encyclopedia of Mathematics and Its Applications, Cambridge University Press, Cambridge, 2018, p. 49-50.
Lorenz Halbeisen, A number-theoretic conjecture and its implication for set theory, Acta Math. Univ. Comenianae 74(2) (2005), 243-254.
Lorenz Halbeisen and Norbert Hungerbuehler, Number theoretic aspects of a combinatorial function, Notes on Number Theory and Discrete Mathematics 5 (1999), 138-150.
Jon Maiga, Upper bound of Fibonacci entry points, 2019.
OEIS Wiki, Shadow transform.
N. J. A. Sloane, Transforms.
Eric Weisstein's World of Mathematics, Unitary Divisor.
Eric Weisstein's World of Mathematics, Unitary Divisor Function.
Wikipedia, Unitary divisor.
Index entries for sequences computed from exponents in factorization of n

Cf. A077610, A048105, A000173, A013928, A000079, A001221, A002110, A034448, A047994, A061142, A277561.
Sum of the k-th powers of the squarefree divisors of n for k=0..10: this sequence (k=0), A048250 (k=1), A351265 (k=2), A351266 (k=3), A351267 (k=4), A351268 (k=5), A351269 (k=6), A351270 (k=7), A351271 (k=8), A351272 (k=9), A351273 (k=10).
Sequences of the form n^k * Product_ {p|n, p prime} (1 + 1/p^k) for k=0..10: this sequence (k=0), A001615 (k=1), A065958 (k=2), A065959 (k=3), A065960 (k=4), A351300 (k=5), A351301 (k=6), A351302 (k=7), A351303 (k=8), A351304 (k=9), this sequence (k=10).
Cf. A020821 (Dgf at s=2), A177057 (Dgf at s=4).

a034444 = length . a077610_row  -- Reinhard Zumkeller, Feb 12 2012

[#[d:d in Divisors(n)|Gcd(d,n div d) eq 1]:n in [1..110]]; // Marius A. Burtea, Jan 11 2020

[&+[Abs(MoebiusMu(d)):d in Divisors(n)]:n in [1..110]]; // Marius A. Burtea, Jan 11 2020

with(numtheory): for n from 1 to 200 do printf(`%d,`,2^nops(ifactors(n)[2])) od:
with(numtheory);
# returns the number of unitary divisors of n and a list of them
f:=proc(n)
local ct,i,t1,ans;
ct:=0; ans:=[];
t1:=divisors(n);
for i from 1 to nops(t1) do
d:=t1[i];
if igcd(d,n/d)=1 then ct:=ct+1; ans:=[op(ans),d]; fi;
od:
RETURN([ct,ans]);
end;
# N. J. A. Sloane, May 01 2013
# alternative Maple program:
a:= n-> 2^nops(ifactors(n)[2]):
seq(a(n), n=1..105);  # Alois P. Heinz, Jan 23 2024
a := n -> 2^NumberTheory:-NumberOfPrimeFactors(n, distinct):  # Peter Luschny, May 13 2025

a[n_] := Count[Divisors[n], d_ /; GCD[d, n/d] == 1]; a /@ Range[105] (* Jean-François Alcover, Apr 05 2011 *)
Table[2^PrimeNu[n],{n,110}] (* Harvey P. Dale, Jul 14 2011 *)

a(n)=1<Charles R Greathouse IV, Feb 11 2011

for(n=1, 100, print1(direuler(p=2, n, (1+X)/(1-X))[n], ", ")) \\ Vaclav Kotesovec, Sep 26 2020

from sympy import divisors, gcd
def a(n):
    return sum(1 for d in divisors(n) if gcd(d, n//d)==1)
# Indranil Ghosh, Apr 16 2017

from sympy import primefactors
def a(n): return 2**len(primefactors(n))
print([a(n) for n in range(1, 101)]) # Indranil Ghosh, Apr 16 2017

(define (A034444 n) (if (= 1 n) n (* 2 (A034444 (A028234 n))))) ;; Antti Karttunen, May 29 2017

a(n) = Sum_{d|n} abs(mu(n)) = 2^(number of different primes dividing n) = 2^A001221(n), with mu(n) = A008683(n). [Added Möbius formula. - Wolfdieter Lang, Jan 11 2020]
a(n) = Product_{ primes p|n } (1 + Legendre(1, p)).
Multiplicative with a(p^k)=2 for p prime and k>0. - Henry Bottomley, Oct 25 2001
a(n) = Sum_{d|n} tau(d^2)*mu(n/d), Dirichlet convolution of A048691 and A008683. - Benoit Cloitre, Oct 03 2002
Dirichlet generating function: zeta(s)^2/zeta(2s). - Franklin T. Adams-Watters, Sep 11 2005
Inverse Mobius transform of A008966. - Franklin T. Adams-Watters, Sep 11 2005
Asymptotically [Finch] the cumulative sum of a(n) = Sum_{n=1..N} a(n) ~ (6/(Pi^2))*N*log(N) + (6/(Pi^2))*(2*gamma - 1 - (12/(Pi^2))*zeta'(2))*N + O(sqrt(N)). - Jonathan Vos Post, May 08 2005 [typo corrected by Vaclav Kotesovec, Sep 13 2018]
a(n) = Sum_{d|n} floor(rad(d)/d), where rad is A007947 and floor(rad(n)/n) = A008966(n). - Enrique Pérez Herrero, Nov 13 2009
a(n) = A000005(n) - A048105(n); number of nonzero terms in row n of table A225817. - Reinhard Zumkeller, Jul 30 2013
G.f.: Sum_{n>0} A008966(n)*x^n/(1-x^n). - Mircea Merca, Feb 25 2014
a(n) = Sum_{d|n} lambda(d)*mu(d), where lambda is A008836. - Enrique Pérez Herrero, Apr 27 2014
a(n) = A277561(A156552(n)). - Antti Karttunen, May 29 2017
a(n) = A005361(n^2)/A005361(n). - Velin Yanev, Jul 26 2017
L.g.f.: -log(Product_{k>=1} (1 - mu(k)^2*x^k)^(1/k)) = Sum_{n>=1} a(n)*x^n/n. - Ilya Gutkovskiy, Jul 30 2018
a(n) = Sum_{d|n} A001615(d) * A023900(n/d). - Torlach Rush, Jan 20 2020
Sum_{d|n, gcd(d, n/d) = 1} a(d) * (-1)^omega(n/d) = 1. - Amiram Eldar, May 29 2020
a(n) = lim_{k->oo} A000005(n^(2*k))/A000005(n^k). - Velin Yanev and Amiram Eldar, Jan 10 2025

More terms from James Sellers, Jun 20 2000

A048250 Sum of the squarefree divisors of n.

Original entry on oeis.org

1, 3, 4, 3, 6, 12, 8, 3, 4, 18, 12, 12, 14, 24, 24, 3, 18, 12, 20, 18, 32, 36, 24, 12, 6, 42, 4, 24, 30, 72, 32, 3, 48, 54, 48, 12, 38, 60, 56, 18, 42, 96, 44, 36, 24, 72, 48, 12, 8, 18, 72, 42, 54, 12, 72, 24, 80, 90, 60, 72, 62, 96, 32, 3, 84, 144, 68, 54, 96, 144, 72, 12, 74

Offset: 1

Labos Elemer

Also sum of divisors of the squarefree kernel of n: a(n) = A000203(A007947(n)). - Reinhard Zumkeller, Jul 19 2002
The absolute values of the Dirichlet inverse of A001615. - R. J. Mathar, Dec 22 2010
Row sums of the triangle in A206778. - Reinhard Zumkeller, Feb 12 2012
Inverse Möbius transform of n * mu(n)^2 = |A055615(n)|. - Wesley Ivan Hurt, Jun 08 2023

			For n=1000, out of the 16 divisors, four are squarefree: {1,2,5,10}. Their sum is 18. Or, 1000 = 2^3*5^3 hence a(1000) = (2+1)*(5+1) = 18.

D. Suryanarayana, On the core of an integer, Indian J. Math. 14 (1972) 65-74.

T. D. Noe, Table of n, a(n) for n = 1..10000
Steven R. Finch, Unitarism and infinitarism.
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
Index entries for sequences related to sums of squares

Cf. A003557, A007913, A007947, A023900, A034448, A206787, A309192.
Sum of the k-th powers of the squarefree divisors of n for k=0..10: A034444 (k=0), this sequence (k=1), A351265 (k=2), A351266 (k=3), A351267 (k=4), A351268 (k=5), A351269 (k=6), A351270 (k=7), A351271 (k=8), A351272 (k=9), A351273 (k=10).
Cf. A240976 (tenth of Dgf at s=3).

a034448 = sum . a206778_row  -- Reinhard Zumkeller, Feb 12 2012

A048250 := proc(n) local ans, i:ans := 1: for i from 1 to nops(ifactors(n)[ 2 ]) do ans := ans*(1+ifactors(n)[ 2 ][ i ] [ 1 ]): od: RETURN(ans) end:
# alternative:
seq(mul(1+p, p = numtheory:-factorset(n)), n=1..1000); # Robert Israel, Mar 18 2015

sumOfSquareFreeDivisors[ n_ ] := Plus @@ Select[ Divisors[ n ], MoebiusMu[ # ] != 0 & ]; Table[ sumOfSquareFreeDivisors[ i ], {i, 85} ]
Table[Total[Select[Divisors[n],SquareFreeQ]],{n,80}] (* Harvey P. Dale, Jan 25 2013 *)
a[1] = 1; a[n_] := Times@@(1 + FactorInteger[n][[;;,1]]); Array[a, 100] (* Amiram Eldar, Dec 19 2018 *)

a(n)=if(n<1,0,sumdiv(n,d,if(core(d)==d,d)))

a(n)=if(n<1,0,direuler(p=2,n,(1+p*X)/(1-X))[n])

a(n)=sumdiv(n,d,moebius(d)^2*d); \\ Joerg Arndt, Jul 06 2011

a(n)=my(f=factor(n)); for(i=1,#f~,f[i,2]=1); sigma(f) \\ Charles R Greathouse IV, Sep 09 2014

from math import prod
from sympy import primefactors
def A048250(n): return prod(p+1 for p in primefactors(n)) # Chai Wah Wu, Apr 20 2023

def A048250(n): return mul(map(lambda p: p+1, prime_divisors(n)))
[A048250(n) for n in (1..73)]  # Peter Luschny, May 23 2013

If n = Product p_i^e_i, a(n) = Product (p_i + 1). - Vladeta Jovovic, Apr 19 2001
Dirichlet g.f.: zeta(s)*zeta(s-1)/zeta(2*s-2). - Michael Somos, Sep 08 2002
a(n) = Sum_{d|n} mu(d)^2*d = Sum_{d|n} |A055615(d)|. - Benoit Cloitre, Dec 09 2002
Pieter Moree (moree(AT)mpim-bonn.mpg.de), Feb 20 2004 can show that Sum_{n <= x} a(n) = x^2/2 + O(x*sqrt{x}) and adds: "As S. R. Finch pointed out to me, in Suryanarayana's paper this is proved under the Riemann hypothesis with error term O(x^{7/5+epsilon})".
a(n) = psi(rad(n)) = A001615(A007947(n)). - Enrique Pérez Herrero, Aug 24 2010
a(n) = rad(n)*psi(n)/n = A001615(n)*A007947(n)/n. - Enrique Pérez Herrero, Aug 31 2010
G.f.: Sum_{k>=1} mu(k)^2*k*x^k/(1 - x^k). - Ilya Gutkovskiy, Jan 03 2017
Lim_{n->oo} (1/n) * Sum_{k=1..n} a(k)/k = 1. - Amiram Eldar, Jun 10 2020
a(n) = Sum_{d divides n} mu(d)^2*core(d), where core(n) = A007913(n). - Peter Bala, Jan 24 2024

A077610 Triangle in which n-th row lists unitary divisors of n.

Original entry on oeis.org

1, 1, 2, 1, 3, 1, 4, 1, 5, 1, 2, 3, 6, 1, 7, 1, 8, 1, 9, 1, 2, 5, 10, 1, 11, 1, 3, 4, 12, 1, 13, 1, 2, 7, 14, 1, 3, 5, 15, 1, 16, 1, 17, 1, 2, 9, 18, 1, 19, 1, 4, 5, 20, 1, 3, 7, 21, 1, 2, 11, 22, 1, 23, 1, 3, 8, 24, 1, 25, 1, 2, 13, 26, 1, 27, 1, 4, 7, 28, 1, 29, 1, 2, 3, 5, 6, 10, 15, 30

Offset: 1

Eric W. Weisstein, Nov 11 2002

n-th row = n-th row of A165430 without repetitions. - Reinhard Zumkeller, Mar 04 2013
Denominators of sequence of all positive rational numbers ordered as follows:  let m = p(i(1))^e(i(1))*...*p(i(k))^e(i(k)) be the prime factorization of m.  Let S(m) be the vector of rationals p(i(k+1-j))^e(i(k+1-j))/p(i(j))^e(i(j)) for j = 1..k.  The sequence (a(n)) is the concatenation of vectors S(m) for m = 1, 2, ...; for numerators see A229994. - Clark Kimberling, Oct 31 2013
The concept of unitary divisors was introduced by the Indian mathematician Ramaswamy S. Vaidyanathaswamy (1894-1960) in 1931. He called them "block factors". The term "unitary divisor" was coined by Cohen (1960). - Amiram Eldar, Mar 09 2024

			1;
1, 2;
1, 3;
1, 4;
1, 5;
1, 2, 3, 6;
1, 7;
1, 8;
1, 9;
1, 2, 5, 10;
1, 11;

Reinhard Zumkeller, Rows n=1..1000 of triangle, flattened
Eckford Cohen, Arithmetical functions associated with the unitary divisors of an integer, Mathematische Zeitschrift, Vol. 74 (1960), pp. 66-80.
R. Vaidyanathaswamy, The theory of multiplicative arithmetic functions, Transactions of the American Mathematical Society, Vol. 33, No. 2 (1931), pp. 579-662.
Eric Weisstein's World of Mathematics, Unitary Divisor.
Index entries for sequences related to enumerating the rationals.

Cf. A037445, A027750, A034444 (row lengths), A034448 (row sums); A206778.

a077610 n k = a077610_row n !! k
a077610_row n = [d | d <- [1..n], let (n',m) = divMod n d,
                     m == 0, gcd d n' == 1]
a077610_tabf = map a077610_row [1..]
-- Reinhard Zumkeller, Feb 12 2012

with(numtheory);
# returns the number of unitary divisors of n and a list of them, from N. J. A. Sloane, May 01 2013
f:=proc(n)
local ct,i,t1,ans;
ct:=0; ans:=[];
t1:=divisors(n);
for i from 1 to nops(t1) do
d:=t1[i];
if igcd(d,n/d)=1 then ct:=ct+1; ans:=[op(ans),d]; fi;
od:
RETURN([ct,ans]);
end;
# Alternatively:
isUnitary := (n, d) -> igcd(n, d) = d and igcd(n/d, d) = 1:
aList := n -> select(k -> isUnitary(n, k), [seq(1..n)]): # Peter Luschny, Jun 13 2025

row[n_] := Select[ Divisors[n], GCD[#, n/#] == 1 &]; Table[row[n], {n, 1, 30}] // Flatten (* Jean-François Alcover, Oct 22 2012 *)

row(n)=my(f=factor(n),k=#f~); Set(vector(2^k,i, prod(j=1,k, if(bittest(i,j-1),1,f[j,1]^f[j,2]))))
v=[];for(n=1,20,v=concat(v,row(n)));v \\ Charles R Greathouse IV, Sep 02 2015

row(n) = {my(d = divisors(n)); select(x->(gcd(x, n/x)==1), d);} \\ Michel Marcus, Oct 11 2015

from math import gcd
def is_unitary(n, d) -> bool: return gcd(n, d) == d and gcd(n//d, d) == 1
def aList(n) -> list[int]: return [k for k in range(1, n+1) if is_unitary(n, k)]
for n in range(1, 31): print(aList(n))  # Peter Luschny, Jun 13 2025

d is unitary divisor of n <=> gcd(n, d) = d and gcd(n/d, d) = 1. - Peter Luschny, Jun 13 2025

A050320 Number of ways n is a product of squarefree numbers > 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 1, 2, 2, 1, 1, 2, 1, 2, 2, 2, 1, 2, 1, 2, 1, 2, 1, 5, 1, 1, 2, 2, 2, 3, 1, 2, 2, 2, 1, 5, 1, 2, 2, 2, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 2, 1, 6, 1, 2, 2, 1, 2, 5, 1, 2, 2, 5, 1, 3, 1, 2, 2, 2, 2, 5, 1, 2, 1, 2, 1, 6, 2, 2, 2, 2, 1, 6, 2, 2, 2, 2, 2, 2, 1, 2, 2, 3, 1, 5, 1

Offset: 1

Christian G. Bower, Sep 15 1999

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24 = 2^3*3 and 375 = 3*5^3 both have prime signature (3,1).
Broughan shows (Theorem 8) that the average value of a(n) is k exp(2*sqrt(log n)/sqrt(zeta(2)))/log(n)^(3/4) where k is about 0.18504.  - Charles R Greathouse IV, May 21 2013
From Gus Wiseman, Aug 20 2020: (Start)
Also the number of set multipartitions (multisets of sets) of the multiset of prime indices of n. For example, the a(n) set multipartitions for n = 2, 6, 36, 60, 360 are:
  {1}  {12}    {12}{12}      {1}{123}      {1}{12}{123}
       {1}{2}  {1}{2}{12}    {12}{13}      {12}{12}{13}
               {1}{1}{2}{2}  {1}{1}{23}    {1}{1}{12}{23}
                             {1}{2}{13}    {1}{1}{2}{123}
                             {1}{3}{12}    {1}{2}{12}{13}
                             {1}{1}{2}{3}  {1}{3}{12}{12}
                                           {1}{1}{1}{2}{23}
                                           {1}{1}{2}{2}{13}
                                           {1}{1}{2}{3}{12}
                                           {1}{1}{1}{2}{2}{3}
(End)

			For n = 36 we have three choices as 36 = 2*2*3*3 = 6*6 = 2*3*6 (but no factorizations with factors 4, 9, 12, 18 or 36 are allowed), thus a(36) = 3. - _Antti Karttunen_, Oct 21 2017

Giovanni Resta, Table of n, a(n) for n = 1..10000
Kevin Broughan, Quadrafree factorisatio numerorum, Rocky Mountain J. Math. 44 (3) (2014) 791-807.
Index entries for sequences computed from exponents in factorization of n

Cf. A001055, A005117, A050325. a(p^k)=1. a(A002110)=A000110.
a(n!)=A103774(n).
Cf. A206778.
Differs from A259936 for the first time at n=36.
A050326 is the strict case.
A045778 counts strict factorizations.
A089259 counts set multipartitions of integer partitions.
A116540 counts normal set multipartitions.
Cf. A008480, A124010, A317829, A318360.

a050320 n = h n $ tail $ a206778_row n where
   h 1 _          = 1
   h _ []         = 0
   h m fs'@(f:fs) =
     if f > m then 0 else if r > 0 then h m fs else h m' fs' + h m fs
     where (m', r) = divMod m f
-- Reinhard Zumkeller, Dec 16 2013

sub[w_, e_] := Block[{v = w}, v[[e]]--; v]; ric[w_, k_] := If[Max[w] == 0, 1, Block[{e, s, p = Flatten@Position[Sign@w, 1]}, s = Select[Prepend[#, First@p] & /@ Subsets[Rest@p], Total[1/2^#] <= k &]; Sum[ric[sub[w, e], Total[1/2^e]], {e, s}]]]; sig[w_] := sig[w] = ric[w, 1];  a[n_] := sig@ Sort[Last /@ FactorInteger[n]]; Array[a, 103] (* Giovanni Resta, May 21 2013 *)
sqfacs[n_]:=If[n<=1,{{}},Join@@Table[(Prepend[#,d]&)/@Select[sqfacs[n/d],Min@@#>=d&],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]]
Table[Length[sqfacs[n]],{n,100}] (* Gus Wiseman, Aug 20 2020 *)

Dirichlet g.f.: Product_{n is squarefree and > 1} (1/(1-1/n^s)).
a(n) = A050325(A101296(n)). - R. J. Mathar, May 26 2017
a(n!) = A103774(n); a(A006939(n)) = A337072(n). - Gus Wiseman, Aug 20 2020

A050326 Number of factorizations of n into distinct squarefree numbers > 1.

Original entry on oeis.org

1, 1, 1, 0, 1, 2, 1, 0, 0, 2, 1, 1, 1, 2, 2, 0, 1, 1, 1, 1, 2, 2, 1, 0, 0, 2, 0, 1, 1, 5, 1, 0, 2, 2, 2, 1, 1, 2, 2, 0, 1, 5, 1, 1, 1, 2, 1, 0, 0, 1, 2, 1, 1, 0, 2, 0, 2, 2, 1, 4, 1, 2, 1, 0, 2, 5, 1, 1, 2, 5, 1, 0, 1, 2, 1, 1, 2, 5, 1, 0, 0, 2, 1, 4, 2, 2, 2, 0, 1, 4, 2, 1, 2, 2, 2, 0, 1, 1, 1, 1, 1, 5, 1

Offset: 1

Christian G. Bower, Oct 15 1999

nonn

a(n) depends only on prime signature of n (cf. A025487). So a(24) = a(375) since 24=2^3*3 and 375=3*5^3 both have prime signature (3,1).
a(A212164(n)) = 0; a(A212166(n)) = 1; a(A006881(n)) = 2; a(A190107(n)) = 3; a(A085987(n)) = 4; a(A225228(n)) = 5; a(A179670(n)) = 7; a(A162143(n)) = 8; a(A190108(n)) = 11; a(A212167(n)) > 0; a(A212168(n)) > 1. - Reinhard Zumkeller, May 03 2013
The comment that a(A212164(n)) = 0 is incorrect. For example, 3600 belongs to A212164 but a(3600) = 1. The positions of zeros in this sequence are A293243. - Gus Wiseman, Oct 10 2017

			The a(30) = 5 factorizations are: 2*3*5, 2*15, 3*10, 5*6, 30. The a(180) = 5 factorizations are: 2*3*5*6, 2*3*30, 2*6*15, 3*6*10, 6*30. - _Gus Wiseman_, Oct 10 2017

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

Cf. A001055, A005117, A045778, A046523, A050320, A050327, a(p^k)=0 (p>1), a(A002110)=A000110, a(n!)=A103775(n), A206778, A293243.

import Data.List (subsequences, genericIndex)
a050326 n = genericIndex a050326_list (n-1)
a050326_list = 1 : f 2 where
   f x = (if x /= s then a050326 s
                    else length $ filter (== x) $ map product $
                         subsequences $ tail $ a206778_row x) : f (x + 1)
         where s = a046523 x
-- Reinhard Zumkeller, May 03 2013

N:= 1000: # to get a(1)..a(N)
A:= Vector(N):
A[1]:= 1:
for n from 2 to N do
  if numtheory:-issqrfree(n) then
     S:= [$1..N/n]; T:= n*S; A[T]:= A[T]+A[S]
    fi;
od:
convert(A,list); # Robert Israel, Oct 10 2017

sqfacs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[sqfacs[n/d],Min@@#>d&]],{d,Select[Rest[Divisors[n]],SquareFreeQ]}]];
Table[Length[sqfacs[n]],{n,100}] (* Gus Wiseman, Oct 10 2017 *)

Dirichlet g.f.: prod{n is squarefree and > 1}(1+1/n^s).

a(n) = A050327(A101296(n)). - R. J. Mathar, May 26 2017

A068068 Number of odd unitary divisors of n. d is a unitary divisor of n if d divides n and gcd(d,n/d)=1.

Original entry on oeis.org

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

Offset: 1

Robert G. Wilson v, Feb 19 2002

Shadow transform of triangular numbers.
a(n) is the number of primitive Pythagorean triangles with inradius n. For the smallest inradius of exactly 2^n primitive Pythagorean triangles see A070826.
Number of primitive Pythagorean triangles with leg 4n. For smallest (even) leg of exactly 2^n PPTs, see A088860. - Lekraj Beedassy, Jul 12 2006
As shown by Chi and Killgrove, a(n) is the total number of primitive Pythagorean triples satisfying area = n * perimeter, or equivalently 2 raised to the power of the number of distinct, odd primes contained in n. - Ant King, Mar 15 2011
This is the case k=0 of the sum over the k-th powers of the odd unitary divisors of n, which is multiplicative with a(2^e)=1 and a(p^e)=1+p^(e*k), p>2, and has Dirichlet g.f. zeta(s)*zeta(s-k)*(1-2^(k-s))/( zeta(2s-k)*(1-2^(k-2*s)) ). - R. J. Mathar, Jun 20 2011
Also the number of odd squarefree divisors of n: a(n) = Sum_{k = 1..A034444(k)} (A077610(n,k) mod 2) = Sum_{k = 1..A034444(k)} (A206778(n,k) mod 2). - Reinhard Zumkeller, Feb 12 2012
a(n) is also the number of even unitary divisors of 2*n. - Amiram Eldar, Jan 28 2023

Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
Henjin Chi and Raymond Killgrove, Problem 1447, Crux Math 15(5), May 1989.
Henjin Chi and Raymond Killgrove, Solution to Problem 1447, Crux Math 16(7), September 1990.
L. J. Gerstein, Pythagorean triples and inner products, Math. Mag. 78 (2005), 205-213.
Lorenz Halbeisen and Norbert Hungerbuehler, Number theoretic aspects of a combinatorial function, Notes on Number Theory and Discrete Mathematics 5 (1999), 138-150. (ps, pdf); see Definition 7 for the shadow transform.
Lorenz Halbeisen, A number-theoretic conjecture and its implication for set theory, Acta Math. Univ. Comenianae 74(2) (2005), 243-254.
R. J. Mathar, Survey of Dirichlet series of multiplicative arithmetic functions, arXiv:1106.4038 [math.NT], 2011-2012.
OEIS Wiki, Shadow transform.
Neville Robbins, On the number of primitive Pythagorean triangles with a given inradius, Fibonacci Quart. 44(4) (2006), 368-369.
N. J. A. Sloane, Transforms.
Eric Weisstein's World of Mathematics, Unitary Divisor.
Wikipedia, Unitary divisor.

Cf. A001221, A001620, A024361, A034444, A056901, A068067, A008683, A323239.

a068068 = length . filter odd . a077610_row
-- Reinhard Zumkeller, Feb 12 2012

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

a[n_] := Length[Select[Divisors[n], OddQ[ # ]&&GCD[ #, n/# ]==1&]]
a[n_] := 2^(PrimeNu[n]+Mod[n, 2]-1); Array[a, 105] (* Jean-François Alcover, Dec 01 2015 *)
f[p_, e_] := If[p == 2, 1, 2]; a[1] = 1; a[n_] := Times @@ f @@@ FactorInteger[n]; Array[a, 100] (* Amiram Eldar, Sep 18 2020 *)

a(n) = sumdiv(n, d, (d%2)*(gcd(d, n/d)==1)); \\ Michel Marcus, May 13 2014

a(n) = 2^omega(n>>valuation(n,2)) \\ Charles R Greathouse IV, May 14 2014

a(n) = A034444(2n)/2. If n is even, a(n) = 2^(omega(n)-1); if n is odd, a(n) = 2^omega(n). Here omega(n) = A001221(n) is the number of distinct prime divisors of n.
Multiplicative with a(2^e) = 1, a(p^e) = 2, p>2. - Christian G. Bower May 18 2005
a(n) = A024361(4n). - Lekraj Beedassy, Jul 12 2006
Dirichlet g.f.: zeta^2(s)/ ( zeta(2*s)*(1+2^(-s)) ). Dirichlet convolution of A034444 and A154269. - R. J. Mathar, Apr 16 2011
a(n) = Sum_{d|n} mu(2*d)^2. - Ridouane Oudra, Aug 11 2019
Sum_{k=1..n} a(k) ~ 4*n*((log(n) + 2*gamma - 1 + log(2)/3) / Pi^2 - 12*zeta'(2) / Pi^4), where gamma is the Euler-Mascheroni constant A001620. - Vaclav Kotesovec, Sep 18 2020
a(n) = Sum_{d divides n, d odd} mu(d)^2. - Peter Bala, Feb 01 2024

Edited by Dean Hickerson, Jun 08 2002

cp's OEIS Frontend

A005117 Squarefree numbers: numbers that are not divisible by a square greater than 1.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Haskell

Magma

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

PARI

PARI

Python

Python

Python

Formula

A034444 a(n) is the number of unitary divisors of n (d such that d divides n, gcd(d, n/d) = 1).

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Magma

Magma

Maple

Mathematica

PARI

PARI

Python

Python

Scheme

Formula

Extensions

A048250 Sum of the squarefree divisors of n.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

Python

Sage

Formula

A077610 Triangle in which n-th row lists unitary divisors of n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI