The increasing use of nonlinear loads in power systems introduces voltage and current components at non-fundamental frequencies, leading to harmonic distortion, which negatively impacts electrical and electronic devices. A common mitigation strategy involves identifying harmonic sources and installing filters nearby. However, due to the high cost of power quality (PQ) meters, comprehensive harmonic level monitoring across the entire power system is impractical. To address this, various methodologies for Harmonic State Estimation (HSE) have been developed, which estimate distortion levels on unmonitored system buses using data from a minimal set of monitored ones. Many HSE techniques rely on optimization algorithms with numerous tuning parameters, complicating their application. This paper proposes a novel methodology for fundamental frequency power flow and harmonic state estimation using the Jaya algorithm, which is characterized by fewer tuning parameters for easier adjustment. It also introduces a strategy to determine the minimal number of buses that need monitoring to achieve system observability. The methodology is validated on the IEEE-14 and IEEE-30 bus systems, demonstrating its effectiveness. The results of the proposed methodology are compared with those obtained using Evolutionary Strategies (ESs), highlighting its enhanced accuracy and computational efficiency.