On large scales, runoff changes are complex due to the regional differences in basin characteristics and interactions among multiple influencing factors. Accurately and quantitatively decomposing the impacts of climate change and human activities on runoff changes in large basins with significant spatial heterogeneity remains a challenge. In this paper, the fuzzy C-means (FCM) algorithm was used to identify the regional differences in hydrological processes of the Yangtze River Basin (YZRB) based on multiple factors. Then, regional Budyko equations with synthetic parameters were constructed for runoff decomposition using the hierarchical Bayesian model, considering the spatial correlations of multiple factors among basins. Results showed that 44 subbasins in the YZRB can be clustered into five clusters based on 10 factors including precipitation, temperature, precipitation seasonality, vegetation coverage, saturated hydraulic conductivity, elevation, slope, topographic wetness index, irrigation water consumption, and proportion of built-up area per unit area. From Cluster 1 to 5, P, T, IRR, and B decreased from about 1400 mm, 17 °C, 57.9 mm, and 1.59 % to about 870 mm, 12 °C, 15.6 mm, and 0.40 %, respectively, while ELEV and SLP increased from about 300 m and 3° to about 3500 m and 10°, respectively. Hierarchical Bayesian models can improve the performance of Budyko equations by better-simulating n, compared with linear regression models. Their R2, RB, and RMSE for the whole YZRB were as follows: 0.9104 versus 0.7591, 4.33 % versus 8.21 %, and 0.06 versus 0.11, respectively. Runoff decomposition from 1982 to 2018 showed that precipitation led to more than 34 % change in runoff in Cluster 2, while the influence of potential evapotranspiration on runoff changes was less than 3.00 % in each cluster. Rising temperature caused runoff to decrease by more than 7.90 % in Cluster 1, 3, and 4. In addition, precipitation seasonality, vegetation, and irrigation water consumption had different influences in different clusters. This study improves the ability of the Budyko equation to quantitatively decompose the impacts of multiple factors on regional runoff changes and provides more reasonable suggestions for water resource management.