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 61-70 of 114 results. Next

A019507 Droll numbers: numbers > 1 whose sum of even prime factors equals the sum of odd prime factors.

Original entry on oeis.org

72, 240, 672, 800, 2240, 4224, 5184, 6272, 9984, 14080, 17280, 33280, 39424, 48384, 52224, 57600, 93184, 116736, 161280, 174080, 192000, 247808, 304128, 373248, 389120, 451584, 487424, 537600, 565248, 585728, 640000, 718848, 1013760, 1089536, 1244160, 1384448
Offset: 1

Views

Author

Mario Velucchi (mathchess(AT)velucchi.it)

Keywords

Examples

			6272 = 2*2*2*2*2*2*2*7*7 is droll since 2+2+2+2+2+2+2 = 14 = 7+7.
		

Crossrefs

For count instead of sum we have A072978.
Partitions of this type are counted by A239261, without zero terms A249914.
For prime indices instead of factors we have A366748, zeros of A366749.
The LHS is A366839 with alternating zeros, for indices A366531, triangle A113686.
The RHS is A366840, for indices A366528, triangle A113685.
A000009 counts partitions into odd parts, ranks A066208.
A035363 counts partitions into even parts, ranks A066207.
A112798 lists prime indices, length A001222, sum A056239.
A257991 counts odd prime indices, even A257992.
A300061 lists numbers with even sum of prime indices, odd A300063.

