d[0] = 1, d[1] = 0,
d[n] = (n-1)(d_[n-1]+d_[n-2]).
answer:
L = [1,0] for n in range(2,10): L.append((n-1)*(L[n-1]+L[n-2])) print L # optional