This research aims to elucidate the organizational patterns of inter-regional economic interdependence to enhance our comprehension of the national economy's structure at a regional scale. Employing a multilayer network model, this study represents economic interdependence among Indonesian regions, utilizing the Inter-Regional Input-Output (IRIO) table. Through the application of various metrics, such as degree and strength distribution, assortativity coefficient, and rich club coefficient, to the multilayer IRIO network, we uncover the organizational patterns of economic exchanges between provinces and economic sectors within Indonesia. Our findings demonstrate that a multilayer network approach reveals the heterogeneous and complex structure of the national economy at the regional level. By analyzing the assortativity pattern and global rich-club coefficient, we illustrate that the IRIO network exhibits a hierarchical organization, where significant provincial-sector nodes are interconnected and form dense rich clubs, extending from a few structural cores to peripheral regions.