How to forecast and fit the optimal model for multiple time series?

Question

I want to do batch forecasting among multiple series, for example, if I want to forecast time series with IDs that end with 1(1,11,21,31...), how can I do that?