A177057 -id:A177057 - cp's OEIS frontend

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

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

A020793 Decimal expansion of 1/6.

Original entry on oeis.org

1, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6

Offset: 0

N. J. A. Sloane

Except for the first term identical to A010722, A040006 and A021019. Except for the first terms the same as A021028, A021100, A021388, A071279, A101272, A168608, A177057,... - M. F. Hasler, Oct 24 2011

Decimal expansion of gamma(1) = 5/3 (with offset 1), where gamma(n) = Cp(n)/Cv(n) = is the n-th Poisson's constant. For the definition of Cp and Cv see A272002. - Natan Arie Consigli, Jul 10 2016

Calvin C. Clawson, Mathematical Mysteries, The Beauty and Magic of Numbers, Springer, 2013, see p. 224.

Wikipedia, Poisson's constant.
Index entries for linear recurrences with constant coefficients, signature (1).

Cf. A010722, A021019, A021028, A021100, A021388, A040006, A070064, A071279, A081822, A101272, A168608, A177057, A272001, A272002.

RealDigits[1/6,10,120][[1]] (* or *) PadRight[{1},120,{6}] (* Harvey P. Dale, Dec 30 2018 *)

a(n)=6-5*!n  \\ M. F. Hasler, Oct 24 2011

a(n) = 6^n mod 10. - Zerinvary Lajos, Nov 26 2009
Equals Sum_{k>=1} 1/7^k. - Bruno Berselli, Jan 03 2014
10 * 1/6 = 5/3 = (5/2 R)/(3/2 R) = Cp(1)/Cv(1) = A272002/A272001, with R = A081822 (or A070064). - Natan Arie Consigli, Jul 10 2016
G.f.: (1 + 5*x)/(1 - x). - Ilya Gutkovskiy, Jul 10 2016
Equals Sum_{k>=1} 1/(k*Pi)^2. - Maciej Kaniewski, Sep 14 2017
Equals Sum_{k>=1} (zeta(2*k)-1)/4^k. - Amiram Eldar, Jun 08 2021
K_{n>=2} 2*n/(2*n - 3) = 5/3. (see Clawson at p. 224). - Stefano Spezia, Jul 01 2024
E.g.f.: 6*exp(x) - 5. - Elmo R. Oliveira, Aug 05 2024

A341029 Twice the area of the convex hull around dragon curve expansion level n.

Original entry on oeis.org

0, 1, 3, 9, 23, 56, 121, 258, 539, 1118, 2273, 4614, 9323, 18806, 37761, 75798, 151979, 304598, 609793, 1220694, 2442923, 4888406, 9779201, 19562838, 39131819, 78273878, 156557313, 313132374, 626289323, 1252619606, 2505277441, 5010625878, 10021350059

Offset: 0

Kevin Ryde, Feb 02 2021

The area of the hull is a half-integer for n=1..4 and even n>=6, so the sequence is a(n) = 2*area to give integers.

Benedek and Panzone determine the vertices of the convex hull around the dragon fractal. The area of that hull is 7/6 (A177057). This is the limit for the finite expansions scaled down to a unit distance start to end: lim_{n->oo} (a(n)/2) / 2^n = 7/6.

			  @   *---@          curve expansion level n=3,
  |   |   |          convex hull vertices marked "@",
  @---*   *---@      area = 4+1/2,
              |      a(3) = 2*area = 9
          @---@

Kevin Ryde, Table of n, a(n) for n = 0..600
Agnes I. Benedek and Rafael Panzone, On Some Notable Plane Sets, II: Dragons, Revista de la Unión Matemática Argentina, volume 39, numbers 1-2, 1994, pages 76-90.
Kevin Ryde, Iterations of the Dragon Curve, see index "HA".
Index entries for linear recurrences with constant coefficients, signature (3,-3,3,2,-12,12,-12,8).

Cf. A177057 (fractal hull area), A341030 (fractal hull perimeter).

Cf. A362566 (bounding box area).

my(h=[22,29,22,31]); a(n) = if(n<2,n, (7<

For n>=2, a(n) = (7/3)*2^n - (h/6)*2^floor(n/2) + c/3, where h = 22,29,22,31 and c = 1,2,3,2 according as n == 0,1,2,3 (mod 4) respectively.
a(n) = 3*a(n-1) - 3*a(n-2) + 3*a(n-3) + 2*a(n-4) - 12*a(n-5) + 12*a(n-6) - 12*a(n-7) + 8*a(n-8) for n>=10.
G.f.: x*(1 + 3*x^2 + 2*x^3 + 3*x^4 + x^5 - 2*x^7 - 4*x^8) /( (1-x) * (1-2*x) * (1+x^2) * (1-2*x^2) * (1+2*x^2) ).
G.f.: 1 + (1/2)*x + (2/3)/(1-x) - (1/3)/(1+x^2) + (1/6)*x/(1+2*x^2) - (11/3 + 5*x)/(1-2*x^2) + (7/3)/(1-2*x).

cp's OEIS Frontend

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

A020793 Decimal expansion of 1/6.

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A341029 Twice the area of the convex hull around dragon curve expansion level n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

Formula

A177022 Decimal expansion of 49/36.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions