We prove that for a one-dimensional infinite lattice, with long-range coupling among sites, the diffusion of an initial delta-like pulse in the bulk, is ballistic at all times. We obtain a closedform expression for the mean square displacement (MSD) as a function of time, and show some cases including finite range coupling, exponentially decreasing coupling and power-law decreasing coupling. For the case of an initial excitation at the edge of the lattice, we find an approximate expression for the MSD that predicts ballistic behavior at long times, in agreement with numerical results.PACS numbers: 63.10.+a, 66.30.-h The physics of discrete systems have been a topic of interest for many years, because it can give rise to completely different phenomenology in comparison with that present in homogeneous continuous systems. In particular, discrete periodic systems are found in many different contexts including condensed matter physics, optics, Bose-Einstein condensates and magnetic metamaterials among others. Under the appropriate approximation, they can all be described by some variant of the discrete Schrödinger (DS) equation [1]. In that way, many of these systems display the phenomenology common to periodic systems such as the presence of a band structure, discrete diffraction, Bloch oscillations, dynamic localization, Zener tunneling, to name a few.Usually, the DS equation is used in the weak-coupling limit, where the interaction among sites includes nearestneighbors only. This is a good approximation in cases where the coupling among sites decays very quickly with distance, like the exponentially-decreasing coupling found in optical waveguide arrays. However, there are cases where it is advisable to go beyond this approximation. An example of that is a split-ring resonator (SRR) array, where the interaction among the basic units is dipolar in nature and therefore, the coupling decreases as the inverse cubic power of the mutual distance.When coupling beyond nearest-neighbors are considered, the number of possible routes of energy exchange increases. The dynamical evolution of excited pulses in finite 1D and 2D lattices with anisotropic couplings and up to second nearest-neighbor couplings, has been explored in [2] by means of the Green function formalism. Experimental observation of the influence of second order coupling in linear and nonlinear optical zig-zag waveguide arrays has been recently carried out [3]. In a different context, a recent work [4], shows that long-range coupling in low dimensional system can induce a phase transition from delocalized to localized modes. Another interesting scenario that can be modeled as a discrete system with long-range couplings is that of complex networks [5], where the distances between nodes are not necessarily physical.In this Letter we carry out an analytical and numerical study on the diffusion of an initially localized pulse propagating in a one-dimensional discrete periodic lattice, in the presence of arbitrary long-range couplings. We focus on two cases of in...