A026422 -id:A026422 - cp's OEIS frontend

A026424 Number of prime divisors (counted with multiplicity) is odd; Liouville function lambda(n) (A008836) is negative.

2, 3, 5, 7, 8, 11, 12, 13, 17, 18, 19, 20, 23, 27, 28, 29, 30, 31, 32, 37, 41, 42, 43, 44, 45, 47, 48, 50, 52, 53, 59, 61, 63, 66, 67, 68, 70, 71, 72, 73, 75, 76, 78, 79, 80, 83, 89, 92, 97, 98, 99, 101, 102, 103, 105, 107, 108, 109, 110, 112

Offset: 1

N. J. A. Sloane

Neither this sequence nor its complement (A028260) contains any infinite arithmetic progression. - Franklin T. Adams-Watters, Sep 05 2008
A066829(a(n)) = 1. - Reinhard Zumkeller, Jun 26 2009
These numbers can be generated by the sieving process described in A066829. - Reinhard Zumkeller, Jul 01 2009
Lexicographically earliest sequence of distinct nonnegative integers with no term being the product of any two not necessarily distinct terms. The equivalent sequence for addition/subtraction is A005408 (the odd numbers), for exponentiation is A259444, and for binary exclusive OR is A000069. - Peter Munn, Mar 16 2018
The equivalent lexicographically earliest sequence with no term being the product of any two distinct terms is A026416. A000028 is similarly the equivalent sequence when A059897 is used as multiplicative operator in place of standard integer multiplication. - Peter Munn, Mar 16 2019

T. D. Noe, Table of n, a(n) for n = 1..10000
S. Ramanujan, Irregular numbers, J. Indian Math. Soc., 5 (1913), 105-106; Coll. Papers 20-21.
Eric Weisstein's World of Mathematics, Prime Sums
Index entries for sequences generated by sieves - _Reinhard Zumkeller_, Jul 01 2009

Cf. A008836, A028260 (complement).
Apart from initial term, same as A026422.
Cf. A026416 and cross-references therein.
Cf. A000028, A000069, A005408, A259444.

a026424 n = a026424_list !! (n-1)
a026424_list = filter (odd . a001222) [1..]
-- Reinhard Zumkeller, Oct 05 2011

isA026424 := proc(n)
    if type(numtheory[bigomega](n) ,'odd') then
        true;
    else
        false;
    end if;
