A064627 Positive integers not in A064494.

2, 3, 5, 7, 10, 11, 13, 17, 19, 21, 23, 25, 27, 29, 30, 31, 37, 41, 42, 43, 45, 47, 50, 53, 54, 59, 61, 63, 65, 66, 67, 69, 71, 73, 74, 78, 79, 83, 85, 86, 89, 93, 97, 99, 101, 103, 105, 107, 109, 110, 112, 113, 115, 117, 119, 123, 126, 127, 129, 131, 135, 137, 138

Offset: 1

Author

Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Oct 16 2001

Keywords

• nonn

Crossrefs

Cf. A064494.

Programs

SageMath

def divsign(s, k):
    if not k.divides(s): return 0
    return (-1)^(s//k)*k
def A(n):
    s = n + 1
    for k in srange(n, 1, -1):
        s -= divsign(s, k)
    return s
# Use with caution: search range must be adjusted as necessary!
def A064627List(size):
    return sorted(Set([A(n) for n in (1..3*size)]))[0:size]
print(A064627List(63)) # Peter Luschny, Sep 16 2019