Programs

  • Maple
    f:= proc(k, m) # numbers whose sum of prime factors >= m is k; m is prime
       local S,p,j;
       option remember;
       if k = 0 then return [1]
       elif m > k then return []
       fi;
       S:= NULL:
       p:= nextprime(m);
       for j from k by -m to 0 do
         S:= S, op(map(`*`,  procname(j,p) , m^((k-j)/m)))
       od;
       [S]
    end proc:
    g:= proc(N) local m,R;
      R:= NULL;
      for m from 1 while 2^m < N do
       R:= R, op(map(`*`,select(`<=`,f(2*m,3), N/2^m),2^m));
      od;
      sort([R])
    end proc:
    g(10^8); # Robert Israel, Feb 20 2025
  • Mathematica
    Select[Range[2, 2*10^6, 2], First[#] == Total[Rest[#]] & [Times @@@ FactorInteger[#]] &] (* Paolo Xausa, Feb 19 2025 *)
  • PARI
    isok(n) = {if (n % 2, return (0)); f = factor(n); return (2*f[1,2] == sum(i=2, #f~, f[i,1]*f[i,2]));} \\ Michel Marcus, Jun 21 2013

Formula

These are even numbers k such that A366839(k/2) = A366840(k). - Gus Wiseman, Oct 25 2023 (corrected Feb 19 2025)

Extensions

Name edited by Paolo Xausa, Feb 19 2025

A324967 Number of distinct even prime indices of n.

Original entry on oeis.org

0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 0, 1, 1, 1, 1, 0, 0, 1, 1, 0, 2, 0, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 1, 0, 1, 1, 1, 1, 2, 0, 0, 2, 1, 0, 1, 0, 0, 1, 1, 0, 1, 1, 1, 1, 0, 1, 2, 1, 0, 1, 1, 0, 2, 0, 1, 1, 0, 0, 1, 1, 1, 1, 0, 1, 1, 1, 1, 2, 1, 0, 1, 0, 0, 2, 0, 1, 2
Offset: 1

Views

Author

Gus Wiseman, Mar 21 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.
If x and y are coprime then a(x*y) = a(x) + a(y). - Robert Israel, Mar 24 2019

Examples

			180180 has prime indices {1,1,2,2,3,4,5,6}, so a(180180) = 3.
		

Crossrefs

Programs

  • Maple
    f:= proc(n) nops(select(type,map(numtheory:-pi,numtheory:-factorset(n)),even)) end proc:
    map(f, [$1..100]); # Robert Israel, Mar 24 2019
  • Mathematica
    Table[Count[If[n==1,{},FactorInteger[n]],{?(EvenQ[PrimePi[#]]&),}],{n,100}]
  • PARI
    a(n) = my(f=factor(n)[,1]); sum(k=1, #f, !(primepi(f[k]) % 2)); \\ Michel Marcus, Mar 22 2019

Formula

a(n) = A001221(n) - A324966(n). - Robert Israel, Mar 24 2019
G.f.: Sum_{k>=1} x^prime(2*k) / (1 - x^prime(2*k)). - Ilya Gutkovskiy, Feb 12 2020
Additive with a(p^e) = 1 if primepi(p) is even and 0 otherwise. - Amiram Eldar, Oct 06 2023

A372588 Numbers k > 1 such that (greatest binary index of k) + (greatest prime index of k) is odd.

Original entry on oeis.org

2, 6, 7, 8, 10, 11, 15, 18, 19, 21, 24, 26, 27, 28, 29, 32, 33, 34, 40, 41, 44, 45, 46, 47, 50, 51, 55, 59, 60, 62, 65, 70, 71, 72, 74, 76, 78, 79, 81, 84, 86, 87, 89, 91, 95, 96, 98, 101, 104, 105, 106, 107, 108, 111, 112, 113, 114, 116, 117, 122, 126, 128
Offset: 1

Views

Author

Gus Wiseman, May 14 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
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 even version is A372589.

Examples

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
        {2}   2  (1)
      {2,3}   6  (2,1)
    {1,2,3}   7  (4)
        {4}   8  (1,1,1)
      {2,4}  10  (3,1)
    {1,2,4}  11  (5)
  {1,2,3,4}  15  (3,2)
      {2,5}  18  (2,2,1)
    {1,2,5}  19  (8)
    {1,3,5}  21  (4,2)
      {4,5}  24  (2,1,1,1)
    {2,4,5}  26  (6,1)
  {1,2,4,5}  27  (2,2,2)
    {3,4,5}  28  (4,1,1)
  {1,3,4,5}  29  (10)
        {6}  32  (1,1,1,1,1)
      {1,6}  33  (5,2)
      {2,6}  34  (7,1)
      {4,6}  40  (3,1,1,1)
    {1,4,6}  41  (13)
    {3,4,6}  44  (5,1,1)
  {1,3,4,6}  45  (3,2,2)
		

Crossrefs

For sum (A372428, zeros A372427) we have A372586.
For minimum (A372437) we have A372439, complement A372440.
For length (A372441, zeros A071814) we have A372590, complement A372591.
Positions of odd terms in A372442, zeros A372436.
The complement is A372589.
For just binary indices:
- length: A000069, complement A001969
- sum: A158705, complement A158704
- minimum: A003159, complement A036554
- maximum: A053738, complement A053754
For just prime indices:
- length: A026424 (count A027193), complement A028260 (count A027187)
- sum: A300063 (count A058695), complement A300061 (count A058696)
- minimum: A340932 (count A026804), complement A340933 (count A026805)
- maximum: A244991 (count A027193), complement A244990 (count A027187)
A005408 lists odd numbers.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031368 lists odd-indexed primes, even A031215.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.

Programs

  • Mathematica
    Select[Range[2,100],OddQ[IntegerLength[#,2]+PrimePi[FactorInteger[#][[-1,1]]]]&]

Formula

Numbers k such that A070939(k) + A061395(k) is odd.

A379306 Number of squarefree prime indices of n.

Original entry on oeis.org

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

Views

Author

Gus Wiseman, Dec 25 2024

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 39 are {2,6}, so a(39) = 2.
The prime indices of 70 are {1,3,4}, so a(70) = 2.
The prime indices of 98 are {1,4,4}, so a(98) = 1.
The prime indices of 294 are {1,2,4,4}, a(294) = 2.
The prime indices of 1911 are {2,4,4,6}, so a(1911) = 2.
The prime indices of 2548 are {1,1,4,4,6}, so a(2548) = 3.
		

Crossrefs

Positions of first appearances are A000079.
Positions of zero are A379307, counted by A114374 (strict A256012).
Positions of one are A379316, counted by A379308 (strict A379309).
A000040 lists the primes, differences A001223.
A005117 lists the squarefree numbers, differences A076259.
A008966 is the characteristic function for the squarefree numbers.
A013929 lists the nonsquarefree numbers, differences A078147.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
A061398 counts squarefree numbers between primes, zeros A068360.
A377038 gives k-th differences of squarefree numbers.
Other counts of prime indices:
- A087436 postpositive, see A038550.
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379310 nonsquarefree, see A302478.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

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

Formula

Totally additive with a(prime(k)) = A008966(k).

A379317 Positive integers with a unique even prime index.

Original entry on oeis.org

3, 6, 7, 12, 13, 14, 15, 19, 24, 26, 28, 29, 30, 33, 35, 37, 38, 43, 48, 51, 52, 53, 56, 58, 60, 61, 65, 66, 69, 70, 71, 74, 75, 76, 77, 79, 86, 89, 93, 95, 96, 101, 102, 104, 106, 107, 112, 113, 116, 119, 120, 122, 123, 130, 131, 132, 138, 139, 140, 141, 142
Offset: 1

Views

Author

Gus Wiseman, Dec 29 2024

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 terms together with their prime indices begin:
   3: {2}
   6: {1,2}
   7: {4}
  12: {1,1,2}
  13: {6}
  14: {1,4}
  15: {2,3}
  19: {8}
  24: {1,1,1,2}
  26: {1,6}
  28: {1,1,4}
  29: {10}
  30: {1,2,3}
  33: {2,5}
  35: {3,4}
  37: {12}
  38: {1,8}
  43: {14}
  48: {1,1,1,1,2}
		

Crossrefs

Partitions of this type are counted by A038348 (strict A096911).
For all even parts we have A066207, counted by A035363 (strict A000700).
For no even parts we have A066208, counted by A000009 (strict A035457).
Positions of 1 in A257992.
A000040 lists the primes, differences A001223.
A055396 gives least prime index, greatest A061395.
A056239 adds up prime indices, row sums of A112798, counted by A001222.
Other counts of prime indices:
- A330944 nonprime, see A000586, A000607, A076610, A330945.
- A379311 old prime, see A204389, A320629, A379312-A379315.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    Select[Range[100],Length[Select[prix[#],EvenQ]]==1&]

A332821 One part of a 3-way classification of the positive integers. Numbers n for which A048675(n) == 1 (mod 3).

Original entry on oeis.org

2, 5, 9, 11, 12, 16, 17, 21, 23, 28, 30, 31, 39, 40, 41, 47, 49, 52, 54, 57, 59, 66, 67, 70, 72, 73, 75, 76, 83, 87, 88, 91, 96, 97, 100, 102, 103, 109, 111, 116, 126, 127, 128, 129, 130, 133, 135, 136, 137, 138, 148, 149, 154, 157, 159, 165, 167, 168, 169, 172, 175, 179, 180, 183, 184, 186, 190, 191, 197, 203, 211, 212
Offset: 1

Views

Author

Antti Karttunen and Peter Munn, Feb 25 2020

Keywords

Comments

The positive integers are partitioned between A332820, this sequence and A332822.
For each prime p, the terms include exactly one of p and p^2. The primes alternate between this sequence and A332822. This sequence has the primes with odd indexes, those in A031368.
The terms are the even numbers in A332822 halved. The terms are also the numbers m such that 5m is in A332822, and so on for alternate primes: 11, 17, 23 etc. Likewise, the terms are the numbers m such that 3m is in A332820, and so on for alternate primes: 7, 13, 19 etc.
The numbers that are half of the even terms of this sequence are in A332820, which consists exactly of those numbers. The numbers that are one third of the terms that are multiples of 3 are in A332822, which consists exactly of those numbers. For larger primes, an alternating pattern applies as described in the previous paragraph.
If we take each odd term of this sequence and replace each prime in its factorization by the next smaller prime, the resulting number is in A332822, which consists entirely of those numbers.
The product of any 2 terms of this sequence is in A332822, the product of any 3 terms is in A332820, and the product of a term of A332820 and a term of this sequence is in this sequence. So if a number k is present, k^2 is in A332822, k^3 is in A332820, and k^4 is in this sequence.
If k is an even number, exactly one of {k/2, k, 2k} is in the sequence (cf. A191257 / A067368 / A213258); and generally if k is a multiple of a prime p, exactly one of {k/p, k, k*p} is in the sequence.

Crossrefs

Positions of ones in A332823; equivalently, numbers in row 3k+1 of A277905 for some k >= 0.
Subsequences: intersection of A026478 and A066208, A031368 (prime terms), A033431\{0}, A052934\{1}, A069486, A099800, A167747\{1}, A244725\{0}, A244728\{0}, A338911 (semiprime terms).

Programs

  • Mathematica
    Select[Range@ 212, Mod[Total@ #, 3] == 1 &@ Map[#[[-1]]*2^(PrimePi@ #[[1]] - 1) &, FactorInteger[#]] &] (* Michael De Vlieger, Mar 15 2020 *)
  • PARI
    isA332821(n) =  { my(f = factor(n)); (1==((sum(k=1, #f~, f[k, 2]*2^primepi(f[k, 1]))/2)%3)); };

Formula

{a(n) : n >= 1} = {2 * A332820(k) : k >= 1} U {A003961(A332822(k)) : k >= 1}.
{a(n) : n >= 1} = {A332822(k)^2 : k >= 1} U {A331590(2, A332820(k)) : k >= 1}.

A340786 Number of factorizations of 4n into an even number of even factors > 1.

Original entry on oeis.org

1, 1, 1, 3, 1, 2, 1, 3, 2, 2, 1, 4, 1, 2, 2, 6, 1, 3, 1, 4, 2, 2, 1, 6, 2, 2, 2, 4, 1, 4, 1, 7, 2, 2, 2, 7, 1, 2, 2, 6, 1, 4, 1, 4, 3, 2, 1, 10, 2, 3, 2, 4, 1, 4, 2, 6, 2, 2, 1, 8, 1, 2, 3, 12, 2, 4, 1, 4, 2, 4, 1, 10, 1, 2, 3, 4, 2, 4, 1, 10, 3, 2, 1, 8, 2, 2
Offset: 1

Views

Author

Gus Wiseman, Jan 31 2021

Keywords

Examples

			The a(n) factorizations for n = 6, 12, 24, 36, 60, 80, 500:
  4*6   6*8      2*48      2*72      4*60      4*80          40*50
  2*12  2*24     4*24      4*36      6*40      8*40          4*500
        4*12     6*16      6*24      8*30      10*32         8*250
        2*2*2*6  8*12      8*18      10*24     16*20         10*200
                 2*2*4*6   12*12     12*20     2*160         20*100
                 2*2*2*12  2*2*6*6   2*120     2*2*2*40      2*1000
                           2*2*2*18  2*2*2*30  2*2*4*20      2*2*10*50
                                     2*2*6*10  2*2*8*10      2*2*2*250
                                               2*4*4*10      2*10*10*10
                                               2*2*2*2*2*10
		

Crossrefs

Note: A-numbers of Heinz-number sequences are in parentheses below.
Positions of ones are 1 and A000040, or A008578.
A version for partitions is A027187 (A028260).
Allowing odd length gives A108501 (bisection of A340785).
Allowing odd factors gives A339846.
An odd version is A340102.
- Factorizations -
A001055 counts factorizations, with strict case A045778.
A316439 counts factorizations by product and length.
A340101 counts factorizations into odd factors.
A340653 counts balanced factorizations.
A340831/A340832 count factorizations with odd maximum/minimum.
- Even -
A027187 counts partitions of even maximum (A244990).
A058696 counts partitions of even numbers (A300061).
A067661 counts strict partitions of even length (A030229).
A236913 counts partitions of even length and sum (A340784).
A340601 counts partitions of even rank (A340602).

Programs

  • Maple
    g:= proc(n, m, p)
     option remember;
     local F,r,x,i;
     # number of factorizations of n into even factors > m with number of factors == p (mod 2)
     if n = 1 then if p = 0 then return 1 else return 0 fi fi;
     if m > n  or n::odd then return 0 fi;
     F:= sort(convert(select(t -> t > m and t::even, numtheory:-divisors(n)),list));
     r:= 0;
     for x in F do
       for i from 1 while n mod x^i = 0 do
         r:= r + procname(n/x^i, x, (p+i) mod 2)
     od od;
     r
    end proc:
    f:= n -> g(4*n, 1, 0):
    map(f, [$1..100]); # Robert Israel, Mar 16 2023
  • Mathematica
    facs[n_]:=If[n<=1,{{}},Join@@Table[Map[Prepend[#,d]&,Select[facs[n/d],Min@@#>=d&]],{d,Rest[Divisors[n]]}]];
    Table[Length[Select[facs[4n],EvenQ[Length[#]]&&Select[#,OddQ]=={}&]],{n,100}]
  • PARI
    A340786aux(n, m=n, p=0) = if(1==n, (0==p), my(s=0); fordiv(n, d, if((d>1)&&(d<=m)&&!(d%2), s += A340786aux(n/d, d, 1-p))); (s));
    A340786(n) = A340786aux(4*n); \\ Antti Karttunen, Apr 14 2022

A352142 Numbers whose prime factorization has all odd indices and all odd exponents.

Original entry on oeis.org

1, 2, 5, 8, 10, 11, 17, 22, 23, 31, 32, 34, 40, 41, 46, 47, 55, 59, 62, 67, 73, 82, 83, 85, 88, 94, 97, 103, 109, 110, 115, 118, 125, 127, 128, 134, 136, 137, 146, 149, 155, 157, 160, 166, 167, 170, 179, 184, 187, 191, 194, 197, 205, 206, 211, 218, 227, 230
Offset: 1

Views

Author

Gus Wiseman, Mar 18 2022

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, sum A056239, length A001222.
A number's prime signature is the sequence of positive exponents in its prime factorization, which is row n of A124010, length A001221, sum A001222.
These are the Heinz numbers of integer partitions with all odd parts and all odd multiplicities, counted by A117958.

Examples

			The terms together with their prime indices begin:
   1 = 1
   2 = prime(1)
   5 = prime(3)
   8 = prime(1)^3
  10 = prime(1) prime(3)
  11 = prime(5)
  17 = prime(7)
  22 = prime(1) prime(5)
  23 = prime(9)
  31 = prime(11)
  32 = prime(1)^5
  34 = prime(1) prime(7)
  40 = prime(1)^3 prime(3)
		

Crossrefs

The restriction to primes is A031368.
The first condition alone is A066208, counted by A000009.
These partitions are counted by A117958.
The squarefree case is A258116, even A258117.
The second condition alone is A268335, counted by A055922.
The even-even version is A352141 counted by A035444.
A000290 = exponents all even, counted by A035363.
A056166 = exponents all prime, counted by A055923.
A066207 = indices all even, counted by A035363 (complement A086543).
A109297 = same indices as exponents, counted by A114640.
A112798 lists prime indices, reverse A296150, length A001222, sum A056239.
A124010 gives prime signature, sorted A118914, length A001221, sum A001222.
A162641 counts even prime exponents, odd A162642.
A257991 counts odd prime indices, even A257992.
A325131 = disjoint indices from exponents, counted by A114639.
A346068 = indices and exponents all prime, counted by A351982.
A351979 = odd indices with even exponents, counted by A035457.
A352140 = even indices with odd exponents, counted by A055922 aerated.
A352143 = odd indices with odd conjugate indices, counted by A053253 aerated.

Programs

  • Mathematica
    Select[Range[100],#==1||And@@OddQ/@PrimePi/@First/@FactorInteger[#]&&And@@OddQ/@Last/@FactorInteger[#]&]
  • Python
    from itertools import count, islice
    from sympy import primepi, factorint
    def A352142_gen(startvalue=1): # generator of terms >= startvalue
        return filter(lambda k:all(map(lambda x:x[1]%2 and primepi(x[0])%2, factorint(k).items())),count(max(startvalue,1)))
    A352142_list = list(islice(A352142_gen(),30)) # Chai Wah Wu, Mar 18 2022

Formula

Intersection of A066208 and A268335.
A257991(a(n)) = A001222(a(n)).
A162642(a(n)) = A001221(a(n)).
A257992(a(n)) = A162641(a(n)) = 0.

A366848 Odd numbers whose odd prime indices are relatively prime.

Original entry on oeis.org

55, 85, 155, 165, 187, 205, 253, 255, 275, 295, 335, 341, 385, 391, 415, 425, 451, 465, 485, 495, 527, 545, 561, 595, 605, 615, 635, 649, 697, 713, 715, 737, 745, 759, 765, 775, 785, 799, 803, 825, 885, 895, 913, 935, 943, 955, 1003, 1005, 1023, 1025, 1045
Offset: 1

Views

Author

Gus Wiseman, Nov 01 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.

Examples

			The odd prime indices of 345 are {3,9}, which are not relatively prime, so 345 is not in the sequence.
The odd prime indices of 825 are {3,3,5}, which are relatively prime, so 825 is in the sequence
The terms together with their prime indices begin:
    55: {3,5}
    85: {3,7}
   155: {3,11}
   165: {2,3,5}
   187: {5,7}
   205: {3,13}
   253: {5,9}
   255: {2,3,7}
   275: {3,3,5}
   295: {3,17}
   335: {3,19}
   341: {5,11}
   385: {3,4,5}
   391: {7,9}
   415: {3,23}
   425: {3,3,7}
   451: {5,13}
   465: {2,3,11}
   485: {3,25}
   495: {2,2,3,5}
		

Crossrefs

Including even terms and prime indices gives A289509, ones of A289508, counted by A000837.
Including even prime indices gives A302697, counted by A302698.
Including even terms gives A366846, counted by A366850.
For halved even instead of odd prime indices we have A366849.
A000041 counts integer partitions, strict A000009 (also into odds).
A066208 lists numbers with all odd prime indices, even A066207.
A112798 lists prime indices, length A001222, sum A056239.
A257991 counts odd prime indices, even A257992.
A366528 adds up odd prime indices, partition triangle A113685.
A366531 = 2*A366533 adds up even prime indices, triangle A113686/A174713.

Programs

  • Mathematica
    Select[Range[1000], OddQ[#]&&GCD@@Select[PrimePi/@First/@FactorInteger[#], OddQ]==1&]

A372586 Numbers k such that (sum of binary indices of k) + (sum of prime indices of k) is odd.

Original entry on oeis.org

1, 2, 3, 4, 5, 8, 9, 12, 15, 16, 17, 20, 21, 29, 32, 36, 42, 43, 45, 46, 47, 48, 51, 53, 54, 55, 59, 60, 61, 63, 64, 65, 66, 67, 68, 71, 73, 78, 79, 80, 81, 84, 89, 91, 93, 94, 95, 97, 99, 101, 105, 110, 111, 113, 114, 115, 116, 118, 119, 121, 122, 125, 127
Offset: 1

Views

Author

Gus Wiseman, May 14 2024

Keywords

Comments

A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793.
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 even version is A372587.

Examples

			The terms (center), their binary indices (left), and their weakly decreasing prime indices (right) begin:
            {1}   1  ()
            {2}   2  (1)
          {1,2}   3  (2)
            {3}   4  (1,1)
          {1,3}   5  (3)
            {4}   8  (1,1,1)
          {1,4}   9  (2,2)
          {3,4}  12  (2,1,1)
      {1,2,3,4}  15  (3,2)
            {5}  16  (1,1,1,1)
          {1,5}  17  (7)
          {3,5}  20  (3,1,1)
        {1,3,5}  21  (4,2)
      {1,3,4,5}  29  (10)
            {6}  32  (1,1,1,1,1)
          {3,6}  36  (2,2,1,1)
        {2,4,6}  42  (4,2,1)
      {1,2,4,6}  43  (14)
      {1,3,4,6}  45  (3,2,2)
      {2,3,4,6}  46  (9,1)
    {1,2,3,4,6}  47  (15)
          {5,6}  48  (2,1,1,1,1)
		

Crossrefs

Positions of odd terms in A372428, zeros A372427.
For minimum (A372437) we have A372439, complement A372440.
For length (A372441, zeros A071814) we have A372590, complement A372591.
For maximum (A372442, zeros A372436) we have A372588, complement A372589.
The complement is A372587.
For just binary indices:
- length: A000069, complement A001969
- sum: A158705, complement A158704
- minimum: A003159, complement A036554
- maximum: A053738, complement A053754
For just prime indices:
- length: A026424 (count A027193), complement A028260 (count A027187)
- sum: A300063 (count A058695), complement A300061 (count A058696)
- minimum: A340932 (count A026804), complement A340933 (count A026805)
- maximum: A244991 (count A027193), complement A244990 (count A027187)
A005408 lists odd numbers.
A019565 gives Heinz number of binary indices, adjoint A048675.
A029837 gives greatest binary index, least A001511.
A031368 lists odd-indexed primes, even A031215.
A048793 lists binary indices, length A000120, reverse A272020, sum A029931.
A061395 gives greatest prime index, least A055396.
A070939 gives length of binary expansion.
A112798 lists prime indices, length A001222, reverse A296150, sum A056239.

Programs

  • Mathematica
    prix[n_]:=If[n==1,{},Flatten[Cases[FactorInteger[n],{p_,k_}:>Table[PrimePi[p],{k}]]]];
    bix[n_]:=Join@@Position[Reverse[IntegerDigits[n,2]],1];
    Select[Range[100],OddQ[Total[bix[#]]+Total[prix[#]]]&]

Formula

Numbers k such that A029931(k) + A056239(k) is odd.
Previous Showing 61-70 of 114 results. Next