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 14 results. Next

A344651 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with alternating sum k, with k ranging from n mod 2 to n in steps of 2.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 2, 2, 1, 4, 2, 1, 3, 5, 2, 1, 7, 5, 2, 1, 5, 9, 5, 2, 1, 12, 10, 5, 2, 1, 7, 17, 10, 5, 2, 1, 19, 19, 10, 5, 2, 1, 11, 28, 20, 10, 5, 2, 1, 30, 33, 20, 10, 5, 2, 1, 15, 47, 35, 20, 10, 5, 2, 1, 45, 57, 36, 20, 10, 5, 2, 1, 22, 73, 62, 36, 20, 10, 5, 2, 1
Offset: 0

Views

Author

Gus Wiseman, Jun 05 2021

Keywords

Comments

The alternating sum of a partition (y_1,...,y_k) is Sum_i (-1)^(i-1) y_i. This is equal to the number of odd parts in the conjugate partition, so T(n,k) is the number of integer partitions of n with k odd parts in the conjugate partition, which is also the number of partitions of n with k odd parts.
Also the number of integer partitions of n with odd-indexed parts (odd bisection) summing to k, ceiling(n/2) <= k <= n. The even-indexed version is A346633. - Gus Wiseman, Nov 29 2021

Examples

			Triangle begins:
   1
   1
   1   1
   2   1
   2   2   1
   4   2   1
   3   5   2   1
   7   5   2   1
   5   9   5   2   1
  12  10   5   2   1
   7  17  10   5   2   1
  19  19  10   5   2   1
  11  28  20  10   5   2   1
  30  33  20  10   5   2   1
  15  47  35  20  10   5   2   1
  45  57  36  20  10   5   2   1
  22  73  62  36  20  10   5   2   1
  67  92  64  36  20  10   5   2   1
  30 114 102  65  36  20  10   5   2   1
  97 147 107  65  36  20  10   5   2   1
Row n = 10 counts the following partitions (A = 10):
  (55)          (64)         (73)       (82)     (91)   (A)
  (3322)        (442)        (433)      (622)    (811)
  (4411)        (541)        (532)      (721)
  (222211)      (3331)       (631)      (7111)
  (331111)      (4222)       (5221)     (61111)
  (22111111)    (4321)       (6211)
  (1111111111)  (5311)       (42211)
                (22222)      (52111)
                (32221)      (511111)
                (33211)      (4111111)
                (43111)
                (322111)
                (421111)
                (2221111)
                (3211111)
                (31111111)
                (211111111)
The conjugate version is:
  (A)      (55)      (3331)     (331111)    (31111111)   (1111111111)
  (64)     (73)      (5311)     (511111)    (211111111)
  (82)     (91)      (7111)     (3211111)
  (442)    (433)     (33211)    (4111111)
  (622)    (532)     (43111)    (22111111)
  (4222)   (541)     (52111)
  (22222)  (631)     (61111)
           (721)     (322111)
           (811)     (421111)
           (3322)    (2221111)
           (4321)
           (4411)
           (5221)
           (6211)
           (32221)
           (42211)
           (222211)
		

Crossrefs

This is A103919 with all zeros removed.
The strict version is A152146 interleaved with A152157.
The rows are those of A239830 interleaved with those of A239829.
The reverse version is the right half of A344612.
The strict reverse version is the right half of A344739.
A000041 counts partitions of 2n with alternating sum 0, ranked by A000290.
A027187 counts partitions with rev-alternating sum <= 0, ranked by A028260.
A124754 lists alternating sums of standard compositions (reverse: A344618).
A316524 is the alternating sum of the prime indices of n (reverse: A344616).
A325534/A325535 count separable/inseparable partitions.
A344607 counts partitions with rev-alternating sum >= 0, ranked by A344609.
A344608 counts partitions with rev-alternating sum < 0, ranked by A119899.
A344610 counts partitions of n by positive rev-alternating sum.
A344611 counts partitions of 2n with rev-alternating sum >= 0.
A345197 counts compositions by sum, length, and alternating sum.
A346697 gives the sum of odd-indexed prime indices (reverse: A346699).
A346702 represents the odd bisection of compositions, sums A209281.

