A304793 -id:A304793 - cp's OEIS frontend

A319695 Number of distinct values obtained when Euler phi (A000010) is applied to proper divisors of n.

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

Offset: 1

Antti Karttunen, Oct 02 2018

nonn

			For n = 6, it has three proper divisors: 1, 2, 3, and applying A000010 to these gives 1, 1 and 2, with just two distinct values, thus a(6) = 2.

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

Cf. A000010, A319696.

Cf. also A304793, A305611, A316555, A316556, A319685 for similarly constructed sequences.

A319695(n) = { my(m=Map(),s,k=0); fordiv(n,d,if((d

A325801 Number of divisors of n minus the number of distinct positive subset-sums of the prime indices of n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, May 23 2019

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, with sum A056239(n). A positive subset-sum of an integer partition is any sum of a nonempty submultiset of it.

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

Positions of 0's are A299702.
Positions of 1's are A325802.
Positions of positive integers are A299729.
Cf. A000005, A002033, A056239, A108917, A112798, A276024, A304793.
Cf. A325694, A325780, A325781, A325792, A325793, A325799.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p] k]];
Table[DivisorSigma[0,n]-Length[Union[hwt/@Divisors[n]]],{n,100}]

A325801(n) = (numdiv(n) - A299701(n));
A299701(n) = { my(f = factor(n), pids = List([])); for(i=1,#f~, while(f[i,2], f[i,2]--; listput(pids,primepi(f[i,1])))); pids = Vec(pids); my(sv=vector(vecsum(pids))); for(b=1,(2^length(pids))-1,sv[sumbybits(pids,b)] = 1); 1+vecsum(sv); }; \\ Not really an optimal way to count these.
sumbybits(v,b) = { my(s=0,i=1); while(b>0,s += (b%2)*v[i]; i++; b >>= 1); (s); }; \\ Antti Karttunen, May 26 2019

a(n) = A000005(n) - A299701(n).

A325986 Heinz numbers of complete strict integer partitions.

Original entry on oeis.org

1, 2, 6, 30, 42, 210, 330, 390, 462, 510, 546, 714, 798, 2310, 2730, 3570, 3990, 4290, 4830, 5610, 6006, 6090, 6270, 6510, 6630, 7410, 7590, 7854, 8778, 8970, 9282, 9570, 9690, 10230, 10374, 10626, 11310, 11730, 12090, 12210, 12558, 13398, 13566, 14322, 14430

Offset: 1

Gus Wiseman, May 30 2019

nonn

Strict partitions are counted by A000009, while complete partitions are counted by A126796.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
An integer partition of n is complete (A126796, A325781) if every number from 0 to n is the sum of some submultiset of the parts.
The enumeration of these partitions by sum is given by A188431.

			The sequence of terms together with their prime indices begins:
      1: {}
      2: {1}
      6: {1,2}
     30: {1,2,3}
     42: {1,2,4}
    210: {1,2,3,4}
    330: {1,2,3,5}
    390: {1,2,3,6}
    462: {1,2,4,5}
    510: {1,2,3,7}
    546: {1,2,4,6}
    714: {1,2,4,7}
    798: {1,2,4,8}
   2310: {1,2,3,4,5}
   2730: {1,2,3,4,6}
   3570: {1,2,3,4,7}
   3990: {1,2,3,4,8}
   4290: {1,2,3,5,6}
   4830: {1,2,3,4,9}
   5610: {1,2,3,5,7}

Cf. A002033, A056239, A103295, A112798, A126796, A188431, A299702, A304793.

Cf. A325780, A325781, A325782, A325788, A325790, A325988.

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p] k]];
Select[Range[1000],SquareFreeQ[#]&&Union[hwt/@Divisors[#]]==Range[0,hwt[#]]&]

Intersection of A005117 (strict partitions) and A325781 (complete partitions).

A326018 Heinz numbers of knapsack partitions such that no addition of one part up to the maximum is knapsack.

Original entry on oeis.org

1925, 12155, 20995, 23375, 37145

Offset: 1

Gus Wiseman, Jun 03 2019

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
An integer partition is knapsack if every submultiset has a different sum.
The enumeration of these partitions by sum is given by A326016.

			The sequence of terms together with their prime indices begins:
   1925: {3,3,4,5}
  12155: {3,5,6,7}
  20995: {3,6,7,8}
  23375: {3,3,3,5,7}
  37145: {3,7,8,9}

Cf. A002033, A108917, A275972, A299702, A299729, A304793.

Cf. A325780, A325782, A325857, A325862, A325878, A325880, A326015, A326016.

ksQ[y_]:=UnsameQ@@Total/@Union[Subsets[y]];
Select[Range[2,200],With[{phm=If[#==1,{},Flatten[Cases[FactorInteger[#],{p_,k_}:>Table[PrimePi[p],{k}]]]]},ksQ[phm]&&Select[Table[Sort[Append[phm,i]],{i,Max@@phm}],ksQ]=={}]&]

A365919 Heinz numbers of integer partitions with the same number of distinct positive subset-sums as distinct non-subset-sums.

Original entry on oeis.org

1, 3, 9, 21, 22, 27, 63, 76, 81, 117, 147, 175, 186, 189, 243, 248, 273, 286, 290, 322, 345, 351, 399, 418, 441, 513, 516, 567, 688, 715, 729, 819, 1029, 1053, 1062, 1156, 1180, 1197, 1323, 1375, 1416, 1484, 1521, 1539, 1701, 1827, 1888, 1911, 2068, 2115, 2130

Offset: 1

Gus Wiseman, Sep 25 2023

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: {}
     3: {2}
     9: {2,2}
    21: {2,4}
    22: {1,5}
    27: {2,2,2}
    63: {2,2,4}
    76: {1,1,8}
    81: {2,2,2,2}
   117: {2,2,6}
   147: {2,4,4}
   175: {3,3,4}
   186: {1,2,11}
   189: {2,2,2,4}
   243: {2,2,2,2,2}

The LHS is A304793, counted by A365658, with empty sets A299701.
The RHS is A325799, counted by A365923 (strict A365545).
A046663 counts partitions without a subset summing to k, strict A365663.
A056239 adds up prime indices, row sums of A112798.
A276024 counts positive subset-sums of partitions, strict A284640.
A325781 ranks complete partitions, counted by A126796.
A365830 ranks incomplete partitions, counted by A365924.
A365918 counts non-subset-sums of partitions, strict A365922.
Cf. A001223, A005117, A006827, A073491, A188431, A304792, A365831.

prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
smu[y_]:=Union[Total/@Rest[Subsets[y]]];
nmz[y_]:=Complement[Range[Total[y]],Total/@Subsets[y]];
Select[Range[100],Length[smu[prix[#]]]==Length[nmz[prix[#]]]&]

Positive integers k such that A304793(k) = A325799(k).

A316555 Number of distinct GCDs of nonempty submultisets of the integer partition with Heinz number n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, Jul 06 2018

nonn

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Number of distinct values obtained when A289508 is applied to all divisors of n larger than one. - Antti Karttunen, Sep 28 2018

			455 is the Heinz number of (6,4,3) which has possible GCDs of nonempty submultisets {1,2,3,4,6} so a(455) = 5.

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

Cf. A056239, A108917, A122768, A275972, A289508, A289509, A296150, A316313, A316314, A316430, A316556, A316557.

Cf. also A304793, A305611, A319685, A319695 for other similarly constructed sequences.

Table[Length[Union[GCD@@@Rest[Subsets[If[n==1,{},Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]]]]]]],{n,100}]

A289508(n) = gcd(apply(p->primepi(p),factor(n)[,1]));
A316555(n) = { my(m=Map(),s,k=0); fordiv(n,d,if((d>1)&&!mapisdefined(m,s=A289508(d)), mapput(m,s,s); k++)); (k); }; \\ Antti Karttunen, Sep 28 2018

More terms from Antti Karttunen, Sep 28 2018

A319353 Filter sequence combining weights of partitions with Heinz-numbers equal to the proper divisors of n.

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 4, 5, 6, 2, 7, 2, 8, 9, 10, 2, 11, 2, 11, 12, 13, 2, 14, 15, 16, 12, 17, 2, 18, 2, 11, 19, 20, 21, 22, 2, 23, 24, 18, 2, 25, 2, 26, 27, 28, 2, 29, 30, 31, 32, 33, 2, 25, 34, 25, 35, 36, 2, 37, 2, 38, 39, 40, 41, 42, 2, 43, 44, 45, 2, 37, 2, 46, 47, 48, 49, 50, 2, 51, 52, 53, 2, 54, 55, 56, 57, 42, 2, 58, 59, 60, 61, 62, 63, 37, 2, 64, 65

Offset: 1

Antti Karttunen, Sep 17 2018

nonn

Restricted growth sequence transform of A319352.
For all i, j:
  a(i) = a(j) => A301855(i) = A301855(j).
  a(i) = a(j) => A304793(i) = A304793(j).

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

Cf. A056239, A319352.

up_to = 65537;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
A319352(n) = { my(m=1); fordiv(n, d, if(dA056239(d)))); (m); };
v319353 = rgs_transform(vector(up_to,n,A319352(n)));
A319353(n) = v319353[n];

A325800 Numbers whose sum of prime indices is equal to the number of distinct subset-sums of their prime indices.

Original entry on oeis.org

3, 10, 28, 66, 88, 156, 208, 306, 340, 408, 544, 570, 684, 760, 912, 966, 1216, 1242, 1288, 1380, 1656, 1840, 2208, 2436, 2610, 2900, 2944, 3132, 3248, 3480, 3906, 4092, 4176, 4340, 4640, 4650, 5022, 5208, 5456, 5568, 5580, 6200, 6696, 6944, 7326, 7424, 7440

Offset: 1

Gus Wiseman, May 23 2019

nonn

First differs from A325793 in lacking 70.
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, with sum A056239(n). A subset-sum of an integer partition is any sum of a submultiset of it.
The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k), so these are Heinz numbers of integer partitions whose sum is equal to their number of distinct subset-sums. The enumeration of these partitions by sum is given by A126796 interlaced with zeros.

			340 has prime indices {1,1,3,7} which sum to 12 and have 12 distinct subset-sums: {0, 1, 2, 3, 4, 5, 7, 8, 9, 10, 11, 12}, so 340 is in the sequence.
The sequence of terms together with their prime indices begins:
     3: {2}
    10: {1,3}
    28: {1,1,4}
    66: {1,2,5}
    88: {1,1,1,5}
   156: {1,1,2,6}
   208: {1,1,1,1,6}
   306: {1,2,2,7}
   340: {1,1,3,7}
   408: {1,1,1,2,7}
   544: {1,1,1,1,1,7}
   570: {1,2,3,8}
   684: {1,1,2,2,8}
   760: {1,1,1,3,8}
   912: {1,1,1,1,2,8}
   966: {1,2,4,9}
  1216: {1,1,1,1,1,1,8}
  1242: {1,2,2,2,9}
  1288: {1,1,1,4,9}
  1380: {1,1,2,3,9}

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

Positions of 1's in A325799.
Includes A239885 except for 1.
Cf. A002033, A056239, A108917, A112798, A126796, A299701, A299702.
Cf. A325694, A325780, A325781, A325792, A325793, A325801, A325802.

filter:= proc(n) local F,t,S,i,r;
  F:= map(t -> [numtheory:-pi(t[1]),t[2]], ifactors(n)[2]);
  S:= {0}:
  for t in F do
   S:= map(s -> seq(s + i*t[1],i=0..t[2]),S);
  od;
  nops(S) = add(t[1]*t[2],t=F)
end proc:
select(filter, [$1..10000]); # Robert Israel, Oct 30 2024

hwt[n_]:=Total[Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]*k]];
Select[Range[1000],hwt[#]==Length[Union[hwt/@Divisors[#]]]&]

A056239(a(n)) = A299701(a(n)) = A304793(a(n)) + 1.

A319352 a(n) = Product_{d|n, dA056239(d)), where A056239(d) gives the weight of the partition whose Heinz-number is d.

Original entry on oeis.org

1, 2, 2, 6, 2, 30, 2, 30, 10, 42, 2, 1050, 2, 66, 70, 210, 2, 2310, 2, 2310, 110, 78, 2, 80850, 14, 102, 110, 4290, 2, 210210, 2, 2310, 130, 114, 154, 1651650, 2, 138, 170, 210210, 2, 510510, 2, 6630, 10010, 174, 2, 11561550, 22, 7854, 190, 9690, 2, 510510, 182, 510510, 230, 186, 2, 2555102550, 2, 222, 20570, 30030, 238, 881790, 2

Offset: 1

Antti Karttunen, Sep 17 2018

nonn

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

Cf. A056239, A319353 (rgs-transform).

A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
A319352(n) = { my(m=1); fordiv(n, d, if(dA056239(d)))); (m); };

a(n) = Product_{d|n, dA056239(d)).
For all n >= 1:
A001221(a(n)) = A304793(n).
A001222(a(n)) = A032741(n).
1+A056169(a(n)) = A301855(n).

A304795 Number of positive special sums of the integer partition with Heinz number n.

Original entry on oeis.org

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

Offset: 1

Gus Wiseman, May 18 2018

nonn

A positive special sum of y is a number n > 0 such that exactly one submultiset of y sums to n. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

			The a(36) = 4 special sums are 1, 3, 5, 6, corresponding to the submultisets (1), (21), (221), (2211), with Heinz numbers 2, 6, 18, 36.

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

Cf. A000712, A056239, A108917, A122768, A276024, A284640, A296150, A299701, A299702, A301854, A301855, A301957, A304793, A304796.

primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
uqsubs[y_]:=Join@@Select[GatherBy[Union[Rest[Subsets[y]]],Total],Length[#]===1&];
Table[Length[uqsubs[primeMS[n]]],{n,100}]

up_to = 65537;
A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
v056239 = vector(up_to,n,A056239(n));
A304795(n) = { my(m=Map(),s,k=0,c); fordiv(n,d,if(!mapisdefined(m,s = v056239[d],&c), mapput(m,s,1), mapput(m,s,c+1))); sumdiv(n,d,(1==mapget(m,v056239[d])))-1; }; \\ Antti Karttunen, Jul 02 2018

More terms from Antti Karttunen, Jul 02 2018

cp's OEIS Frontend

A319695 Number of distinct values obtained when Euler phi (A000010) is applied to proper divisors of n.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

PARI

A325801 Number of divisors of n minus the number of distinct positive subset-sums of the prime indices of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A325986 Heinz numbers of complete strict integer partitions.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A326018 Heinz numbers of knapsack partitions such that no addition of one part up to the maximum is knapsack.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

A365919 Heinz numbers of integer partitions with the same number of distinct positive subset-sums as distinct non-subset-sums.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Formula

A316555 Number of distinct GCDs of nonempty submultisets of the integer partition with Heinz number n.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A319353 Filter sequence combining weights of partitions with Heinz-numbers equal to the proper divisors of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A325800 Numbers whose sum of prime indices is equal to the number of distinct subset-sums of their prime indices.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs