Разбиваем все точки на 4 группы (по четвертям 1, 2, 3, 4). Точки на осях просто не учитываются. А дальше выражаем ответ через эти 4 значения.
def one_to_n(n):
n -= 1
res = n * (n + 1) // 2
return res
s = 0
if a1 + a2 > 2:
s += one_to_n(a1) * a2
s += one_to_n(a2) * a1
if a3 + a4 > 2:
s += one_to_n(a3) * a4
s += one_to_n(a4) * a3
print(s)