Простыми преобразованиями эту задачу не решить, будем использовать арифметику остатков.
1-ое свойство, которое понадобится
То есть мы спокойно можем заменить каждое слагаемое сравнимым с ним по модулю m. То есть каждое слагаемое в нашей сумме будем рассматривать отдельно.
2-ое свойство, которое нам понадобится:
То есть довольно аналогичная вещь в произведении
На нашем примере все увидим
Находим остатки по модулю 31
Рассматриваем первое слагаемое. Просто двойка не годится, нам нужно найти ближайшее к 31 число, превосходящее его (иногда там в отрицательные числа залезаем, например, 
, но сейчас это не нужно), нам повезло, это 32
Учитываем, что 
, получаем
То есть остаток от деления первого слагаемое на 31 получился равным 10. Прекрасно, аналогично со вторым
Остаток 21, чудесно. Выполняем последний шаг.
То есть остаток от деления исходного числа на 31 равен 0, следовательно, исходное число делится на 31, что и требовалось доказать.
S(16) = (2a1 + 15d)*16/2 = (2a1 + 15d)*8 = 144
Получаем систему
{ a1 + 6d = 6
{ 2a1 + 15d = 18
Умножаем 1 уравнение на -2
{ -2a1 - 12d = -12
{ 2a1 + 15d = 18
Складываем уравнения
3d = 6; d = 2
a1 = 6 - 6d = 6 - 12 = -6
Сумма 40 членов
S(40) = (2a1 + 39d)*40/2 = (-2*6 + 39*2)*20 = 66*20 = 1320