Быстрое преобразование Фурье в задаче с кратными узлами

Спектральная плотность функции, заданной в $n$ узлах равномерной сетки вместе со своими производными первого порядка, определяется за $O(n\ln n)$ операций.