I thought we don't do hmwk questions on this board... and this is basic e&m to boot.
Anyway, I'll be nice. Do you know how to calculate the electric potential from 1 point charge? (A point in 3 dimensions that is.) You do, since it's on your formula sheet: V(r) = k*q/r in spherical coordinates. Where does this come from? The electric potential satisfies poisson's equation: laplacian(V) = rho/epsilon_0, where rho is the charge density. For the case of point charges, rho is linear combination of dirac-delta function. If you've had some ODEs, try to solve: 1/r^2 * diff(r^2 * diff(V(r),r),r) = delta(r)/epsilon_0 (the left-hand side is the nonzero part of the laplacian in spherical coords); the answer is the same V(r) i wrote earlier.
The thing to understand here is that Poisson's Eqn is linear. So solutions are additive. Thus if I have 3 point charges, I don't have to solve another differential equation to find the potential. I can just add up 3 independent solutions for 1 point charge! (linear superposition might ring a bell.) At some level, you knew this b/c you had to exploit linearity to solve the other parts of your problem.
The steps needed to find the electric potential at the origin should be obvious from here.
Also I should point out that this question is actually extremely vague. They don't specify the dimensionality of the problem. A point in 2D does not have the same potential as a point in 3D (a 2D "point" being equivalent to a line of charge in 3D). Basic e&m classes omitting this detail have always irked me.
To see this, find the potential for a point in 2D: solve 1/r*diff(r*diff(V(r),r),r) = delta(r)/epsilon_0 (left-hand side is the nonzero part of the laplacian in cylindrical or polar coords).
Edit: in response to your last comment... the potential is a scalar quantity. It is valid when curl(E) = 0, i.e. E is irrotational. Then you can write E = -gradient(V). And using Gauss's Law, we can write divergence(E) = -laplacian(V) = rho/epsilon_0.
So if you take my math on faith, then it should be obvious that E is a directional while V is scalar, since E = -gradient(V).
But why does that make sense? Well, the electric potential comes from: V = -integral(E dot dL), where E and dL are vector quantities. It's the energy needed to move charge through an E-field starting at point A, going to point B, and following some path. E is the electric field. dL is an incremental length. At the most general level, the integral there depends on the path taken through space.
But that should ring some alarm bells for you. I'm sure you learned in class that E-fields are conservative, path-independent, etc. What happened? Since E is irrotational (and the domain is "simply-connected"), the integral defining V is path-independent; E is conservative.
So if E were NOT conservative, then you could not define a single "global" potential! In fact the potential would then depend on the path you took. Luckily for us, E is conservative, and V is path-independent--no direction needed.
NOTE: E is only conservative under some pretty strict assumptions. But these are the assumptions you work with in intro e&m. Specifically, if you allow magnetic fields (i.e. move the charges "quickly"), then curl(E) != 0 and you need to appeal to something much more complex than Poisson's equation.