end proc:
A026424 := proc(n)
    option remember;
    if n =1 then
        2;
    else
        for a from procname(n-1)+1 do
            if isA026424(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, May 25 2017

Select[Range[2, 112], OddQ[Total[FactorInteger[#]][[2]]] &] (* T. D. Noe, May 07 2011 *)
(* From version 7 on *) Select[Range[2, 112], LiouvilleLambda[#] == -1 &] (* Jean-François Alcover, Aug 19 2013 *)
Select[Range[150],OddQ[PrimeOmega[#]]&] (* Harvey P. Dale, Oct 04 2024 *)

is(n)=bigomega(n)%2 \\ Charles R Greathouse IV, Sep 16 2015

from math import isqrt, prod
from sympy import primerange, integer_nthroot, primepi
def A026424(n):
    def g(x,a,b,c,m): yield from (((d,) for d in enumerate(primerange(b,isqrt(x//c)+1),a)) if m==2 else (((a2,b2),)+d for a2,b2 in enumerate(primerange(b,integer_nthroot(x//c,m)[0]+1),a) for d in g(x,a2,b2,c*b2,m-1)))
    def f(x): return int(n+1+sum(sum(primepi(x//prod(c[1] for c in a))-a[-1][0] for a in g(x,0,1,1,m)) for m in range(2,x.bit_length()+1,2)))
    m, k = n, f(n)
    while m != k: m, k = k, f(k)
    return m # Chai Wah Wu, Apr 10 2025

Sum 1/a(n)^m = (zeta(m)^2-zeta(2m))/(2*zeta(m)), Dirichlet g.f. of A066829. - Ramanujan.
n>=2 is in sequence if n is not the product of two smaller elements. - David W. Wilson, May 06 2005
A001222(a(n)) mod 2 = 1. - Reinhard Zumkeller, Oct 05 2011
Union of A000040, A014612, A014614, A046308 etc. - R. J. Mathar, Jul 09 2012

A026416 A 2-way classification of integers: a(1) = 1, a(2) = 2 and for n > 2, a(n) is the smallest number not of the form a(i)*a(j) for 1 <= i < j < n.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 9, 11, 13, 16, 17, 19, 23, 24, 25, 29, 30, 31, 37, 40, 41, 42, 43, 47, 49, 53, 54, 56, 59, 61, 66, 67, 70, 71, 73, 78, 79, 81, 83, 88, 89, 97, 101, 102, 103, 104, 105, 107, 109, 110, 113, 114, 121, 127, 128, 130, 131, 135, 136, 137, 138, 139

Offset: 1

Clark Kimberling

An equivalent definition is: a(1) = 1, a(2) = 2; and for n > 2, a(n) = least positive integer > a(n-1) and not of the form a(i)*a(j) for 1 <= i < j < n.
a(2) to a(29) match the initial terms of A000028. [corrected by Peter Munn, Mar 15 2019]
This has a simpler definition than A000028, but the resulting pair lacks the crucial property of the A000028/A000379 pair (see the comment in A000028). - N. J. A. Sloane, Sep 28 2007
Contains (for example) 180, so is different from A123193. - Max Alekseyev, Sep 20 2007
From Vladimir Shevelev, Apr 05 2013: (Start)
1) The sequence does not contain (for example) 140, so is different from A000028.
2) Representation of numbers which are absent in the sequence as a product of two different terms of the sequence is, generally speaking, not unique. For example, 210 = 2*105 = 3*70 = 5*42 = 7*30.
(End)
Excluding a(1) = 1, the lexicographically earliest sequence of distinct nonnegative integers such that no term is a product of 2 distinct terms. Removing the latter distinctness requirement, the sequence becomes A026424; and the equivalent sequence where the product is of 2 or more distinct terms is A050376. A000028 is similarly the equivalent sequence when A059897 is used as multiplicative operator in place of standard integer multiplication. - Peter Munn, Mar 15 2019

			a(8) is not 10 because we already have 10 = 2*5. Of course all primes appear. 16 appears because 16 is not a product of earlier terms.

Michael S. Branicky, Table of n, a(n) for n = 1..10000 (terms 1..400 from Vincenzo Librandi)

Complement of A131181. Cf. A000028, A059897.
Cf. A066724, A026477, A050376, A084400.
Similar sequences with different starting conditions: A026417 (1,3), A026419 (1,4), A026420 (2,4), A026421 (3,4).
Related sequences with definition using any products (not necessarily distinct) and with various starting conditions: A026422 (1,2),A026423 (1,3), A026424 (2,3), A026425 (1,4), A026426 (2,4), A026427 (3,4).
See also families of related sequences: A026431 (excluding product-1), A026443 (excluding product+2), A026453 (excluding product-2) and references therein.

a[1]=1; a[2]=2; a[n_] := a[n] = For[k = a[n-1] + 1, True, k++, If[ FreeQ[ Table[ a[i]*a[j], {i, 1, n-2}, {j, i+1, n-1}], k], Return[k]]]; Table[a[n], {n, 1, 101}] (* Jean-François Alcover, May 16 2013 *)

from itertools import count, islice
def agen(): # generator of terms
    a, products = [1, 2], {2}
    yield from a
    for k in count(3):
        if k not in products:
            yield k
            products.update(k*a[i] for i in range(len(a)))
            a.append(k)
        products.discard(k)
print(list(islice(agen(), 62))) # Michael S. Branicky, Jun 09 2025

More terms from Max Alekseyev, Sep 23 2007

Edited by N. J. A. Sloane, Jul 13 2008 at the suggestion of R. J. Mathar and Max Alekseyev

cp's OEIS Frontend

A026424 Number of prime divisors (counted with multiplicity) is odd; Liouville function lambda(n) (A008836) is negative.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Python

Formula