The finite Dirichlet series of the title are defined by the condition that they vanish at as many initial zeros of the zeta function as possible. It turns out that such series can produce extremely good approximations to the values of Riemanns zeta function inside the critical strip. In addition, the coefficients of these series have remarkable number-theoretical properties discovered in large-scale high-precision numerical experiments. So far, we have found no theoretical explanation for the observed phenomena.