  The blazar 4C 31.03 recently underwent a major gamma-ray outburst at the
beginning of 2023 after a prolonged quiescent phase. Fermi-LAT reported a daily
average flux of 5x10^-6 phs cm^-2 s^-1, which is about 60 times its average
value. We investigated this extraordinary outbreak through temporal and
multi-wavelength analysis. From the statistical analysis of the gamma-ray
lightcurves using Bayesian blocks, we identified 3 epochs of prominent flares.
The fastest flux decay during this major outburst was observed within 5.5 +/-
0.7 hours. The highest energy of gamma-ray photons found from the source during
the active phase is ~ 82 GeV. Using the transparency of gamma-rays against pair
production and light crossing time argument, we could obtain the minimum jet
Doppler factor as 17 corresponding to the flaring state.The broadband spectral
energy distribution study performed using synchrotron, SSC and EC emission
processes supports the external Compton scattering of IR photons as the likely
mechanism for the gamma-ray emission from the source. The results of this study
suggest the scenario of the emission region in 4C 31.03, being located beyond
the Broad-line region from the central blackhole. Long-term gamma-ray flux
distribution of 4C 31.03 depicts a double log-normal variability, indicating
that two distinct flux states are active in this energy band. The index
distribution also reveals a two distinct variability patterns, suggesting that
the gamma-ray spectrum can be more precisely described by two photon indices.
