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 31-40 of 155 results. Next

A325248 Heinz number of the omega-sequence of n.

Original entry on oeis.org

1, 2, 2, 6, 2, 18, 2, 10, 6, 18, 2, 90, 2, 18, 18, 14, 2, 90, 2, 90, 18, 18, 2, 126, 6, 18, 10, 90, 2, 50, 2, 22, 18, 18, 18, 42, 2, 18, 18, 126, 2, 50, 2, 90, 90, 18, 2, 198, 6, 90, 18, 90, 2, 126, 18, 126, 18, 18, 2, 630, 2, 18, 90, 26, 18, 50, 2, 90, 18, 50
Offset: 1

Views

Author

Gus Wiseman, Apr 16 2019

Keywords

Comments

We define the omega-sequence of n (row n of A323023) to have length A323014(n) = adjusted frequency depth of n, and the k-th term is Omega(red^{k-1}(n)), where Omega = A001222 and red^{k} is the k-th functional iteration of red = A181819, defined by red(n = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of n. For example, we have 180 -> 18 -> 6 -> 4 -> 3, so the omega-sequence of 180 is (5,3,2,2,1).
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The omega-sequence of 180 is (5,3,2,2,1) with Heinz number 990, so a(180) = 990.
		

Crossrefs

Positions of squarefree terms are A325247.
Positions of normal numbers (A055932) are A325251.
First positions of each distinct term are A325238.
Omega-sequence statistics: A001222 (first omega), A001221 (second omega), A071625 (third omega), A323022 (fourth omega), A304465 (second-to-last omega), A182850 or A323014 (length/frequency depth), A325248 (Heinz number).

Programs

  • Mathematica
    omseq[n_Integer]:=If[n<=1,{},Total/@NestWhileList[Sort[Length/@Split[#]]&,Sort[Last/@FactorInteger[n]],Total[#]>1&]];
    Table[Times@@Prime/@omseq[n],{n,100}]

Formula

A001222(a(n)) = A323014(n).
A061395(a(n)) = A001222(n).
A304465(n) = A055396(a(n)/2).
A325249(n) = A056239(a(n)).
a(n!) = A325275(n).

A336424 Number of factorizations of n where each factor belongs to A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 2, 1, 1, 1, 3, 2, 1, 1, 3, 1, 1, 1, 5, 1, 3, 1, 3, 1, 1, 1, 5, 2, 1, 3, 3, 1, 1, 1, 7, 1, 1, 1, 6, 1, 1, 1, 5, 1, 1, 1, 3, 3, 1, 1, 9, 2, 3, 1, 3, 1, 5, 1, 5, 1, 1, 1, 4, 1, 1, 3, 11, 1, 1, 1, 3, 1, 1, 1, 11, 1, 1, 3, 3, 1, 1, 1, 9, 5, 1, 1, 4, 1, 1
Offset: 1

Views

Author

Gus Wiseman, Aug 03 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) factorizations for n = 2, 4, 8, 60, 16, 36, 32, 48:
  2  4    8      5*12     16       4*9      32         48
     2*2  2*4    3*20     4*4      3*12     4*8        4*12
          2*2*2  3*4*5    2*8      3*3*4    2*16       3*16
                 2*2*3*5  2*2*4    2*18     2*4*4      3*4*4
                          2*2*2*2  2*2*9    2*2*8      2*24
                                   2*2*3*3  2*2*2*4    2*3*8
                                            2*2*2*2*2  2*2*12
                                                       2*2*3*4
                                                       2*2*2*2*3
		

Crossrefs

A327523 is the case when n is restricted to belong to A130091 also.
A001055 counts factorizations.
A007425 counts divisors of divisors.
A045778 counts strict factorizations.
A074206 counts ordered factorizations.
A130091 lists numbers with distinct prime multiplicities.
A181796 counts divisors with distinct prime multiplicities.
A253249 counts nonempty chains of divisors.
A281116 counts factorizations with no common divisor.
A302696 lists numbers whose prime indices are pairwise coprime.
A305149 counts stable factorizations.
A320439 counts factorizations using A289509.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336500 counts divisors of n in A130091 with quotient also in A130091.
A336568 = not a product of two numbers with distinct prime multiplicities.
A336569 counts maximal chains of elements of A130091.
A337256 counts chains of divisors.

Programs

  • Mathematica
    facsusing[s_,n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facsusing[Select[s,Divisible[n/d,#]&],n/d],Min@@#>=d&]],{d,Select[s,Divisible[n,#]&]}]];
    Table[Length[facsusing[Select[Range[2,n],UnsameQ@@Last/@FactorInteger[#]&],n]],{n,100}]

A353834 Nonprime numbers whose prime indices have all equal run-sums.

Original entry on oeis.org

1, 4, 8, 9, 12, 16, 25, 27, 32, 40, 49, 63, 64, 81, 112, 121, 125, 128, 144, 169, 243, 256, 289, 325, 343, 351, 352, 361, 512, 529, 625, 675, 729, 832, 841, 931, 961, 1008, 1024, 1331, 1369, 1539, 1600, 1681, 1728, 1849, 2048, 2176, 2187, 2197, 2209, 2401
Offset: 1

Views

Author

Gus Wiseman, May 26 2022

Keywords

Comments

A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
The sequence of runs of a sequence consists of its maximal consecutive constant subsequences when read left-to-right. For example, the runs of (2,2,1,1,1,3,2,2) are (2,2), (1,1,1), (3), (2,2), with sums (4,3,3,4).

Examples

			The terms together with their prime indices begin:
     1: {}
     4: {1,1}
     8: {1,1,1}
     9: {2,2}
    12: {1,1,2}
    16: {1,1,1,1}
    25: {3,3}
    27: {2,2,2}
    32: {1,1,1,1,1}
    40: {1,1,1,3}
    49: {4,4}
    63: {2,2,4}
    64: {1,1,1,1,1,1}
    81: {2,2,2,2}
   112: {1,1,1,1,4}
   121: {5,5}
   125: {3,3,3}
   128: {1,1,1,1,1,1,1}
For example, 675 is in the sequence because its prime indices {2,2,2,3,3} have run-sums (6,6).
		

Crossrefs

For equal run-lengths we have A072774\A000040, counted by A047966(n)-1.
These partitions are counted by A304442(n) - 1.
These are the nonprime positions of prime powers in A353832.
Including the primes gives A353833.
For distinct run-sums we have A353838\A000040, counted by A353837(n)-1.
For compositions we have A353848\A000079, counted by A353851(n)-1.
A001222 counts prime factors, distinct A001221.
A005811 counts runs in binary expansion, distinct run-lengths A165413.
A056239 adds up prime indices, row sums of A112798 and A296150.
A124010 gives prime signature, sorted A118914.
A300273 ranks collapsible partitions, counted by A275870.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353840-A353846 pertain to partition run-sum trajectory.
A353862 gives greatest run-sum of prime indices, least A353931.
A353866 ranks rucksack partitions, counted by A353864.

Programs

  • Mathematica
    Select[Range[100],!PrimeQ[#]&&SameQ@@Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]*k]&]
  • Python
    from itertools import count, islice
    from sympy import factorint, primepi
    def A353848_gen(startvalue=1): # generator of terms >= startvalue
        return filter(lambda n: n == 1 or (sum((f:=factorint(n)).values()) > 1 and len(set(primepi(p)*e for p, e in f.items())) <= 1), count(max(startvalue,1)))
    A353848_list = list(islice(A353848_gen(),30)) # Chai Wah Wu, May 27 2022

A336571 Number of sets of divisors d|n, 1 < d < n, all belonging to A130091 (numbers with distinct prime multiplicities) and forming a divisibility chain.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 4, 2, 3, 1, 5, 1, 3, 3, 8, 1, 5, 1, 5, 3, 3, 1, 14, 2, 3, 4, 5, 1, 4, 1, 16, 3, 3, 3, 17, 1, 3, 3, 14, 1, 4, 1, 5, 5, 3, 1, 36, 2, 5, 3, 5, 1, 14, 3, 14, 3, 3, 1, 16, 1, 3, 5, 32, 3, 4, 1, 5, 3, 4, 1, 35, 1, 3, 5, 5, 3, 4, 1, 36, 8, 3, 1
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) sets for n = 4, 6, 12, 16, 24, 84, 36:
  {}   {}   {}     {}       {}        {}        {}
  {2}  {2}  {2}    {2}      {2}       {2}       {2}
       {3}  {3}    {4}      {3}       {3}       {3}
            {4}    {8}      {4}       {4}       {4}
            {2,4}  {2,4}    {8}       {7}       {9}
                   {2,8}    {12}      {12}      {12}
                   {4,8}    {2,4}     {28}      {18}
                   {2,4,8}  {2,8}     {2,4}     {2,4}
                            {4,8}     {2,12}    {3,9}
                            {2,12}    {2,28}    {2,12}
                            {3,12}    {3,12}    {2,18}
                            {4,12}    {4,12}    {3,12}
                            {2,4,8}   {4,28}    {3,18}
                            {2,4,12}  {7,28}    {4,12}
                                      {2,4,12}  {9,18}
                                      {2,4,28}  {2,4,12}
                                                {3,9,18}
		

Crossrefs

A336423 is the version for chains containing n.
A336570 is the maximal version.
A000005 counts divisors.
A001055 counts factorizations.
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.
A336500 counts divisors of n in A130091 with quotient also in A130091.

Programs

  • Mathematica
    strchns[n_]:=If[n==1,1,Sum[strchns[d],{d,Select[Most[Divisors[n]],UnsameQ@@Last/@FactorInteger[#]&]}]];
    Table[strchns[n],{n,100}]

A353744 Numbers k such that the k-th composition in standard order has all equal run-lengths.

Original entry on oeis.org

0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 15, 16, 17, 18, 20, 22, 24, 25, 31, 32, 33, 34, 36, 37, 38, 40, 41, 42, 43, 44, 45, 48, 49, 50, 52, 54, 58, 63, 64, 65, 66, 68, 69, 70, 72, 76, 77, 80, 81, 82, 88, 89, 96, 97, 98, 101, 102, 104, 105, 108, 109, 127, 128
Offset: 1

Views

Author

Gus Wiseman, Jun 11 2022

Keywords

Comments

The k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. This gives a bijective correspondence between nonnegative integers and integer compositions.

Examples

			Composition 2362 in standard order is (3,3,1,1,2,2), with run-lengths (2,2,2), so 2362 is in the sequence.
		

Crossrefs

Standard compositions are listed by A066099.
The version for partitions is A072774, counted by A047966.
These compositions are counted by A329738.
For distinct instead of equal run-lengths we have A351596.
For run-sums instead of lengths we have A353848, counted by A353851.
For distinct run-sums we have A353852, counted by A353850.
A003242 counts anti-run compositions, ranked by A333489.
A005811 counts runs in binary expansion.
A300273 ranks collapsible partitions, counted by A275870.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353847 represents the composition run-sum transformation.
A353853-A353859 pertain to composition run-sum trajectory.
A353860 counts collapsible compositions.
A353833 ranks partitions with all equal run-sums, counted by A304442.

Programs

  • Mathematica
    stc[n_]:=Differences[Prepend[Join@@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
    Select[Range[0,100],SameQ@@Length/@Split[stc[#]]&]

A182853 Squarefree composite integers and powers of squarefree composite integers.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 26, 30, 33, 34, 35, 36, 38, 39, 42, 46, 51, 55, 57, 58, 62, 65, 66, 69, 70, 74, 77, 78, 82, 85, 86, 87, 91, 93, 94, 95, 100, 102, 105, 106, 110, 111, 114, 115, 118, 119, 122, 123, 129, 130, 133, 134, 138, 141, 142, 143, 145, 146, 154, 155, 158, 159, 161
Offset: 1

Views

Author

Matthew Vandermast, Jan 04 2011

Keywords

Comments

Numbers that require exactly three iterations to reach a fixed point under the x -> A181819(x) map. In each case, 2 is the fixed point that is reached. (1 is the other fixed point of the x -> A181819(x) map.) Cf. A182850.
Numbers such that A001221(n) > 1 and A071625(n) = 1.

Crossrefs

Numbers n such that A182850(n) = 3. See also A182854, A182855.
Subsequence of A072774 and A182851.
Cf. A120944.

Programs

  • PARI
    isoka(n) = (omega(n) > 1) && issquarefree(n); \\ A120944
    isok(n) = isoka(n) || (ispower(n,,&k) && isoka(k)); \\ Michel Marcus, Jun 24 2017
    
  • Python
    from math import isqrt
    from sympy import mobius, primepi, integer_nthroot
    def A182853(n):
        def g(x): return int(sum(mobius(k)*(x//k**2) for k in range(1, isqrt(x)+1))-primepi(x))
        def f(x): return n-2+x+(y:=x.bit_length())-sum(g(integer_nthroot(x,k)[0]) for k in range(1,y))
        kmin, kmax = 1,2
        while f(kmax) >= kmax:
            kmax <<= 1
        while True:
            kmid = kmax+kmin>>1
            if f(kmid) < kmid:
                kmax = kmid
            else:
                kmin = kmid
            if kmax-kmin <= 1:
                break
        return kmax # Chai Wah Wu, Aug 19 2024
  • Scheme
    (define A182853 (MATCHING-POS 1 1 (lambda (n) (= 3 (A182850 n))))) ;; After the alternative definition of the sequence given by the original author. Requires also MATCHING-POS macro from my IntSeq-library - Antti Karttunen, Feb 05 2016
    

A280292 a(n) = sopfr(n) - sopf(n).

Original entry on oeis.org

0, 0, 0, 2, 0, 0, 0, 4, 3, 0, 0, 2, 0, 0, 0, 6, 0, 3, 0, 2, 0, 0, 0, 4, 5, 0, 6, 2, 0, 0, 0, 8, 0, 0, 0, 5, 0, 0, 0, 4, 0, 0, 0, 2, 3, 0, 0, 6, 7, 5, 0, 2, 0, 6, 0, 4, 0, 0, 0, 2, 0, 0, 3, 10, 0, 0, 0, 2, 0, 0, 0, 7, 0, 0, 5, 2, 0, 0, 0, 6, 9, 0, 0, 2, 0, 0, 0, 4, 0, 3, 0, 2, 0, 0, 0, 8, 0, 7, 3, 7, 0, 0, 0, 4, 0
Offset: 1

Views

Author

Michel Marcus, Dec 31 2016

Keywords

Comments

Alladi and Erdős (1977) proved that for all numbers m>=0, m!=1, the sequence of numbers k such that a(k) = m has a positive asymptotic density which is equal to a rational multiple of 1/zeta(2) = 6/Pi^2 (A059956). For example, when m=0, the sequence is the squarefree numbers (A005117), whose density is 6/Pi^2, and when m=2 the sequence is A081770, whose density is 1/Pi^2. - Amiram Eldar, Nov 02 2020
Sum of prime factors minus sum of distinct prime factors. Counting partitions by this statistic (sum minus sum of distinct parts) gives A364916. - Gus Wiseman, Feb 21 2025

References

  • Jean-Marie De Koninck and Aleksandar Ivić, Topics in Arithmetical Functions: Asymptotic Formulae for Sums of Reciprocals of Arithmetical Functions and Related Fields, Amsterdam, Netherlands: North-Holland, 1980. See pp. 164-166.
  • Steven R. Finch, Mathematical Constants II, Encyclopedia of Mathematics and Its Applications, Cambridge University Press, Cambridge, 2018, p. 165.

Crossrefs

A multiplicative version is A003557, firsts A064549 (sorted A001694).
For length instead of sum we have A046660.
For product instead of sum we have A066503, firsts A381076.
Positions of first appearances are A280286 (sorted A381075).
For indices instead of factors we have A380955, firsts A380956 (sorted A380957).
For exponents instead of factors we have A380958, firsts A380989.
A000040 lists the primes, differences A001223.
A001222 counts prime factors (distinct A001221).
A003963 gives product of prime indices, distinct A156061, excess A380986.
A005117 lists squarefree numbers, complement A013929.
A007947 gives squarefree kernel.
A020639 gives least prime factor (index A055396), greatest A061395 (index A006530).
A027746 lists prime factors, distinct A027748.
A112798 lists prime indices (sum A056239), distinct A304038 (sum A066328).

Programs

  • Mathematica
    Array[Total@ # - Total@ Union@ # &@ Flatten[ConstantArray[#1, #2] & @@@ FactorInteger@ #] &, 105] (* Michael De Vlieger, Feb 25 2019 *)
  • PARI
    sopfr(n) = my(f=factor(n)); sum(j=1, #f~, f[j, 1]*f[j, 2]);
    sopf(n) = my(f=factor(n)); sum(j=1, #f~, f[j, 1]);
    a(n) = sopfr(n) - sopf(n);

Formula

a(n) = A001414(n) - A008472(n).
a(A005117(n)) = 0.
a(n) = A001414(A003557(n)). - Antti Karttunen, Oct 07 2017
Additive with a(1) = 0 and a(p^e) = p*(e-1) for prime p and e > 0. - Werner Schulte, Feb 24 2019
From Amiram Eldar, Nov 02 2020: (Start)
a(n) = a(A057521(n)).
Sum_{n<=x} a(n) ~ x*log(log(x)) + O(x) (Alladi and Erdős, 1977).
Sum_{n<=x, n nonsquarefree} 1/a(n) ~ c*x + O(sqrt(x)*log(x)), where c = Integral_{t=0..1} (F(t)-6/Pi^2)/t dt, and F(t) = Product_{p prime} (1-1/p)*(1-1/(t^p - p)) (De Koninck et al., 1981; Finch, 2018), or, equivalently c = Sum_{k>=2} d(k)/k = 0.1039..., where d(k) = (6/Pi^2)*A338559(k)/A338560(k) is the asymptotic density of the numbers m with a(m) = k (Alladi and Erdős, 1977; Ivić, 2003). (End)

Extensions

More terms from Antti Karttunen, Oct 07 2017

A336423 Number of strict chains of divisors from n to 1 using terms of A130091 (numbers with distinct prime multiplicities).

Original entry on oeis.org

1, 1, 1, 2, 1, 0, 1, 4, 2, 0, 1, 5, 1, 0, 0, 8, 1, 5, 1, 5, 0, 0, 1, 14, 2, 0, 4, 5, 1, 0, 1, 16, 0, 0, 0, 0, 1, 0, 0, 14, 1, 0, 1, 5, 5, 0, 1, 36, 2, 5, 0, 5, 1, 14, 0, 14, 0, 0, 1, 0, 1, 0, 5, 32, 0, 0, 1, 5, 0, 0, 1, 35, 1, 0, 5, 5, 0, 0, 1, 36, 8, 0, 1, 0
Offset: 1

Views

Author

Gus Wiseman, Jul 27 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 = 4, 8, 12, 16, 24, 32:
  4/1    8/1      12/1      16/1        24/1         32/1
  4/2/1  8/2/1    12/2/1    16/2/1      24/2/1       32/2/1
         8/4/1    12/3/1    16/4/1      24/3/1       32/4/1
         8/4/2/1  12/4/1    16/8/1      24/4/1       32/8/1
                  12/4/2/1  16/4/2/1    24/8/1       32/16/1
                            16/8/2/1    24/12/1      32/4/2/1
                            16/8/4/1    24/4/2/1     32/8/2/1
                            16/8/4/2/1  24/8/2/1     32/8/4/1
                                        24/8/4/1     32/16/2/1
                                        24/12/2/1    32/16/4/1
                                        24/12/3/1    32/16/8/1
                                        24/12/4/1    32/8/4/2/1
                                        24/8/4/2/1   32/16/4/2/1
                                        24/12/4/2/1  32/16/8/2/1
                                                     32/16/8/4/1
                                                     32/16/8/4/2/1
		

Crossrefs

A336569 is the maximal case.
A336571 does not require n itself to have distinct prime multiplicities.
A000005 counts divisors.
A007425 counts divisors of divisors.
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 nonempty strict chains of divisors.
A327498 gives the maximum divisor with distinct prime multiplicities.
A336422 counts divisible pairs of divisors, both in A130091.
A336424 counts factorizations using A130091.
A336500 counts divisors of n in A130091 with quotient also in A130091.
A337256 counts strict chains of divisors.

Programs

  • Mathematica
    strchns[n_]:=If[n==1,1,If[!UnsameQ@@Last/@FactorInteger[n],0,Sum[strchns[d],{d,Select[Most[Divisors[n]],UnsameQ@@Last/@FactorInteger[#]&]}]]];
    Table[strchns[n],{n,100}]

A136565 a(n) = sum of the distinct values making up the exponents in the prime-factorization of n.

Original entry on oeis.org

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

Views

Author

Leroy Quet, Jan 07 2008

Keywords

Comments

The sums of the first 10^k terms, for k = 1, 2, ..., are 13, 192, 2089, 21405, 215730, 2162136, 21636277, 216410510, 2164253043, 21642998932, ... . Apparently, the asymptotic mean of this sequence is 2.164... . - Amiram Eldar, Jun 30 2025

Examples

			120 = 2^3 * 3^1 * 5^1. The exponents of the prime factorization are therefore 3,1,1. The distinct values which equal these exponents are 1 and 3. So a(120) = 1+3 = 4.
		

Crossrefs

Programs

  • Mathematica
    Join[{0},Table[Total[Union[Transpose[FactorInteger[n]][[2]]]],{n,2,110}]] (* Harvey P. Dale, Jun 23 2013 *)
  • PARI
    A136565(n) = vecsum(apply(primepi,factor(factorback(apply(e->prime(e),(factor(n)[,2]))))[,1])); \\ Antti Karttunen, Sep 06 2018

Formula

a(n) = A088529(n) = A181591(n) for n: 2 <= n < 24. - Reinhard Zumkeller, Nov 01 2010
a(n) = A066328(A181819(n)). - Antti Karttunen, Sep 06 2018

Extensions

More terms from Diana L. Mecum, Jul 17 2008

A336420 Irregular triangle read by rows where T(n,k) is the number of divisors of the n-th superprimorial A006939(n) with distinct prime multiplicities and k prime factors counted with multiplicity.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 1, 3, 2, 5, 2, 1, 1, 1, 4, 3, 11, 7, 7, 10, 5, 2, 1, 1, 1, 5, 4, 19, 14, 18, 37, 25, 23, 15, 23, 10, 5, 2, 1, 1, 1, 6, 5, 29, 23, 33, 87, 70, 78, 74, 129, 84, 81, 49, 39, 47, 23, 10, 5, 2, 1, 1, 1, 7, 6, 41, 34, 52, 165, 144, 183, 196, 424, 317, 376, 325, 299, 431, 304, 261, 172, 129, 81, 103, 47, 23, 10, 5, 2, 1, 1
Offset: 0

Views

Author

Gus Wiseman, Jul 25 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.
The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).
T(n,k) is also the number of length-n vectors 0 <= v_i <= i summing to k whose nonzero values are all distinct.

Examples

			Triangle begins:
  1
  1  1
  1  2  1  1
  1  3  2  5  2  1  1
  1  4  3 11  7  7 10  5  2  1  1
  1  5  4 19 14 18 37 25 23 15 23 10  5  2  1  1
The divisors counted in row n = 4 are:
  1  2  4     8   16   48   144   432  2160  10800  75600
     3  9    12   24   72   360   720  3024
     5  25   18   40   80   400  1008
     7       20   54  108   504  1200
             27   56  112   540  2800
             28  135  200   600
             45  189  675   756
             50            1350
             63            1400
             75            4725
            175
		

Crossrefs

A000110 gives row sums.
A000124 gives row lengths.
A000142 counts divisors of superprimorials.
A006939 lists superprimorials or Chernoff numbers.
A008278 is the version counting only distinct prime factors.
A008302 counts divisors of superprimorials by bigomega.
A022915 counts permutations of prime indices of superprimorials.
A076954 can be used instead of A006939.
A130091 lists numbers with distinct prime multiplicities.
A146291 counts divisors by bigomega.
A181796 counts divisors with distinct prime multiplicities.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.
A336498 counts divisors of factorials by bigomega.
A336499 uses factorials instead superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    Table[Length[Select[Divisors[chern[n]],PrimeOmega[#]==k&&UnsameQ@@Last/@FactorInteger[#]&]],{n,0,5},{k,0,n*(n+1)/2}]
Previous Showing 31-40 of 155 results. Next