We present a method to construct number-conserving Hamiltonians whose ground states exactly reproduce an arbitrarily chosen BCS-type mean-field state. Such parent Hamiltonians can be constructed not only for the usual s -wave BCS state, but also for more exotic states of this form, including the ground states of Kitaev wires and two-dimensional topological superconductors. This method leads to infinite families of locally interacting fermion models with exact topological superconducting ground states. After explaining the general technique, we apply this method to construct two specific classes of models. The first one is a one-dimensional double wire lattice model with Majorana-like degenerate ground states. The second one is a two-dimensional p x + i p y superconducting model, where we also obtain analytic expressions for topologically degenerate ground states in the presence of vortices. Our models may provide a deeper conceptual understanding of how Majorana zero modes could emerge in condensed matter systems, as well as inspire novel routes to realize them in experiment.