Under certain conditions, two samples of fluid at different initial temperatures present a counterintuitive behavior known as the Mpemba effect: it is the hotter system that cools sooner. Here, we show that the Mpemba effect is present in granular fluids, both in uniformly heated and in freely cooling systems. In both cases, the system remains homogeneous, and no phase transition is present. Analytical quantitative predictions are given for how differently the system must be initially prepared to observe the Mpemba effect, the theoretical predictions being confirmed by both molecular dynamics and Monte Carlo simulations. Possible implications of our analysis for other systems are also discussed. DOI: 10.1103/PhysRevLett.119.148001 Let us consider two identical beakers of water, initially at two different temperatures, put in contact with a thermal reservoir at subzero (on the Celsius scale) temperature. While one may intuitively expect that the initially cooler sample would freeze first, it has been observed that this is not always the case [1]. This paradoxical behavior named the Mpemba effect (ME) has been known since antiquity and discussed by philosophers like Aristotle, Roger Bacon, Francis Bacon, and Descartes [2,3]. Nevertheless, physicists only started to analyze it in the second part of the past century, mainly in popular science or education journals .There is no consensus on the underlying physical mechanisms that bring about the ME. Specifically, water evaporation [4,5,9,24], differences in the gas composition of water [11,17,25], natural convection [6,23,26], or the influence of supercooling, either alone [14,27] or combined with other causes [28][29][30][31], have been claimed to have an impact on the ME. Conversely, the own existence of the ME in water has been recently put in question [32]. Notwithstanding, Mpemba-like effects have also been observed in different physical systems, such as carbon nanotube resonators [33] or clathrate hydrates [34].The ME requires the evolution equation for the temperature to involve other variables, which may facilitate or hinder the temperature relaxation rate. The initial values of those additional variables depend on the way the system has been prepared, i.e., "aged," before starting the relaxation process. Typically, aging and memory effects are associated with slowly evolving systems with a complex energy landscape, such as glassy [35][36][37][38][39][40][41][42][43] or dense granular systems [44][45][46]. However, these effects have also been observed in simpler systems, like granular gases [47][48][49][50] or, very recently, crumpled thin sheets and elastic foams [51].In a general physical system, the study of the ME implies finding those additional variables that control the temperature relaxation and determining how different they have to be initially in order to facilitate its emergence. In order to quantify the effect with the tools of nonequilibrium statistical mechanics, a precise definition thereof is mandatory. An option is to look at the relaxat...