Проделаем сначала численный эксперимент. Поделите оцениваемый остаток суммы

на

и гляньте на
график функции, хотя бы при

. Ничего не напоминает?

Это гауссиан. И сразу видно, что главный вклад в сумму

вносит окрестность

, что и указывал Sonic86. В итоге общая идея - свести эту сумму к гауссиану (ну а тогда, можно надеяться, что вместо суммирования в итоге можно будет перейти к интегрированию). Практическая реализация этого плана:
Полагаем

. (Здесь используется тот факт, что характерная ширина нашего гауссиана имеет порядок

. Этот

не обязательно угадывать - можно работать и по схеме Sonic86:

, просто с

результат лучше выглядит). После манипуляций со Стирлингом и т.п., получим в главном порядке по

:
Интересующая нас сумма по

с шагом 1 в окрестности

переписывается как интегральная сумма от функции

с шагом

. При

пределы интегрирования по

можно заменить на

и мы получаем

Таким, образом, как и было обещано, два первых члена асимптотики вашей суммы имеют вид

Любопытно, что главный член определяется первым членом ряда, а второй член - суммой многих членов центрального участка. Это правильный ответ, но для его совершенно строгого обоснования здесь не хватает некоторой оценки суммы начального участка ряда от

(и конечного участка тоже - но там все проще, т.к. члены ряда меньше из-за множителя

). Следующая поправка имеет порядок

.