for のループ回数と掛け算の回数を減らして、代わりに set の操作を使ってみました。
どちらが効率が良いかは計測してないのでわかりません。
import math
X = 10
Y = 15
Z = 20
n = 10
m = 4
points = set()
# 0<=x<=y<=z の座標を列挙
vals = list(map(lambda i: i * m, range(n // m + 1)))
for xi, x in enumerate(vals):
for yi, y in enumerate(vals[xi:]):
for zi, z in enumerate(vals[yi:]):
if x * x + y * y + z * z <= n * n:
points.add((x, y, z))
# x, y, z の座標を入れ替えても球の中にある
#points |= set(map(lambda p: (p[0], p[1], p[2]), points))
points |= set(map(lambda p: (p[0], p[2], p[1]), points))
points |= set(map(lambda p: (p[1], p[2], p[0]), points))
points |= set(map(lambda p: (p[1], p[0], p[2]), points))
points |= set(map(lambda p: (p[2], p[0], p[1]), points))
points |= set(map(lambda p: (p[2], p[1], p[0]), points))
# x, y, z の符号を反転させても球の中にある
points |= set(map(lambda p: ( p[0], p[1], -p[2]), points))
points |= set(map(lambda p: ( p[0], -p[1], p[2]), points))
points |= set(map(lambda p: (-p[0], p[1], p[2]), points))
# 中心を移動
points = [ ( p[0] + X, p[1] + Y, p[2] + Z ) for p in points ]
print(points)