A345957 -id:A345957 - cp's OEIS frontend

A329888 a(n) = A329900(A329602(n)); Heinz number of the even bisection (even-indexed parts) of the integer partition with Heinz number n.

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

Offset: 1

Antti Karttunen, Dec 22 2019

nonn

From Gus Wiseman, Aug 05 2021 and Antti Karttunen, Oct 13 2021: (Start)
Also the product of primes at even positions in the weakly decreasing list (with multiplicity) of prime factors of n. For example, the prime factors of 108 are (3,3,3,2,2), with even bisection (3,2), with product 6, so a(108) = 6.
Proof: A108951(n) gives a number with the same largest prime factor (A006530) and its exponent (A071178) as in n, and with each smaller prime p = 2, 3, 5, 7, ... < A006530(n) having as its exponent the partial sum of the exponents of all prime factors >= p present in n (with primes not present in n having the exponent 0). Then applying A000188 replaces each such "partial sum exponent" k with floor(k/2). Finally, A319626 replaces those halved exponents with their first differences (here the exponent of the largest prime present stays intact, because the next larger prime's exponent is 0 in n). It should be easy to see that if prime q is not present in n (i.e., does not divide it), then neither it is present in a(n). Moreover, if the partial sum exponent of q is odd and only one larger than the partial sum exponent of the next larger prime factor of n, then q will not be present in a(n), while in all other cases q is present in a(n). See also the last example.
(End)

			From _Gus Wiseman_, Aug 15 2021: (Start)
The list of all numbers with image 12 and their corresponding prime factors begins:
  144: (3,3,2,2,2,2)
  216: (3,3,3,2,2,2)
  240: (5,3,2,2,2,2)
  288: (3,3,2,2,2,2,2)
  336: (7,3,2,2,2,2)
  360: (5,3,3,2,2,2)
(End)
The positions from the left are indexed as 1, 2, 3, ..., etc, so e.g., for 240 we pick the second, the fourth and the sixth prime factor, 3, 2 and 2, to obtain a(240) = 3*2*2 = 12. For 288, we similarly pick the second (3), the fourth (2) and the sixth (2) to obtain a(288) = 3*2*2 = 12. - _Antti Karttunen_, Oct 13 2021
Consider n = 11945934 = 2*3*3*3*7*11*13*13*17. Its primorial inflation is A108951(11945934) = 96478365991115908800000 = 2^9 * 3^8 * 5^5 * 7^5 * 11^4 * 13^3 * 17^1. Applying A000188 to this halves each exponent (floored down if the exponent is odd), leaving the factors 2^4 * 3^4 * 5^2 * 7^2 * 11^2 * 13^1 = 2497294800. Then applying A319626 to this number retains the largest prime factor (and its exponent), and subtracts from the exponent of each of the rest of primes the exponent of the next larger prime, so from 2^4 * 3^4 * 5^2 * 7^2 * 11^2 * 13^1 we get 2^(4-4) * 3^(4-2) * 5^(2-2) * 7^(2-2) * 11^(2-1) * 13^1 = 3^2 * 11^1 * 13^1 = 1287 = a(11945934), which is obtained also by selecting every second prime from the list [17, 13, 13, 11, 7, 3, 3, 3, 2] and taking their product. - _Antti Karttunen_, Oct 15 2021

Antti Karttunen, Table of n, a(n) for n = 1..16384
Antti Karttunen, Data supplement: n, a(n) computed for n = 1..65537
Index entries for sequences related to Heinz numbers
Index entries for sequences related to primorial numbers

Cf. A000188, A108951, A319626, A329602, A329900.
A left inverse of A000290.
Positions of 1's are A008578.
Positions of primes are A168645.
The sum of prime indices of a(n) is A346700(n).
The odd version is A346701.
The odd non-reverse version is A346703.
The non-reverse version is A346704.
The version for standard compositions is A346705, odd A346702.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A001414 adds up prime factors, row sums of A027746.
A027187 counts partitions of even length, ranked by A028260.
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A346633 adds up the even bisection of standard compositions.
A346698 adds up the even bisection of prime indices.
Cf. A000097, A035363, A102750, A236913, A253553, A344606, A344617, A344653, A345957, A345958, A345959.

Table[Times@@Last/@Partition[Reverse[Flatten[Apply[ConstantArray,FactorInteger[n],{1}]]],2],{n,100}] (* Gus Wiseman, Oct 13 2021 *)

PARI
```
A329888(n) = A329900(A329602(n));
```

A329888(n) = if(1==n,n,my(f=factor(n),m=1,p=0); forstep(k=#f~,1,-1,while(f[k,2], m *= f[k,1]^(p%2); f[k,2]--; p++)); (m)); \\ (After Wiseman's new interpretation) - Antti Karttunen, Sep 21 2021

a(n) = A329900(A329602(n)) = A329900(A000188(A108951(n))).
A108951(a(n)) = A329602(n).
a(n^2) = n for all n >= 1.
a(n) * A346701(n) = n. - Gus Wiseman, Aug 07 2021
A056239(a(n)) = A346700(n). - Gus Wiseman, Aug 07 2021
Antti Karttunen, Sep 21 2021
From Antti Karttunen, Oct 13 2021: (Start)
a(n) = A319626(A329602(n)) = A319626(A000188(A108951(n))).
For all x in A102750, a(x) = a(A253553(x)). (End)

Name amended with Gus Wiseman's new interpretation - Antti Karttunen, Oct 13 2021

A346701 Heinz number of the odd bisection (odd-indexed parts) of the integer partition with Heinz number n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 03 2021

nonn

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.

			The partition (2,2,2,1,1) has Heinz number 108 and odd bisection (2,2,1) with Heinz number 18, so a(108) = 18.
The partitions (3,2,2,1,1), (3,2,2,2,1), (3,3,2,1,1) have Heinz numbers 180, 270, 300 and all have odd bisection (3,2,1) with Heinz number 30, so a(180) = a(270) = a(300) = 30.

Positions of last appearances are A000290 without the first term 0.
Positions of primes are A037143 (complement: A033942).
The even version is A329888.
Positions of first appearances are A342768.
The sum of prime indices of a(n) is A346699(n), non-reverse: A346697.
The non-reverse version is A346703.
The even non-reverse version is A346704.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum, reverse A344612.
A209281 (shifted) adds up the odd bisection of standard compositions.
A316524 gives the alternating sum of prime indices, reverse A344616.
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.
A344606 counts alternating permutations of prime indices.
A344617 gives the sign of the alternating sum of prime indices.
A346700 gives the sum of the even bisection of reversed prime indices.
Cf. A025047, A027187, A027193, A053738, A106356, A277103, A341446, A344653, A345957, A345958, A345959, A346698, A346702.

primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
Table[Times@@Prime/@First/@Partition[Append[Reverse[primeMS[n]],0],2],{n,100}]

a(n) * A329888(n) = n.

A056239(a(n)) = A346699(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

A347466 Number of factorizations of n^2.

Original entry on oeis.org

1, 2, 2, 5, 2, 9, 2, 11, 5, 9, 2, 29, 2, 9, 9, 22, 2, 29, 2, 29, 9, 9, 2, 77, 5, 9, 11, 29, 2, 66, 2, 42, 9, 9, 9, 109, 2, 9, 9, 77, 2, 66, 2, 29, 29, 9, 2, 181, 5, 29, 9, 29, 2, 77, 9, 77, 9, 9, 2, 269, 2, 9, 29, 77, 9, 66, 2, 29, 9, 66, 2, 323, 2, 9, 29, 29

Offset: 1

Gus Wiseman, Sep 23 2021

nonn

A factorization of n is a weakly increasing sequence of positive integers > 1 with product n.

			The a(1) = 1 through a(8) = 11 factorizations:
  ()  (4)    (9)    (16)       (25)   (36)       (49)   (64)
      (2*2)  (3*3)  (2*8)      (5*5)  (4*9)      (7*7)  (8*8)
                    (4*4)             (6*6)             (2*32)
                    (2*2*4)           (2*18)            (4*16)
                    (2*2*2*2)         (3*12)            (2*4*8)
                                      (2*2*9)           (4*4*4)
                                      (2*3*6)           (2*2*16)
                                      (3*3*4)           (2*2*2*8)
                                      (2*2*3*3)         (2*2*4*4)
                                                        (2*2*2*2*4)
                                                        (2*2*2*2*2*2)

Antti Karttunen, Table of n, a(n) for n = 1..16384
Index entries for sequences computed from exponents in factorization of n

Positions of 2's are the primes (A000040), which have squares A001248.
The restriction to powers of 2 is A058696.
The additive version (partitions) is A072213.
The case of integer alternating product is A347459, nonsquared A347439.
A000290 lists squares, complement A000037.
A001055 counts factorizations.
A339846 counts even-length factorizations.
A339890 counts odd-length factorizations.
A347050 = factorizations with alternating permutation, complement A347706.
Cf. A000041, A062312, A120452, A144338, A273013, A330972, A345957, A346635, A347437, A347438, A347457, A347460, A347464.

b:= proc(n, k) option remember; `if`(n>k, 0, 1)+`if`(isprime(n), 0,
      add(`if`(d>k, 0, b(n/d, d)), d=numtheory[divisors](n) minus {1, n}))
    end:
a:= proc(n) option remember; b((l-> mul(ithprime(i)^l[i], i=1..nops(l)))(
      sort(map(i-> i[2], ifactors(n^2)[2]), `>`))$2)
    end:
seq(a(n), n=1..76);  # Alois P. Heinz, Oct 14 2021

facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
Table[Length[facs[n^2]],{n,25}]

A001055(n, m=n) = if(1==n, 1, my(s=0); fordiv(n, d, if((d>1)&&(d<=m), s += A001055(n/d, d))); (s));
A347466(n) = A001055(n^2); \\ Antti Karttunen, Oct 13 2021

a(n) = A001055(A000290(n)).

A342768 a(n) = A342767(n, n).

Original entry on oeis.org

1, 2, 3, 8, 5, 12, 7, 32, 27, 20, 11, 48, 13, 28, 45, 128, 17, 108, 19, 80, 63, 44, 23, 192, 125, 52, 243, 112, 29, 180, 31, 512, 99, 68, 175, 432, 37, 76, 117, 320, 41, 252, 43, 176, 405, 92, 47, 768, 343, 500, 153, 208, 53, 972, 275, 448, 171, 116, 59, 720

Offset: 1

Rémy Sigrist, Apr 02 2021

nonn

This sequence has similarities with A087019.

These are the positions of first appearances of each positive integer in A346701, and also in A346703. - Gus Wiseman, Aug 09 2021

			For n = 42:
- 42 = 2 * 3 * 7, so:
          2 3 7
        x 2 3 7
        -------
          2 3 7
        2 3 3
    + 2 2 2
    -----------
      2 2 3 3 7
- hence a(42) = 2 * 2 * 3 * 3 * 7 = 252.

Rémy Sigrist, Table of n, a(n) for n = 1..10000
Rémy Sigrist, PARI program for A342768
Index entries for sequences related to dismal (or lunar) arithmetic

Cf. A007947, A087019, A342767.
The sum of prime indices of a(n) is 2*A056239(n) - A061395(n).
The version for even indices is A129597(n) = 2*a(n) for n > 1.
The sorted version is A346635.
These are the positions of first appearances in A346701 and in A346703.
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A027193 counts partitions of odd length, ranked by A026424.
A209281 adds up the odd bisection of standard compositions (even: A346633).
A346697 adds up the odd bisection of prime indices (reverse: A346699).
Cf. A000290, A006530, A033942, A037143, A344653, A345957, A346698, A346700, A346704.

Table[n^2/FactorInteger[n][[-1,1]],{n,100}] (* Gus Wiseman, Aug 09 2021 *)

PARI
```
See Links section.
```

a(n) = n iff n = 1 or n is a prime number.
a(p^k) = p^(2*k-1) for any k > 0 and any prime number p.
A007947(a(n)) = A007947(n).
A001222(a(n)) = 2*A001222(n) - 1 for any n > 1.
From Gus Wiseman, Aug 09 2021: (Start)
A001221(a(n)) = A001221(n).
If g = A006530(n) is the greatest prime factor of n, then a(n) = n^2/g.
a(n) = A129597(n)/2.
(End)

A346635 Numbers whose division (or multiplication) by their greatest prime factor yields a perfect square. Numbers k such that k*A006530(k) is a perfect square.

Original entry on oeis.org

1, 2, 3, 5, 7, 8, 11, 12, 13, 17, 19, 20, 23, 27, 28, 29, 31, 32, 37, 41, 43, 44, 45, 47, 48, 52, 53, 59, 61, 63, 67, 68, 71, 73, 76, 79, 80, 83, 89, 92, 97, 99, 101, 103, 107, 108, 109, 112, 113, 116, 117, 124, 125, 127, 128, 131, 137, 139, 148, 149, 151, 153

Offset: 1

Gus Wiseman, Aug 10 2021

nonn

This is the sorted version of A342768(n) = position of first appearance of n in A346701 (but A346703 works also).

			The terms together with their prime indices begin:
     1: {}          31: {11}            71: {20}
     2: {1}         32: {1,1,1,1,1}     73: {21}
     3: {2}         37: {12}            76: {1,1,8}
     5: {3}         41: {13}            79: {22}
     7: {4}         43: {14}            80: {1,1,1,1,3}
     8: {1,1,1}     44: {1,1,5}         83: {23}
    11: {5}         45: {2,2,3}         89: {24}
    12: {1,1,2}     47: {15}            92: {1,1,9}
    13: {6}         48: {1,1,1,1,2}     97: {25}
    17: {7}         52: {1,1,6}         99: {2,2,5}
    19: {8}         53: {16}           101: {26}
    20: {1,1,3}     59: {17}           103: {27}
    23: {9}         61: {18}           107: {28}
    27: {2,2,2}     63: {2,2,4}        108: {1,1,2,2,2}
    28: {1,1,4}     67: {19}           109: {29}
    29: {10}        68: {1,1,7}        112: {1,1,1,1,4}

Removing 1 gives a subset of A026424.
The unsorted even version is A129597.
The unsorted version is A342768(n) = A342767(n,n).
Except the first term, the even version is 2*a(n).
A000290 lists squares.
A001221 counts distinct prime factors.
A001222 counts all prime factors.
A006530 gives the greatest prime factor.
A061395 gives the greatest prime index.
A027193 counts partitions of odd length.
A056239 adds up prime indices, row sums of A112798.
A209281 = odd bisection sum of standard compositions (even: A346633).
A316524 = alternating sum of prime indices (sign: A344617, rev.: A344616).
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.
A344606 counts alternating permutations of prime indices.
A346697 = odd bisection sum of prime indices (weights of A346703).
A346699 = odd bisection sum of reversed prime indices (weights of A346701).
Cf. A028260, A033942, A035363, A037143, A341446, A344653, A345957, A345958, A345959, A346698, A346700, A346704.

filter:= proc(n) issqr(n/max(numtheory:-factorset(n))) end proc:
filter(1):= true:
select(filter, [$1..200]); # Robert Israel, Nov 26 2022

sqrQ[n_]:=IntegerQ[Sqrt[n]];
Select[Range[100],sqrQ[#*FactorInteger[#][[-1,1]]]&]

isok(m) = (m==1) || issquare(m/vecmax(factor(m)[,1])); \\ Michel Marcus, Aug 12 2021

a(n) = A129597(n)/2 for n > 1.

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

A347042 Number of divisors d > 1 of n such that bigomega(d) divides bigomega(n), where bigomega = A001222.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Aug 17 2021

nonn

			The a(n) divisors for selected n:
  n = 1:   2:   4:   6:   24:  30:  36:  60:  96:  144: 210: 216: 240: 360:
      ---------------------------------------------------------------------
      {}   2    2    2    2    2    2    2    2    2    2    2    2    2
                4    3    3    3    3    3    3    3    3    3    3    3
                     6    4    5    4    4    4    4    5    4    4    4
                          6    30   6    5    6    6    6    6    5    5
                          24        9    6    8    8    7    8    6    6
                                    36   10   12   9    10   9    8    8
                                         15   96   12   14   12   10   9
                                         60        18   15   18   12   10
                                                   144  21   27   15   12
                                                        35   216  20   15
                                                        210       30   18
                                                                  240  20
                                                                       30
                                                                       45
                                                                       360

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

Positions of 1's are A000040.
The smallest of these divisors is A020639
The case of divisors with half bigomega is A345957 (rounded: A096825).
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.
Cf. A026424, A033676, A033677, A335433, A335448, A347045, A347046.

Table[Length[Select[Rest[Divisors[n]],IntegerQ[PrimeOmega[n]/PrimeOmega[#]]&]],{n,100}]

a(n) = my(bn=bigomega(n)); sumdiv(n, d, if (d>1, !(bn % bigomega(d)))); \\ Michel Marcus, Aug 18 2021

from sympy import divisors, primeomega
def a(n):
    bigomegan = primeomega(n)
    return sum(bigomegan%primeomega(d) == 0 for d in divisors(n)[1:])
print([a(n) for n in range(1, 88)]) # Michael S. Branicky, Aug 18 2021

from sympy import factorint, divisors
from sympy.utilities.iterables import multiset_combinations
def A347042(n):
    fs = factorint(n,multiple=True)
    return sum(len(list(multiset_combinations(fs,d))) for d in divisors(len(fs),generator=True)) # Chai Wah Wu, Aug 21 2021

A129597 Central diagonal of array A129595.

Original entry on oeis.org

1, 4, 6, 16, 10, 24, 14, 64, 54, 40, 22, 96, 26, 56, 90, 256, 34, 216, 38, 160, 126, 88, 46, 384, 250, 104, 486, 224, 58, 360, 62, 1024, 198, 136, 350, 864, 74, 152, 234, 640, 82, 504, 86, 352, 810, 184, 94, 1536, 686, 1000, 306, 416, 106, 1944, 550, 896, 342

Offset: 1

Antti Karttunen, May 01 2007, based on Marc LeBrun's Jan 11 2006 message on SeqFan mailing list

nonn

These are the positions of first appearances of each positive integer in A346704. - Gus Wiseman, Oct 16 2021

Antti Karttunen, Table of n, a(n) for n = 1..10000

a(n) = A129595(n,n).
The sum of prime indices of a(n) is 2*A056239(n) - A061395(n) + 1 for n > 1.
The version for odd indices is A342768(n) = a(n)/2 for n > 1.
Except the first term, the sorted version is 2*A346635.
These are the positions of first appearances in A346704.
A001221 counts distinct prime factors.
A001222 counts prime factors with multiplicity.
A027187 counts partitions of even length, ranked by A028260.
A346633 adds up the even bisection of standard compositions (odd: A209281).
A346698 adds up the even bisection of prime indices (reverse: A346699).
Cf. A000290, A006530, A037143, A329888, A344606, A345957, A346697, A346700, A346701.

Table[If[n==1,1,2*n^2/FactorInteger[n][[-1,1]]],{n,100}] (* Gus Wiseman, Aug 10 2021 *)

A129597(n) = if(1==n, n, my(f=factor(n)); (2*n*n)/f[#f~, 1]); \\ Antti Karttunen, Oct 16 2021

From Gus Wiseman, Aug 10 2021: (Start)
For n > 1, A001221(a(n)) = A099812(n).
If g = A006530(n) is the greatest prime factor of n > 1, then a(n) = 2n^2/g.
a(n) = A100484(A000720(n)) = 2n iff n is prime.
a(n > 1) = 2*A342768(n).
(End)

cp's OEIS Frontend

A329888 a(n) = A329900(A329602(n)); Heinz number of the even bisection (even-indexed parts) of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

PARI

Formula

Extensions

A346701 Heinz number of the odd bisection (odd-indexed parts) of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

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

Programs

Mathematica

Python

Python

Formula

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Python

Python

Formula

A347466 Number of factorizations of n^2.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A342768 a(n) = A342767(n, n).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A346635 Numbers whose division (or multiplication) by their greatest prime factor yields a perfect square. Numbers k such that k*A006530(k) is a perfect square.

Views

Author

Keywords

Comments