The numerical integration of the differential equations describing chemical kinetics consumes the majority of computational time in combustion simulations that involve direct coupling of chemistry and flow, such as transported probability density function (PDF) methods, direct numerical simulation (DNS), conditional moment closure (CMC), unsteady flamelet, multiple mapping closure (MMC), thickened flame model, linear eddy model (LEM), partially stirred reactor (PaSR) as in OpenFOAM and laminar flame computation. This step can be accelerated by tabulation, and artificial neural networks (ANNs) have recently emerged as a powerful technique in this domain. To be applicable to a wide family of problems, an ANN tabulation approach must be based on data generated by an abstract process, rather than from the turbulent flame to be simulated. In the present work, the hybrid flamelet/random data and multiple multilayer perceptrons (HFRD-MMLP) method (Ding et al., Combust. Flame 231, 111493, 2021) for non-premixed flames is taken as a basis to develop a thermochemistry tabulation method for premixed flames. In the spirit of maintaining an essentially random data set that still originates in meaningful composition states, a set of one-dimensional premixed flame simulations is employed to generate data that are used as starting points for a random data generation process and subsequently discarded. The approach is applied to large eddy simulations (LES) of the Cambridge/Sandia swirl burner in configurations five and six, with the transported PDF method employed to provide closure for the filtered reaction source terms and the stochastic fields method used for numerical solution. Very good agreement in both major and minor species is observed between the LES-PDF simulations using direct integration of the reaction source term and the ones with the ANNs. Furthermore, the average time taken for reaction source term computations is reduced by fourteen times, while memory requirements constitute only 1.4 MB.