Michal P Krompiec
1) Look up PCM in the manual. Moreover, there are quite a few recent (and less recent) reviews on this topic.
2) You normally do not include the solvent molecule in the calculation, the solvent is modelled as a "polarizable continuum" surrounding your molecule (solute). PCM method(s) are relatively easy to use (black-box). It is possible to include some solvent molecules around your solute, this is called explicit solvation and is considerably more difficult to do. You probably want to optimize the geometry of your solute in the simulated solvent environment, though.
Best wishes and good luck,