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

A366740 Positive integers whose semiprime divisors do not all have different Heinz weights (sum of prime indices, A056239).

Original entry on oeis.org

90, 180, 210, 270, 360, 420, 450, 462, 525, 540, 550, 630, 720, 810, 840, 858, 900, 910, 924, 990, 1050, 1080, 1100, 1155, 1170, 1260, 1326, 1350, 1386, 1440, 1470, 1530, 1575, 1620, 1650, 1666, 1680, 1710, 1716, 1800, 1820, 1848, 1870, 1890, 1911, 1938, 1980
Offset: 1

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.
From Robert Israel, Nov 06 2023: (Start)
Positive integers divisible by the product of four primes, prime(i)*prime(j)*prime(k)*prime(l), i < j <= k < l, with i + l = j + k.
All positive multiples of terms are terms. (End)

Examples

			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), which are not all different, so 90 is in the sequence.
The terms together with their prime indices begin:
    90: {1,2,2,3}
   180: {1,1,2,2,3}
   210: {1,2,3,4}
   270: {1,2,2,2,3}
   360: {1,1,1,2,2,3}
   420: {1,1,2,3,4}
   450: {1,2,2,3,3}
   462: {1,2,4,5}
   525: {2,3,3,4}
   540: {1,1,2,2,2,3}
   550: {1,3,3,5}
   630: {1,2,2,3,4}
   720: {1,1,1,1,2,2,3}
		

Crossrefs

The complement is too dense.
For all divisors instead of just semiprimes we have A299729, strict A316402.
Distinct semi-sums of prime indices are counted by A366739.
Partitions of this type are counted by A366753, non-binary A366754.
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 distinct subset-sums of prime indices, positive A304793.
A299702 ranks knapsack partitions, counted by A108917, strict A275972.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264

