Polycystic Ovary Syndrome (PCOS) is one of the most prevalent endocrine disorders affecting women of reproductive age and is a leading cause of infertility. Ultrasound imaging is widely used for PCOS diagnosis; however, visual assessment of ovarian morphology is highly subjective, time-consuming, and dependent on clinical expertise. Quality differences in ultrasound images, very near to similar visual patterns among PCOS and NOT PCOS images, and noise in the images increase the threat of improper diagnosis. These problems suggest a need for an accurate, automatic, and computer-assisted PCOS diagnostic system. This research aims to create a deep learning-assisted automatic PCOS diagnostic system which can detect and classify the Polycystic Ovary Syndrome from the gray-scale ultrasound ovarian images. In addition to high classification accuracy, the proposed framework incorporates an explicit explainability pipeline that highlights diagnostically relevant ovarian regions, such as follicular distributions and stromal patterns, thereby supporting clinically interpretable decision making. The proposed HAREN framework addresses the limitations of single backbone models, and attention augmented variants, such as vanilla ResNet50 and ResNet50 with hybrid attention by leveraging ensemble learning and residual feature fusion. HAREN combines three architecturally diverse and complementary pretrained CNN backbones (ResNet50, DenseNet121, and EfficientNetB0) to enhance feature diversity. In addition, a novel hybrid attention mechanism combining channel, spatial, and cross-scale attention is introduced to emphasize diagnostically relevant ovarian regions. A residual fusion strategy is employed to preserve discriminative features and stabilize training, and an explicit explainability pipeline is incorporated to support Grad CAM-based visual interpretation. This network first converts the ultrasound grayscale ovarian images to RGB , followed by the extraction of important features applying backbones, which are augmented with attention mechanisms. The network, trained with categorical crossentropy loss, was evaluated using comprehensive performance metrics on 11,784 ultrasound images (6,784 PCOS and 5,000 NOT PCOS). HAREN achieved 99.33% accuracy, 98.96% precision, 98.97% recall, 98.96% F1 score, and an AUC of 99.93%, outperforming conventional models. Overall, it delivers an accurate, reliable, and interpretable solution for automated PCOS detection, demonstrating strong potential for clinical decision support systems