In this letter we present, in a number conserving framework, a model of interacting fermions in a two-wire geometry supporting non-local zero-energy Majorana-like edge excitations. The model has an exactly solvable line, on varying the density of fermions, described by a topologically non-trivial ground state wave-function. Away from the exactly solvable line we study the system by means of the numerical density matrix renormalization group. We characterize its topological properties through the explicit calculation of a degenerate entanglement spectrum and of the braiding operators which are exponentially localized at the edges. Furthermore, we establish the presence of a gap in its single particle spectrum while the Hamiltonian is gapless, and compute the correlations between the edge modes as well as the superfluid correlations. The topological phase covers a sizeable portion of the phase diagram, the solvable line being one of its boundaries.Introduction -Large part of the enormous attention devoted in the last years to topological superconductors owes to the exotic quasiparticles such as Majorana modes, which localize at their boundaries (edges, vortices, . . . ) [1,2] and play a key role in several robust quantum information protocols [3]. Kitaev's p-wave superconducting quantum wire [4] provides a minimal setting showcasing all the key aspects of topological states of matter in fermionic systems. The existence of a socalled "sweet point" supporting an exact and easy-tohandle analytical solution puts this model at the heart of our understanding of systems supporting Majorana modes. Various implementations in solid state [5,6] and ultracold atoms [7,8] via proximity to superconducting or superfluid reservoirs have been proposed, and experimental signatures of edge modes were reported [9]. Kitaev's model is an effective mean-field model and its Hamiltonian does not commute with the particle number operator. Considerable activity has been devoted to understanding models supporting Majorana edge modes in a number-conserving setting [10][11][12][13][14], as in various experimental platforms (e.g. solid state [10,11] or ultracold atoms [12,13]) this property is naturally present. It was realised that a simple way to promote particle number conservation to a symmetry of the model, while keeping the edge state physics intact, was to consider at least two coupled wires rather than a single one [10][11][12]. However, since attractive interactions are pivotal to generate superconducting order in the canonical ensemble, one usually faces a complex interacting many-body problem. Therefore, approximations such as bosonization [10][11][12], or numerical approaches [13] were invoked. An exactly solvable model of a topological superconductor in a number conserving setting, which would directly complement Kitaev's scenario, is missing (see however [14]).In this letter we present an exactly solvable model of a topological superconductor which supports exotic Majorana-like quasiparticles at its ends and retains the fermionic ...