Laser-induced breakdown spectroscopy (LIBS) of pharmaceutical drugs that contain paracetamol was investigated in air and argon atmospheres. The characteristic neutral and ionic spectral lines of various elements and molecular signatures of CN violet and C₂ Swan band systems were observed. The relative hardness of all drug samples was measured as well. Principal component analysis, a multivariate method, was applied in the data analysis for demarcation purposes of the drug samples. The CN violet and C₂ Swan spectral radiances were investigated for evaluation of a possible correlation of the chemical and molecular structures of the pharmaceuticals. Complementary Raman and Fourier-transform-infrared spectroscopies were used to record the molecular spectra of the drug samples. The application of the above techniques for drug screening are important for the identification and mitigation of drugs that contain additives that may cause adverse side-effects.