|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,1
|
|
LINKS
|
|
|
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]
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
Klaus Strassburger (strass(AT)ddfi.uni-duesseldorf.de), Oct 16 2001
|
|
STATUS
|
approved
|
|
|
|