Programs

  • Mathematica
    ats[y_]:=Sum[(-1)^(i-1)*y[[i]],{i,Length[y]}];
    Table[Length[Select[IntegerPartitions[n],ats[#]==k&]],{n,0,15},{k,Mod[n,2],n,2}]

A346697 Sum of the odd-indexed parts (odd bisection) of the multiset of prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 01 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.

Examples

			The prime indices of 1100 are {1,1,3,3,5}, so a(1100) = 1 + 3 + 5 = 9.
The prime indices of 2100 are {1,1,2,3,3,4}, so a(2100) = 1 + 2 + 3 = 6.
		

Crossrefs

The version for standard compositions is A209281(n+1) (even: A346633).
Subtracting the even version gives A316524 (reverse: A344616).
The even version is A346698.
The reverse version is A346699.
The even reverse version is A346700.
A000120 and A080791 count binary digits 1 and 0, with difference A145037.
A000302 counts compositions with odd alternating sum, ranked by A053738.
A001414 adds up prime factors, row sums of A027746.
A029837 adds up parts of standard compositions (alternating: A124754).
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[First/@Partition[Append[primeMS[n],0],2]],{n,100}]

Formula

a(n) = A056239(n) - A346698(n).
a(n) = A316524(n) + A346698(n).
a(n odd omega) = A346699(n).
a(n even omega) = A346700(n).
A344616(n) = A346699(n) - A346700(n).

A346698 Sum of the even-indexed parts (even bisection) of the multiset of prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 01 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.

Examples

			The prime indices of 1100 are {1,1,3,3,5}, so a(1100) = 1 + 3 = 4.
The prime indices of 2100 are {1,1,2,3,3,4}, so a(2100) = 1 + 3 + 4 = 8.
		

Crossrefs

Subtracting from the odd version gives A316524 (reverse: A344616).
The version for standard compositions is A346633 (odd: A209281(n+1)).
The odd version is A346697.
The even reverse version is A346699.
The reverse version is A346700.
A000120 and A080791 count binary digits 1 and 0, with difference A145037.
A001414 adds up prime factors, row-sums of A027746.
A029837 adds up parts of standard compositions (alternating: A124754).
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[Append[primeMS[n],0],2]],{n,100}]
  • PARI
    A346698(n) = if(1==n,0,my(f=factor(n),s=0,p=0); for(k=1,#f~,while(f[k,2], s += (p%2)*primepi(f[k,1]); f[k,2]--; p++)); (s)); \\ Antti Karttunen, Nov 30 2021

Formula

a(n) = A056239(n) - A346697(n).
a(n) = A346697(n) - A316524(n).
a(n even omega) = A346699(n).
a(n odd omega) = A346700(n).
A344616(n) = A346699(n) - A346700(n).

Extensions

Data section extended up to 105 terms by Antti Karttunen, Nov 30 2021

A209281 Start with first run [0,1] then, for n >= 2, the n-th run has length 2^n and is the concatenation of [a(1),a(2),...,a(2^n/2)] and [n-a(1),n-a(2),...,n-a(2^n/2)].

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 1, 2, 4, 3, 2, 3, 1, 2, 3, 2, 5, 4, 3, 4, 2, 3, 4, 3, 1, 2, 3, 2, 4, 3, 2, 3, 6, 5, 4, 5, 3, 4, 5, 4, 2, 3, 4, 3, 5, 4, 3, 4, 1, 2, 3, 2, 4, 3, 2, 3, 5, 4, 3, 4, 2, 3, 4, 3, 7, 6, 5, 6, 4, 5, 6, 5, 3, 4, 5, 4, 6, 5, 4, 5, 2, 3, 4, 3, 5, 4, 3
Offset: 1

Views

Author

Benoit Cloitre, Jan 16 2013

Keywords

Comments

Also the sum of the odd bisection (odd-indexed parts) of the (n-1)-th composition in standard order, where the k-th composition in standard order (graded reverse-lexicographic, A066099) is obtained by taking the set of positions of 1's in the reversed binary expansion of k, prepending 0, taking first differences, and reversing again. Note that this sequence counts {} as composition number 1 (instead of the usual 0). For example, composition number 741 in standard order is (2,1,1,3,2,1), with odd bisection (2,1,2), so a(742) = 2 + 1 + 2 = 5. - Gus Wiseman, Aug 24 2021

Examples

			[0,1] -> [0,1] U [2-0,2-1] =
[0,1,2,1] -> [0,1,2,1] U [3-0,3-1,3-2,3-1] =
[0,1,2,1,3,2,1,2] etc.
From _Gus Wiseman_, Aug 08 2021: (Start)
As a triangle without the initial 0, row-lengths A000079:
  1
  2 1
  3 2 1 2
  4 3 2 3 1 2 3 2
  5 4 3 4 2 3 4 3 1 2 3 2 4 3 2 3
  6 5 4 5 3 4 5 4 2 3 4 3 5 4 3 4 1 2 3 2 4 3 2 3 5 4 3 4 2 3 4 3
(End)
		

Crossrefs

Cf. A010060 (Thue-Morse), A103204 (indices of 1's).
Cf. A029837 (binary order), A000120 (binary weight), A006068 (inverse Gray), A272020 (bit positions).
Cf. A089215.
As a triangle: A000079 (row lengths), A001792 (row sums).
Other composition part sums: A124754. A346633.
Also the sum of row A346702(n-1) of A066099.
Cf. A346697 (on prime indices).

Programs

  • Mathematica
    Table[Total[First/@Partition[Append[Differences[Prepend[Join@@Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse,0],2]],{n,0,100}] (* Gus Wiseman, Aug 08 2021 *)
  • PARI
    /* compute 2^15 terms */ v=[0,1];for(n=2,15,v=concat(v,vector(2^n/2,i,n-v[i]));a(n)=v[n];)
    
  • PARI
    a(n) = n--; my(s=1,ns); while((ns=n>>s), n=bitxor(n,ns); s<<=1); hammingweight(n); \\ Kevin Ryde, May 14 2022

Formula

Let T(n)=A010060(n) then for n>=1 a(2n)=a(n)+1-T(n-1) and a(2n+1)=a(n+1)+T(n).
For n>=2 a(n) = a(ceiling(n/2))+T(n-1) hence:
a(n) = Sum_{k=0..ceiling(log(n-1)/log(2))} T(floor((n-1)/2^k)).
For k>=0 a(3*2^k+1)=1 (more precisely a(n)=1 iff n is in A103204), a(2^k+1)=k+1, a(5*2^k+1)=2, a(7*2^k+1)=k+2 etc.
From Gus Wiseman, Aug 18 2021: (Start)
a(n + 1) = (A029837(n) + A124754(n))/2.
a(n + 1) = A029837(n) - A346633(n).
a(n + 1) = A346633(n) - A124754(n).
a(n + 1) = A029837(A346702(n)).
(End)
From Kevin Ryde, May 14 2022: (Start)
a(n) = A000120(A006068(n-1)), binary weight of inverse binary Gray code.
a(n) = Sum_{k=1..A000120(n-1)} (-1)^(k-1) * A272020(n-1,k), alternating sum of 1-bit positions.
a(n) = A089215(n) - 1.
(End)

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

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 08 2021

Keywords

Examples

			The prime factors of 108 are (2,2,3,3,3), with even bisection (2,3), with product 6, so a(108) = 6.
The prime factors of 720 are (2,2,2,2,3,3,5), with even bisection (2,2,3), with product 12, so a(720) = 12.
		

Crossrefs

Positions of first appearances are A129597.
Positions of 1's are A008578.
Positions of primes are A168645.
The sum of prime indices of a(n) is A346698(n).
The odd version is A346703 (sum: A346697).
The odd reverse version is A346701 (sum: A346699).
The reverse version appears to be A329888 (sum: A346700).
A001221 counts distinct prime factors.
A001222 counts all prime factors.
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).
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.

Programs

  • Maple
    f:= proc(n) local F,i;
      F:= ifactors(n)[2];
      F:= sort(map(t -> t[1]$t[2],F));
      mul(F[i],i=2..nops(F),2)
    end proc:
    map(f, [$1..100]); # Robert Israel, Aug 12 2024
  • Mathematica
    Table[Times@@Last/@Partition[Flatten[Apply[ConstantArray,FactorInteger[n],{1}]],2],{n,100}]

Formula

a(n) * A346703(n) = n.
A056239(a(n)) = A346698(n).

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

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

Original entry on oeis.org

0, 1, 2, 1, 3, 2, 4, 2, 2, 3, 5, 3, 6, 4, 3, 2, 7, 3, 8, 4, 4, 5, 9, 3, 3, 6, 4, 5, 10, 4, 11, 3, 5, 7, 4, 3, 12, 8, 6, 4, 13, 5, 14, 6, 5, 9, 15, 4, 4, 4, 7, 7, 16, 4, 5, 5, 8, 10, 17, 4, 18, 11, 6, 3, 6, 6, 19, 8, 9, 5, 20, 4, 21, 12, 5, 9, 5, 7, 22, 5, 4
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 with Heinz number 1100 is (5,3,3,1,1), so a(1100) = 5 + 3 + 1 = 9.
The partition with Heinz number 2100 is (4,3,3,2,1,1), so a(2100) = 4 + 3 + 1 = 8.
		

Crossrefs

The version for standard compositions is A209281(n+1) (even: A346633).
Subtracting the even version gives A344616 (non-reverse: A316524).
The even version is A346700.
The non-reverse version (multisets instead of partitions) is A346697.
The even non-reverse version is A346698.
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[First/@Partition[Append[Reverse[primeMS[n]],0],2]],{n,100}]

Formula

a(n) = A056239(A346701(n)).
a(n) = A056239(n) - A346700(n).
a(n) = A344616(n) + A346700(n).
a(n odd omega) = A346697(n).
a(n even omega) = A346698(n).
A316524(n) = A346697(n) - A346698(n).

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

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

Views

Author

Rémy Sigrist, Apr 02 2021

Keywords

Comments

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

Examples

			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.
		

Crossrefs

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).

Programs

  • Mathematica
    Table[n^2/FactorInteger[n][[-1,1]],{n,100}] (* Gus Wiseman, Aug 09 2021 *)
  • PARI
    See Links section.

Formula

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)
Showing 1-10 of 14 results. Next