A358136 -id:A358136 - cp's OEIS frontend

A048675 If n = p_i^e_i * ... * p_k^e_k, p_i < ... < p_k primes (with p_i = prime(i)), then a(n) = (1/2) * (e_i * 2^i + ... + e_k * 2^k).

0, 1, 2, 2, 4, 3, 8, 3, 4, 5, 16, 4, 32, 9, 6, 4, 64, 5, 128, 6, 10, 17, 256, 5, 8, 33, 6, 10, 512, 7, 1024, 5, 18, 65, 12, 6, 2048, 129, 34, 7, 4096, 11, 8192, 18, 8, 257, 16384, 6, 16, 9, 66, 34, 32768, 7, 20, 11, 130, 513, 65536, 8, 131072, 1025, 12, 6, 36, 19

Offset: 1

Antti Karttunen, Jul 14 1999

nonn

The original motivation for this sequence was to encode the prime factorization of n in the binary representation of a(n), each such representation being unique as long as this map is restricted to A005117 (squarefree numbers, resulting a permutation of nonnegative integers A048672) or any of its subsequence, resulting an injective function like A048623 and A048639.
However, also the restriction to A260443 (not all terms of which are squarefree) results a permutation of nonnegative integers, namely A001477, the identity permutation.
When a polynomial with nonnegative integer coefficients is encoded with the prime factorization of n (e.g., as in A206296, A260443), then a(n) gives the evaluation of that polynomial at x=2.
The primitive completely additive integer sequence that satisfies a(n) = a(A225546(n)), n >= 1. By primitive, we mean that if b is another such sequence, then there is an integer k such that b(n) = k * a(n) for all n >= 1. - Peter Munn, Feb 03 2020
If the binary rank of an integer partition y is given by Sum_i 2^(y_i-1), and the Heinz number is Product_i prime(y_i), then a(n) is the binary rank of the integer partition with Heinz number n. Note the function taking a set s to Sum_i 2^(s_i-1) is the inverse of A048793 (binary indices), and the function taking a multiset m to Product_i prime(m_i) is the inverse of A112798 (prime indices). - Gus Wiseman, May 22 2024

			From _Gus Wiseman_, May 22 2024: (Start)
The A018819(7) = 6 cases of binary rank 7 are the following, together with their prime indices:
   30: {1,2,3}
   40: {1,1,1,3}
   54: {1,2,2,2}
   72: {1,1,1,2,2}
   96: {1,1,1,1,1,2}
  128: {1,1,1,1,1,1,1}
(End)

Antti Karttunen, Table of n, a(n) for n = 1..1024
Hans Havermann, Factorization of the first 10000 terms, in format [[primes], [exponents]]

Row 2 of A104244.
Similar logarithmic functions: A001414, A056239, A090880, A289506, A293447.
Left inverse of the following sequences: A000079, A019565, A038754, A068911, A134683, A260443, A332824.
A003961, A028234, A032742, A055396, A064989, A067029, A225546, A297845 are used to express relationship between terms of this sequence.
Cf. also A048623, A048676, A099884, A277896 and tables A277905, A285325.
Cf. A297108 (Möbius transform), A332813 and A332823 [= a(n) mod 3].
Pairs of sequences (f,g) that satisfy a(f(n)) = g(n), possibly with offset change: (A000203,A331750), (A005940,A087808), (A007913,A248663), (A007947,A087207), (A097248,A048675), (A206296,A000129), (A248692,A056239), (A283477,A005187), (A284003,A006068), (A285101,A028362), (A285102,A068052), (A293214,A001065), (A318834,A051953), (A319991,A293897), (A319992,A293898), (A320017,A318674), (A329352,A069359), (A332461,A156552), (A332462,A156552), (A332825,A000010) and apparently (A163511,A135529).
See comments/formulas in A277333, A331591, A331740 giving their relationship to this sequence.
The formula section details how the sequence maps the terms of A329050, A329332.
A277892, A322812, A322869, A324573, A324575 give properties of the n-th term of this sequence.
The term k appears A018819(k) times.
The inverse transformation is A019565 (Heinz number of binary indices).
The version for distinct prime indices is A087207.
Numbers k such that a(k) is prime are A277319, counts A372688.
Grouping by image gives A277905.
A014499 lists binary indices of prime numbers.
A061395 gives greatest prime index, least A055396.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.
Binary indices:
- listed A048793, sum A029931
- reversed A272020
- opposite A371572, sum A230877
- length A000120, complement A023416
- min A001511, opposite A000012
- max A070939, opposite A070940
- complement A368494, sum A359400
- opposite complement A371571, sum A359359
Cf. A000720, A035100, A071814, A096111, A225620, A243055, A304818, A355536, A358136, A372890.
Cf. A087207, A097248.

