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.

Previous Showing 31-40 of 68 results. Next

A301855 Number of divisors d|n such that no other divisor of n has the same Heinz weight A056239(d).

Original entry on oeis.org

1, 2, 2, 3, 2, 4, 2, 4, 3, 4, 2, 4, 2, 4, 4, 5, 2, 6, 2, 6, 4, 4, 2, 4, 3, 4, 4, 6, 2, 6, 2, 6, 4, 4, 4, 5, 2, 4, 4, 6, 2, 8, 2, 6, 6, 4, 2, 4, 3, 6, 4, 6, 2, 8, 4, 8, 4, 4, 2, 4, 2, 4, 4, 7, 4, 8, 2, 6, 4, 6, 2, 4, 2, 4, 6, 6, 4, 8, 2, 6, 5, 4, 2, 6, 4, 4, 4, 8, 2, 6, 4, 6, 4, 4, 4, 4, 2, 6, 6, 9, 2, 8, 2, 8, 8
Offset: 1

Views

Author

Gus Wiseman, Mar 27 2018

Keywords

Examples

			The a(24) = 4 special divisors are 1, 2, 12, 24.
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    uqsubs[y_]:=Join@@Select[GatherBy[Union[Subsets[y]],Total],Length[#]===1&];
    Table[Length[uqsubs[primeMS[n]]],{n,100}]
  • PARI
    A056239(n) = { my(f); if(1==n, 0, f=factor(n); sum(i=1, #f~, f[i,2] * primepi(f[i,1]))); }
    A301855(n) = if(1==n,n,my(m=Map(),w,s); fordiv(n,d,w = A056239(d); if(!mapisdefined(m, w, &s), mapput(m,w,Set([d])), mapput(m,w,setunion(Set([d]),s)))); sumdiv(n,d,(1==length(mapget(m,A056239(d)))))); \\ Antti Karttunen, Jul 01 2018

Extensions

More terms from Antti Karttunen, Jul 01 2018

A301900 Heinz numbers of strict non-knapsack partitions. Squarefree numbers such that more than one divisor has the same Heinz weight A056239(d).

Original entry on oeis.org

30, 70, 154, 165, 210, 273, 286, 330, 390, 442, 462, 510, 546, 561, 570, 595, 646, 690, 714, 741, 770, 858, 870, 874, 910, 930, 1045, 1110, 1122, 1155, 1173, 1190, 1230, 1254, 1290, 1326, 1330, 1334, 1365, 1410, 1430, 1482, 1495, 1590, 1610, 1653, 1770
Offset: 1

Views

Author

Gus Wiseman, Mar 28 2018

Keywords

Comments

An integer partition is knapsack if every distinct submultiset has a different sum. The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).

Examples

			Sequence of strict non-knapsack partitions begins: (321), (431), (541), (532), (4321), (642), (651), (5321), (6321), (761), (5421), (7321), (6421), (752), (8321), (743), (871), (9321), (7421), (862), (5431), (6521).
		

Crossrefs

Programs

  • Mathematica
    wt[n_]:=If[n===1,0,Total[Cases[FactorInteger[n],{p_,k_}:>k*PrimePi[p]]]];
    Select[Range[1000],SquareFreeQ[#]&&!UnsameQ@@wt/@Divisors[#]&]

Formula

Complement of A005117 in A299702.

A301988 Nonprime Heinz numbers of integer partitions whose product is equal to their sum.

Original entry on oeis.org

9, 30, 84, 108, 200, 264, 624, 1120, 1440, 1632, 3648, 7040, 8832, 12544, 16128, 20736, 22272, 33280, 47616, 76800, 113664, 157696, 174080, 202752, 251904, 528384, 778240, 1155072, 1490944, 1916928, 2605056, 3440640, 3768320, 3964928, 4423680, 5799936
Offset: 1

Views

Author

Gus Wiseman, Mar 30 2018

Keywords

Comments

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

Examples

			Sequence of reversed integer partitions begins: (22), (123), (1124), (11222), (11133), (11125), (111126), (1111134), (11111223), (1111127), (11111128), (111111135), (111111129), (1111111144), (11111111224), (111111112222).
		

Crossrefs

Programs

  • Mathematica
    primeMS[n_]:=If[n===1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[200000],!PrimeQ[#]&&Total[primeMS[#]]===Times@@primeMS[#]&]

A305611 Number of distinct positive subset-sums of the multiset of prime factors of n.

Original entry on oeis.org

0, 1, 1, 2, 1, 3, 1, 3, 2, 3, 1, 5, 1, 3, 3, 4, 1, 5, 1, 5, 3, 3, 1, 7, 2, 3, 3, 5, 1, 6, 1, 5, 3, 3, 3, 8, 1, 3, 3, 7, 1, 7, 1, 5, 5, 3, 1, 9, 2, 5, 3, 5, 1, 7, 3, 7, 3, 3, 1, 9, 1, 3, 5, 6, 3, 7, 1, 5, 3, 6, 1, 10, 1, 3, 5, 5, 3, 7, 1, 9, 4, 3, 1, 10, 3, 3
Offset: 1

Views

Author

Gus Wiseman, Jun 06 2018

Keywords

Comments

An integer n is a positive subset-sum of a multiset y if there exists a nonempty submultiset of y with sum n.
One less than the number of distinct values obtained when A001414 is applied to all divisors of n. - Antti Karttunen, Jun 13 2018

Examples

			The a(12) = 5 positive subset-sums of {2, 2, 3} are 2, 3, 4, 5, and 7.
		

Crossrefs

Programs

  • Mathematica
    Table[Length[Union[Total/@Rest[Subsets[Join@@Cases[FactorInteger[n],{p_,k_}:>Table[p,{k}]]]]]],{n,100}]
  • PARI
    up_to = 65537;
    A001414(n) = ((n=factor(n))[, 1]~*n[, 2]); \\ From A001414.
    v001414 = vector(up_to,n,A001414(n));
    A305611(n) = { my(m=Map(),s,k=0); fordiv(n,d,if(!mapisdefined(m,s = v001414[d]), mapput(m,s,s); k++)); (k-1); }; \\ Antti Karttunen, Jun 13 2018
    
  • Python
    from sympy import factorint
    from sympy.utilities.iterables import multiset_combinations
    def A305611(n):
        fs = factorint(n)
        return len(set(sum(d) for i in range(1,sum(fs.values())+1) for d in multiset_combinations(fs,i))) # Chai Wah Wu, Aug 23 2021

A335549 Number of normal patterns matched by the multiset of prime indices of n in weakly increasing order.

Original entry on oeis.org

1, 2, 2, 3, 2, 3, 2, 4, 3, 3, 2, 5, 2, 3, 3, 5, 2, 5, 2, 5, 3, 3, 2, 7, 3, 3, 4, 5, 2, 4, 2, 6, 3, 3, 3, 7, 2, 3, 3, 7, 2, 4, 2, 5, 5, 3, 2, 9, 3, 5, 3, 5, 2, 7, 3, 7, 3, 3, 2, 7, 2, 3, 5, 7, 3, 4, 2, 5, 3, 4, 2, 10, 2, 3, 5, 5, 3, 4, 2, 9, 5, 3, 2, 7, 3, 3, 3
Offset: 1

Views

Author

Gus Wiseman, Jun 21 2020

Keywords

Comments

First differs from A181796 at a(90) = 8 A181796(90) = 7.
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 a (normal) pattern to be a finite sequence covering an initial interval of positive integers. Patterns are counted by A000670 and ranked by A333217. A sequence S is said to match a pattern P if there is a not necessarily contiguous subsequence of S whose parts have the same relative order as P. For example, (3,1,1,3) matches (1,1,2), (2,1,1), and (2,1,2), but avoids (1,2,1), (1,2,2), and (2,2,1).

Examples

			The Heinz number of (1,2,2,3) is 90 and it matches 8 patterns: (), (1), (11), (12), (112), (122), (123), (1223); so a(90) = 8.
		

Crossrefs

The version for standard compositions instead of prime indices is A335454.
Permutations of prime indices are counted by A008480.
Permutations are counted by A000142 and ranked by A333218.
Patterns are counted by A000670 and ranked by A333217.
Subset-sums are counted by A304792 and ranked by A299701.
Patterns matched by compositions of n are counted by A335456(n).
Minimal patterns avoided by a standard composition are counted by A335465.

Programs

  • Mathematica
    primeMS[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    mstype[q_]:=q/.Table[Union[q][[i]]->i,{i,Length[Union[q]]}];
    Table[Length[Union[mstype/@Subsets[primeMS[n]]]],{n,100}]

A366739 Number of distinct semi-sums of the multiset of prime indices of n. Number of distinct sums of prime indices of semiprime divisors of n (counted by A086971).

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Nov 04 2023

Keywords

Comments

First differs from A086971 at a(90) = 3, A086971(90) = 4.
A prime index of n is a number m such that prime(m) divides n. The multiset of prime indices of n is row n of A112798.
We define a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.

Examples

			The prime indices of 90 are {1,2,2,3}, with semi-sums
  3 = 1+2
  4 = 1+3 (or 2+2)
  5 = 2+3
so a(90) = 3.
Alternatively, the semiprime divisors of 90 are (6,9,10,15), with prime indices ({1,2},{2,2},{1,3},{2,3}) with sums (3,4,4,5) so a(90) = 3.
		

Crossrefs

The non-binary version is A299701.
Summing over partitions gives A366738, strict A366741.
For all sums of pairs of elements we have A367095.
Positions of first appearances are A367097.
A001222 counts prime factors (or prime indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, conjugate A065119.
A056239 adds up prime indices, row sums of A112798.
A299702 ranks knapsack partitions, counted by A108917.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Union[Total/@Subsets[prix[n],{2}]]],{n,100}]
  • PARI
    A366739(n) = #Set(apply(d->((f)->sum(i=1,#f~,f[i,2]*primepi(f[i,1])))(factor(d)), select(d->2==bigomega(d), divisors(n)))); \\ Antti Karttunen, Jan 20 2025

Formula

a(n) <= A086971(n). - Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A367095 Number of distinct sums of pairs (repeats allowed) of prime indices of n.

Original entry on oeis.org

0, 1, 1, 1, 1, 3, 1, 1, 1, 3, 1, 3, 1, 3, 3, 1, 1, 3, 1, 3, 3, 3, 1, 3, 1, 3, 1, 3, 1, 5, 1, 1, 3, 3, 3, 3, 1, 3, 3, 3, 1, 6, 1, 3, 3, 3, 1, 3, 1, 3, 3, 3, 1, 3, 3, 3, 3, 3, 1, 5, 1, 3, 3, 1, 3, 6, 1, 3, 3, 6, 1, 3, 1, 3, 3, 3, 3, 6, 1, 3, 1, 3, 1, 6, 3, 3, 3, 3, 1, 5, 3, 3, 3, 3, 3, 3, 1, 3, 3, 3, 1, 6, 1, 3, 5
Offset: 1

Views

Author

Gus Wiseman, Nov 06 2023

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.
Is the image missing only 2 and 4?

Examples

			The prime indices of 15 are {2,3}, with sums of pairs:
  2+2 = 4
  2+3 = 5
  3+3 = 6
so a(15) = 3.
The prime indices of 180 are {1,1,2,2,3}, with sums of pairs:
  1+1 = 2
  1+2 = 3
  1+3 = 4
  2+2 = 4
  2+3 = 5
  3+3 = 6
so a(180) = 5.
		

Crossrefs

Depends only on squarefree kernel A007947. (Even more exactly, on A322591 - Antti Karttunen, Jan 20 2025)
Positions of first appearances appear to be a subset of A325986.
For 2-element submultisets we have A366739, for all submultisets A299701.
A001222 counts prime factors (also indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, conjugate A065119.
A056239 adds up prime indices, row sums of A112798.
A304793 counts positive subset-sums of prime indices.
A367096 lists semiprime divisors, count A086971.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Length[Union[Total/@Tuples[prix[n],2]]],{n,100}]
  • PARI
    A367095(n) = if(1==n, 0, my(pis=apply(primepi,factor(n)[,1]), pairsums = vector(binomial(1+#pis,2)), k=0); for(i=1,#pis,for(j=i,#pis,k++; pairsums[k] = pis[i]+pis[j])); #Set(pairsums)); \\ Antti Karttunen, Jan 20 2025

Extensions

Data section extended to a(105) by Antti Karttunen, Jan 20 2025

A367093 Least positive integer with n more semiprime divisors than semi-sums of prime indices.

Original entry on oeis.org

1, 90, 630, 2310, 6930, 34650, 30030, 90090, 450450, 570570, 510510, 1531530, 7657650, 14804790, 11741730, 9699690, 29099070, 145495350
Offset: 0

Views

Author

Gus Wiseman, Nov 05 2023

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 a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.
Are all primorials after 210 included?

Examples

			The terms together with their prime indices begin:
       1: {}
      90: {1,2,2,3}
     630: {1,2,2,3,4}
    2310: {1,2,3,4,5}
    6930: {1,2,2,3,4,5}
   34650: {1,2,2,3,3,4,5}
   30030: {1,2,3,4,5,6}
   90090: {1,2,2,3,4,5,6}
  450450: {1,2,2,3,3,4,5,6}
  570570: {1,2,3,4,5,6,8}
  510510: {1,2,3,4,5,6,7}
		

Crossrefs

The first part (semiprime divisors) is A086971, firsts A220264.
The second part (semi-sums of prime indices) is A366739, firsts A367097.
All sums of pairs of prime indices are counted by A367095.
The non-binary version is A367105.
A001222 counts prime factors (or prime indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, conjugate A065119.
A056239 adds up prime indices, row sums of A112798.
A299701 counts subset-sums of prime indices, positive A304793.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264

Programs

  • Mathematica
    nn=10000;
    w=Table[Length[Union[Subsets[prix[n],{2}]]]-Length[Union[Total/@Subsets[prix[n],{2}]]],{n,nn}];
    spnm[y_]:=Max@@NestWhile[Most,y,Union[#]!=Range[0,Max@@#]&];
    Table[Position[w,k][[1,1]],{k,0,spnm[w]}]
  • Python
    from itertools import count
    from sympy import factorint, primepi
    from sympy.utilities.iterables import multiset_combinations
    def A367093(n):
        for k in count(1):
            c, a = 0, set()
            for s in (sum(p) for p in multiset_combinations({primepi(i):j for i,j in factorint(k).items()},2)):
                if s not in a:
                    a.add(s)
                else:
                    c += 1
                if c > n:
                    break
            if c == n:
                return k # Chai Wah Wu, Nov 13 2023

Formula

a(n) is the least positive integer such that A086971(a(n)) - A366739(a(n)) = n.

Extensions

a(12)-a(16) from Chai Wah Wu, Nov 13 2023
a(17) from Chai Wah Wu, Nov 18 2023

A367404 Triangle read by rows where T(n,k) is the number of integer partitions of n with a semi-sum k.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 3, 2, 2, 2, 5, 3, 4, 2, 3, 7, 5, 6, 4, 3, 3, 11, 7, 9, 6, 6, 3, 4, 15, 11, 13, 10, 9, 6, 4, 4, 22, 15, 20, 13, 15, 9, 8, 4, 5, 30, 22, 27, 21, 21, 15, 12, 8, 5, 5, 42, 30, 39, 28, 30, 21, 20, 12, 10, 5, 6, 56, 42, 53, 41, 42, 33, 28, 20, 15, 10, 6, 6
Offset: 2

Views

Author

Gus Wiseman, Nov 17 2023

Keywords

Comments

We define a semi-sum of a multiset to be any sum of a 2-element submultiset. This is different from sums of pairs of elements. For example, 2 is the sum of a pair of elements of {1}, but there are no semi-sums.

Examples

			The partition y = (3,2,1,1) has semi-sum 3 = 2+1, but no semi-sum 6, so y is counted under T(7,3) but not under T(7,6).
Triangle begins:
   1
   1   1
   2   1   2
   3   2   2   2
   5   3   4   2   3
   7   5   6   4   3   3
  11   7   9   6   6   3   4
  15  11  13  10   9   6   4   4
  22  15  20  13  15   9   8   4   5
  30  22  27  21  21  15  12   8   5   5
  42  30  39  28  30  21  20  12  10   5   6
  56  42  53  41  42  33  28  20  15  10   6   6
  77  56  73  55  60  42  44  28  25  15  12   6   7
Row n = 7 counts the following partitions:
  (511)      (421)     (331)    (421)   (511)  (61)
  (4111)     (3211)    (322)    (4111)  (421)  (52)
  (3211)     (2221)    (3211)   (322)   (331)  (43)
  (31111)    (22111)   (31111)  (3211)
  (22111)    (211111)  (2221)
  (211111)             (22111)
  (1111111)
		

Crossrefs

Column k = 0 is A000041.
Column n = k is A004526.
The complement for all submultisets is A046663, strict A365663.
For subsets instead of partitions we have A365541, non-binary A365381.
The non-binary version is A365543, strict A365661.
Row sums are A366738.
The strict case is A367405.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n], MemberQ[Total/@Subsets[#, {2}],k]&]], {n,2,10}, {k,2,n}]

A365832 Triangle read by rows where T(n,k) is the number of strict integer partitions of n with k distinct sums of nonempty subsets.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Sep 28 2023

Keywords

Examples

			The partition (7,6,1) has sums 1, 6, 7, 8, 13, 14, so is counted under T(14,6).
Triangle begins:
  1
  0  1
  0  1  0
  0  1  0  1
  0  1  0  1  0
  0  1  0  2  0  0
  0  1  0  2  0  0  1
  0  1  0  3  0  0  0  1
  0  1  0  3  0  0  1  1  0
  0  1  0  4  0  0  0  3  0  0
  0  1  0  4  0  0  2  2  0  0  1
  0  1  0  5  0  0  0  5  0  0  0  1
  0  1  0  5  0  0  2  5  0  0  0  0  2
  0  1  0  6  0  0  0  8  0  0  0  1  0  2
  0  1  0  6  0  0  3  7  0  0  0  0  3  1  1
  0  1  0  7  0  0  0 12  0  0  0  1  0  4  0  2
  0  1  0  7  0  0  3 11  0  0  0  1  3  2  2  1  1
  0  1  0  8  0  0  0 16  0  0  0  1  0  7  0  3  0  2
  0  1  0  8  0  0  4 15  0  0  0  1  3  3  6  2  0  0  3
  0  1  0  9  0  0  0 21  0  0  0  2  0  9  0  7  0  1  0  4
  0  1  0  9  0  0  4 20  0  0  1  0  4  8  5  5  0  0  2  0  5
Row n = 14 counts the following partitions (A..E = 10..14):
  (E)  .  (D1)  .  .  (761)  (B21)  .  .  .  .  (6521)  (8321)  (7421)
          (C2)        (752)  (A31)              (6431)
          (B3)        (743)  (941)              (5432)
          (A4)               (932)
          (95)               (851)
          (86)               (842)
                             (653)
		

Crossrefs

Row sums are A000009.
Rightmost column n = k is A188431, non-strict A126796.
The one-based weighted row sums are A284640.
The corresponding rank statistic is A299701.
The non-strict version is A365658.
Central column n = 2k in the non-strict case is A365660.
Reverse-weighted row-sums are A365922, non-strict A276024.
A000041 counts integer partitions.
A000124 counts distinct sums of subsets of {1..n}.
A365543 counts partitions with a submultiset summing to k, strict A365661.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[n],UnsameQ@@#&&Length[Union[Total/@Rest[Subsets[#]]]]==k&]],{n,0,15},{k,0,n}]
Previous Showing 31-40 of 68 results. Next