Modeling of HeN+ clusters I. Ab initio and analytical potential energy surfaces for the ground state and two lowest excited states of He3+