A096825 -id:A096825 - cp's OEIS frontend

A096827 Number of antichains in divisor lattice D(n).

2, 3, 3, 4, 3, 6, 3, 5, 4, 6, 3, 10, 3, 6, 6, 6, 3, 10, 3, 10, 6, 6, 3, 15, 4, 6, 5, 10, 3, 20, 3, 7, 6, 6, 6, 20, 3, 6, 6, 15, 3, 20, 3, 10, 10, 6, 3, 21, 4, 10, 6, 10, 3, 15, 6, 15, 6, 6, 3, 50, 3, 6, 10, 8, 6, 20, 3, 10, 6, 20, 3, 35, 3, 6, 10, 10, 6, 20, 3, 21, 6, 6, 3, 50, 6, 6, 6, 15, 3, 50, 6

Offset: 1

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

nonn

The divisor lattice D(n) is the lattice of the divisors of the natural number n.
The empty set is counted as an antichain in D(n).
a(n) = gamma(n+1) where gamma is degree of cardinal completeness of Łukasiewicz n-valued logic. - Artur Jasinski, Mar 01 2010

Alexander S. Karpenko, Lukasiewicz's Logics and Prime Numbers, Luniver Press, Beckington, 2006. See Table I p. 113.

Arkadiusz Wesolowski, Table of n, a(n) for n = 1..990
Index entries for sequences related to Łukasiewicz

Cf. A096825, A096826, A097699.
Cf. A175177, A175178. - Artur Jasinski, Mar 01 2010
Cf. A000005, A008480, A253249, A285572 A285573.

nn=200;
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]],{n,nn}] (* Gus Wiseman, Aug 24 2018 *)

a(n) = A285573(n) + 1. - Gus Wiseman, Aug 24 2018

More terms from John W. Layman, Aug 20 2004

A345957 Number of divisors of n with exactly half as many prime factors as n, counting multiplicity.

Original entry on oeis.org

1, 0, 0, 1, 0, 2, 0, 0, 1, 2, 0, 0, 0, 2, 2, 1, 0, 0, 0, 0, 2, 2, 0, 2, 1, 2, 0, 0, 0, 0, 0, 0, 2, 2, 2, 3, 0, 2, 2, 2, 0, 0, 0, 0, 0, 2, 0, 0, 1, 0, 2, 0, 0, 2, 2, 2, 2, 2, 0, 4, 0, 2, 0, 1, 2, 0, 0, 0, 2, 0, 0, 0, 0, 2, 0, 0, 2, 0, 0, 0, 1, 2, 0, 4, 2, 2, 2

Offset: 1

Gus Wiseman, Aug 16 2021

nonn

These divisors do not necessarily include the central divisors (A207375), and may not themselves be central.

			The a(n) divisors for selected n:
  n = 1:  6:  36:  60:  210:  840:  900:  1260:  1296:  3600:
     --------------------------------------------------------
      1   2    4    4     6     8    12     12     16     16
          3    6    6    10    12    18     18     24     24
               9   10    14    20    20     20     36     36
                   15    15    28    30     28     54     40
                         21    30    45     30     81     60
                         35    42    50     42            90
                               70    75     45           100
                              105           63           150
                                            70           225
                                           105

The case of powers of 2 is A000035.
Positions of even terms are A000037.
Positions of odd terms are A000290.
Positions of 0's are A026424.
Positions of 1's are A056798.
The rounded version is A096825.
The case of all divisors (not just 2) is A347042.
The smallest of these divisors is A347045 (rounded: A347043).
The greatest of these divisors is A347046 (rounded: A347044).
A000005 counts divisors.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A056239 adds up prime indices, row sums of A112798.
A207375 lists central divisors.
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.
A334997 counts chains of divisors of n by length.
Cf. A001227, A001414, A028260, A033676, A033677, A073093, A074206, A217581, A344653, A346697-A346704.