nthprime := proc(n) local i; if(isprime(n)) then for i from 1 to 1000000 do if(ithprime(i) = n) then RETURN(i); fi; od; else RETURN(0); fi; end; # nthprime(2) = 1, nthprime(3) = 2, nthprime(5) = 3, etc. - this is also A049084.
A048675 := proc(n) local s,d; s := 0; for d in ifactors(n)[ 2 ] do s := s + d[ 2 ]*(2^(nthprime(d[ 1 ])-1)); od; RETURN(s); end;
# simpler alternative
f:= n -> add(2^(numtheory:-pi(t[1])-1)*t[2], t=ifactors(n)[2]):
map(f, [$1..100]); # Robert Israel, Oct 10 2016

a[1] = 0; a[n_] := Total[ #[[2]]*2^(PrimePi[#[[1]]]-1)& /@ FactorInteger[n] ]; Array[a, 100] (* Jean-François Alcover, Mar 15 2016 *)

a(n) = my(f = factor(n)); sum(k=1, #f~, f[k,2]*2^primepi(f[k,1]))/2; \\ Michel Marcus, Oct 10 2016

\\ The following program reconstructs terms (e.g. for checking purposes) from the factorization file prepared by Hans Havermann:
v048675sigs = readvec("a048675.txt");
A048675(n) = if(n<=2,n-1,my(prsig=v048675sigs[n],ps=prsig[1],es=prsig[2]); prod(i=1,#ps,ps[i]^es[i])); \\ Antti Karttunen, Feb 02 2020

from sympy import factorint, primepi
def a(n):
    if n==1: return 0
    f=factorint(n)
    return sum([f[i]*2**(primepi(i) - 1) for i in f])
print([a(n) for n in range(1, 51)]) # Indranil Ghosh, Jun 19 2017

a(1) = 0, a(n) = 1/2 * (e1*2^i1 + e2*2^i2 + ... + ez*2^iz) if n = p_{i1}^e1*p_{i2}^e2*...*p_{iz}^ez, where p_i is the i-th prime. (e.g. p_1 = 2, p_2 = 3).
Totally additive with a(p^e) = e * 2^(PrimePi(p)-1), where PrimePi(n) = A000720(n). [Missing factor e added to the comment by Antti Karttunen, Jul 29 2015]
From Antti Karttunen, Jul 29 2015: (Start)
a(1) = 0; for n > 1, a(n) = 2^(A055396(n)-1) + a(A032742(n)). [Where A055396(n) gives the index of the smallest prime dividing n and A032742(n) gives the largest proper divisor of n.]
a(1) = 0; for n > 1, a(n) = (A067029(n) * (2^(A055396(n)-1))) + a(A028234(n)).
Other identities. For all n >= 0:
a(A019565(n)) = n.
a(A260443(n)) = n.
a(A206296(n)) = A000129(n).
a(A005940(n+1)) = A087808(n).
a(A007913(n)) = A248663(n).
a(A007947(n)) = A087207(n).
a(A283477(n)) = A005187(n).
a(A284003(n)) = A006068(n).
a(A285101(n)) = A028362(1+n).
a(A285102(n)) = A068052(n).
Also, it seems that a(A163511(n)) = A135529(n) for n >= 1. (End)
a(1) = 0, a(2n) = 1+a(n), a(2n+1) = 2*a(A064989(2n+1)). - Antti Karttunen, Oct 11 2016
From Peter Munn, Jan 31 2020: (Start)
a(n^2) = a(A003961(n)) = 2 * a(n).
a(A297845(n,k)) = a(n) * a(k).
a(n) = a(A225546(n)).
a(A329332(n,k)) = n * k.
a(A329050(n,k)) = 2^(n+k).
(End)
From Antti Karttunen, Feb 02-25 2020, Feb 01 2021: (Start)
a(n) = Sum_{d|n} A297108(d) = Sum_{d|A225546(n)} A297108(d).
a(n) = a(A097248(n)).
For n >= 2:
A001221(a(n)) = A322812(n), A001222(a(n)) = A277892(n).
A000203(a(n)) = A324573(n), A033879(a(n)) = A324575(n).
For n >= 1, A331750(n) = a(A000203(n)).
For n >= 1, the following chains hold:
A293447(n) >= a(n) >= A331740(n) >= A331591(n).
a(n) >= A087207(n) >= A248663(n).
(End)
a(n) = A087207(A097248(n)). - Flávio V. Fernandes, Jul 16 2025

Entry revised by Antti Karttunen, Jul 29 2015

More linking formulas added by Antti Karttunen, Apr 18 2017

A358194 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with partial sums summing to k, where k ranges from n to n(n+1)/2.

Original entry on oeis.org

1, 1, 1, 1, 1, 0, 1, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 2, 1, 0, 1, 1, 1, 1, 1, 0, 1, 1, 0, 1, 1, 0, 1, 1

Offset: 0

Gus Wiseman, Dec 31 2022

The partial sums of a sequence (a, b, c, ...) are (a, a+b, a+b+c, ...).

			Triangle begins:
  1
  1
  1 1
  1 0 1 1
  1 0 1 1 0 1 1
  1 0 0 1 1 0 1 1 0 1 1
  1 0 0 1 1 1 1 0 1 1 0 1 1 0 1 1
  1 0 0 0 1 1 1 0 1 1 1 1 1 0 1 1 0 1 1 0 1 1
  1 0 0 0 1 1 1 1 0 1 1 1 2 1 0 1 1 1 1 1 0 1 1 0 1 1 0 1 1
For example, the T(15,59) = 5 partitions are: (8,2,2,2,1), (7,3,3,1,1), (6,5,2,1,1), (4,3,2,2,2,2), (3,3,3,3,2,1).

Row sums are A000041.
The version for compositions is A053632.
Row lengths are A152947.
The version for reversed partitions is A264034.
A048793 = partial sums of reversed standard compositions, sum A029931.
A358134 = partial sums of standard compositions, sum A359042.
A358136 = partial sums of prime indices, sum A318283.
A359361 = partial sums of reversed prime indices, sum A304818.
Cf. A000009, A325362, A358137, A359397.

Table[Length[Select[IntegerPartitions[n],Total[Accumulate[#]]==k&]],{n,0,8},{k,n,n*(n+1)/2}]

A358137 Heinz number of the partial sums of the prime indices of n.

Original entry on oeis.org

1, 2, 3, 6, 5, 10, 7, 30, 21, 14, 11, 42, 13, 22, 33, 210, 17, 110, 19, 66, 39, 26, 23, 330, 65, 34, 273, 78, 29, 130, 31, 2310, 51, 38, 85, 546, 37, 46, 57, 390, 41, 170, 43, 102, 357, 58, 47, 2730, 133, 238, 69, 114, 53, 1870, 95, 510, 87, 62, 59, 714, 61

Offset: 1

Gus Wiseman, Oct 31 2022

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 terms together with their prime indices begin:
      1: {}
      2: {1}
      3: {2}
      6: {1,2}
      5: {3}
     10: {1,3}
      7: {4}
     30: {1,2,3}
     21: {2,4}
     14: {1,4}
     11: {5}
     42: {1,2,4}
     13: {6}
     22: {1,5}
     33: {2,5}
    210: {1,2,3,4}
     17: {7}
    110: {1,3,5}

The sorted version is A325362.
The prime indices are rows of A358136, partial sums of rows of A112798.
A000040 lists the primes.
A000041 counts partitions, strict A000009.
A003963 multiplies prime indices.
A056239 adds up prime indices.
Cf. A000720, A001221, A001222, A215366, A296150, A318283, A355536, A358134.

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

A001222(a(n)) = A001222(n).

A359361 Irregular triangle read by rows whose n-th row lists the partial sums of the integer partition with Heinz number n.

Original entry on oeis.org

1, 2, 1, 2, 3, 2, 3, 4, 1, 2, 3, 2, 4, 3, 4, 5, 2, 3, 4, 6, 4, 5, 3, 5, 1, 2, 3, 4, 7, 2, 4, 5, 8, 3, 4, 5, 4, 6, 5, 6, 9, 2, 3, 4, 5, 3, 6, 6, 7, 2, 4, 6, 4, 5, 6, 10, 3, 5, 6, 11, 1, 2, 3, 4, 5, 5, 7, 7, 8, 4, 7, 2, 4, 5, 6, 12, 8, 9, 6, 8, 3, 4, 5, 6, 13

Offset: 2

Gus Wiseman, Dec 30 2022

The partial sums of a sequence (a, b, c, ...) are (a, a+b, a+b+c, ...).

The Heinz number of a partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k). The partition with Heinz number n is the reversed n-th row of A112798.

			Triangle begins:
   2: 1
   3: 2
   4: 1 2
   5: 3
   6: 2 3
   7: 4
   8: 1 2 3
   9: 2 4
  10: 3 4
  11: 5
  12: 2 3 4
  13: 6
  14: 4 5
  15: 3 5
  16: 1 2 3 4
For example, the integer partition with Heinz number 90 is (3,2,2,1), so row n = 90 is (3,5,7,8).

Row-lengths are A001222.
The version for standard compositions is A048793, non-reversed A358134.
Last element in each row is A056239.
First element in each row is A061395
Rows are the partial sums of rows of A296150.
Row-sums are A304818.
A reverse version is A358136, row sums A318283, Heinz numbers A358137.
The sorted Heinz numbers of rows are A359397.
A000041 counts partitions, strict A000009.
A112798 lists prime indices, product A003963.
A355536 lists differences of prime indices.
Cf. A000720, A001221, A055396, A261079, A325362.

T:= n-> ListTools[PartialSums](sort([seq(numtheory
       [pi](i[1])$i[2], i=ifactors(n)[2])], `>`))[]:
seq(T(n), n=2..50);  # Alois P. Heinz, Jan 01 2023

Table[Accumulate[Reverse[Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]]],{n,2,30}]

A359755 Positions of first appearances in the sequence of weighted sums of prime indices (A304818).

Original entry on oeis.org

1, 2, 3, 4, 6, 7, 8, 10, 12, 15, 16, 18, 20, 24, 26, 28, 36, 40, 46, 48, 50, 52, 56, 62, 68, 74, 76, 86, 88, 92, 94, 106, 107, 118, 122, 124, 131, 134, 136, 142, 146, 152, 158, 164, 166, 173, 178, 188, 193, 194, 199, 202, 206, 214, 218, 226, 229, 236, 239, 254

Offset: 1

Gus Wiseman, Jan 15 2023

nonn

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 weighted sum of a sequence (y_1,...,y_k) is Sum_{i=1..k} i*y_i.

			The terms together with their prime indices begin:
    1: {}
    2: {1}
    3: {2}
    4: {1,1}
    6: {1,2}
    7: {4}
    8: {1,1,1}
   10: {1,3}
   12: {1,1,2}
   15: {2,3}
   16: {1,1,1,1}
   18: {1,2,2}
   20: {1,1,3}
   24: {1,1,1,2}

The version for standard compositions is A089633, zero-based A359756.
Positions of first appearances in A304818, reverse A318283.
The zero-based version is A359675, unsorted A359676.
The reverse zero-based version is A359680, unsorted A359681.
This is the sorted version of A359682, reverse A359679.
The reverse version is A359754.
A053632 counts compositions by weighted sum.
A112798 lists prime indices, length A001222, sum A056239.
A320387 counts multisets by weighted sum, zero-based A359678.
A358136 lists partial sums of prime indices, ranked by A358137, rev A359361.
Cf. A029931, A124757, A243055, A358194, A359497, A359674, A359683.

nn=1000;
primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
ots[y_]:=Sum[i*y[[i]],{i,Length[y]}];
seq=Table[ots[primeMS[n]],{n,1,nn}];
Select[Range[nn],FreeQ[seq[[Range[#-1]]],seq[[#]]]&]

A242628 Irregular table enumerating partitions; n-th row has partitions in previous row with each part incremented, followed by partitions in previous row with an additional part of size 1.

Original entry on oeis.org

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

Offset: 1

Franklin T. Adams-Watters, May 19 2014

This can be calculated using the binary expansion of n; see the PARI program.
The n-th row consists of all partitions with hook size (maximum + number of parts - 1) equal to n.
The partitions in row n of this sequence are the conjugates of the partitions in row n of A125106 taken in reverse order.
Row n is also the reversed partial sums plus one of the n-th composition in standard order (A066099) minus one. - Gus Wiseman, Nov 07 2022

			The table starts:
  1;
  2; 1,1;
  3; 2,2; 2,1; 1,1,1;
  4; 3,3; 3,2; 2,2,2; 3,1 2,2,1 2,1,1 1,1,1,1;
  ...

Alois P. Heinz, Rows n = 1..12, flattened
Gus Wiseman, Statistics, classes, and transformations of standard compositions

Cf. A241596 (another version of this list of partitions), A125106, A240837, A112531, A241597 (compositions).
For other schemes to list integer partitions, please see for example A227739, A112798, A241918, A114994.
First element in each row is A008687.
Last element in each row is A065120.
Heinz numbers of rows are A253565.
Another version is A358134.
Cf. A000120, A029837, A029931, A030303, A066099, A070939, A253566, A355536, A358135, A358136.

b:= proc(n) option remember; `if`(n=1, [[1]],
      [map(x-> map(y-> y+1, x), b(n-1))[],
       map(x-> [x[], 1], b(n-1))[]])
    end:
T:= n-> map(x-> x[], b(n))[]:
seq(T(n), n=1..7);  # Alois P. Heinz, Sep 25 2015

T[1] = {{1}};
T[n_] := T[n] = Join[T[n-1]+1, Append[#, 1]& /@ T[n-1]];
Array[T, 7] // Flatten (* Jean-François Alcover, Jan 25 2021 *)

apart(n) = local(r=[1]); while(n>1,if(n%2==0,for(k=1,#r,r[k]++),r=concat(r,[1]));n\=2);r \\ Generates n-th partition.

A359042 Sum of partial sums of the n-th composition in standard order (A066099).

Original entry on oeis.org

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

Offset: 0

Gus Wiseman, Dec 20 2022

nonn

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. This gives a bijective correspondence between nonnegative integers and integer compositions.

			The 29th composition in standard order is (1,1,2,1), with partial sums (1,2,4,5), with sum 12, so a(29) = 12.

Gus Wiseman, Statistics, classes, and transformations of standard compositions

See link for sequences related to standard compositions.
Each n appears A000009(n) times.
The reverse version is A029931.
Comps counted by this statistic are A053632, ptns A264034, rev ptns A358194.
This is the sum of partial sums of rows of A066099.
The version for Heinz numbers of partitions is A318283, row sums of A358136.
Row sums of A358134.
A011782 counts compositions.
A065120 gives first part of standard compositions, last A001511.
A242628 lists adjusted partial sums, ranked by A253565, row sums A359043.
A358135 gives last minus first of standard compositions.
Cf. A000120, A029837, A070939, A133494, A253566, A358133, A358137.

stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
Table[Total[Accumulate[stc[n]]],{n,0,100}]

A358134 Triangle read by rows whose n-th row lists the partial sums of the n-th composition in standard order (row n of A066099).

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Oct 31 2022

			Triangle begins:
  1
  2
  1 2
  3
  2 3
  1 3
  1 2 3
  4
  3 4
  2 4
  2 3 4
  1 4
  1 3 4
  1 2 4
  1 2 3 4

Gus Wiseman, Statistics, classes, and transformations of standard compositions

See link for sequences related to standard compositions.
First element in each row is A065120.
Rows are the partial sums of rows of A066099.
Last element in each row is A070939.
An adjusted version is A242628, ranked by A253565.
The first differences instead of partial sums are A358133.
The version for Heinz numbers of partitions is A358136, ranked by A358137.
Row sums are A359042.
A011782 counts compositions.
A351014 counts distinct runs in standard compositions.
A358135 gives last minus first of standard compositions.
Cf. A000120, A001511, A029837, A029931, A048896, A070939, A133494, A253566, A355536, A357135.

stc[n_]:=Differences[Prepend[Join @@ Position[Reverse[IntegerDigits[n,2]],1],0]]//Reverse;
Join@@Table[Accumulate[stc[n]],{n,100}]

A359674 Zero-based weighted sum of the prime indices of n in weakly increasing order.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jan 13 2023

nonn

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 zero-based weighted sum of a sequence (y_1,...,y_k) is Sum_{i=1..k} (i-1)*y_i.

			The prime indices of 12 are {1,1,2}, so a(12) = 0*1 + 1*1 + 2*2 = 5.

Positions of last appearances (except 0) are A001248.
Positions of 0's are A008578.
The version for standard compositions is A124757, reverse A231204.
The one-based version is A304818, reverse A318283.
Positions of first appearances are A359675, reverse A359680.
First position of n is A359676(n), reverse A359681.
The reverse version is A359677, firsts A359679.
Number of appearances of positive n is A359678(n).
A053632 counts compositions by zero-based weighted sum.
A112798 lists prime indices, length A001222, sum A056239.
A358136 lists partial sums of prime indices, ranked by A358137, rev A359361.
Cf. A029931, A055932, A243055, A359043, A358194, A359360, A359497, A359681, A359682, A359755.

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

A261079 Sum of index differences between prime factors of n, summed over all unordered pairs of primes present (with multiplicity) in the prime factorization of n.

Original entry on oeis.org

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

Offset: 1

Antti Karttunen, Sep 23 2015

nonn

			For n = 1 the prime factorization is empty, thus there is nothing to sum, so a(1) = 0.
For n = 6 = 2*3 = prime(1) * prime(2), a(6) = 1 because the (absolute value of) difference between prime indices of 2 and 3 is 1.
For n = 10 = 2*5 = prime(1) * prime(3), a(10) = 2 because the difference between prime indices of 2 and 5 is 2.
For n = 12 = 2*2*3 = prime(1) * prime(1) * prime(2), a(12) = 2 because the difference between prime indices of 2 and 3 is 1, and the pair (2,3) occurs twice as one can pick either one of the two 2's present in the prime factorization to be a pair of a single 3. Note that the index difference between 2 and 2 is 0, thus the pair (2,2) of prime divisors does not contribute to the sum.
For n = 36 = 2*2*3*3, a(36) = 4 because the index difference between 2 and 3 is 1, and the prime factor pair (2,3) occurs 2^2 = four times in total. As the index difference is zero between 2 and 2 as well as between 3 and 3, the pairs (2,2) and (3,3) do not contribute to the sum.

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

Cf. A000720.
Cf. A000961 (positions of zeros), A006094 (positions of ones).
Cf. A243055, A242411.
Cf. also A260737.
A055396 gives minimum prime index, maximum A061395.
A112798 list prime indices, length A001222, sum A056239.
A304818 adds up partial sums of reversed prime indices, row sums of A359361.
A318283 adds up partial sums of prime indices, row sums of A358136.
Cf. A029931, A316413, A325362, A326836, A326844, A355536, A358137, A359362.

Table[Function[p, Total@ Map[Function[b, Times @@ {First@ Differences@ PrimePi@ b, Count[Subsets[p, {2}], c_ /; SameQ[c, b]]}], Subsets[Union@ p, {2}]]][Flatten@ Replace[FactorInteger@ n, {p_, e_} :> ConstantArray[p, e], 2]], {n, 120}] (* Michael De Vlieger, Mar 08 2017 *)

a(n) = A304818(n) - A318283(n). - Gus Wiseman, Jan 09 2023

a(n) = 2*A304818(n) - A359362(n). - Gus Wiseman, Jan 09 2023

cp's OEIS Frontend

A048675 If n = p_i^e_i * ... * p_k^e_k, p_i < ... < p_k primes (with p_i = prime(i)), then a(n) = (1/2) * (e_i * 2^i + ... + e_k * 2^k).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

PARI

Python

Formula

Extensions

A358194 Irregular triangle read by rows where T(n,k) is the number of integer partitions of n with partial sums summing to k, where k ranges from n to n(n+1)/2.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A358137 Heinz number of the partial sums of the prime indices of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A359361 Irregular triangle read by rows whose n-th row lists the partial sums of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Maple

Mathematica

A359755 Positions of first appearances in the sequence of weighted sums of prime indices (A304818).

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A242628 Irregular table enumerating partitions; n-th row has partitions in previous row with each part incremented, followed by partitions in previous row with an additional part of size 1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

A359042 Sum of partial sums of the n-th composition in standard order (A066099).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

A358134 Triangle read by rows whose n-th row lists the partial sums of the n-th composition in standard order (row n of A066099).

Views

Author