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 11-20 of 41 results. Next

A367223 Number of subsets of {1..n} whose cardinality cannot be written as a nonnegative linear combination of the elements.

Original entry on oeis.org

0, 0, 1, 2, 4, 8, 15, 27, 49, 90, 165, 301, 548, 998, 1819, 3316, 6040, 10986, 19959, 36253, 65904, 119986, 218796, 399461, 729752, 1333162, 2434411, 4441954, 8097478, 14746715, 26830230, 48773790, 88605927, 160900978, 292140427, 530487359, 963610200, 1751171679, 3183997509
Offset: 0

Views

Author

Gus Wiseman, Nov 14 2023

Keywords

Examples

			3 cannot be written as a nonnegative linear combination of 2, 4, and 5, so {2,4,5} is counted under a(6).
The a(2) = 1 through a(6) = 15 subsets:
  {2}  {2}  {2}    {2}      {2}
       {3}  {3}    {3}      {3}
            {4}    {4}      {4}
            {3,4}  {5}      {5}
                   {3,4}    {6}
                   {3,5}    {3,4}
                   {4,5}    {3,5}
                   {2,4,5}  {3,6}
                            {4,5}
                            {4,6}
                            {5,6}
                            {2,4,5}
                            {2,4,6}
                            {2,5,6}
                            {4,5,6}
		

Crossrefs

