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.

Showing 1-10 of 16 results. Next

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

Views

Author

Gus Wiseman, Aug 16 2021

Keywords

Comments

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

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    Table[Length[Select[Divisors[n],PrimeOmega[#]==PrimeOmega[n]/2&]],{n,100}]
  • PARI
    a(n) = my(nb=bigomega(n)); sumdiv(n, d, 2*bigomega(d) == nb); \\ Michel Marcus, Aug 16 2021
    
  • Python
    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
    
  • Python
    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

A347457 Heinz numbers of integer partitions with integer alternating product.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 31, 32, 34, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 52, 53, 54, 56, 57, 58, 59, 61, 62, 63, 64, 65, 67, 68, 71, 72, 73, 74, 75, 76, 78
Offset: 1

Views

Author

Gus Wiseman, Sep 26 2021

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.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
Also numbers whose multiset of prime indices has integer reverse-alternating product.

Examples

			The prime indices of 525 are {2,3,3,4}, with reverse-alternating product 2, so 525 is in the sequence
The prime indices of 135 are {2,2,2,3}, with reverse-alternating product 3/2, so 135 is not in the sequence.
		

Crossrefs

The reciprocal version is A028982.
Allowing any alternating product > 1 gives A028983, reverse A347465.
Factorizations of this type are counted by A347437.
These partitions are counted by A347446.
The reverse reciprocal version A347451.
The odd-length case is A347453.
The reverse version is A347454.
The complement is A347455.
A056239 adds up prime indices, row sums of A112798.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A335433 lists numbers whose prime indices are separable, complement A335448.
A347461 counts possible alternating products of partitions, reverse A347462.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Select[Range[100],IntegerQ[altprod[Reverse[primeMS[#]]]]&]

A347450 Numbers whose multiset of prime indices has alternating product <= 1.

Original entry on oeis.org

1, 2, 4, 6, 8, 9, 10, 14, 15, 16, 18, 21, 22, 24, 25, 26, 32, 33, 34, 35, 36, 38, 39, 40, 46, 49, 50, 51, 54, 55, 56, 57, 58, 60, 62, 64, 65, 69, 72, 74, 77, 81, 82, 84, 85, 86, 87, 88, 90, 91, 93, 94, 95, 96, 98, 100, 104, 106, 111, 115, 118, 119, 121, 122
Offset: 1

Views

Author

Gus Wiseman, Sep 24 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.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
Also Heinz numbers integer partitions with reverse-alternating product <= 1, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Also numbers whose multiset of prime indices has alternating sum <= 1.

Examples

			The initial terms and their prime indices:
      1: {}            26: {1,6}           56: {1,1,1,4}
      2: {1}           32: {1,1,1,1,1}     57: {2,8}
      4: {1,1}         33: {2,5}           58: {1,10}
      6: {1,2}         34: {1,7}           60: {1,1,2,3}
      8: {1,1,1}       35: {3,4}           62: {1,11}
      9: {2,2}         36: {1,1,2,2}       64: {1,1,1,1,1,1}
     10: {1,3}         38: {1,8}           65: {3,6}
     14: {1,4}         39: {2,6}           69: {2,9}
     15: {2,3}         40: {1,1,1,3}       72: {1,1,1,2,2}
     16: {1,1,1,1}     46: {1,9}           74: {1,12}
     18: {1,2,2}       49: {4,4}           77: {4,5}
     21: {2,4}         50: {1,3,3}         81: {2,2,2,2}
     22: {1,5}         51: {2,7}           82: {1,13}
     24: {1,1,1,2}     54: {1,2,2,2}       84: {1,1,2,4}
     25: {3,3}         55: {3,5}           85: {3,7}
		

Crossrefs

The additive version (alternating sum <= 0) is A028260.
The reverse version is A028982, counted by A119620.
Allowing any alternating product < 1 gives A119899.
Factorizations of this type are counted by A339846, complement A339890.
Allowing any alternating product >= 1 gives A344609, multiplicative A347456.
Partitions of this type are counted by A347443.
Allowing any integer alternating product gives A347454, reciprocal A347451.
The complement is A347465, reverse A028983, counted by A347448.
A056239 adds up prime indices, row sums of A112798.
A236913 counts partitions of 2n with reverse-alternating sum <= 0.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A335433 lists numbers whose prime indices are separable, complement A335448.
A344606 counts alternating permutations of prime indices.
A347457 lists Heinz numbers of partitions with integer alternating product.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Select[Range[100],altprod[primeMS[#]]<=1&]

Formula

Union of A028982 and A119899.
Union of A028260 and A001105.

A347454 Numbers whose multiset of prime indices has integer alternating product.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 12, 13, 16, 17, 18, 19, 20, 23, 25, 27, 28, 29, 31, 32, 36, 37, 41, 42, 43, 44, 45, 47, 48, 49, 50, 52, 53, 59, 61, 63, 64, 67, 68, 71, 72, 73, 75, 76, 78, 79, 80, 81, 83, 89, 92, 97, 98, 99, 100, 101, 103, 107, 108, 109, 112, 113
Offset: 1

Views

Author

Gus Wiseman, Sep 26 2021

Keywords

Comments

First differs from A265640 in having 42.
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 alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
Also Heinz numbers of partitions with integer reverse-alternating product, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The terms and their prime indices begin:
      1: {}            20: {1,1,3}         47: {15}
      2: {1}           23: {9}             48: {1,1,1,1,2}
      3: {2}           25: {3,3}           49: {4,4}
      4: {1,1}         27: {2,2,2}         50: {1,3,3}
      5: {3}           28: {1,1,4}         52: {1,1,6}
      7: {4}           29: {10}            53: {16}
      8: {1,1,1}       31: {11}            59: {17}
      9: {2,2}         32: {1,1,1,1,1}     61: {18}
     11: {5}           36: {1,1,2,2}       63: {2,2,4}
     12: {1,1,2}       37: {12}            64: {1,1,1,1,1,1}
     13: {6}           41: {13}            67: {19}
     16: {1,1,1,1}     42: {1,2,4}         68: {1,1,7}
     17: {7}           43: {14}            71: {20}
     18: {1,2,2}       44: {1,1,5}         72: {1,1,1,2,2}
     19: {8}           45: {2,2,3}         73: {21}
		

Crossrefs

The even-length case is A000290.
The additive version is A026424.
Allowing any alternating product < 1 gives A119899, strict A028260.
Allowing any alternating product >= 1 gives A344609, multiplicative A347456.
Factorizations of this type are counted by A347437.
These partitions are counted by A347445, reverse A347446.
Allowing any alternating product <= 1 gives A347450.
The reciprocal version is A347451.
The odd-length case is A347453.
The version for reversed prime indices is A347457, complement A347455.
Allowing any alternating product > 1 gives A347465, reverse A028983.
A056239 adds up prime indices, row sums of A112798.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A335433 lists numbers whose prime indices are separable, complement A335448.
A344606 counts alternating permutations of prime indices.
A347461 counts possible alternating products of partitions.
A347462 counts possible reverse-alternating products of partitions.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Select[Range[100],IntegerQ[altprod[primeMS[#]]]&]

A346703 Product of primes at odd positions in the weakly increasing list (with multiplicity) of prime factors of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 08 2021

Keywords

Examples

			The prime factors of 108 are (2,2,3,3,3), with odd bisection (2,3,3), with product 18, so a(108) = 18.
The prime factors of 720 are (2,2,2,2,3,3,5), with odd bisection (2,2,3,5), with product 60, so a(720) = 60.
		

Crossrefs

Positions of 2's are A001747.
Positions of primes are A037143 (complement: A033942).
The even reverse version appears to be A329888.
Positions of first appearances are A342768.
The sum of prime indices of a(n) is A346697(n), reverse: A346699.
The reverse version is A346701.
The even 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).
A335433/A335448 rank separable/inseparable partitions.
A344606 counts alternating permutations of prime indices.
A344617 gives the sign of the alternating sum of prime indices.
A346633 adds up the even bisection of standard compositions.
A346698 gives the sum of the even bisection of prime indices.
A346700 gives the sum of the even bisection of reversed prime indices.

Programs

  • Mathematica
    Table[Times@@First/@Partition[Append[Flatten[Apply[ConstantArray,FactorInteger[n],{1}]],0],2],{n,100}]

Formula

a(n) * A346704(n) = n.
A056239(a(n)) = A346697(n).

A346700 Sum of the even bisection (even-indexed parts) of the integer partition with Heinz number n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 03 2021

Keywords

Comments

First differs from A334107 at a(64) = 3, A334107(64) = 2.
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.

Examples

			The partition with Heinz number 1100 is (5,3,3,1,1), so a(1100) = 3 + 1 = 4.
The partition with Heinz number 2100 is (4,3,3,2,1,1), so a(2100) = 3 + 2 + 1 = 6.
		

Crossrefs

Sum of prime indices of A329888(n).
Subtracting from the odd version gives A344616 (non-reverse: A316524).
The unreversed version for standard compositions is A346633.
The odd non-reverse version is A346697.
The non-reverse version (multisets instead of partitions) is A346698.
The odd version is A346699.
A001414 adds up prime factors, row sums of A027746.
A056239 adds up prime indices, row sums of A112798.
A103919 counts partitions by sum and alternating sum (reverse: A344612).
A325534 counts separable partitions, ranked by A335433.
A325535 counts inseparable partitions, ranked by A335448.
A344606 counts alternating permutations of prime indices.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Total[Last/@Partition[Reverse[primeMS[n]],2]],{n,100}]
  • PARI
    A346700(n) = if(1==n,0,my(f=factor(n),s=0,p=0); forstep(k=#f~,1,-1,while(f[k,2], s += (p%2)*primepi(f[k,1]); f[k,2]--; p++)); (s)); \\ Antti Karttunen, Sep 21 2021

Formula

a(n) = A056239(n) - A346699(n).
a(n) = A346699(n) - A344616(n).
a(n even omega) = A346697(n).
a(n odd omega) = A346698(n).
A316524(n) = A346697(n) - A346698(n).
a(n) = A056239(A329888(n)). - Gus Wiseman and Antti Karttunen, Oct 13 2021

Extensions

Data section extended up to 105 terms by Antti Karttunen, Sep 21 2021

A347453 Heinz numbers of odd-length integer partitions with integer alternating (or reverse-alternating) product.

Original entry on oeis.org

2, 3, 5, 7, 8, 11, 12, 13, 17, 18, 19, 20, 23, 27, 28, 29, 31, 32, 37, 41, 42, 43, 44, 45, 47, 48, 50, 52, 53, 59, 61, 63, 67, 68, 71, 72, 73, 75, 76, 78, 79, 80, 83, 89, 92, 97, 98, 99, 101, 103, 107, 108, 109, 112, 113, 114, 116, 117, 124, 125, 127, 128, 130
Offset: 1

Views

Author

Gus Wiseman, Sep 24 2021

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.
We define the alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^(i-1)).
Also numbers whose multiset of prime indices has odd length and integer alternating product, where a prime index of n is a number m such that prime(m) divides n.

Examples

			The terms and their prime indices begin:
      2: {1}         29: {10}            61: {18}
      3: {2}         31: {11}            63: {2,2,4}
      5: {3}         32: {1,1,1,1,1}     67: {19}
      7: {4}         37: {12}            68: {1,1,7}
      8: {1,1,1}     41: {13}            71: {20}
     11: {5}         42: {1,2,4}         72: {1,1,1,2,2}
     12: {1,1,2}     43: {14}            73: {21}
     13: {6}         44: {1,1,5}         75: {2,3,3}
     17: {7}         45: {2,2,3}         76: {1,1,8}
     18: {1,2,2}     47: {15}            78: {1,2,6}
     19: {8}         48: {1,1,1,1,2}     79: {22}
     20: {1,1,3}     50: {1,3,3}         80: {1,1,1,1,3}
     23: {9}         52: {1,1,6}         83: {23}
     27: {2,2,2}     53: {16}            89: {24}
     28: {1,1,4}     59: {17}            92: {1,1,9}
		

Crossrefs

The reciprocal version is A000290.
Allowing any alternating product <= 1 gives A001105.
Allowing any alternating product gives A026424.
Factorizations of this type are counted by A347441.
These partitions are counted by A347444.
Allowing any length gives A347454.
Allowing any alternating product > 1 gives A347465.
A027193 counts odd-length partitions.
A056239 adds up prime indices, row sums of A112798.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A335433 lists numbers whose prime indices are separable, complement A335448.
A344606 counts alternating permutations of prime indices.
A347446 counts partitions with integer alternating product.
A347457 ranks partitions with integer alt product, complement A347455.
A347461 counts possible alternating products of partitions.
A347462 counts possible reverse-alternating products of partitions.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Select[Range[100],OddQ[PrimeOmega[#]]&&IntegerQ[altprod[primeMS[#]]]&]

A347451 Numbers whose multiset of prime indices has integer reciprocal alternating product.

Original entry on oeis.org

1, 2, 4, 6, 8, 9, 10, 14, 16, 18, 21, 22, 24, 25, 26, 32, 34, 36, 38, 39, 40, 46, 49, 50, 54, 56, 57, 58, 62, 64, 65, 72, 74, 81, 82, 84, 86, 87, 88, 90, 94, 96, 98, 100, 104, 106, 111, 115, 118, 121, 122, 126, 128, 129, 133, 134, 136, 142, 144, 146, 150, 152
Offset: 1

Views

Author

Gus Wiseman, Sep 24 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.
We define the reciprocal alternating product of a sequence (y_1,...,y_k) to be Product_i y_i^((-1)^i).
Also Heinz numbers integer partitions with integer reverse-reciprocal alternating product, where the Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			The terms and their prime indices begin:
      1: {}            32: {1,1,1,1,1}       65: {3,6}
      2: {1}           34: {1,7}             72: {1,1,1,2,2}
      4: {1,1}         36: {1,1,2,2}         74: {1,12}
      6: {1,2}         38: {1,8}             81: {2,2,2,2}
      8: {1,1,1}       39: {2,6}             82: {1,13}
      9: {2,2}         40: {1,1,1,3}         84: {1,1,2,4}
     10: {1,3}         46: {1,9}             86: {1,14}
     14: {1,4}         49: {4,4}             87: {2,10}
     16: {1,1,1,1}     50: {1,3,3}           88: {1,1,1,5}
     18: {1,2,2}       54: {1,2,2,2}         90: {1,2,2,3}
     21: {2,4}         56: {1,1,1,4}         94: {1,15}
     22: {1,5}         57: {2,8}             96: {1,1,1,1,1,2}
     24: {1,1,1,2}     58: {1,10}            98: {1,4,4}
     25: {3,3}         62: {1,11}           100: {1,1,3,3}
     26: {1,6}         64: {1,1,1,1,1,1}    104: {1,1,1,6}
		

Crossrefs

The version for reversed prime indices is A028982, counted by A119620.
The additive version is A119899, strict A028260.
Allowing any alternating product >= 1 gives A344609.
Factorizations of this type are counted by A347439.
Allowing any alternating product <= 1 gives A347450.
The non-reciprocal version is A347454.
Allowing any alternating product > 1 gives A347465, reverse A028983.
A056239 adds up prime indices, row sums of A112798.
A316524 gives the alternating sum of prime indices (reverse: A344616).
A335433 lists numbers whose prime indices are separable, complement A335448.
A344606 counts alternating permutations of prime indices.
A347457 ranks partitions with integer alternating product.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    altprod[q_]:=Product[q[[i]]^(-1)^(i-1),{i,Length[q]}];
    Select[Range[100],IntegerQ[1/altprod[primeMS[#]]]&]

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

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Dec 22 2019

Keywords

Comments

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)

Examples

			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
		

Crossrefs

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.

Programs

  • Mathematica
    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));
    
  • PARI
    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

Formula

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)
For all x in A102750, a(x) = a(A253553(x)). (End)

Extensions

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

Views

Author

Gus Wiseman, Aug 03 2021

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.

Examples

			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.
		

Crossrefs

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.

Programs

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

Formula

a(n) * A329888(n) = n.
A056239(a(n)) = A346699(n).
Showing 1-10 of 16 results. Next