Earthquake-induced landslides pose significant threats to human life and property. Assessing their hazard reveals the probability of regional landslide occurrence under seismic action, providing a basis for pre-earthquake planning and post-earthquake rescue operations. This study integrates landslide data from 26 historical earthquakes through efficient data preprocessing. A Fault Direction Effect (FDE) factor is introduced, quantitatively characterizing the angle between natural slope aspect and fault rupture direction, establishing a seismic geological effect evaluation index with defined physical significance. We constructed a multi-factor database comprising 18 disaster-inducing factors and conducted statistical analyses of the spatial distribution of earthquake-induced landslides. Based on the Newmark model, we developed eighteen coupled Newmark and machine learning models suitable for shallow translational failures, hereafter referred to as the Newmarki-X dual-drive models. Using metrics such as AUC, precision, recall, F1 score, and Kappa coefficient, supplemented by leave-one-out cross-validation, and taking the 2008 Wenchuan, 2015 Nepal, and 2010 Haiti earthquake data, the best performing model (improved Newmark coupled with Backpropagation Neural Network(BPNN) incorporating peak ground acceleration (PGA)) showed significant performance improvements: pre-earthquake prediction AUC reached 0.88 and post-earthquake inversion AUC 0.94, outperforming traditional single-method models in both accuracy and stability. By employing big data analysis, machine learning, and model integration techniques, this study develops a multi-technique earthquake-induced landslide hazard evaluation system, offering quantitative support for regional disaster risk management.