The following sequences count and rank integer partitions and finite sets according to whether their length is a subset-sum or linear combination of the parts. The current sequence is starred.
sum-full sum-free comb-full comb-free
-------------------------------------------
A007865/A085489/A151897 count certain types of sum-free subsets.
A088809/A093971/A364534 count certain types of sum-full subsets.
A124506 appears to count combination-free subsets, differences of A326083.
A365046 counts combination-full subsets, differences of A364914.
Triangles:
A116861 counts positive linear combinations of strict partitions of k.
A364916 counts linear combinations of strict partitions of k.
A366320 counts subsets without a subset summing to k, with A365381.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y}, {i,0,Floor[n/k]}]}, Select[Tuples[s], Total[Times@@@#]==n&]];
    Table[Length[Select[Subsets[Range[n]], combs[Length[#],Union[#]]=={}&]], {n,0,10}]
  • Python
    from itertools import combinations
    from sympy.utilities.iterables import partitions
    def A367223(n):
        c, mlist = 0, []
        for m in range(1,n+1):
            t = set()
            for p in partitions(m):
                t.add(tuple(sorted(p.keys())))
            mlist.append([set(d) for d in t])
        for k in range(1,n+1):
            for w in combinations(range(1,n+1),k):
                ws = set(w)
                for s in mlist[k-1]:
                    if s <= ws:
                        break
                else:
                    c += 1
        return c # Chai Wah Wu, Nov 16 2023

Formula

a(n) = 2^n - A367222(n).

Extensions

a(14)-a(33) from Chai Wah Wu, Nov 15 2023
a(34)-a(38) from Max Alekseyev, Feb 25 2025

A364915 Number of integer partitions of n such that no distinct part can be written as a nonnegative linear combination of other distinct parts.

Original entry on oeis.org

1, 1, 2, 2, 3, 3, 4, 5, 6, 7, 8, 12, 10, 16, 16, 19, 21, 29, 25, 37, 35, 44, 46, 60, 55, 75, 71, 90, 90, 114, 110, 140, 138, 167, 163, 217, 201, 248, 241, 298, 303, 359, 355, 425, 422, 520, 496, 594, 603, 715, 706, 834, 826, 968, 972, 1153, 1147, 1334, 1315, 1530
Offset: 0

Views

Author

Gus Wiseman, Aug 22 2023

Keywords

Examples

			The a(1) = 1 through a(10) = 8 partitions (A=10):
  1  2   3    4     5      6       7        8         9          A
     11  111  22    32     33      43       44        54         55
              1111  11111  222     52       53        72         64
                           111111  322      332       333        73
                                   1111111  2222      522        433
                                            11111111  3222       3322
                                                      111111111  22222
                                                                 1111111111
The partition (5,4,3) has no part that can be written as a nonnegative linear combination of the others, so is counted under a(12).
The partition (6,4,3,2) has 6=4+2, or 6=3+3, or 6=2+2+2, or 4=2+2, so is not counted under a(15).
		

Crossrefs

For sums instead of combinations we have A237667, binary A236912.
For subsets instead of partitions we have A326083, complement A364914.
The strict case is A364350.
The complement is A365068, strict A364839.
The positive case is A365072, strict A365006.
A000041 counts integer partitions, strict A000009.
A007865 counts binary sum-free sets w/ re-usable parts, complement A093971.
A008284 counts partitions by length, strict A008289.
A116861 and A364916 count linear combinations of strict partitions.
A364912 counts linear combinations of partitions of k.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y}, {i,0,Floor[n/k]}]}, Select[Tuples[s], Total[Times@@@#]==n&]];
    Table[Length[Select[IntegerPartitions[n], Function[ptn,!Or@@Table[combs[ptn[[k]],Delete[ptn,k]]!={}, {k,Length[ptn]}]]@*Union]], {n,0,15}]
  • Python
    from sympy.utilities.iterables import partitions
    def A364915(n):
        if n <= 1: return 1
        alist, c = [set(tuple(sorted(set(p))) for p in partitions(i)) for i in range(n)], 1
        for p in partitions(n,k=n-1):
            s = set(p)
            if not any(set(t).issubset(s-{q}) for q in s for t in alist[q]):
                c += 1
        return c # Chai Wah Wu, Sep 23 2023

Formula

a(n) = A000041(n) - A365068(n).

Extensions

a(37)-a(59) from Chai Wah Wu, Sep 25 2023

A280292 a(n) = sopfr(n) - sopf(n).

Original entry on oeis.org

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

Views

Author

Michel Marcus, Dec 31 2016

Keywords

Comments

Alladi and Erdős (1977) proved that for all numbers m>=0, m!=1, the sequence of numbers k such that a(k) = m has a positive asymptotic density which is equal to a rational multiple of 1/zeta(2) = 6/Pi^2 (A059956). For example, when m=0, the sequence is the squarefree numbers (A005117), whose density is 6/Pi^2, and when m=2 the sequence is A081770, whose density is 1/Pi^2. - Amiram Eldar, Nov 02 2020
Sum of prime factors minus sum of distinct prime factors. Counting partitions by this statistic (sum minus sum of distinct parts) gives A364916. - Gus Wiseman, Feb 21 2025

References

  • Jean-Marie De Koninck and Aleksandar Ivić, Topics in Arithmetical Functions: Asymptotic Formulae for Sums of Reciprocals of Arithmetical Functions and Related Fields, Amsterdam, Netherlands: North-Holland, 1980. See pp. 164-166.
  • Steven R. Finch, Mathematical Constants II, Encyclopedia of Mathematics and Its Applications, Cambridge University Press, Cambridge, 2018, p. 165.

Crossrefs

A multiplicative version is A003557, firsts A064549 (sorted A001694).
For length instead of sum we have A046660.
For product instead of sum we have A066503, firsts A381076.
Positions of first appearances are A280286 (sorted A381075).
For indices instead of factors we have A380955, firsts A380956 (sorted A380957).
For exponents instead of factors we have A380958, firsts A380989.
A000040 lists the primes, differences A001223.
A001222 counts prime factors (distinct A001221).
A003963 gives product of prime indices, distinct A156061, excess A380986.
A005117 lists squarefree numbers, complement A013929.
A007947 gives squarefree kernel.
A020639 gives least prime factor (index A055396), greatest A061395 (index A006530).
A027746 lists prime factors, distinct A027748.
A112798 lists prime indices (sum A056239), distinct A304038 (sum A066328).

Programs

  • Mathematica
    Array[Total@ # - Total@ Union@ # &@ Flatten[ConstantArray[#1, #2] & @@@ FactorInteger@ #] &, 105] (* Michael De Vlieger, Feb 25 2019 *)
  • PARI
    sopfr(n) = my(f=factor(n)); sum(j=1, #f~, f[j, 1]*f[j, 2]);
    sopf(n) = my(f=factor(n)); sum(j=1, #f~, f[j, 1]);
    a(n) = sopfr(n) - sopf(n);

Formula

a(n) = A001414(n) - A008472(n).
a(A005117(n)) = 0.
a(n) = A001414(A003557(n)). - Antti Karttunen, Oct 07 2017
Additive with a(1) = 0 and a(p^e) = p*(e-1) for prime p and e > 0. - Werner Schulte, Feb 24 2019
From Amiram Eldar, Nov 02 2020: (Start)
a(n) = a(A057521(n)).
Sum_{n<=x} a(n) ~ x*log(log(x)) + O(x) (Alladi and Erdős, 1977).
Sum_{n<=x, n nonsquarefree} 1/a(n) ~ c*x + O(sqrt(x)*log(x)), where c = Integral_{t=0..1} (F(t)-6/Pi^2)/t dt, and F(t) = Product_{p prime} (1-1/p)*(1-1/(t^p - p)) (De Koninck et al., 1981; Finch, 2018), or, equivalently c = Sum_{k>=2} d(k)/k = 0.1039..., where d(k) = (6/Pi^2)*A338559(k)/A338560(k) is the asymptotic density of the numbers m with a(m) = k (Alladi and Erdős, 1977; Ivić, 2003). (End)

Extensions

More terms from Antti Karttunen, Oct 07 2017

A364911 Triangle read by rows where T(n,k) is the number of integer partitions with sum <= n and with distinct parts summing to k.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Aug 27 2023

Keywords

Comments

Also the number of ways to write any number up to n as a positive linear combination of a strict integer partition of k.

Examples

			Triangle begins:
  1
  1  1
  1  2  1
  1  3  1  2
  1  4  2  3  2
  1  5  2  5  3  3
  1  6  3  8  4  4  4
  1  7  3 11  6  6  6  5
  1  8  4 14  9  8 10  7  6
  1  9  4 19 11 11 14 11  9  8
  1 10  5 23 14 15 21 15 14 11 10
  1 11  5 28 17 19 28 22 20 17 15 12
  1 12  6 34 21 22 40 28 28 24 24 17 15
  1 13  6 40 25 27 50 38 37 34 35 27 22 18
  1 14  7 46 29 32 65 49 50 43 51 38 35 26 22
  1 15  7 54 33 38 79 62 63 59 68 55 50 41 32 27
Row n = 5 counts the following partitions:
    .    1           2     3         4       5
         1+1         2+2   1+2       1+3     1+4
         1+1+1             1+1+2     1+1+3   2+3
         1+1+1+1           1+1+1+2
         1+1+1+1+1         1+2+2
Row n = 5 counts the following positive linear combinations:
  .  1*1  1*2  1*3      1*4      1*5
     2*1  2*2  1*2+1*1  1*3+1*1  1*3+1*2
     3*1       1*2+2*1  1*3+2*1  1*4+1*1
     4*1       1*2+3*1
     5*1       2*2+1*1
		

Crossrefs

Column n = k is A000009.
Column k = 0 is A000012.
Column k = 1 is A000027.
Row sums are A000070.
Column k = 2 is A008619.
Columns are partial sums of columns of A116861.
Column k = 3 appears to be the partial sums of A137719.
Diagonal n = 2k is A364910.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A114638 counts partitions where (length) = (sum of distinct parts).
A116608 counts partitions by number of distinct parts.
A364350 counts combination-free strict partitions, complement A364839.

Programs

  • Mathematica
    Table[Length[Select[Array[IntegerPartitions,n+1,0,Join],Total[Union[#]]==k&]],{n,0,9},{k,0,n}]
  • PARI
    T(n)={[Vecrev(p) | p<-Vec(prod(k=1, n, 1 - y^k + y^k/(1 - x^k), 1/(1 - x) + O(x*x^n)))]}
    { my(A=T(10)); for(n=1, #A, print(A[n])) } \\ Andrew Howroyd, Jan 11 2024

Formula

G.f.: A(x,y) = (1/(1 - x)) * Product_{k>=1} (1 - y^k + y^k/(1 - x^k)). - Andrew Howroyd, Jan 11 2024

A380955 Sum of prime indices of n (with multiplicity) minus sum of distinct prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Feb 11 2025

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 96 are {1,1,1,1,1,2}, with sum 7, and with distinct prime indices {1,2}, with sum 3, so a(96) = 7 - 3 = 4.
		

Crossrefs

Positions of 0's are A005117, complement A013929.
For length instead of sum we have A046660.
Positions of 1's are A081770.
For factors instead of indices we have A280292, firsts A280286 (sorted A381075).
A multiplicative version is A290106.
Counting partitions by this statistic gives A364916.
Dominates A374248.
Positions of first appearances are A380956, sorted A380957.
For prime multiplicities instead of prime indices we have A380958.
For product instead of sum we have A380986.
A000040 lists the primes, differences A001223.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, length A001222.
A304038 lists distinct prime indices, sum A066328, length A001221.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Table[Total[prix[n]]-Total[Union[prix[n]]],{n,100}]

Formula

a(n) = A056239(n) - A066328(n).
Additive: a(m*n) = a(m) + a(n) if gcd(m,n) = 1.

A364910 Number of integer partitions of 2n whose distinct parts sum to n.

Original entry on oeis.org

1, 1, 1, 3, 3, 4, 12, 11, 19, 23, 54, 55, 103, 115, 178, 289, 389, 507, 757, 970, 1343, 2033, 2579, 3481, 4840, 6312, 8317, 10998, 15459, 19334, 26368, 33480, 44709, 56838, 74878, 93369, 128109, 157024, 206471, 258357, 338085, 417530, 544263, 669388, 859570, 1082758, 1367068
Offset: 0

Views

Author

Gus Wiseman, Aug 16 2023

Keywords

Comments

Also the number of ways to write n as a nonnegative linear combination of the parts of a strict integer partition of n.

Examples

			The a(0) = 1 through a(7) = 11 partitions:
  ()  (11)  (22)  (33)     (44)      (55)       (66)         (77)
                  (2211)   (3311)    (3322)     (4422)       (4433)
                  (21111)  (311111)  (4411)     (5511)       (5522)
                                     (4111111)  (33321)      (6611)
                                                (42222)      (442211)
                                                (322221)     (4222211)
                                                (332211)     (4421111)
                                                (3222111)    (42221111)
                                                (3321111)    (422111111)
                                                (32211111)   (611111111)
                                                (51111111)   (4211111111)
                                                (321111111)
The a(0) = 1 through a(7) = 11 linear combinations:
  0  1*1  1*2  1*3      1*4      1*5      1*6          1*7
               0*2+3*1  0*3+4*1  0*4+5*1  0*4+3*2      0*6+7*1
               1*2+1*1  1*3+1*1  1*3+1*2  0*5+6*1      1*4+1*3
                                 1*4+1*1  1*4+1*2      1*5+1*2
                                          1*5+1*1      1*6+1*1
                                          0*3+0*2+6*1  0*4+0*2+7*1
                                          0*3+1*2+4*1  0*4+1*2+5*1
                                          0*3+2*2+2*1  0*4+2*2+3*1
                                          0*3+3*2+0*1  0*4+3*2+1*1
                                          1*3+0*2+3*1  1*4+0*2+3*1
                                          1*3+1*2+1*1  1*4+1*2+1*1
                                          2*3+0*2+0*1
		

Crossrefs

The case with no zero coefficients is A000009.
Central diagonal of A116861.
A version based on Heinz numbers is A364906.
Using all partitions (not just strict) we get A364907.
The version for compositions is A364908, strict A364909.
Main diagonal of A364916.
Using strict partitions of any number from 1 to n gives A365002.
These partitions have ranks A365003.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A323092 counts double-free partitions, ranks A320340.

Programs

  • Mathematica
    Table[Length[Select[IntegerPartitions[2n],Total[Union[#]]==n&]],{n,0,15}]
  • PARI
    a(n) = {my(res = 0); forpart(p = 2*n,s = Set(p); if(vecsum(s) == n, res++)); res} \\ David A. Corneth, Aug 20 2023
    
  • Python
    from sympy.utilities.iterables import partitions
    def A364910(n): return sum(1 for d in partitions(n<<1,k=n) if sum(set(d))==n) # Chai Wah Wu, Sep 13 2023

Formula

a(n) = A116861(2n,n).
a(n) = A364916(n,n).

Extensions

More terms from David A. Corneth, Aug 20 2023

A365312 Number of strict integer partitions with sum <= n that cannot be linearly combined using nonnegative coefficients to obtain n.

Original entry on oeis.org

0, 0, 0, 1, 1, 3, 2, 6, 4, 8, 7, 16, 6, 24, 17, 24, 20, 46, 22, 62, 31, 63, 57, 106, 35, 122, 90, 137, 88, 212, 74, 262, 134, 267, 206, 345, 121, 476, 294, 484, 232, 698, 242, 837, 389, 763, 571, 1185, 318, 1327, 634, 1392, 727, 1927, 640, 2056, 827, 2233, 1328
Offset: 0

Views

Author

Gus Wiseman, Sep 05 2023

Keywords

Examples

			The strict partition (7,3,2) has 19 = 1*7 + 2*3 + 3*2 so is not counted under a(19).
The strict partition (9,6,3) cannot be linearly combined to obtain 19, so is counted under a(19).
The a(0) = 0 through a(11) = 16 strict partitions:
  .  .  .  (2)  (3)  (2)  (4)  (2)    (3)  (2)    (3)    (2)
                     (3)  (5)  (3)    (5)  (4)    (4)    (3)
                     (4)       (4)    (6)  (5)    (6)    (4)
                               (5)    (7)  (6)    (7)    (5)
                               (6)         (7)    (8)    (6)
                               (4,2)       (8)    (9)    (7)
                                           (4,2)  (6,3)  (8)
                                           (6,2)         (9)
                                                         (10)
                                                         (4,2)
                                                         (5,4)
                                                         (6,2)
                                                         (6,3)
                                                         (6,4)
                                                         (7,3)
                                                         (8,2)
		

Crossrefs

The complement for positive coefficients is counted by A088314.
For positive coefficients we have A088528.
The complement is counted by A365311.
For non-strict partitions we have A365378, complement A365379.
The version for subsets is A365380, complement A365073.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A116861 and A364916 count linear combinations of strict partitions.
A364350 counts combination-free strict partitions, non-strict A364915.
A364839 counts combination-full strict partitions, non-strict A364913.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Length[Select[Select[Join@@Array[IntegerPartitions,n], UnsameQ@@#&],combs[n,#]=={}&]],{n,0,10}]
  • Python
    from math import isqrt
    from sympy.utilities.iterables import partitions
    def A365312(n):
        a = {tuple(sorted(set(p))) for p in partitions(n)}
        return sum(1 for m in range(1,n+1) for b in partitions(m,m=isqrt(1+(n<<3))>>1) if max(b.values()) == 1 and not any(set(d).issubset(set(b)) for d in a)) # Chai Wah Wu, Sep 13 2023

Extensions

a(26)-a(58) from Chai Wah Wu, Sep 13 2023

A365311 Number of strict integer partitions with sum <= n that can be linearly combined using nonnegative coefficients to obtain n.

Original entry on oeis.org

0, 1, 2, 3, 5, 6, 11, 12, 20, 24, 35, 38, 63, 63, 92, 112, 148, 160, 230, 244, 339, 383, 478, 533, 726, 781, 978, 1123, 1394, 1526, 1960, 2112, 2630, 2945, 3518, 3964, 4856, 5261, 6307, 7099, 8464, 9258, 11140, 12155, 14419, 16093, 18589, 20565, 24342, 26597, 30948
Offset: 0

Views

Author

Gus Wiseman, Sep 04 2023

Keywords

Examples

			The strict partition (6,3) cannot be linearly combined to obtain 10, so is not counted under a(10).
The strict partition (4,2) has 6 = 1*4 + 1*2 so is counted under a(6), but (4,2) cannot be linearly combined to obtain 7 so is not counted under a(7).
The a(1) = 1 through a(7) = 12 strict partitions:
  (1)  (1)  (1)    (1)    (1)    (1)      (1)
       (2)  (3)    (2)    (5)    (2)      (7)
            (2,1)  (4)    (2,1)  (3)      (2,1)
                   (2,1)  (3,1)  (6)      (3,1)
                   (3,1)  (3,2)  (2,1)    (3,2)
                          (4,1)  (3,1)    (4,1)
                                 (3,2)    (4,3)
                                 (4,1)    (5,1)
                                 (4,2)    (5,2)
                                 (5,1)    (6,1)
                                 (3,2,1)  (3,2,1)
                                          (4,2,1)
		

Crossrefs

For positive coefficients we have A088314.
The positive complement is counted by A088528.
The version for subsets is A365073.
The complement is counted by A365312.
For non-strict partitions we have A365379.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A116861 and A364916 count linear combinations of strict partitions.
A364350 counts combination-free strict partitions, non-strict A364915.
A364839 counts combination-full strict partitions, non-strict A364913.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]},Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Length[Select[Select[Join@@Array[IntegerPartitions,n],UnsameQ@@#&],combs[n,#]!={}&]],{n,10}]
  • Python
    from math import isqrt
    from sympy.utilities.iterables import partitions
    def A365311(n):
        a = {tuple(sorted(set(p))) for p in partitions(n)}
        return sum(1 for m in range(1,n+1) for b in partitions(m,m=isqrt(1+(n<<3))>>1) if max(b.values()) == 1 and any(set(d).issubset(set(b)) for d in a)) # Chai Wah Wu, Sep 13 2023

Extensions

a(26)-a(50) from Chai Wah Wu, Sep 13 2023

A280286 a(n) is the least k such that sopfr(k) - sopf(k) = n.

Original entry on oeis.org

4, 9, 8, 25, 16, 49, 32, 81, 64, 121, 128, 169, 256, 625, 512, 289, 1024, 361, 2048, 1444, 1331, 529, 5324, 2116, 2197, 4232, 8788, 841, 17576, 961, 7569, 3844, 4913, 7688, 19652, 1369, 6859, 5476, 12321, 1681, 34225, 1849, 15129, 7396, 12167, 2209, 46225, 8836, 19881
Offset: 2

Views

Author

Michel Marcus, Dec 31 2016

Keywords

Crossrefs

Cf. A001414 (sopfr), A008472 (sopf), A001248, A280163.
A multiplicative version is A064549 (sorted A001694), firsts of A003557.
For length instead of sum we have A151821.
These are the positions of first appearances in A280292 = A001414 - A008472.
For indices instead of factors we have A380956 (sorted A380957), firsts of A380955.
A multiplicative version for indices is A380987 (sorted A380988), firsts of A290106.
For prime exponents instead of factors we have A380989, firsts of A380958.
The sorted version is A381075.
For product instead of sum see A381076, sorted firsts of A066503.
A000040 lists the primes, differences A001223.
A005117 lists squarefree numbers, complement A013929.
A020639 gives least prime factor (index A055396), greatest A061395 (index A006530).
A056239 adds up prime indices, row sums of A112798.

Programs

  • Mathematica
    prifacs[n_]:=If[n==1,{},Flatten[Apply[ConstantArray,FactorInteger[n],{1}]]];
    q=Table[Total[prifacs[n]]-Total[Union[prifacs[n]]],{n,1000}];
    mnrm[s_]:=If[Min@@s==1,mnrm[DeleteCases[s-1,0]]+1,0];
    Table[Position[q,k][[1,1]],{k,2,mnrm[q/.(0->1)]}] (* Gus Wiseman, Feb 20 2025 *)
  • PARI
    sopfr(n) = my(f=factor(n)); sum(j=1, #f~, f[j,1]*f[j,2]);
    sopf(n) = my(f=factor(n)); sum(j=1, #f~, f[j,1]);
    a(n) = {my(k = 2); while (sopfr(k) - sopf(k) != n, k++); k;}

Formula

For p prime, a(p) = p^2 (see A001248).

A365002 Number of ways to write n as a nonnegative linear combination of a strict integer partition.

Original entry on oeis.org

1, 1, 2, 4, 8, 10, 26, 32, 63, 84, 157, 207, 383, 477, 768, 1108, 1710, 2261, 3536, 4605, 6869, 9339, 13343, 17653, 25785, 33463, 46752, 61549, 85614, 110861, 153719, 197345, 268623, 346845, 463513, 593363, 797082, 1011403, 1335625, 1703143, 2232161, 2820539
Offset: 0

Views

Author

Gus Wiseman, Aug 22 2023

Keywords

Comments

A way of writing n as a (nonnegative) linear combination of a finite sequence y is any sequence of pairs (k_i,y_i) such that k_i >= 0 and Sum k_i*y_i = n. For example, the pairs ((3,1),(1,1),(1,1),(0,2)) are a way of writing 5 as a linear combination of (1,1,1,2), namely 5 = 3*1 + 1*1 + 1*1 + 0*2. Of course, there are A000041(n) ways to write n as a linear combination of (1..n).

Examples

			The a(1) = 1 through a(5) = 10 ways:
  1*1  1*2  1*3      1*4      1*5
       2*1  3*1      2*2      5*1
            0*2+3*1  4*1      0*2+5*1
            1*2+1*1  0*2+4*1  0*3+5*1
                     0*3+4*1  0*4+5*1
                     1*2+2*1  1*2+3*1
                     1*3+1*1  1*3+1*2
                     2*2+0*1  1*3+2*1
                              1*4+1*1
                              2*2+1*1
		

Crossrefs

Row sums of lower-left half of A364916 as an array.
Column sums of right half of A364916 as a triangle.
For all positive coefficients we have A000041, non-strict A006951.
A000041 counts integer partitions, strict A000009.
A008284 counts partitions by length, strict A008289.
A364350 counts combination-free strict partitions, complement A364839.
A364913 counts combination-full partitions.

Programs

  • Mathematica
    combs[n_,y_]:=With[{s=Table[{k,i},{k,y},{i,0,Floor[n/k]}]}, Select[Tuples[s],Total[Times@@@#]==n&]];
    Table[Sum[Length[combs[n,y]], {y,Select[Join@@IntegerPartitions/@Range[n], UnsameQ@@#&]}],{n,0,15}]
  • Python
    from itertools import combinations
    from collections import Counter
    from sympy.utilities.iterables import partitions
    def A365002(n):
        aset = Counter(tuple(sorted(set(p))) for p in partitions(n))
        return sum(sum(aset[t] for t in aset if set(t).issubset(set(q))) for l in range(1,n+1) for q in combinations(range(1,n+1),l) if sum(q)<=n) # Chai Wah Wu, Sep 20 2023

Extensions

a(16)-a(34) from Chai Wah Wu, Sep 20 2023
a(35)-a(38) from Chai Wah Wu, Sep 21 2023
a(0)=1 and a(39)-a(41) from Alois P. Heinz, Jan 11 2024
Previous Showing 11-20 of 41 results. Next