login
A064627
Positive integers not in A064494.
4
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,1
PROG
(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
CROSSREFS
Cf. A064494.
Sequence in context: A338174 A058590 A187711 * A324813 A235050 A117286
KEYWORD
nonn
AUTHOR
Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Oct 16 2001
STATUS
approved