We introduce in this paper an exact solvable BCS-Hubbard model in arbitrary dimensions. The model describes a p-wave BCS superconductor with equal spin pairing moving on a bipartite (cubic, square etc.) lattice with on site Hubbard interaction $U$. We show that the model becomes exactly solvable for arbitrary $U$ when the BCS pairing amplitude $Delta$ equals the hopping amplitude $t$. The nature of the solution is described in detail in this paper. The construction of the exact solution is parallel to the exactly solvable Kitaev honeycomb model for $S=1/2$ quantum spins and can be viewed as a generalization of Kitaevs construction to $S=1/2$ interacting lattice fermions. The BCS-Hubbard model discussed in this paper is just an example of a large class of exactly solvable lattice fermion models that can be constructed similarly.