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.

Previous Showing 41-50 of 146 results. Next

A034695 Tau_6 (the 6th Piltz divisor function), the number of ordered 6-factorizations of n; Dirichlet convolution of number-of-divisors function (A000005) with A007426.

Original entry on oeis.org

1, 6, 6, 21, 6, 36, 6, 56, 21, 36, 6, 126, 6, 36, 36, 126, 6, 126, 6, 126, 36, 36, 6, 336, 21, 36, 56, 126, 6, 216, 6, 252, 36, 36, 36, 441, 6, 36, 36, 336, 6, 216, 6, 126, 126, 36, 6, 756, 21, 126, 36, 126, 6, 336, 36, 336, 36, 36, 6, 756, 6, 36, 126, 462, 36, 216, 6, 126
Offset: 1

Views

Author

Keywords

References

  • T. M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, pages 29 and 38
  • Leveque, William J., Fundamentals of Number Theory. New York:Dover Publications, 1996, ISBN 9780486689067, p .167-Exercise 5.b.

Crossrefs

Cf. A000005 (tau_2), A007425 (tau_3), A007426 (tau_4), A061200 (tau_5).
Cf. A061204.
Column k=6 of A077592.

Programs

  • Mathematica
    tau[n_, 1] = 1; tau[n_, k_] := tau[n, k] = Plus @@ (tau[ #, k - 1] & /@ Divisors[n]); Table[ tau[n, 6], {n, 68}] (* Robert G. Wilson v, Nov 02 2005 *)
    tau[1, k_] := 1; tau[n_, k_] := Times @@ (Binomial[Last[#]+k-1, k-1]& /@ FactorInteger[n]); Table[tau[n, 6], {n, 1, 100}] (* Amiram Eldar, Sep 13 2020 *)
  • PARI
    a(n) = my(f=factor(n)); for (i=1, #f~, f[i,1] = binomial(f[i,2] + 5, f[i,2]); f[i,2]=1); factorback(f); \\ Michel Marcus, Jun 09 2014
    
  • PARI
    for(n=1, 100, print1(numerator(direuler(p=2, n, 1/(1-X)^6)[n]), ", ")) \\ Vaclav Kotesovec, May 06 2025
    
  • Python
    from math import prod, comb
    from sympy import factorint
    def A034695(n): return prod(comb(5+e,5) for e in factorint(n).values()) # Chai Wah Wu, Dec 22 2024

Formula

Dirichlet g.f.: zeta^6(s).
Multiplicative with a(p^e) = binomial(e+5, e). - Mitch Harris, Jun 27 2005
The Piltz divisor functions hold for tau_j(*)tau_k = tau_{j+k}, where (*) means Dirichlet convolution.
G.f.: Sum_{k>=1} tau_5(k)*x^k/(1 - x^k). - Ilya Gutkovskiy, Oct 30 2018

Extensions

More terms from Robert G. Wilson v, Nov 02 2005

A212171 Prime signature of n (nonincreasing version): row n of table lists positive exponents in canonical prime factorization of n, in nonincreasing order.

Original entry on oeis.org

1, 1, 2, 1, 1, 1, 1, 3, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 4, 1, 2, 1, 1, 2, 1, 1, 1, 1, 1, 1, 3, 1, 2, 1, 1, 3, 2, 1, 1, 1, 1, 1, 1, 5, 1, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 3, 1, 1, 1, 1, 1, 1, 2, 1, 2, 1, 1, 1, 1, 4, 1, 2, 2, 1, 1, 1, 2, 1, 1, 3, 1, 1, 1, 3, 1
Offset: 2

Views

Author

Matthew Vandermast, Jun 03 2012

Keywords

Comments

Length of row n equals A001221(n).
The multiset of positive exponents in n's prime factorization completely determines a(n) for a host of OEIS sequences, including several "core" sequences. Of those not cross-referenced here or in A212172, many can be found by searching the database for A025487.
(Note: Differing opinions may exist about whether the prime signature of n should be defined as this multiset itself, or as a symbol or collection of symbols that identify or "signify" this multiset. The definition of this sequence is designed to be compatible with either view, as are the original comments. When n >= 2, the customary ways to signify the multiset of exponents in n's prime factorization are to list the constituent exponents in either nonincreasing or nondecreasing order; this table gives the nonincreasing version.)
Table lists exponents in the order in which they appear in the prime factorization of a member of A025487. This ordering is common in database comments (e.g., A008966).
Each possible multiset of an integer's positive prime factorization exponents corresponds to a unique partition that contains the same elements (cf. A000041). This includes the multiset of 1's positive exponents, { } (the empty multiset), which corresponds to the partition of 0.
Differs from A124010 from a(23) on, corresponding to the factorization of 18 = 2^1*3^2 which is here listed as row 18 = [2, 1], but as [1, 2] (in the order of the prime factors) in A124010 and also in A118914 which lists the prime signatures in nondecreasing order (so that row 12 = 2^2*3^1 is also [1, 2]). - M. F. Hasler, Apr 08 2022

Examples

			First rows of table read:
  1;
  1;
  2;
  1;
  1,1;
  1;
  3;
  2;
  1,1;
  1;
  2,1;
  ...
The multiset of positive exponents in the prime factorization of 6 = 2*3 is {1,1} (1s are often left implicit as exponents). The prime signature of 6 is therefore {1,1}.
12 = 2^2*3 has positive exponents 2 and 1 in its prime factorization, as does 18 = 2*3^2. Rows 12 and 18 of the table both read {2,1}.
		

Crossrefs

Cf. A025487, A001221 (row lengths), A001222 (row sums). A118914 gives the nondecreasing version. A124010 lists exponents in n's prime factorization in natural order, with A124010(1) = 0.
A212172 cross-references over 20 sequences that depend solely on n's prime exponents >= 2, including the "core" sequence A000688. Other sequences determined by the exponents in the prime factorization of n include:
Additive: A001221, A001222, A056169.
A highly incomplete selection of sequences, each definable by the set of prime signatures possessed by its members: A000040, A000290, A000578, A000583, A000961, A001248, A001358, A001597, A001694, A002808, A004709, A005117, A006881, A013929, A030059, A030229, A052486.

Programs

  • Magma
    &cat[Reverse(Sort([pe[2]:pe in Factorisation(n)])):n in[1..76]]; // Jason Kimberley, Jun 13 2012
    
  • PARI
    apply( {A212171_row(n)=vecsort(factor(n)[,2]~,,4)}, [1..40])\\ M. F. Hasler, Apr 19 2022

Formula

Row n of A118914, reversed.
Row n of A124010 for n > 1, with exponents sorted in nonincreasing order. Equivalently, row A046523(n) of A124010 for n > 1.

A330570 Partial sums of A097988 (d_3(n)^2).

Original entry on oeis.org

1, 10, 19, 55, 64, 145, 154, 254, 290, 371, 380, 704, 713, 794, 875, 1100, 1109, 1433, 1442, 1766, 1847, 1928, 1937, 2837, 2873, 2954, 3054, 3378, 3387, 4116, 4125, 4566, 4647, 4728, 4809, 6105, 6114, 6195, 6276, 7176, 7185, 7914, 7923, 8247, 8571, 8652, 8661, 10686, 10722, 11046, 11127, 11451, 11460, 12360
Offset: 1

Views

Author

N. J. A. Sloane, Jan 08 2020

Keywords

Comments

This and the following sequences (and continuing in A331071) were inspired by the papers of Hooley, Indlekofer, Motohashi, Redmond, Titchmarsh, etc.

Crossrefs

Programs

  • Mathematica
    Accumulate[a[n_]:=DivisorSum[n, DivisorSigma[0, #]&]^2; Array[a, 60]] (* Vincenzo Librandi, Jan 11 2020 *)
  • PARI
    lista(nmax) = {my(s = 0); for(n = 1, nmax, s += vecprod(apply(e -> (e+1)*(e+2)/2, factor(n)[,2]))^2; print1(s, ", "));} \\ Amiram Eldar, Apr 19 2024

Formula

a(n) ~ c * n * log(n)^8 /8!, where c = Product_{p prime} ((1-1/p)^4 * (1 + 4/p + 1/p^2)) = 0.049321673579400091761... (Titchmarsh, 1942). - Amiram Eldar, Apr 19 2024

A077049 Left summatory matrix, T, by antidiagonals upwards.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 0, 0, 0, 1, 1, 1, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0
Offset: 1

Views

Author

Clark Kimberling, Oct 22 2002

Keywords

Comments

If S = (s(1), s(2), ...) is a sequence written as a column vector, then T*S is the summatory sequence of S; i.e., its n-th term is Sum_{k|n} s(k). T is the inverse of the left Moebius transformation matrix, A077050. Except for the first term in some cases, column 1 of T^(-2) is A007427, column 1 of T^(-1) is A008683, Column c of T^2 is A000005, column 1 of T^3 is A007425.
This is essentially the same as A051731, which includes only the triangle. Note that the standard in the OEIS is left to right antidiagonals, which would make this the right summatory matrix, and A077051 the left one. - Franklin T. Adams-Watters, Apr 08 2009
From Gary W. Adamson, Apr 28 2010: (Start)
As defined with antidiagonals of the array = the triangle shown in the example section. Row sums of this triangle = A032741 (with a different offset): 1, 1, 2, 1, 3, 1, 3, ...
Let the triangle = M. Then lim_{n->inf} M^n = A002033, the left-shifted vector considered as a sequence: (1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 8, ...). (End)

Examples

			T(4,2) = 1 since 2 divides 4. Northwest corner:
  1 0 0 0 0 0
  1 1 0 0 0 0
  1 0 1 0 0 0
  1 1 0 1 0 0
  1 0 0 0 1 0
  1 1 1 0 0 1
From _Gary W. Adamson_, Apr 28 2010: (Start)
First few rows of the triangle (when T is read by antidiagonals upwards):
  1;
  1, 0;
  1, 1, 0;
  1, 0, 0, 0;
  1, 1, 1, 0, 0;
  1, 0, 0, 0, 0, 0;
  1, 1, 0, 1, 0, 0, 0;
  1, 0, 1, 0, 0, 0, 0, 0;
  1, 1, 0, 0, 1, 0, 0, 0, 0;
  ... (End)
		

Crossrefs

Cf. A051731, A077050, A077051, A077052, A000005 (row sums).
Cf. A032741, A002033. - Gary W. Adamson, Apr 28 2010

Programs

  • Maple
    A077049 := proc(n,k)
        if modp(n,k) = 0 then
            1;
        else
            0 ;
        end if;
    end proc:
    for d from 2 to 10 do
        for k from 1 to d-1 do
            n := d-k ;
            printf("%d,",A077049(n,k)) ;
        end do:
    end do: # R. J. Mathar, Jul 22 2017
  • Mathematica
    With[{nn = 14}, DeleteCases[#, -1] & /@ Transpose@ Table[Take[#, nn] &@ Flatten@ Join[ConstantArray[-1, k - 1], ConstantArray[Reverse@ IntegerDigits[2^(k - 1), 2], Ceiling[(nn - k + 1)/k]]], {k, nn}]] // Flatten (* Michael De Vlieger, Jul 22 2017 *)
  • PARI
    nn=10; matrix(nn, nn, n, k, if (n % k, 0, 1)) \\ Michel Marcus, May 21 2015
    
  • Python
    def T(n, k):
        return 1 if n%k==0 else 0
    for n in range(1, 11): print([T(n - k + 1, k) for k in range(1, n + 1)]) # Indranil Ghosh, Jul 22 2017

Formula

T(n,k)=1 if k|n, otherwise T(n,k)=0, k >= 1, n >= 1.
From Boris Putievskiy, May 08 2013: (Start)
As table T(n,k) = floor(k/n) - floor((k-1)/n).
As linear sequence a(n) = floor(A004736(n)/A002260(n)) - floor((A004736(n)-1)/A002260(n)); a(n) = floor(j/i)-floor((j-1)/i), where i=n-t*(t+1)/2, j=(t*t+3*t+4)/2-n, t=floor((-1+sqrt(8*n-7))/2). (End)

Extensions

Name edited by Petros Hadjicostas, Jul 27 2019

A163767 a(n) = tau_{n}(n) = number of ordered n-factorizations of n.

Original entry on oeis.org

1, 2, 3, 10, 5, 36, 7, 120, 45, 100, 11, 936, 13, 196, 225, 3876, 17, 3078, 19, 4200, 441, 484, 23, 62400, 325, 676, 3654, 11368, 29, 27000, 31, 376992, 1089, 1156, 1225, 443556, 37, 1444, 1521, 459200, 41, 74088, 43, 43560, 46575, 2116, 47, 11995200, 1225
Offset: 1

Views

Author

Paul D. Hanna, Aug 04 2009

Keywords

Comments

Also the number of length n - 1 chains of divisors of n. - Gus Wiseman, May 07 2021

Examples

			Successive Dirichlet self-convolutions of the all 1's sequence begin:
(1),1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,... (A000012)
1,(2),2,3,2,4,2,4,3,4,2,6,2,4,4,5,... (A000005)
1,3,(3),6,3,9,3,10,6,9,3,18,3,9,9,15,... (A007425)
1,4,4,(10),4,16,4,20,10,16,4,40,4,16,16,35,... (A007426)
1,5,5,15,(5),25,5,35,15,25,5,75,5,25,25,70,... (A061200)
1,6,6,21,6,(36),6,56,21,36,6,126,6,36,36,126,... (A034695)
1,7,7,28,7,49,(7),84,28,49,7,196,7,49,49,210,... (A111217)
1,8,8,36,8,64,8,(120),36,64,8,288,8,64,64,330,... (A111218)
1,9,9,45,9,81,9,165,(45),81,9,405,9,81,81,495,... (A111219)
1,10,10,55,10,100,10,220,55,(100),10,550,10,100,... (A111220)
1,11,11,66,11,121,11,286,66,121,(11),726,11,121,... (A111221)
1,12,12,78,12,144,12,364,78,144,12,(936),12,144,... (A111306)
...
where the main diagonal forms this sequence.
From _Gus Wiseman_, May 07 2021: (Start)
The a(1) = 1 through a(5) = 5 chains of divisors:
  ()  (1)  (1/1)  (1/1/1)  (1/1/1/1)
      (2)  (3/1)  (2/1/1)  (5/1/1/1)
           (3/3)  (2/2/1)  (5/5/1/1)
                  (2/2/2)  (5/5/5/1)
                  (4/1/1)  (5/5/5/5)
                  (4/2/1)
                  (4/2/2)
                  (4/4/1)
                  (4/4/2)
                  (4/4/4)
(End)
		

Crossrefs

Main diagonal of A077592.
Diagonal n = k + 1 of the array A334997.
The version counting all multisets of divisors (not just chains) is A343935.
A000005 counts divisors.
A001055 counts factorizations (strict: A045778, ordered: A074206).
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A067824 counts strict chains of divisors starting with n.
A122651 counts strict chains of divisors summing to n.
A146291 counts divisors of n with k prime factors (with multiplicity).
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts nonempty strict chains of divisors of n.
A251683/A334996 count strict nonempty length-k divisor chains from n to 1.
A337255 counts strict length-k chains of divisors starting with n.
A339564 counts factorizations with a selected factor.
A343662 counts strict length-k chains of divisors (row sums: A337256).
Cf. A060690.

Programs

  • Mathematica
    Table[Times@@(Binomial[#+n-1,n-1]&/@FactorInteger[n][[All,2]]),{n,1,50}] (* Enrique Pérez Herrero, Dec 25 2013 *)
  • PARI
    {a(n,m=n)=if(n==1,1,if(m==1,1,sumdiv(n,d,a(d,1)*a(n/d,m-1))))}
    
  • Python
    from math import prod, comb
    from sympy import factorint
    def A163767(n): return prod(comb(n+e-1,e) for e in factorint(n).values()) # Chai Wah Wu, Jul 05 2024

Formula

a(p) = p for prime p.
a(n) = n^k when n is the product of k distinct primes (conjecture).
a(n) = n-th term of the n-th Dirichlet self-convolution of the all 1's sequence.
a(2^n) = A060690(n). - Alois P. Heinz, Jun 12 2024

A074816 a(n) = 3^A001221(n) = 3^omega(n).

Original entry on oeis.org

1, 3, 3, 3, 3, 9, 3, 3, 3, 9, 3, 9, 3, 9, 9, 3, 3, 9, 3, 9, 9, 9, 3, 9, 3, 9, 3, 9, 3, 27, 3, 3, 9, 9, 9, 9, 3, 9, 9, 9, 3, 27, 3, 9, 9, 9, 3, 9, 3, 9, 9, 9, 3, 9, 9, 9, 9, 9, 3, 27, 3, 9, 9, 3, 9, 27, 3, 9, 9, 27, 3, 9, 3, 9, 9, 9, 9, 27, 3, 9, 3, 9, 3, 27, 9, 9, 9, 9, 3, 27, 9, 9, 9, 9, 9, 9, 3, 9, 9, 9
Offset: 1

Views

Author

Benoit Cloitre, Sep 08 2002

Keywords

Comments

Old name was: a(n) = sum(d|n, tau(d)*mu(d)^2 ).
Terms are powers of 3.
The inverse Mobius transform of A074823, as the Dirichlet g.f. is product_{primes p} (1+2*p^(-s)) and the Dirichlet g.f. of A074816 is product_{primes p} (1+2*p^(-s))/(1-p^(-s)). - R. J. Mathar, Feb 09 2011
If n is squarefree, then a(n) = #{(x, y) : x, y positive integers, lcm (x, y) = n}. See Crandall & Pomerance. - Michel Marcus, Mar 23 2016

References

  • Richard Crandall and Carl Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2001; see Exercise 2.3 p. 108.

Crossrefs

Programs

Formula

a(n) = 3^m if n is divisible by m distinct primes. i.e., a(n)=3 if n is in A000961; a(n)=9 if n is in A007774 ...
a(n) = 3^A001221(n) = 3^omega(n). Multiplicative with a(p^e)=3. - Vladeta Jovovic, Sep 09 2002.
a(n) = tau_3(rad(n)) = A007425(A007947(n)). - Enrique Pérez Herrero, Jun 24 2010
a(n) = abs(Sum_{d|n} A000005(d^3)*mu(d)). - Enrique Pérez Herrero, Jun 28 2010
a(n) = Sum_{d|n, gcd(d, n/d) = 1} 2^omega(d) (The total number of unitary divisors of the unitary divisors of n). - Amiram Eldar, May 29 2020, Dec 27 2024
a(n) = Sum_{d1|n, d2|n} mu(d1*d2)^2. - Wesley Ivan Hurt, Feb 04 2022
Dirichlet g.f.: zeta(s)^3 * Product_{primes p} (1 - 3/p^(2*s) + 2/p^(3*s)). - Vaclav Kotesovec, Feb 16 2022

Extensions

Simpler definition at the suggestion of Michel Marcus. - N. J. A. Sloane, Mar 25 2016

A174465 G.f.: exp( Sum_{n>=1} A174466(n)*x^n/n ) where A174466(n) = Sum_{d|n} d*sigma(n/d)*tau(d).

Original entry on oeis.org

1, 1, 4, 7, 19, 31, 74, 122, 258, 430, 835, 1378, 2557, 4162, 7382, 11932, 20471, 32676, 54634, 86251, 141001, 220371, 353413, 546783, 863043, 1322425, 2057525, 3125092, 4801297, 7230393, 10984924, 16410474, 24679719, 36593278, 54526145, 80272501
Offset: 0

Views

Author

Paul D. Hanna, Apr 04 2010

Keywords

Comments

Compare to the g.f. of the number of planar partitions of n (A000219):
exp( Sum_{n>=1} sigma_2(n)*x^n/n ) where sigma_2(n) = Sum_{d|n} d*sigma(n/d)*phi(d).
tau(n) = A000005(n) = the number of divisors of n,
and sigma(n) = A000203(n) = sum of divisors of n.

Crossrefs

Programs

  • Mathematica
    nmax = 50; CoefficientList[Series[Product[1/(1-x^(i*j*k)), {i, 1, nmax}, {j, 1, nmax/i}, {k, 1, nmax/i/j}], {x, 0, nmax}], x] (* Vaclav Kotesovec, Jan 04 2017 *)
    nmax = 50; A007425 = Table[Sum[DivisorSigma[0, d], {d, Divisors[n]}], {n, 1, nmax}]; s = 1 - x; Do[s *= Sum[Binomial[A007425[[k]], j]*(-1)^j*x^(j*k), {j, 0, nmax/k}]; s = Expand[s]; s = Take[s, Min[nmax + 1, Exponent[s, x] + 1, Length[s]]];, {k, 2, nmax}]; CoefficientList[Series[1/s, {x, 0, nmax}], x] (* Vaclav Kotesovec, Aug 30 2018 *)
  • PARI
    {a(n)=polcoeff(exp(sum(m=1,n,x^m/m*sumdiv(m,d,d*sigma(m/d)*sigma(d,0)))+x*O(x^n)),n)}

Formula

G.f.: Product_{i>=1, j>=1, k>=1} 1/(1 - x^(i*j*k)). - Vaclav Kotesovec, Jan 04 2017
G.f.: Product_{k>=1} 1/(1 - x^k)^tau_3(k), where tau_3() = A007425. - Ilya Gutkovskiy, May 22 2018

A336422 Number of ways to choose a divisor of a divisor of n, both having distinct prime exponents.

Original entry on oeis.org

1, 3, 3, 6, 3, 5, 3, 10, 6, 5, 3, 13, 3, 5, 5, 15, 3, 13, 3, 13, 5, 5, 3, 24, 6, 5, 10, 13, 3, 7, 3, 21, 5, 5, 5, 21, 3, 5, 5, 24, 3, 7, 3, 13, 13, 5, 3, 38, 6, 13, 5, 13, 3, 24, 5, 24, 5, 5, 3, 20, 3, 5, 13, 28, 5, 7, 3, 13, 5, 7, 3, 42, 3, 5, 13, 13, 5, 7, 3
Offset: 1

Views

Author

Gus Wiseman, Jul 26 2020

Keywords

Comments

A number has distinct prime exponents iff its prime signature is strict.

Examples

			The a(n) ways for n = 1, 2, 4, 6, 8, 12, 30, 210:
  1/1/1  2/1/1  4/1/1  6/1/1  8/1/1  12/1/1    30/1/1  210/1/1
         2/2/1  4/2/1  6/2/1  8/2/1  12/2/1    30/2/1  210/2/1
         2/2/2  4/2/2  6/2/2  8/2/2  12/2/2    30/2/2  210/2/2
                4/4/1  6/3/1  8/4/1  12/3/1    30/3/1  210/3/1
                4/4/2  6/3/3  8/4/2  12/3/3    30/3/3  210/3/3
                4/4/4         8/4/4  12/4/1    30/5/1  210/5/1
                              8/8/1  12/4/2    30/5/5  210/5/5
                              8/8/2  12/4/4            210/7/1
                              8/8/4  12/12/1           210/7/7
                              8/8/8  12/12/2
                                     12/12/3
                                     12/12/4
                                     12/12/12
		

Crossrefs

A336421 is the case of superprimorials.
A007425 counts divisors of divisors.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A327498 gives the maximum divisor with distinct prime exponents.
A336500 counts divisors with quotient also having distinct prime exponents.
A336568 = not a product of two numbers with distinct prime exponents.

Programs

  • Mathematica
    strdivs[n_]:=Select[Divisors[n],UnsameQ@@Last/@FactorInteger[#]&];
    Table[Sum[Length[strdivs[d]],{d,strdivs[n]}],{n,30}]

A336569 Number of maximal strict chains of divisors from n to 1 using elements of A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 2, 1, 0, 0, 1, 1, 2, 1, 2, 0, 0, 1, 3, 1, 0, 1, 2, 1, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 3, 1, 0, 1, 2, 2, 0, 1, 4, 1, 2, 0, 2, 1, 3, 0, 3, 0, 0, 1, 0, 1, 0, 2, 1, 0, 0, 1, 2, 0, 0, 1, 5, 1, 0, 2, 2, 0, 0, 1, 4, 1, 0, 1, 0, 0, 0, 0
Offset: 1

Views

Author

Gus Wiseman, Jul 29 2020

Keywords

Comments

A number's prime signature (row n of A124010) is the sequence of positive exponents in its prime factorization, so a number has distinct prime multiplicities iff all the exponents in its prime signature are distinct.

Examples

			The a(n) chains for n = 12, 72, 144, 192 (ones not shown):
  12/3    72/18/2       144/72/18/2       192/96/48/24/12/3
  12/4/2  72/18/9/3     144/72/18/9/3     192/64/32/16/8/4/2
          72/24/12/3    144/48/24/12/3    192/96/32/16/8/4/2
          72/24/8/4/2   144/72/24/12/3    192/96/48/16/8/4/2
          72/24/12/4/2  144/48/16/8/4/2   192/96/48/24/8/4/2
                        144/48/24/8/4/2   192/96/48/24/12/4/2
                        144/72/24/8/4/2
                        144/48/24/12/4/2
                        144/72/24/12/4/2
		

Crossrefs

A336423 is the non-maximal version.
A336570 is the version for chains not necessarily containing n.
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A007425 counts divisors of divisors.
A032741 counts proper divisors.
A045778 counts strict factorizations.
A071625 counts distinct prime multiplicities.
A074206 counts strict chains of divisors from n to 1.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts chains of divisors.
A336422 counts divisible pairs of divisors, both in A130091.
A336424 counts factorizations using A130091.
A336571 counts divisor sets of elements of A130091.

Programs

  • Mathematica
    strsigQ[n_]:=UnsameQ@@Last/@FactorInteger[n];
    fasmax[y_]:=Complement[y,Union@@(Most[Subsets[#]]&/@y)];
    strchs[n_]:=If[n==1,{{}},If[!strsigQ[n],{},Join@@Table[Prepend[#,d]&/@strchs[d],{d,Select[Most[Divisors[n]],strsigQ]}]]];
    Table[Length[fasmax[strchs[n]]],{n,100}]

A211776 a(n) = Product_{d | n} tau(d).

Original entry on oeis.org

1, 2, 2, 6, 2, 16, 2, 24, 6, 16, 2, 288, 2, 16, 16, 120, 2, 288, 2, 288, 16, 16, 2, 9216, 6, 16, 24, 288, 2, 4096, 2, 720, 16, 16, 16, 46656, 2, 16, 16, 9216, 2, 4096, 2, 288, 288, 16, 2, 460800, 6, 288, 16, 288, 2, 9216, 16, 9216, 16, 16, 2, 5308416, 2
Offset: 1

Views

Author

Jaroslav Krizek, Apr 20 2012

Keywords

Examples

			For n = 6: divisors of 6: 1, 2, 3, 6; tau(d): 1, 2, 2, 4; product _{d | n} tau(d) = 1*2*2*4 = 16, where tau = A000005.
		

Crossrefs

Cf. A000005, A001221, A007425 (Sum_{d | n} tau(d)).

Programs

  • Maple
    A211776 := proc(n)
        mul( A000005(d),d=numtheory[divisors](n)) ;
    end proc:
    seq(A211776(n),n=1..20) ; # R. J. Mathar, Feb 13 2019
  • Mathematica
    Table[Product[DivisorSigma[0, i], {i, Divisors[n]}], {n, 100}] (* T. D. Noe, Apr 26 2012 *)
    a[1] = 1; a[n_] := Module[{e = FactorInteger[n][[;; , 2]]}, d = Times @@ (e + 1); Times @@ ((e + 1)!^(d/(e + 1)))]; Array[a, 100] (* using the Formula section,  Amiram Eldar, Aug 04 2020 *)
  • PARI
    A211776(n) = { my(m=1); fordiv(n, d, m *= numdiv(d)); m };
    A211776(n) = prod(d=1, n, if((n % d), 1, numdiv(d)));
    \\ Antti Karttunen, May 19 2017

Formula

a(n) = Product_{i=1..omega(n)} (b_i+1)!^(tau(n)/(b_i+1)), where omega(n) is the number of distinct prime factors of n, tau(n) is the number of divisors of n, and n = p_1^(b_1)*p_2^(b_2)* ... *p_{omega(n)}^(b_{omega(n)}). - Anand Rao Tadipatri, Aug 04 2020
Previous Showing 41-50 of 146 results. Next