Table[Length[Select[Divisors[n],PrimeOmega[#]==PrimeOmega[n]/2&]],{n,100}]

a(n) = my(nb=bigomega(n)); sumdiv(n, d, 2*bigomega(d) == nb); \\ Michel Marcus, Aug 16 2021

from sympy import divisors, factorint
def a(n):
    npf = len(factorint(n, multiple=True))
    divs = divisors(n)
    return sum(2*len(factorint(d, multiple=True)) == npf for d in divs)
print([a(n) for n in range(1, 88)]) # Michael S. Branicky, Aug 17 2021
(Python 3.8+)
from itertools import combinations
from math import prod, comb
from sympy import factorint
def A345957(n):
    if n == 1:
        return 1
    fs = factorint(n)
    elist = list(fs.values())
    q, r = divmod(sum(elist),2)
    k = len(elist)
    if r:
        return 0
    c = 0
    for i in range(k+1):
        m = (-1)**i
        for d in combinations(range(k),i):
            t = k+q-sum(elist[j] for j in d)-i-1
            if t >= 0:
                c += m*comb(t,k-1)
    return c # Chai Wah Wu, Aug 20 2021

from sympy import factorint
from sympy.utilities.iterables import multiset_combinations
def A345957(n):
    if n == 1:
        return 1
    fs = factorint(n,multiple=True)
    q, r = divmod(len(fs),2)
    return 0 if r else len(list(multiset_combinations(fs,q))) # Chai Wah Wu, Aug 20 2021

A347044 Greatest divisor of n with half (rounded up) as many prime factors (counting multiplicity) as n.

Original entry on oeis.org

1, 2, 3, 2, 5, 3, 7, 4, 3, 5, 11, 6, 13, 7, 5, 4, 17, 9, 19, 10, 7, 11, 23, 6, 5, 13, 9, 14, 29, 15, 31, 8, 11, 17, 7, 9, 37, 19, 13, 10, 41, 21, 43, 22, 15, 23, 47, 12, 7, 25, 17, 26, 53, 9, 11, 14, 19, 29, 59, 15, 61, 31, 21, 8, 13, 33, 67, 34, 23, 35, 71

Offset: 1

Gus Wiseman, Aug 16 2021

nonn

Appears to contain each positive integer at least once, but only a finite number of times.

			The divisors of 123456 with half bigomega are: 16, 24, 5144, 7716, so a(123456) = 7716.

Amiram Eldar, Table of n, a(n) for n = 1..10000

The greatest divisor without the condition is A006530 (smallest: A020639).
Divisors of this type are counted by A096825 (exact: A345957).
The case of powers of 2 is A163403.
The smallest divisor of this type is given by A347043 (exact: A347045).
The exact version is A347046.
A000005 counts divisors.
A001221 counts distinct prime factors.
A001222 counts all prime factors (also called bigomega).
A038548 counts inferior (or superior) divisors (strict: A056924).
A056239 adds up prime indices, row sums of A112798.
A207375 lists central divisors (min: A033676, max: A033677).
A340387 lists numbers whose sum of prime indices is twice bigomega.
A340609 lists numbers whose maximum prime index divides bigomega.
A340610 lists numbers whose maximum prime index is divisible by bigomega.
A347042 counts divisors d|n such that bigomega(d) divides bigomega(n).
Cf. A000720, A027746, A060775, A106529, A217581, A244990, A335433, A335448.

Table[Max[Select[Divisors[n],PrimeOmega[#]==Ceiling[PrimeOmega[n]/2]&]],{n,100}]
a[n_] := Module[{p = Flatten[Table[#[[1]], {#[[2]]}] & /@ FactorInteger[n]], np}, np = Length[p]; Times @@ p[[Floor[np/2] + 1;; np]]]; Array[a, 100] (* Amiram Eldar, Nov 02 2024 *)

from sympy import divisors, factorint
def a(n):
    npf = len(factorint(n, multiple=True))
    for d in divisors(n)[::-1]:
        if len(factorint(d, multiple=True)) == (npf+1)//2: return d
    return 1
print([a(n) for n in range(1, 72)]) # Michael S. Branicky, Aug 18 2021

from math import prod
from sympy import factorint
def A347044(n):
    fs = factorint(n,multiple=True)
    l = len(fs)
    return prod(fs[l//2:]) # Chai Wah Wu, Aug 20 2021

a(n) = Product_{k=floor(A001222(n)/2)+1..A001222(n)} A027746(n,k). - Amiram Eldar, Nov 02 2024

A347043 Smallest divisor of n with half (rounded up) as many prime factors (counting multiplicity) as n.

Original entry on oeis.org

1, 2, 3, 2, 5, 2, 7, 4, 3, 2, 11, 4, 13, 2, 3, 4, 17, 6, 19, 4, 3, 2, 23, 4, 5, 2, 9, 4, 29, 6, 31, 8, 3, 2, 5, 4, 37, 2, 3, 4, 41, 6, 43, 4, 9, 2, 47, 8, 7, 10, 3, 4, 53, 6, 5, 4, 3, 2, 59, 4, 61, 2, 9, 8, 5, 6, 67, 4, 3, 10, 71, 8, 73, 2, 15, 4, 7, 6, 79, 8

Offset: 1

Gus Wiseman, Aug 15 2021

nonn

Appears to contain every positive integer at least once.

This is correct. For any integer m, let p be any prime > m. Then a(m*p^A001222(m)) = m. - Sebastian Karlsson, Oct 11 2022

			The divisors of 123456 with half bigomega are: 16, 24, 5144, 7716, so a(123456) = 16.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Positions of 2's are A001747.
Positions of odd terms are A005408.
Positions of even terms are A005843.
The case of powers of 2 is A016116.
The smallest divisor without the condition is A020639 (greatest: A006530).
These divisors are counted by A096825 (exact: A345957).
The greatest of these divisors is A347044 (exact: A347046).
The exact version is A347045.
A000005 counts divisors.
A001221 counts distinct prime factors.
A001222 counts all prime factors (also called bigomega).
A056239 adds up prime indices, row sums of A112798.
A207375 lists central divisors (min: A033676, max: A033677).
A340387 lists numbers whose sum of prime indices is twice bigomega.
A340609 lists numbers whose maximum prime index divides bigomega.
A340610 lists numbers whose maximum prime index is divisible by bigomega.
A347042 counts divisors d|n such that bigomega(d) divides bigomega(n).
Cf. A000720, A027746, A038548, A106529, A140271, A244991, A324522, A335433, A335448.

Table[Min[Select[Divisors[n],PrimeOmega[#]==Ceiling[PrimeOmega[n]/2]&]],{n,100}]
a[n_] := Module[{p = Flatten[Table[#[[1]], {#[[2]]}] & /@ FactorInteger[n]]}, Times @@ p[[1 ;; Ceiling[Length[p]/2]]]]; Array[a, 100] (* Amiram Eldar, Nov 02 2024 *)

a(n) = my(bn=ceil(bigomega(n)/2)); fordiv(n, d, if (bigomega(d)==bn, return (d))); \\ Michel Marcus, Aug 18 2021

from sympy import divisors, factorint
def a(n):
    npf = len(factorint(n, multiple=True))
    for d in divisors(n):
        if len(factorint(d, multiple=True)) == (npf+1)//2: return d
    return 1
print([a(n) for n in range(1, 81)]) # Michael S. Branicky, Aug 18 2021

from math import prod
from sympy import factorint
def A347043(n):
    fs = factorint(n,multiple=True)
    l = len(fs)
    return prod(fs[:(l+1)//2]) # Chai Wah Wu, Aug 20 2021

a(n) = Product_{k=1..ceiling(A001222(n)/2)} A027746(n,k). - Amiram Eldar, Nov 02 2024

A361200 Product of the left half (exclusive) of the multiset of prime factors of n; a(1) = 0.

Original entry on oeis.org

0, 1, 1, 2, 1, 2, 1, 2, 3, 2, 1, 2, 1, 2, 3, 4, 1, 2, 1, 2, 3, 2, 1, 4, 5, 2, 3, 2, 1, 2, 1, 4, 3, 2, 5, 4, 1, 2, 3, 4, 1, 2, 1, 2, 3, 2, 1, 4, 7, 2, 3, 2, 1, 6, 5, 4, 3, 2, 1, 4, 1, 2, 3, 8, 5, 2, 1, 2, 3, 2, 1, 4, 1, 2, 3, 2, 7, 2, 1, 4, 9, 2, 1, 4, 5, 2, 3

Offset: 1

Gus Wiseman, Mar 10 2023

nonn

			The prime factors of 250 are {2,5,5,5}, with left half (exclusive) {2,5}, with product 10, so a(250) = 10.

Amiram Eldar, Table of n, a(n) for n = 1..10000

Positions of 1's are A000040.
Positions of 2's are A037143.
The inclusive version is A347043.
The right inclusive version A347044.
The right version is A361201.
A000005 counts divisors.
A001221 counts distinct prime factors.
A006530 gives greatest prime factor.
A112798 lists prime indices, length A001222, sum A056239.
A360616 gives half of bigomega (exclusive), inclusive A360617.
A360673 counts multisets by right sum (exclusive), inclusive A360671.
First for prime indices, second for partitions, third for prime factors:
- A360676 gives left sum (exclusive), counted by A360672, product A361200.
- A360677 gives right sum (exclusive), counted by A360675, product A361201.
- A360678 gives left sum (inclusive), counted by A360675, product A347043.
- A360679 gives right sum (inclusive), counted by A360672, product A347044.
Cf. A001248, A026424, A027746, A096825, A347045, A347046, A360005.

Table[If[n==1,0,Times@@Take[Join@@ConstantArray@@@FactorInteger[n],Floor[PrimeOmega[n]/2]]],{n,100}]
a[n_] := Module[{p = Flatten[Table[#[[1]], {#[[2]]}] & /@ FactorInteger[n]]}, Times @@ p[[1 ;; Floor[Length[p]/2]]]]; a[1] = 0; Array[a, 100] (* Amiram Eldar, Nov 02 2024 *)

a(n) * A347044(n) = n.
A361201(n) * A347043(n) = n.
a(n) = Product_{k=1..floor(A001222(n)/2)} A027746(n,k) for n >= 2. - Amiram Eldar, Nov 02 2024

A361201 Product of the right half (exclusive) of the multiset of prime factors of n; a(1) = 0.

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 2, 3, 5, 1, 3, 1, 7, 5, 4, 1, 3, 1, 5, 7, 11, 1, 6, 5, 13, 3, 7, 1, 5, 1, 4, 11, 17, 7, 9, 1, 19, 13, 10, 1, 7, 1, 11, 5, 23, 1, 6, 7, 5, 17, 13, 1, 9, 11, 14, 19, 29, 1, 15, 1, 31, 7, 8, 13, 11, 1, 17, 23, 7, 1, 9, 1, 37, 5, 19, 11, 13, 1

Offset: 1

Gus Wiseman, Mar 10 2023

			The prime factors of 250 are {2,5,5,5}, with right half (exclusive) {5,5}, with product 25, so a(250) = 25.

Robert Israel, Table of n, a(n) for n = 1..10000

Positions of 1's are A000040.
Positions of first appearances are A123666.
The left inclusive version A347043.
The inclusive version is A347044.
The left version is A361200.
A000005 counts divisors.
A001221 counts distinct prime factors.
A006530 gives greatest prime factor.
A112798 lists prime indices, length A001222, sum A056239.
A360616 gives half of bigomega (exclusive), inclusive A360617.
A360673 counts multisets by right sum (exclusive), inclusive A360671.
First for prime indices, second for partitions, third for prime factors:
- A360676 gives left sum (exclusive), counted by A360672, product A361200.
- A360677 gives right sum (exclusive), counted by A360675, product A361201.
- A360678 gives left sum (inclusive), counted by A360675, product A347043.
- A360679 gives right sum (inclusive), counted by A360672, product A347044.
Cf. A001248, A026424, A096825, A347045, A347046, A360005.

f:= proc(n) local F;
  F:= ifactors(n)[2];
  F:= sort(map(t -> t[1]$t[2],F));
  convert(F[ceil(nops(F)/2)+1 ..-1],`*`)
end proc:
f(1):= 0:
map(f, [$1..100]); # Robert Israel, Aug 12 2024

Table[If[n==1,0,Times@@Take[Join@@ConstantArray@@@FactorInteger[n],-Floor[PrimeOmega[n]/2]]],{n,100}]

A361200(n) * A347044(n) = n.

A361201(n) * A347043(n) = n.

A347045 Smallest divisor of n with exactly half as many prime factors (counting multiplicity) as n, or 1 if there are none.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 16 2021

nonn

			The divisors of 90 with half bigomega are: 6, 9, 10, 15, so a(90) = 6.

Amiram Eldar, Table of n, a(n) for n = 1..10000

The smallest divisor without the condition is A020639 (greatest: A006530).
Positions of 1's are A026424.
Positions of even terms are A063745 = 2*A026424.
The case of powers of 2 is A072345.
Positions of 2's are A100484.
Divisors of this type are counted by A345957 (rounded: A096825).
The rounded version is A347043.
The greatest divisor of this type is A347046 (rounded: A347044).
A000005 counts divisors.
A001221 counts distinct prime factors.
A001222 counts all prime factors (also called bigomega).
A056239 adds up prime indices, row sums of A112798.
A207375 lists central divisors (min: A033676, max: A033677).
A340387 lists numbers whose sum of prime indices is twice bigomega.
A340609 lists numbers whose maximum prime index divides bigomega.
A340610 lists numbers whose maximum prime index is divisible by bigomega.
A347042 counts divisors d|n such that bigomega(d) divides bigomega(n).
Cf. A000720, A001747, A027746, A038548, A056924, A063962, A106529, A140271, A244991, A324522, A335433, A335448.

Table[If[#=={},1,Min[#]]&@Select[Divisors[n], PrimeOmega[#]==PrimeOmega[n]/2&],{n,100}]
a[n_] := Module[{p = Flatten[Table[#[[1]], {#[[2]]}] & /@ FactorInteger[n]], np}, np = Length[p]; If[OddQ[np], 1, Times @@ p[[1 ;; np/2]]]]; Array[a, 100] (* Amiram Eldar, Nov 02 2024 *)

from sympy import divisors, factorint
def a(n):
    npf = len(factorint(n, multiple=True))
    for d in divisors(n)[1:-1]:
        if 2*len(factorint(d, multiple=True)) == npf: return d
    return 1
print([a(n) for n in range(1, 88)]) # Michael S. Branicky, Aug 18 2021

from math import prod
from sympy import factorint
def A347045(n):
    fs = factorint(n,multiple=True)
    q, r = divmod(len(fs),2)
    return 1 if r else prod(fs[:q]) # Chai Wah Wu, Aug 20 2021

a(n) = Product_{k=1..A001222(n)/2} A027746(n,k) if A001222(n) is even, and 1 otherwise. - Amiram Eldar, Nov 02 2024

A347046 Greatest divisor of n with exactly half as many prime factors (counting multiplicity) as n, or 1 if there are none.

Original entry on oeis.org

1, 1, 1, 2, 1, 3, 1, 1, 3, 5, 1, 1, 1, 7, 5, 4, 1, 1, 1, 1, 7, 11, 1, 6, 5, 13, 1, 1, 1, 1, 1, 1, 11, 17, 7, 9, 1, 19, 13, 10, 1, 1, 1, 1, 1, 23, 1, 1, 7, 1, 17, 1, 1, 9, 11, 14, 19, 29, 1, 15, 1, 31, 1, 8, 13, 1, 1, 1, 23, 1, 1, 1, 1, 37, 1, 1, 11, 1, 1, 1, 9

Offset: 1

Gus Wiseman, Aug 16 2021

nonn

Problem: What are the positions of last appearances > 1?

			The divisors of 90 with half bigomega are: 6, 9, 10, 15, so a(90) = 15.

Amiram Eldar, Table of n, a(n) for n = 1..10000

The greatest divisor without the condition is A006530 (smallest: A020639).
Positions of 1's are A026424.
The case of powers of 2 is A072345.
Positions of first appearances are A123667 (sorted: A123666).
Divisors of this type are counted by A345957 (rounded: A096825).
The rounded version is A347044.
The smallest divisor of this is A347045 (rounded: A347043).
A000005 counts divisors.
A001221 counts distinct prime factors.
A001222 counts all prime factors (also called bigomega).
A056239 adds up prime indices, row sums of A112798.
A207375 lists central divisors (min: A033676, max: A033677).
A340387 lists numbers whose sum of prime indices is twice bigomega.
A340609 lists numbers whose maximum prime index divides bigomega.
A340610 lists numbers whose maximum prime index is divisible by bigomega.
A347042 counts divisors d|n such that bigomega(d) divides bigomega(n).
Cf. A000720, A027746, A038548, A060775, A106529, A244991, A335433, A335448.

Table[If[#=={},1,Max[#]]&@Select[Divisors[n], PrimeOmega[#]==PrimeOmega[n]/2&],{n,100}]
a[n_] := Module[{p = Flatten[Table[#[[1]], {#[[2]]}] & /@ FactorInteger[n]], np}, np = Length[p]; If[OddQ[np], 1, Times @@ p[[np/2+1 ;; np]]]]; Array[a, 100] (* Amiram Eldar, Nov 02 2024 *)

from sympy import divisors, factorint
def a(n):
    npf = len(factorint(n, multiple=True))
    for d in divisors(n)[-1:0:-1]:
        if 2*len(factorint(d, multiple=True)) == npf: return d
    return 1
print([a(n) for n in range(1, 82)]) # Michael S. Branicky, Aug 18 2021

from math import prod
from sympy import factorint
def A347046(n):
    fs = factorint(n,multiple=True)
    q, r = divmod(len(fs),2)
    return 1 if r else prod(fs[q:]) # Chai Wah Wu, Aug 20 2021

a(n) = Product_{k=A001222(n)/2+1..A001222(n)} A027746(n,k) if A001222(n) is even, and 1 otherwise. - Amiram Eldar, Nov 02 2024

A334144 Consider the mapping k -> (k - (k/p)), where prime p | k. a(n) = maximum distinct terms at any position j among the various paths to 1.

Original entry on oeis.org

1, 1, 1, 1, 1, 2, 2, 1, 2, 2, 2, 2, 2, 2, 3, 1, 1, 2, 2, 2, 3, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 1, 4, 2, 4, 3, 3, 3, 3, 2, 2, 4, 4, 3, 4, 3, 3, 2, 4, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 1, 5, 5, 5, 2, 5, 5, 5, 3, 3, 3, 4, 3, 6, 4, 4, 2, 3, 2, 2, 4, 3, 4, 4, 3, 3, 5, 5, 3, 5, 3, 5, 2, 2, 4, 6, 3, 3, 3, 3, 3, 6, 3

Offset: 1

Michael De Vlieger, Peter Kagey, Antti Karttunen, Apr 15 2020

nonn

Let i = A064097(n) be the common path length and let 1 <= j <= i. Given a path P, we find for any j relatively few distinct values. Regarding a common path length i, see A333123 comment 2, and proof at A064097.

Maximum term in row n of A334184.

			For n=15, the paths are shown vertically at left, and the graph obtained appears at right:
  15   15   15   15   15  =>         15
   |    |    |    |    |            _/ \_
   |    |    |    |    |           /     \
  10   10   12   12   12  =>     10       12
   |    |    |    |    |         | \_   _/ |
   |    |    |    |    |         |   \ /   |
   5    8    6    6    8  =>     5    8    6
   |    |    |    |    |          \_  |  _/|
   |    |    |    |    |            \_|_/  |
   4    4    3    4    4  =>          4    3
   |    |    |    |    |              |  _/
   |    |    |    |    |              |_/
   2    2    2    2    2  =>          2
   |    |    |    |    |              |
   |    |    |    |    |              |
   1    1    1    1    1  =>          1
Because the maximum number of distinct terms in any row is 3, a(15) = 3.

Peter Kagey, Table of n, a(n) for n = 1..10000
Peter Kagey, Math.StackExchange, Does a graded poset on the positive integers generated from subtracting factors define a lattice?
Richard Stanley, MIT OpenCourseWare, Lecture Notes in Algebraic Combinatorics: The Sperner Property
Wikipedia, Sperner property of a partially ordered set

Cf. A064097, A332992, A332999, A333123, A334111, A334184.

Cf. also A096825.

Max[Length@ Union@ # & /@ Transpose@ #] & /@ Nest[Function[{a, n}, Append[a, Join @@ Table[Flatten@ Prepend[#, n] & /@ a[[n - n/p]], {p, FactorInteger[n][[All, 1]]}]]] @@ {#, Length@ # + 1} &, {{{1}}}, 105]
(* Second program: *)
g[n_] := Block[{lst = {{n}}}, While[lst[[-1]] != {1}, lst = Join[lst, {Union@ Flatten[# - #/(First@ # & /@ FactorInteger@ #) & /@ lst[[-1]] ]}]]; Max[Length /@ lst]]; Array[g, 105] (* Robert G. Wilson v, May 08 2020 *)

A345926 Number of distinct possible alternating sums of permutations of the multiset of prime indices of 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, 3, 2, 2, 2, 2, 2, 2, 1, 2, 2, 3

Offset: 1

Gus Wiseman, Jul 14 2021

nonn

First differs from A096825 at a(90) = 3, A096825(90) = 4.
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 alternating sum of a sequence (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. Of course, the alternating sum of prime indices is also the reverse-alternating sum of reversed prime indices.
Also the number of possible values of A056239(d) where d is a divisor of n with half as many prime factors (rounded up) as n.

			Grouping the 12 permutations of {1,2,2,3} by alternating sum k gives:
  k = -2: (1223) (1322) (2213) (2312)
  k =  0: (1232) (2123) (2321) (3212)
  k =  2: (2132) (2231) (3122) (3221)
so a(90) = 3.

The version for prime factors instead of indices is A343943.
A000005 counts divisors.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A001414 adds up prime factors, row sums of A027746.
A056239 adds up prime indices,  row sums of A112798.
A097805 counts compositions by alternating (or reverse-alternating) sum.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A316524 gives the alternating sum of prime indices (reverse: A344616).
A345197 counts compositions by length and alternating sum.
A344610 counts partitions by sum and positive reverse-alternating sum.
Cf. A008549, A032443, A083399, A096825, A239830, A344607, A344609, A344651, A345957, A345960, A345961.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
Table[Length[Union[ats/@Permutations[primeMS[n]]]],{n,100}]

from sympy import factorint, primepi
from sympy.utilities.iterables import multiset_combinations
def A345926(n):
    fs = dict((primepi(a),b) for (a,b) in factorint(n).items())
    return len(set(sum(d) for d in multiset_combinations(fs, (sum(fs.values())+1)//2))) # Chai Wah Wu, Aug 23 2021

cp's OEIS Frontend

A096827 Number of antichains in divisor lattice D(n).

Views

Author

Keywords

Comments

References

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A345957 Number of divisors of n with exactly half as many prime factors as n, counting multiplicity.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Python

Python

A347044 Greatest divisor of n with half (rounded up) as many prime factors (counting multiplicity) as n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Python

Formula

A347043 Smallest divisor of n with half (rounded up) as many prime factors (counting multiplicity) as n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Python

Formula

A361200 Product of the left half (exclusive) of the multiset of prime factors of n; a(1) = 0.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A361201 Product of the right half (exclusive) of the multiset of prime factors of n; a(1) = 0.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Maple

Mathematica

Formula

A347045 Smallest divisor of n with exactly half as many prime factors (counting multiplicity) as n, or 1 if there are none.

Views

Author

Keywords

Examples

Links

Crossrefs