Пане /KIT\, ну не тупо ж мініфікувати (до речі, ви дволітерні ідентифікатори використовуєте, коли ще не вичерпалися однолітерні). Принаймні, не до основних оптимізацій.
Наприклад, зараз я бачу, що parallelogram_matches завжди має хоча б пару однакових параметрів (а використовує вона їх симетрично). Отже, можна зробити

def parallelogram_matches(x,z):
    return x*(x+1)*(z+1)+(x+1)*x*(z+1)+(x+1)*(x+1)*z

а це далі скорочується до

def parallelogram_matches(x,z):
    return 2*x*(x+1)*(z+1)+(x+1)*(x+1)*z

можна ще спробувати

def parallelogram_matches(x,z):
    return (x+1)*(2*x*(z+1)+(x+1)*z)

А якщо дужки в зовнішних розкрити, то матимемо

def parallelogram_matches(x,z):
    return (x+1)*(3*x*z+2*x+z)

поки що виходить з 85 символів 63 проти ваших 57 - але ж це ще без жодної мініфікації, та й ще з десяток символів далі по коду вилетить, бо тепер замість parallelogram_matches(s+1,s+1,s) пишемо parallelogram_matches(s+1,s).
До речі, я якогось біса тут parallelogram замість parallelepiped написав. А насправді ж треба писати box, і тоді

def box_matches(x,z):
    return (x+1)*(3*x*z+2*x+z)

53 символи без мініфікації - проти ваших 57.
Ще можна подумати, як ефективно отримувати найбільше число, 2-ий чи 3-ій степіть котрого не перевищує задане. Оскільки 1/3 трохи менше за третину, int(125**(1/3)) дає 4; через це мені довелося дописати оте

    side = int(n**(1/3))
    if (side+1)**3==n:
        side+=1

Але, гадаю, це можна вирішити якось простіше (і коротше); скажімо,

side = int((n+.5)**(1/3))

і т.д.

Я зрозумів, що зарано взявся за цю задачку... Вона занадто складна.
Але все-одно - дуже дякую) Повернусь до неї пізніше...

А ви гадаєте, str створюється без циклу та ділення?
Ну і так, схоже, на Python задача не розв'язується - 2*10^9 дає не менше 2*10^7 операцій, а ділення обчислюється не дуже швидко. Навіть divmod не рятує. Можна, звісно, погратися у віднімання від масиву, але навряд чи там вийде краще, а той факт, що всі числа, крім першого n, діляться на 9, не дуже допомагає.

ну от наприклад - варіант скорочення:
Ви якось ту суму цифр рахуєте складно
якщо вже робити це через  літерал числа, то так простіше:

>>> number = 123456
>>> print(sum(map(int, list(str(number)))))
21
>>> 

поміряємо швидкість:

from time import time
number = ''.join([str(i) for i in range(10**3)])
number = int(number)
print(number)
number_copy = number

Ваш код:

t_start = time()
m = 0
t = len(str(number))
for i in range(t):
    m = m + (number % 10)
    number = number // 10
print(m)
print(time() - t_start)

13500
0.015638351440429688

короткий варіант:

t_start = time()
print(sum(map(int, list(str(number_copy)))))
print(time() - t_start)

13500
0.0008478164672851562

ну,тут ще можна трохи оптимізувати, вернуся до компа  - наишу.але, нефакт що поможе
типу
x, y = y, x+1

ну, принаймі хоч код спроститься:

inp, outp = int(input()), 0
while inp > 0:
    inp -= sum(map(int, list(str(inp))))
    outp += 1    
print(outp)

Та там просто неможливо цю задачку зробити на python і щоб вона вклалася в 2 секунди. Як не крути...

без мапів та конвертації в стрінги:

def get_sum(number):
    summ = 0
    while number:
        number, r = divmod(number,10)
        summ += r
    return summ

inp, outp = int(input()), 0
while inp > 0:
    inp -= get_sum(inp)
    outp += 1    
    
print(outp)

але  while - повільна штука, тому це не поможе

Виключні ситуації повільніші за if.