We find all self-duality functions of the form D(ξ,η)=∏xd(ξx,ηx)\\documentclass[12pt]{minimal}\t\t\t\t\\usepackage{amsmath}\t\t\t\t\\usepackage{wasysym}\t\t\t\t\\usepackage{amsfonts}\t\t\t\t\\usepackage{amssymb}\t\t\t\t\\usepackage{amsbsy}\t\t\t\t\\usepackage{mathrsfs}\t\t\t\t\\usepackage{upgreek}\t\t\t\t\\setlength{\\oddsidemargin}{-69pt}\t\t\t\t\\begin{document}$$\\begin{aligned} D(\\xi , \\eta )= \\prod _{x} d(\\xi _x, \\eta _x) \\end{aligned}$$\\end{document}for a class of interacting particle systems. We call these duality functions of simple factorized form. The functions we recover are self-duality functions for interacting particle systems such as zero-range processes, symmetric inclusion and exclusion processes, as well as duality and self-duality functions for their continuous counterparts. The approach is based on, firstly, a general relation between factorized duality functions and stationary product measures and, secondly, an intertwining relation provided by generating functions. For the interacting particle systems, these self-duality and duality functions turn out to be generalizations of those previously obtained in Giardinà et al. (J Stat Phys 135:25–55, 2009) and, more recently, in Franceschini and Giardinà (Preprint, arXiv:1701.09115, 2016) . Thus, we discover that only these two families of dualities cover all possible cases. Moreover, the same method discloses all simple factorized self-duality functions for interacting diffusion systems such as the Brownian energy process, where both the process and its dual are in continuous variables.