Задача дискретного логарифмирования (ДЛ) скорее всего является трудной, но скорее всего не является NP-полной. Но насколько далеко ей до полноты, а если ее усложнить?
Возмем систему целочисленных уравнений от ненулевых неизвестных

с ненулевыми коэффициэнтами

и

"башенного" вида:

,
где

простое и

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

и порядок (неперестановочный) использования неизвестных

, уравнение можно записать как:

, где

,

при

,

. Обратная функция к

по

это дискретный логарифм. Нелинейность

по

делает систему труднорешаемой.
При этом спутанность неизвестных

друг с другом в одном уравнении имеет "последовательный" характер, только за счет нелинейности

, при линейной

мы получили бы линейное уравнение.
В то время как в квадратичном диофантовом уравнении (КДУ) спутанность неизвестных

друг с другом в одном уравнении имеет "паралельный" характер - гораздо сильнее:

.
Видимо поэтому их система (СКДУ) дает NP-полную задачу (поиск одного из решений, если оно есть). Системы нелинейных диофантовом уравнении без такой спутанности, т.е. с уравнениями

, скорее всего не дают NP-полноты.
Предложенная задача СБДЛ при

не решается с помощью (не сводится полиномно к) ДЛ, она более общая (для задачи факторизации такого обобщения не сделать), значит можно предположить ее NP-полноту. Или, наоборот, "последовательная" (слабая) спутанность неизвестных

делает такие задачи не NP-полными, для любой не NP-полной функции, обратной к

по

. Какие результаты или гипотезы об этом известны?
Более упрощенная задача - с одним уравнением и одним, заменяющем остальные, неизвестным:

. Она имеет промежуточную сложность между ДЛ и СБДЛ, но совсем специфична.