-(x^-n ExpIntegralE[-n, -((1 + n) Log[x])] Log[x] (x Log[x])^n)/n!