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 51-60 of 155 results. Next

A325276 Irregular triangle read by rows where row n is the omega-sequence of n!.

Original entry on oeis.org

1, 2, 2, 1, 4, 2, 2, 1, 5, 3, 2, 2, 1, 7, 3, 3, 1, 8, 4, 3, 2, 2, 1, 11, 4, 3, 2, 2, 1, 13, 4, 3, 2, 2, 1, 15, 4, 4, 1, 16, 5, 4, 2, 2, 1, 19, 5, 4, 2, 2, 1, 20, 6, 4, 2, 2, 1, 22, 6, 4, 2, 1, 24, 6, 5, 2, 2, 1, 28, 6, 5, 2, 2, 1, 29, 7, 5, 2, 2, 1
Offset: 0

Views

Author

Gus Wiseman, Apr 18 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).

Examples

			Triangle begins:
  {}
  {}
   1
   2  2  1
   4  2  2  1
   5  3  2  2  1
   7  3  3  1
   8  4  3  2  2  1
  11  4  3  2  2  1
  13  4  3  2  2  1
  15  4  4  1
  16  5  4  2  2  1
  19  5  4  2  2  1
  20  6  4  2  2  1
  22  6  4  2  1
  24  6  5  2  2  1
  28  6  5  2  2  1
  29  7  5  2  2  1
  32  7  5  2  2  1
  33  8  5  2  2  1
  36  8  5  2  2  1
  38  8  5  2  2  1
  40  8  6  2  2  1
  41  9  6  2  2  1
  45  9  6  2  2  1
  47  9  6  2  2  1
  49  9  6  3  2  2  1
  52  9  6  3  2  2  1
  55  9  6  3  2  2  1
  56 10  6  3  2  2  1
  59 10  6  3  2  2  1
		

Crossrefs

Row lengths are A325272. Row sums are A325274. Row n is row A325275(n) of A112798. Second-to-last column is A325273. Column k = 1 is A022559. Column k = 2 is A000720. Column k = 3 is A071626.
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), A325249 (sum).

Programs

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

A336419 Number of divisors d of the n-th superprimorial A006939(n) with distinct prime exponents such that the quotient A006939(n)/d also has distinct prime exponents.

Original entry on oeis.org

1, 2, 4, 10, 24, 64, 184, 536, 1608, 5104, 16448, 55136, 187136, 658624, 2339648, 8618208, 31884640, 121733120, 468209408, 1849540416, 7342849216
Offset: 0

Views

Author

Gus Wiseman, Jul 25 2020

Keywords

Comments

A number has distinct prime exponents iff its prime signature is strict.
The n-th superprimorial or Chernoff number is A006939(n) = Product_{i = 1..n} prime(i)^(n - i + 1).

Examples

			The a(0) = 1 through a(3) = 10 divisors:
  1  2  12  360
-----------------
  1  1   1    1
     2   3    5
         4    8
        12    9
             18
             20
             40
             45
             72
            360
		

Crossrefs

A000110 shifted once to the left dominates this sequence.
A006939 lists superprimorials or Chernoff numbers.
A022915 counts permutations of prime indices of superprimorials.
A130091 lists numbers with distinct prime exponents.
A181796 counts divisors with distinct prime exponents.
A181818 gives products of superprimorials.
A317829 counts factorizations of superprimorials.
A336417 counts perfect-power divisors of superprimorials.

