OFFSET
1,2
COMMENTS
The list contains all primes p (A000040) because Sum_{k = 1..p} gcd(p, floor(p / k)) = 2*p - 1 and Sum_{k = 1..p} gcd(p, ceiling(p / k)) = 2*p - 1.
EXAMPLE
m = 5: Sum_{k = 1..5} gcd(5, floor(5 / k)) = 9, Sum_{k = 1..5} gcd(5, ceiling(5 / k)) = 9, 9 = 9, thus m = 5 is a term.
m = 35: Sum_{k = 1..35} gcd(35, floor(35 / k)) = 83, Sum_{k = 1..35} gcd(35, ceiling(35 / k)) = 83, 83 = 83, thus m = 35 is a term.
PROG
(PARI) isok(m) = sum(k=1, m, gcd(m, floor(m/k))) == sum(k=1, m, gcd(m, ceil(m/k))); \\ Michel Marcus, Jun 28 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
Ctibor O. Zizka, Jun 27 2025
STATUS
approved