Programs

  • Maple
    N:= 10^4: # for terms <= N
    P:= select(isprime, [$1..N]): nP:= nops(P):
    R:= {}:
    for i from 1 while P[i]*P[i+1]^2*P[i+2] < N do
      for j from i+1 while P[i]*P[j]^2 * P[j+1] < N do
        for k from j do
          l:= j+k-i;
          if l <= k or l > nP then break fi;
          v:= P[i]*P[j]*P[k]*P[l];
          if v <= N then
            R:= R union {seq(t,t=v..N,v)};
          fi
    od od od:
    sort(convert(R,list)); # Robert Israel, Nov 06 2023
  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[1000],!UnsameQ@@Total/@Union[Subsets[prix[#],{2}]]&]

Formula

These are numbers k such that A086971(k) > A366739(k).

A325793 Positive integers whose number of divisors is equal to their sum of prime indices.

Original entry on oeis.org

3, 10, 28, 66, 70, 88, 208, 228, 306, 340, 364, 490, 495, 525, 544, 550, 675, 744, 870, 966, 1160, 1216, 1242, 1254, 1288, 1326, 1330, 1332, 1672, 1768, 1785, 1870, 2002, 2064, 2145, 2295, 2457, 2900, 2944, 3250, 3280, 3430, 3468, 3540, 3724, 4125, 4144, 4248
Offset: 1

Views

Author

Gus Wiseman, May 23 2019

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, with sum A056239(n).

Examples

			The term 70 is in the sequence because it has 8 divisors {1, 2, 5, 7, 10, 14, 35, 70} and its sum of prime indices is also 1 + 3 + 4 = 8.
The sequence of terms together with their prime indices begins:
     3: {2}
    10: {1,3}
    28: {1,1,4}
    66: {1,2,5}
    70: {1,3,4}
    88: {1,1,1,5}
   208: {1,1,1,1,6}
   228: {1,1,2,8}
   306: {1,2,2,7}
   340: {1,1,3,7}
   364: {1,1,4,6}
   490: {1,3,4,4}
   495: {2,2,3,5}
   525: {2,3,3,4}
   544: {1,1,1,1,1,7}
   550: {1,3,3,5}
   675: {2,2,2,3,3}
   744: {1,1,1,2,11}
   870: {1,2,3,10}
   966: {1,2,4,9}
		

Crossrefs

Positions of 0's in A325794.
Contains A239885 except for 1.

Programs

  • Maple
    filter:= proc(n) local F,t;
      F:= ifactors(n)[2];
      add(numtheory:-pi(t[1])*t[2],t=F) = mul(t[2]+1,t=F)
    end proc:
    select(filter, [$1..10000]); # Robert Israel, Oct 16 2023
  • Mathematica
    Select[Range[100],DivisorSigma[0,#]==Total[Cases[FactorInteger[#],{p_,k_}:>PrimePi[p]*k]]&]

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

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

A319685 Number of distinct values obtained when arithmetic derivative (A003415) is applied to proper divisors of n.

Original entry on oeis.org

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

Views

Author

Antti Karttunen, Oct 02 2018

Keywords

Comments

Because for all d|n, dA003415(d) < A003415(n), it follows that the terms here are one less than in A319686.
Differs from A033273(n) = A000005(n) - A001221(n) at n = 1, 112, 156, 224, 280, 312, 336, 342, 380, 448, 468, 525, 560, 608, 624, 660, 672, 684, 756, 760, 780, 784, 840, 870, 896, 936, 984, 1008, ...

Examples

			The proper divisors of 112 are [1, 2, 4, 7, 8, 14, 16, 28, 56]. Applying arithmetic derivative A003415 to these, we obtain values [0, 1, 4, 1, 12, 9, 32, 32, 92], of which only 7 are distinct: 0, 1, 4, 9, 12, 32, and 92, thus a(112) = 7.
		

Crossrefs

One less than A319686.
Cf. A003415.
Cf. also A304793, A305611, A316555, A316556, A319695 for similarly constructed sequences.

Programs

  • Mathematica
    d[0] = d[1] = 0; d[n_] := d[n] = n * Plus @@ ((Last[#]/First[#]) & /@ FactorInteger[n]); a[n_] := CountDistinct[d /@ Most[Divisors[n]]]; Array[a, 100] (* Amiram Eldar, Apr 17 2024 *)
  • PARI
    A003415(n) = {my(fac); if(n<1, 0, fac=factor(n); sum(i=1, matsize(fac)[1], n*fac[i, 2]/fac[i, 1]))}; \\ From A003415
    A319685(n) = { my(m=Map(),s,k=0); fordiv(n,d,if((dA003415(d)), mapput(m,s,s); k++)); (k); };

Formula

a(n) = A319686(n)-1.

A325791 Number of necklace permutations of {1..n} such that every positive integer from 1 to n * (n + 1)/2 is the sum of some circular subsequence.

Original entry on oeis.org

1, 1, 1, 2, 4, 20, 82, 252, 1074, 7912, 39552, 152680, 776094, 5550310, 30026848, 108376910
Offset: 0

Views

Author

Gus Wiseman, May 23 2019

Keywords

Comments

A necklace permutation is a permutation that is either empty or whose first part is the minimum. A circular subsequence is a sequence of consecutive terms where the last and first parts are also considered consecutive. The only circular subsequence of maximum length is the sequence itself, not any rotation of it. For example, the circular subsequences of (1,3,2) are: (), (1), (2), (3), (1,3), (2,1), (3,2), (1,3,2).

Examples

			The a(1) = 1 through a(5) = 20 permutations:
  (1)  (1,2)  (1,2,3)  (1,2,3,4)  (1,2,3,4,5)
              (1,3,2)  (1,3,2,4)  (1,2,3,5,4)
                       (1,4,2,3)  (1,2,4,3,5)
                       (1,4,3,2)  (1,2,4,5,3)
                                  (1,2,5,4,3)
                                  (1,3,2,5,4)
                                  (1,3,4,2,5)
                                  (1,3,4,5,2)
                                  (1,3,5,2,4)
                                  (1,3,5,4,2)
                                  (1,4,2,3,5)
                                  (1,4,2,5,3)
                                  (1,4,3,2,5)
                                  (1,4,5,2,3)
                                  (1,4,5,3,2)
                                  (1,5,2,3,4)
                                  (1,5,2,4,3)
                                  (1,5,3,2,4)
                                  (1,5,3,4,2)
                                  (1,5,4,3,2)
		

Crossrefs

Programs

  • Mathematica
    subalt[q_]:=Union[ReplaceList[q,{_,s__,_}:>{s}],DeleteCases[ReplaceList[q,{t___,,u___}:>{u,t}],{}]];
    Table[Length[Select[Permutations[Range[n]],#=={}||First[#]==1&&Range[n*(n+1)/2]==Union[Total/@subalt[#]]&]],{n,0,5}]

Extensions

a(11)-a(15) from Bert Dobbelaere, Nov 01 2020

A316556 Number of distinct LCMs 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, 4, 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, 4, 1, 2, 2, 1, 2, 4, 1, 2, 3, 4, 1, 2, 1, 2, 3, 2, 3, 3, 1, 2, 1, 2, 1, 3, 3, 2, 2, 2, 1, 4, 3, 2, 3, 2, 3, 2, 1, 2, 3, 2, 1, 4, 1, 2, 5
Offset: 1

Views

Author

Gus Wiseman, Jul 06 2018

Keywords

Comments

The Heinz number of an integer partition (y_1,...,y_k) is prime(y_1)*...*prime(y_k).
Number of distinct values obtained when A290103 is applied to all divisors of n larger than one. - Antti Karttunen, Sep 25 2018

Examples

			462 is the Heinz number of (5,4,2,1) which has possible LCMs of nonempty submultisets {1,2,4,5,10,20} so a(462) = 6.
		

Crossrefs

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

Programs

  • Mathematica
    Table[Length[Union[LCM@@@Rest[Subsets[If[n==1,{},Cases[FactorInteger[n],{p_,k_}:>PrimePi[p]]]]]]],{n,100}]
  • PARI
    A290103(n) = lcm(apply(p->primepi(p),factor(n)[,1]));
    A316556(n) = { my(m=Map(),s,k=0); fordiv(n,d,if((d>1)&&!mapisdefined(m,s=A290103(d)), mapput(m,s,s); k++)); (k); }; \\ Antti Karttunen, Sep 25 2018

Extensions

More terms from Antti Karttunen, Sep 25 2018

A367097 Least positive integer whose multiset of prime indices has exactly n distinct semi-sums.

Original entry on oeis.org

1, 4, 12, 30, 60, 210, 330, 660, 2730, 3570, 6270, 12540, 53130, 79170, 110670, 221340, 514140, 1799490, 2284590, 4196010, 6750870, 13501740, 37532220, 97350330, 131362770, 189620970, 379241940, 735844830, 1471689660
Offset: 0

Views

Author

Gus Wiseman, Nov 09 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.
From David A. Corneth, Nov 15 2023: (Start)
Terms are cubefree.
bigomega(a(n)) = A001222(a(n)) >= A002024(n) + 1 = floor(sqrt(2n) + 1/2) + 1 for n > 0. (End)

Examples

			The prime indices of 60 are {1,1,2,3}, with four semi-sums {2,3,4,5}, and 60 is the first number whose prime indices have four semi-sums, so a(4) = 60.
The terms together with their prime indices begin:
       1: {}
       4: {1,1}
      12: {1,1,2}
      30: {1,2,3}
      60: {1,1,2,3}
     210: {1,2,3,4}
     330: {1,2,3,5}
     660: {1,1,2,3,5}
    2730: {1,2,3,4,6}
    3570: {1,2,3,4,7}
    6270: {1,2,3,5,8}
   12540: {1,1,2,3,5,8}
   53130: {1,2,3,4,5,9}
   79170: {1,2,3,4,6,10}
  110670: {1,2,3,4,7,11}
  221340: {1,1,2,3,4,7,11}
  514140: {1,1,2,3,5,8,13}
		

Crossrefs

The non-binary version is A259941, firsts of A299701.
These are the positions of first appearances in A366739.
A001222 counts prime factors (or prime indices), distinct A001221.
A001358 lists semiprimes, squarefree A006881, complement A100959.
A056239 adds up prime indices, row sums of A112798.
A299702 ranks knapsack partitions, counted by A108917.
A366738 counts semi-sums of partitions, strict A366741.
Semiprime divisors are listed by A367096 and have:
- square count: A056170
- sum: A076290
- squarefree count: A079275
- count: A086971
- firsts: A220264

Programs

  • Mathematica
    nn=1000;
    w=Table[Length[Union[Total/@Subsets[prix[n],{2}]]],{n,nn}];
    spnm[y_]:=Max@@NestWhile[Most,y,Union[#]!=Range[0,Max@@#]&];
    v=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 A367097(n): return next(k for k in count(1) if len({sum(s) for s in multiset_combinations({primepi(i):j for i,j in factorint(k).items()},2)}) == n) # Chai Wah Wu, Nov 13 2023

Formula

2 | a(n) for n > 0. - David A. Corneth, Nov 13 2023

Extensions

a(17)-a(22) from Chai Wah Wu, Nov 13 2023
a(23)-a(28) from David A. Corneth, Nov 13 2023
Previous Showing 11-20 of 34 results. Next