cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

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

Views

Author

N. J. A. Sloane

Keywords

Comments

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

Examples

			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.

References

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

Links

Crossrefs

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

Programs

  • Haskell
    a034444 = length . a077610_row  -- Reinhard Zumkeller, Feb 12 2012
  • Magma
    [#[d:d in Divisors(n)|Gcd(d,n div d) eq 1]:n in [1..110]]; // Marius A. Burtea, Jan 11 2020
  • Magma
    [&+[Abs(MoebiusMu(d)):d in Divisors(n)]:n in [1..110]]; // Marius A. Burtea, Jan 11 2020
  • Maple
    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
  • Mathematica
    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 *)
  • PARI
    a(n)=1<Charles R Greathouse IV, Feb 11 2011
  • PARI
    for(n=1, 100, print1(direuler(p=2, n, (1+X)/(1-X))[n], ", ")) \\ Vaclav Kotesovec, Sep 26 2020
  • Python
    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
  • Python
    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
  • Scheme
    (define (A034444 n) (if (= 1 n) n (* 2 (A034444 (A028234 n))))) ;; Antti Karttunen, May 29 2017

Formula

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

Extensions

More terms from James Sellers, Jun 20 2000

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.