Programs

  • Mathematica
    chern[n_]:=Product[Prime[i]^(n-i+1),{i,n}];
    Table[Length[Select[Divisors[chern[n]],UnsameQ@@Last/@FactorInteger[#]&&UnsameQ@@Last/@FactorInteger[chern[n]/#]&]],{n,0,6}]
  • PARI
    recurse(n,k,b,d)={if(k>n, 1, sum(i=0, k, if((i==0||!bittest(b,i)) && (i==k||!bittest(d,k-i)), self()(n, k+1, bitor(b, 1<Andrew Howroyd, Aug 30 2020

Extensions

a(10)-a(20) from Andrew Howroyd, Aug 31 2020

A353507 Product of multiplicities of the prime exponents (signature) of n; a(1) = 0.

Original entry on oeis.org

0, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 3, 1, 1, 2, 2, 2, 2, 1, 2, 2, 1, 1, 3, 1, 1, 1, 2, 1, 1, 1, 1, 2, 1, 1, 1, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 3, 1, 1, 2, 3, 1, 1, 1, 2, 1, 1, 2, 3, 1, 1, 1, 2, 1, 2, 2, 2, 2
Offset: 1

Views

Author

Gus Wiseman, May 19 2022

Keywords

Comments

Warning: If the prime multiplicities of n are a multiset y, this sequence gives the product of multiplicities in y, not the product of y.
Differs from A351946 at A351946(1260) = 4, a(1260) = 2.
Differs from A327500 at A327500(450) = 3, a(450) = 2.
We set a(1) = 0 so that the positions of first appearances are the primorials A002110.
Also the product of the prime metasignature of n (row n of A238747).

Examples

			The prime signature of 13860 is (2,2,1,1,1), with multiplicities (2,3), so a(13860) = 6.
		

Crossrefs

Positions of first appearances are A002110.
The prime indices themselves have product A003963, counted by A339095.
The prime signature itself has product A005361, counted by A266477.
A001222 counts prime factors with multiplicity, distinct A001221.
A056239 adds up prime indices, row sums of A112798 and A296150.
A071625 counts distinct prime exponents (third omega).
A124010 gives prime signature, sorted A118914.
A130091 lists numbers with distinct prime exponents, counted by A098859.
A181819 gives prime shadow, with an inverse A181821.
A238747 gives prime metasignature, sorted A353742.
A323022 gives fourth omega.

Programs

  • Maple
    f:= proc(n) local M,s;
      M:= ifactors(n)[2][..,2];
      mul(numboccur(s,M),s=convert(M,set));
    end proc:
    f(1):= 0:
    map(f, [$1..100]); # Robert Israel, May 19 2023
  • Mathematica
    Table[If[n==1,0,Times@@Length/@Split[Sort[Last/@FactorInteger[n]]]],{n,100}]
    Join[{0},Table[Times@@(Length/@Split[FactorInteger[n][[;;,2]]]),{n,2,100}]] (* Harvey P. Dale, Oct 20 2024 *)
  • Python
    from math import prod
    from itertools import groupby
    from sympy import factorint
    def A353507(n): return 0 if n == 1 else prod(len(list(g)) for k, g in groupby(factorint(n).values())) # Chai Wah Wu, May 20 2022

Formula

A353836 Triangle read by rows where T(n,k) is the number of integer partitions of n with k distinct run-sums.

Original entry on oeis.org

1, 0, 1, 0, 2, 0, 0, 2, 1, 0, 0, 4, 1, 0, 0, 0, 2, 5, 0, 0, 0, 0, 5, 5, 1, 0, 0, 0, 0, 2, 12, 1, 0, 0, 0, 0, 0, 7, 12, 3, 0, 0, 0, 0, 0, 0, 3, 19, 8, 0, 0, 0, 0, 0, 0, 0, 5, 27, 9, 1, 0, 0, 0, 0, 0, 0, 0, 2, 33, 20, 1, 0, 0, 0, 0, 0, 0, 0
Offset: 0

Views

Author

Gus Wiseman, May 26 2022

Keywords

Comments

The run-sums of a sequence are the sums of its maximal consecutive constant subsequences (runs). For example, the run-sums of (2,2,1,1,1,3,2,2) are (4,3,3,4).

Examples

			Triangle begins:
  1
  0  1
  0  2  0
  0  2  1  0
  0  4  1  0  0
  0  2  5  0  0  0
  0  5  5  1  0  0  0
  0  2 12  1  0  0  0  0
  0  7 12  3  0  0  0  0  0
  0  3 19  8  0  0  0  0  0  0
  0  5 27  9  1  0  0  0  0  0  0
  0  2 33 20  1  0  0  0  0  0  0  0
  0 13 28 34  2  0  0  0  0  0  0  0  0
  0  2 48 46  5  0  0  0  0  0  0  0  0  0
  0  5 65 51 14  0  0  0  0  0  0  0  0  0  0
  0  4 57 99 15  1  0  0  0  0  0  0  0  0  0  0
For example, row n = 8 counts the following partitions:
  (8)         (53)       (431)
  (44)        (62)       (521)
  (422)       (71)       (3221)
  (2222)      (332)
  (41111)     (611)
  (221111)    (3311)
  (11111111)  (4211)
              (5111)
              (22211)
              (32111)
              (311111)
              (2111111)
		

Crossrefs

Row sums are A000041.
Counting distinct parts instead of run-sums gives A116608.
Column k = 1 is A304442, ranked by A353833 (nonprime A353834).
The rank statistic is A353835, weak A353861, for compositions A353849.
A275870 counts collapsible partitions, ranked by A300273.
A351014 counts distinct runs in standard compositions.
A353832 represents the operation of taking run-sums of a partition.
A353837 counts partitions with all distinct run-sums, ranked by A353838.
A353840-A353846 pertain to partition run-sum trajectory.
A353864 counts rucksack partitions, ranked by A353866.
A353865 counts perfect rucksack partitions, ranked by A353867.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n], Length[Union[Total/@Split[#]]]==k&]],{n,0,15},{k,0,n}]

A367579 Irregular triangle read by rows where row n is the multiset multiplicity kernel (MMK) of the multiset of prime indices of n.

Original entry on oeis.org

1, 2, 1, 3, 1, 1, 4, 1, 2, 1, 1, 5, 1, 2, 6, 1, 1, 2, 2, 1, 7, 1, 2, 8, 1, 3, 2, 2, 1, 1, 9, 1, 2, 3, 1, 1, 2, 1, 4, 10, 1, 1, 1, 11, 1, 2, 2, 1, 1, 3, 3, 1, 1, 12, 1, 1, 2, 2, 1, 3, 13, 1, 1, 1, 14, 1, 5, 2, 3, 1, 1, 15, 1, 2, 4, 1, 3, 2, 2, 1, 6, 16, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Nov 25 2023

Keywords

Comments

Row n = 1 is empty.
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.
We define the multiset multiplicity kernel MMK(m) of a multiset m by the following property, holding for all distinct multiplicities k >= 1. If S is the set of elements of multiplicity k in m, then min(S) has multiplicity |S| in MMK(m). For example, MMK({1,1,2,2,3,4,5}) = {1,1,3,3,3}, and MMK({1,2,3,4,5,5,5,5}) = {1,1,1,1,5}.
Note: I chose the word 'kernel' because, as with A007947 and A304038, MMK(m) is constructed using the same underlying elements as m and has length equal to the number of distinct elements of m. However, it is not necessarily a submultiset of m.

Examples

			The first 45 rows:
     1: {}      16: {1}       31: {11}
     2: {1}     17: {7}       32: {1}
     3: {2}     18: {1,2}     33: {2,2}
     4: {1}     19: {8}       34: {1,1}
     5: {3}     20: {1,3}     35: {3,3}
     6: {1,1}   21: {2,2}     36: {1,1}
     7: {4}     22: {1,1}     37: {12}
     8: {1}     23: {9}       38: {1,1}
     9: {2}     24: {1,2}     39: {2,2}
    10: {1,1}   25: {3}       40: {1,3}
    11: {5}     26: {1,1}     41: {13}
    12: {1,2}   27: {2}       42: {1,1,1}
    13: {6}     28: {1,4}     43: {14}
    14: {1,1}   29: {10}      44: {1,5}
    15: {2,2}   30: {1,1,1}   45: {2,3}
		

Crossrefs

Indices of empty and singleton rows are A000961.
Row lengths are A001221.
Depends only on rootless base A052410, see A007916.
Row minima are A055396.
Rows have A071625 distinct elements.
Indices of constant rows are A072774.
Indices of strict rows are A130091.
Rows have Heinz numbers A367580.
Row sums are A367581.
Row maxima are A367583, opposite A367587.
Index of first row with Heinz number n is A367584.
Sorted row indices of first appearances are A367585.
Indices of rows of the form {1,1,...} are A367586.
Agrees with sorted prime signature at A367683, counted by A367682.
A submultiset of prime indices at A367685, counted by A367684.
A007947 gives squarefree kernel.
A112798 lists prime indices, length A001222, sum A056239, reverse A296150.
A124010 lists prime multiplicities (prime signature), sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A238747 gives prime metasignature, reversed A353742.
A304038 lists distinct prime indices, length A001221, sum A066328.
A367582 counts partitions by sum of multiset multiplicity kernel.

Programs

  • Mathematica
    mmk[q_]:=With[{mts=Length/@Split[q]}, Sort[Table[Min@@Select[q,Count[q,#]==i&], {i,mts}]]];
    Table[mmk[PrimePi/@Join@@ConstantArray@@@If[n==1, {},FactorInteger[n]]], {n,100}]

Formula

For all positive integers n and k, row n^k is the same as row n.

A353841 Length of the trajectory of the partition run-sum transformation of n, using Heinz numbers; a(1) = 0.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, May 25 2022

Keywords

Comments

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). This gives a bijective correspondence between positive integers and integer partitions.
Starting with n, this is one plus the number of times one must apply A353832 to reach a squarefree number.
Also Kimberling's depth statistic (defined in A237685 and A237750) plus one.

Examples

			The trajectory for a(1080) = 4 is the following, with prime indices shown on the right:
  1080: {1,1,1,2,2,2,3}
   325: {3,3,6}
   169: {6,6}
    37: {12}
The trajectory for a(87780) = 5 is the following, with prime indices shown on the right:
  87780: {1,1,2,3,4,5,8}
  65835: {2,2,3,4,5,8}
  51205: {3,4,4,5,8}
  19855: {3,5,8,8}
   2915: {3,5,16}
The trajectory for a(39960) = 5 is the following, with prime indices shown on the right:
  39960: {1,1,1,2,2,2,3,12}
  12025: {3,3,6,12}
   6253: {6,6,12}
   1369: {12,12}
     89: {24}
		

Crossrefs

Positions of 1's are A005117.
The version for run-lengths instead of sums is A182850 or A323014.
Positions of first appearances are A353743.
These are the row-lengths of A353840.
Other sequences pertaining to this trajectory are A353842-A353845.
Counting partitions by this statistic gives A353846.
The version for compositions is A353854, run-lengths of A353853.
A001222 counts prime factors, distinct A001221.
A005811 counts runs in binary expansion.
A056239 adds up prime indices, row sums of A112798 and A296150.
A300273 ranks collapsible partitions, counted by A275870.
A318928 gives runs-resistance of binary expansion.
A353832 represents the operation of taking run-sums of a partition.
A353833 ranks partitions with all equal run-sums, counted by A304442.
A353835 counts distinct run-sums of prime indices, weak A353861.
A353838 ranks partitions with all distinct run-sums, counted by A353837.
A353866 ranks rucksack partitions, counted by A353864.

Programs

  • Mathematica
    Table[If[n==1,0,Length[NestWhileList[Times@@Prime/@Cases[If[#==1,{},FactorInteger[#]],{p_,k_}:>PrimePi[p]*k]&,n,!SquareFreeQ[#]&]]],{n,100}]
  • PARI
    pis_to_runs(n) = { my(runs=List([]), f=factor(n)); for(i=1,#f~,while(f[i,2], listput(runs,primepi(f[i,1])); f[i,2]--)); (runs); };
    A353832(n) = if(1==n,n,my(pruns = pis_to_runs(n), m=1, runsum=pruns[1]); for(i=2,#pruns,if(pruns[i] == pruns[i-1], runsum += pruns[i], m *= prime(runsum); runsum = pruns[i])); (m*prime(runsum)));
    A353841(n) = if(1==n,0,for(i=1,oo,if(issquarefree(n), return(i), n = A353832(n)))); \\ Antti Karttunen, Jan 20 2025

Formula

a(1) = 0, and for n > 1, if A008966(n) = 1 [n is in A005117], a(n) = 1, otherwise a(n) = 1+a(A353832(n)). [See comments] - Antti Karttunen, Jan 20 2025

Extensions

More terms from Antti Karttunen, Jan 20 2025

A367581 Sum of the multiset multiplicity kernel (in which each multiplicity becomes the least element of that multiplicity) of the prime indices of n.

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 4, 1, 2, 2, 5, 3, 6, 2, 4, 1, 7, 3, 8, 4, 4, 2, 9, 3, 3, 2, 2, 5, 10, 3, 11, 1, 4, 2, 6, 2, 12, 2, 4, 4, 13, 3, 14, 6, 5, 2, 15, 3, 4, 4, 4, 7, 16, 3, 6, 5, 4, 2, 17, 5, 18, 2, 6, 1, 6, 3, 19, 8, 4, 3, 20, 3, 21, 2, 5, 9, 8, 3, 22, 4, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Nov 28 2023

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.
We define the multiset multiplicity kernel MMK(m) of a multiset m by the following property, holding for all distinct multiplicities k >= 1. If S is the set of elements of multiplicity k in m, then min(S) has multiplicity |S| in MMK(m). For example, MMK({1,1,2,2,3,4,5}) = {1,1,3,3,3}, and MMK({1,2,3,4,5,5,5,5}) = {1,1,1,1,5}. As an operation on multisets, MMK is represented by A367579, and as an operation on their Heinz numbers, it is represented by A367580.

Examples

			The multiset multiplicity kernel of {1,2,2,3} is {1,1,2}, so a(90) = 4.
		

Crossrefs

Positions of 1's are A000079 without 1.
Positions of first appearances are A008578.
Depends only on rootless base A052410, see A007916, A052409.
The triangle A367579 has these as row sums, ranks A367580.
The triangle for this rank statistic is A367582.
For maximum instead of sum we have A367583, opposite A367587.
A007947 gives squarefree kernel.
A112798 lists prime indices, length A001222, sum A056239, reverse A296150.
A124010 gives prime signature, sorted A118914.
A181819 gives prime shadow, with an inverse A181821.
A238747 gives prime metasignature, reverse A353742.
A304038 lists distinct prime indices, length A001221, sum A066328.

Programs

  • Mathematica
    mmk[q_]:=With[{mts=Length/@Split[q]}, Sort[Table[Min@@Select[q,Count[q,#]==i&], {i,mts}]]];
    Table[Total[mmk[PrimePi/@Join@@ConstantArray@@@FactorInteger[n]]], {n,100}]

Formula

a(n^k) = a(n) for all positive integers n and k.
a(n) = A056239(A367580(n)).
If n is squarefree, a(n) = A055396(n)*A001222(n).

A325278 Smallest number with adjusted frequency depth n.

Original entry on oeis.org

1, 2, 4, 6, 12, 60, 2520, 1286485200, 35933692027611398678865941374040400000
Offset: 0

Views

Author

Gus Wiseman, Apr 17 2019

Keywords

Comments

The adjusted frequency depth of a positive integer n is 0 if n = 1, and otherwise it is 1 plus the number of times one must apply A181819 to reach a prime number, where A181819(k = p^i*...*q^j) = prime(i)*...*prime(j) = product of primes indexed by the prime exponents of k. For example, 180 has adjusted frequency depth 5 because we have: 180 -> 18 -> 6 -> 4 -> 3.
Differs from A182857 in having 2 instead of 3.

Crossrefs

A subsequence of 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
    nn=10000;
    fd[n_]:=Switch[n,1,0,?PrimeQ,1,,1+fd[Times@@Prime/@Last/@FactorInteger[n]]];
    fds=fd/@Range[nn];
    Sort[Table[Position[fds,x][[1,1]],{x,Union[fds]}]]

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}]
Previous Showing 51-60 of 155 results. Next