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 21-30 of 75 results. Next

A096825 Maximal size of an antichain in divisor lattice D(n).

Original entry on oeis.org

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

Views

Author

Yuval Dekel (dekelyuval(AT)hotmail.com) and Vladeta Jovovic, Aug 17 2004

Keywords

Comments

The divisor lattice D(n) is the lattice of the divisors of the natural number n.
Also the number of divisors of n with half (rounded either way) as many prime factors (counting multiplicity) as n. - Gus Wiseman, Aug 24 2018

Examples

			There are two maximal size antichains of divisors of 180, namely {12, 18, 20, 30, 45} and {4, 6, 9, 10, 15}. Both have length 5 so a(180) = 5. - _Gus Wiseman_, Aug 24 2018
		

Crossrefs

Programs

  • Maple
    a:=proc(n) local klist,x; klist:=ifactors(n)[2,1..-1,2]; coeff(normal(mul((1-x^(k+1))/(1-x),k=klist)),x,floor(add(k,k=klist)/2)) end: seq(a(n), n=1..100);
  • Mathematica
    a[n_] := Module[{pp, kk, x}, {pp, kk} = Transpose[FactorInteger[n]]; Coefficient[ Product[ Total[x^Range[0, k]], {k, kk}], x, Quotient[ Total[ kk], 2] ] ]; Array[a, 100] (* Jean-François Alcover, Nov 20 2017 *)
    Table[Length[Select[Divisors[n],PrimeOmega[#]==Round[PrimeOmega[n]/2]&]],{n,50}] (* Gus Wiseman, Aug 24 2018 *)
  • PARI
    a(n)=if(n<6||isprimepower(n), return(1)); my(d=divisors(n),r=1,u); d=d[2..#d-1];for(k=0,2^#d-1,if(hammingweight(k)<=r,next); u=vecextract(d,k); for(i=1,#u, for(j=i+1,#u, if(u[j]%u[i]==0, next(3))));r=#u);r \\ Charles R Greathouse IV, May 14 2013
    
  • Python
    from sympy import factorint
    from sympy.utilities.iterables import multiset_combinations
    def A096825(n):
        fs = factorint(n)
        return len(list(multiset_combinations(fs,sum(fs.values())//2))) # Chai Wah Wu, Aug 23 2021
  • Sage
    def A096825(n) :
        if n==1 : return 1
        R. = QQ[]; mults = [x[1] for x in factor(n)]
        return prod((t^(m+1)-1)//(t-1) for m in mults)[sum(mults)//2]
    # Eric M. Schmidt, May 11 2013
    

Formula

a(n) is the coefficient at x^k in (1+x+...+x^k_1)*...*(1+x+...+x^k_q) where n=p_1^k_1*...*p_q^k_q is the prime factorization of n and k=floor((k_1+...+k_q)/2). - Alec Mihailovs (alec(AT)mihailovs.com), Aug 22 2004

Extensions

More terms from Alec Mihailovs (alec(AT)mihailovs.com), Aug 22 2004

A342193 Numbers with no prime index dividing all the other prime indices.

Original entry on oeis.org

1, 15, 33, 35, 45, 51, 55, 69, 75, 77, 85, 91, 93, 95, 99, 105, 119, 123, 135, 141, 143, 145, 153, 155, 161, 165, 175, 177, 187, 195, 201, 203, 205, 207, 209, 215, 217, 219, 221, 225, 231, 245, 247, 249, 253, 255, 265, 275, 279, 285, 287, 291, 295, 297, 299
Offset: 1

Views

Author

Gus Wiseman, Apr 11 2021

Keywords

Comments

Alternative name: 1 and numbers with smallest prime index not dividing all the other prime indices.
First differs from A339562 in having 45.
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.
Also 1 and Heinz numbers of integer partitions with smallest part not dividing all the others (counted by A338470). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}         105: {2,3,4}      201: {2,19}
     15: {2,3}      119: {4,7}        203: {4,10}
     33: {2,5}      123: {2,13}       205: {3,13}
     35: {3,4}      135: {2,2,2,3}    207: {2,2,9}
     45: {2,2,3}    141: {2,15}       209: {5,8}
     51: {2,7}      143: {5,6}        215: {3,14}
     55: {3,5}      145: {3,10}       217: {4,11}
     69: {2,9}      153: {2,2,7}      219: {2,21}
     75: {2,3,3}    155: {3,11}       221: {6,7}
     77: {4,5}      161: {4,9}        225: {2,2,3,3}
     85: {3,7}      165: {2,3,5}      231: {2,4,5}
     91: {4,6}      175: {3,3,4}      245: {3,4,4}
     93: {2,11}     177: {2,17}       247: {6,8}
     95: {3,8}      187: {5,7}        249: {2,23}
     99: {2,2,5}    195: {2,3,6}      253: {5,9}
		

Crossrefs

The complement is counted by A083710 (strict: A097986).
The complement with no 1's is A083711 (strict: A098965).
These partitions are counted by A338470 (strict: A341450).
The squarefree case is A339562, with squarefree complement A339563.
The case with maximum prime index not divisible by all others is A343338.
The case with maximum prime index divisible by all others is A343339.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A001221 counts distinct prime factors.
A006128 counts partitions with a selected position (strict: A015723).
A056239 adds up prime indices, row sums of A112798.
A299702 lists Heinz numbers of knapsack partitions.
A339564 counts factorizations with a selected factor.

Programs

  • Mathematica
    Select[Range[100],#==1||With[{p=PrimePi/@First/@FactorInteger[#]},!And@@IntegerQ/@(p/Min@@p)]&]

A343337 Numbers with no prime index divisible by all the other prime indices.

Original entry on oeis.org

1, 15, 30, 33, 35, 45, 51, 55, 60, 66, 69, 70, 75, 77, 85, 90, 91, 93, 95, 99, 102, 105, 110, 119, 120, 123, 132, 135, 138, 140, 141, 143, 145, 150, 153, 154, 155, 161, 165, 170, 175, 177, 180, 182, 186, 187, 190, 198, 201, 203, 204, 205, 207, 209, 210, 215
Offset: 1

Views

Author

Gus Wiseman, Apr 13 2021

Keywords

Comments

Alternative name: 1 and numbers whose greatest prime index is not divisible by all the other prime indices.
First differs from A318992 in lacking 195.
First differs from A343343 in lacking 195.
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.
Also Heinz numbers of partitions with greatest part not divisible by all the others (counted by A343341). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}            90: {1,2,2,3}      141: {2,15}
     15: {2,3}         91: {4,6}          143: {5,6}
     30: {1,2,3}       93: {2,11}         145: {3,10}
     33: {2,5}         95: {3,8}          150: {1,2,3,3}
     35: {3,4}         99: {2,2,5}        153: {2,2,7}
     45: {2,2,3}      102: {1,2,7}        154: {1,4,5}
     51: {2,7}        105: {2,3,4}        155: {3,11}
     55: {3,5}        110: {1,3,5}        161: {4,9}
     60: {1,1,2,3}    119: {4,7}          165: {2,3,5}
     66: {1,2,5}      120: {1,1,1,2,3}    170: {1,3,7}
     69: {2,9}        123: {2,13}         175: {3,3,4}
     70: {1,3,4}      132: {1,1,2,5}      177: {2,17}
     75: {2,3,3}      135: {2,2,2,3}      180: {1,1,2,2,3}
     77: {4,5}        138: {1,2,9}        182: {1,4,6}
     85: {3,7}        140: {1,1,3,4}      186: {1,2,11}
For example, 195 has prime indices {2,3,6}, and 6 is divisible by both 2 and 3, so 195 does not belong to the sequence.
		

Crossrefs

The complement is counted by A130689.
The dual version is A342193.
The case with smallest prime index not dividing all the others is A343338.
The case with smallest prime index dividing by all the others is A343340.
These are the Heinz numbers of the partitions counted by A343341.
Including the dual version gives A343343.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A006128 counts partitions with a selected position.
A056239 adds up prime indices, row sums of A112798.
A067824 counts strict chains of divisors starting with n.
A253249 counts strict chains of divisors.
A339564 counts factorizations with a selected factor.

Programs

  • Mathematica
    Select[Range[1000],#==1||With[{p=PrimePi/@First/@FactorInteger[#]},!And@@IntegerQ/@(Max@@p/p)]&]

A339563 Squarefree numbers > 1 whose smallest prime index divides all the other prime indices.

Original entry on oeis.org

2, 3, 5, 6, 7, 10, 11, 13, 14, 17, 19, 21, 22, 23, 26, 29, 30, 31, 34, 37, 38, 39, 41, 42, 43, 46, 47, 53, 57, 58, 59, 61, 62, 65, 66, 67, 70, 71, 73, 74, 78, 79, 82, 83, 86, 87, 89, 94, 97, 101, 102, 103, 106, 107, 109, 110, 111, 113, 114, 115, 118, 122, 127
Offset: 1

Views

Author

Gus Wiseman, Apr 10 2021

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.
Also Heinz numbers of strict integer partitions whose smallest part divides all the others (counted by A097986). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
      2: {1}       29: {10}        59: {17}
      3: {2}       30: {1,2,3}     61: {18}
      5: {3}       31: {11}        62: {1,11}
      6: {1,2}     34: {1,7}       65: {3,6}
      7: {4}       37: {12}        66: {1,2,5}
     10: {1,3}     38: {1,8}       67: {19}
     11: {5}       39: {2,6}       70: {1,3,4}
     13: {6}       41: {13}        71: {20}
     14: {1,4}     42: {1,2,4}     73: {21}
     17: {7}       43: {14}        74: {1,12}
     19: {8}       46: {1,9}       78: {1,2,6}
     21: {2,4}     47: {15}        79: {22}
     22: {1,5}     53: {16}        82: {1,13}
     23: {9}       57: {2,8}       83: {23}
     26: {1,6}     58: {1,10}      86: {1,14}
		

Crossrefs

These partitions are counted by A097986 (non-strict: A083710).
The case with no 1's is counted by A098965 (non-strict: A083711).
The squarefree complement is A339562, ranked by A341450.
The complement of the not necessarily squarefree version is A342193.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A001055 counts factorizations.
A001221 counts distinct prime factors.
A005117 lists squarefree numbers.
A006128 counts partitions with a selected position (strict: A015723).
A056239 adds up prime indices, row sums of A112798.
A338470 counts partitions with no dividing part.

Programs

  • Mathematica
    Select[Range[2,100],SquareFreeQ[#]&&With[{p=PrimePi/@First/@FactorInteger[#]},And@@IntegerQ/@(p/Min@@p)]&]

A337135 a(1) = 1; for n > 1, a(n) = Sum_{d|n, d <= sqrt(n)} a(d).

Original entry on oeis.org

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

Views

Author

Ilya Gutkovskiy, Nov 21 2020

Keywords

Comments

From Gus Wiseman, Mar 05 2021: (Start)
This sequence counts all of the following essentially equivalent things:
1. Chains of distinct inferior divisors from n to 1, where a divisor d|n is inferior if d <= n/d. Inferior divisors are counted by A038548 and listed by A161906.
2. Chains of divisors from n to 1 whose first-quotients (in analogy with first-differences) are term-wise greater than or equal to their decapitation (maximum element removed). For example, the divisor chain q = 60/4/2/1 has first-quotients (15,2,2), which are >= (4,2,1), so q is counted under a(60).
3. Chains of divisors from n to 1 such that x >= y^2 for all adjacent x, y.
4. Factorizations of n where each factor is greater than or equal to the product of all previous factors.
(End)

Examples

			From _Gus Wiseman_, Mar 05 2021: (Start)
The a(n) chains for n = 1, 2, 4, 12, 16, 24, 36, 60:
  1  2/1  4/1    12/1    16/1      24/1      36/1      60/1
          4/2/1  12/2/1  16/2/1    24/2/1    36/2/1    60/2/1
                 12/3/1  16/4/1    24/3/1    36/3/1    60/3/1
                         16/4/2/1  24/4/1    36/4/1    60/4/1
                                   24/4/2/1  36/6/1    60/5/1
                                             36/4/2/1  60/6/1
                                             36/6/2/1  60/4/2/1
                                                       60/6/2/1
The a(n) factorizations for n = 2, 4, 12, 16, 24, 36, 60:
    2  4    12   16     24     36     60
       2*2  2*6  2*8    3*8    4*9    2*30
            3*4  4*4    4*6    6*6    3*20
                 2*2*4  2*12   2*18   4*15
                        2*2*6  3*12   5*12
                               2*2*9  6*10
                               2*3*6  2*2*15
                                      2*3*10
(End)
		

Crossrefs

Cf. A002033, A008578 (positions of 1's), A068108.
The restriction to powers of 2 is A018819.
Not requiring inferiority gives A074206 (ordered factorizations).
The strictly inferior version is A342083.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A000929, or A342098 forbidding equality.
A000005 counts divisors, with sum A000203.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A167865 counts strict chains of divisors > 1 summing to n.
A207375 lists central divisors.
A253249 counts strict chains of divisors.
A334996 counts ordered factorizations by product and length.
A334997 counts chains of divisors of n by length.
A342086 counts strict factorizations of divisors.
- Inferior: A033676, A066839, A072499, A161906.
- Superior: A033677, A070038, A161908.
- Strictly Inferior: A060775, A070039, A333806, A341674.
- Strictly Superior: A048098, A064052, A140271, A238535, A341673.

Programs

  • Maple
    a:= proc(n) option remember; `if`(n=1, 1, add(
          `if`(d<=n/d, a(d), 0), d=numtheory[divisors](n)))
        end:
    seq(a(n), n=1..128);  # Alois P. Heinz, Jun 24 2021
  • Mathematica
    a[1] = 1; a[n_] := a[n] = DivisorSum[n, a[#] &, # <= Sqrt[n] &]; Table[a[n], {n, 95}]
    (* second program *)
    asc[n_]:=Prepend[#,n]&/@Prepend[Join@@Table[asc[d],{d,Select[Divisors[n],#Gus Wiseman, Mar 05 2021 *)

Formula

G.f.: Sum_{k>=1} a(k) * x^(k^2) / (1 - x^k).
a(2^n) = A018819(n). - Gus Wiseman, Mar 08 2021

A339562 Squarefree numbers with no prime index dividing all the other prime indices.

Original entry on oeis.org

1, 15, 33, 35, 51, 55, 69, 77, 85, 91, 93, 95, 105, 119, 123, 141, 143, 145, 155, 161, 165, 177, 187, 195, 201, 203, 205, 209, 215, 217, 219, 221, 231, 247, 249, 253, 255, 265, 285, 287, 291, 295, 299, 301, 309, 323, 327, 329, 335, 341, 345, 355, 357, 377, 381
Offset: 1

Views

Author

Gus Wiseman, Apr 10 2021

Keywords

Comments

First differs from A342193 in lacking 45.
Alternative name: 1 and squarefree numbers with smallest prime index not dividing all the other prime indices.
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.
Also 1 and Heinz numbers of strict integer partitions with smallest part not dividing all the others (counted by A341450). The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), giving a bijective correspondence between positive integers and integer partitions.

Examples

			The sequence of terms together with their prime indices begins:
      1: {}         141: {2,15}     219: {2,21}
     15: {2,3}      143: {5,6}      221: {6,7}
     33: {2,5}      145: {3,10}     231: {2,4,5}
     35: {3,4}      155: {3,11}     247: {6,8}
     51: {2,7}      161: {4,9}      249: {2,23}
     55: {3,5}      165: {2,3,5}    253: {5,9}
     69: {2,9}      177: {2,17}     255: {2,3,7}
     77: {4,5}      187: {5,7}      265: {3,16}
     85: {3,7}      195: {2,3,6}    285: {2,3,8}
     91: {4,6}      201: {2,19}     287: {4,13}
     93: {2,11}     203: {4,10}     291: {2,25}
     95: {3,8}      205: {3,13}     295: {3,17}
    105: {2,3,4}    209: {5,8}      299: {6,9}
    119: {4,7}      215: {3,14}     301: {4,14}
    123: {2,13}     217: {4,11}     309: {2,27}
		

Crossrefs

The squarefree complement is A339563.
These partitions are counted by A341450.
The not necessarily squarefree version is A342193.
A000005 counts divisors.
A000070 counts partitions with a selected part.
A001221 counts distinct prime factors.
A005117 lists squarefree numbers.
A006128 counts partitions with a selected position (strict: A015723).
A056239 adds up prime indices (row sums of A112798).
A083710 counts partitions with a dividing part (strict: A097986).

Programs

  • Mathematica
    Select[Range[100],#==1||SquareFreeQ[#]&&With[{p=PrimePi/@First/@FactorInteger[#]},!And@@IntegerQ/@(p/Min@@p)]&]

A342087 Number of chains of divisors starting with n and having no adjacent parts x <= y^2.

Original entry on oeis.org

1, 2, 2, 2, 2, 4, 2, 4, 2, 4, 2, 6, 2, 4, 4, 4, 2, 6, 2, 6, 4, 4, 2, 8, 2, 4, 4, 6, 2, 8, 2, 6, 4, 4, 4, 8, 2, 4, 4, 8, 2, 10, 2, 6, 6, 4, 2, 12, 2, 6, 4, 6, 2, 10, 4, 8, 4, 4, 2, 14, 2, 4, 6, 6, 4, 10, 2, 6, 4, 8, 2, 16, 2, 4, 6, 6, 4, 10, 2, 12, 4, 4, 2, 14
Offset: 1

Views

Author

Gus Wiseman, Mar 05 2021

Keywords

Comments

An alternative wording: Number of chains of divisors starting with n and having all adjacent parts x > y^2.

Examples

			The chains for n = 1, 2, 6, 12, 24, 42, 48:
   1    2      6        12        24        42          48
        2/1    6/1      12/1      24/1      42/1        48/1
               6/2      12/2      24/2      42/2        48/2
               6/2/1    12/3      24/3      42/3        48/3
                        12/2/1    24/4      42/6        48/4
                        12/3/1    24/2/1    42/2/1      48/6
                                  24/3/1    42/3/1      48/2/1
                                  24/4/1    42/6/1      48/3/1
                                            42/6/2      48/4/1
                                            42/6/2/1    48/6/1
                                                        48/6/2
                                                        48/6/2/1
		

Crossrefs

The restriction to powers of 2 is A018819.
Not requiring strict inferiority gives A067824.
The weakly inferior version is twice A337135.
The case ending with 1 is counted by A342083.
The strictly superior version is A342084.
The weakly superior version is A342085.
The additive version is A342098, or A000929 allowing equality.
A000005 counts divisors, with sum A000203.
A001055 counts factorizations.
A003238 counts chains of divisors summing to n-1, with strict case A122651.
A038548 counts inferior (or superior) divisors.
A056924 counts strictly inferior (or strictly superior) divisors.
A067824 counts strict chains of divisors starting with n.
A074206 counts ordered factorizations.
A167865 counts strict chains of divisors > 1 summing to n.
A253249 counts strict chains of divisors.
A334997 counts chains of divisors of n by length.

Programs

  • Mathematica
    cem[n_]:=Prepend[Prepend[#,n]&/@Join@@cem/@Most[Divisors[n]],{n}];
    Table[Length[Select[cem[n],And@@Thread[Divide@@@Partition[#,2,1]>Rest[#]]&]],{n,30}]

Formula

For n > 1, a(n) = 2*A342083(n).

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

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}]

A337256 Number of strict chains of divisors of n.

Original entry on oeis.org

2, 4, 4, 8, 4, 12, 4, 16, 8, 12, 4, 32, 4, 12, 12, 32, 4, 32, 4, 32, 12, 12, 4, 80, 8, 12, 16, 32, 4, 52, 4, 64, 12, 12, 12, 104, 4, 12, 12, 80, 4, 52, 4, 32, 32, 12, 4, 192, 8, 32, 12, 32, 4, 80, 12, 80, 12, 12, 4, 176, 4, 12, 32, 128, 12, 52, 4, 32, 12, 52
Offset: 1

Views

Author

Gus Wiseman, Aug 23 2020

Keywords

Examples

			The a(n) chains for n = 1, 2, 4, 6, 8 (empty chains shown as 0):
  0  0    0      0      0
  1  1    1      1      1
     2    2      2      2
     2/1  4      3      4
          2/1    6      8
          4/1    2/1    2/1
          4/2    3/1    4/1
          4/2/1  6/1    4/2
                 6/2    8/1
                 6/3    8/2
                 6/2/1  8/4
                 6/3/1  4/2/1
                        8/2/1
                        8/4/1
                        8/4/2
                        8/4/2/1
		

Crossrefs

A067824 is the case of chains starting with n (or ending with 1).
A074206 is the case of chains from n to 1.
A253249 is the nonempty case.
A000005 counts divisors.
A001055 counts factorizations.
A001222 counts prime factors with multiplicity.
A074206 counts chains of divisors from n to 1.
A122651 counts chains of divisors summing to n.
A167865 counts chains of divisors > 1 summing to n.
A334996 appears to count chains of divisors from n to 1 by length.
A337070 counts chains of divisors starting with A006939(n).
A337071 counts chains of divisors starting with n!.
A337255 counts chains of divisors starting with n by length.

Programs

  • Mathematica
    stableSets[u_,Q_]:=If[Length[u]==0,{{}},With[{w=First[u]},Join[stableSets[DeleteCases[u,w],Q],Prepend[#,w]&/@stableSets[DeleteCases[u,r_/;r==w||Q[r,w]||Q[w,r]],Q]]]];
    Table[Length[stableSets[Divisors[n],!(Divisible[#1,#2]||Divisible[#2,#1])&]],{n,10}]

Formula

a(n) = A253249(n) + 1.
Previous Showing 21-30 of 75 results. Next