<p>We propose a bootstrap program for the <i>form factor squared</i> with operator tr(<i>ϕ</i>²) in maximally supersymmetric Yang-Mills theory in the planar limit, which plays a central role for perturbative calculations of important physical observables such as energy correlators. The tree-level <i>N</i> -point form factor (FF) squared can be obtained by cutting <i>N</i> propagators of a collection of two-point “master diagrams” at (<i>N</i> – 1) loops: for <i>N</i> = 3, 4, 5, 6 there are merely 1, 2, 4, 13 topologies of such diagrams respectively, and their numerators are strongly constrained by power-counting (including “no triangle” property) and other constraints such as the “rung rule”. Moreover, these two-point diagrams provide a “unification” of FF squared at different numbers of loops and legs, which is similar to extracting (planar) amplitude squared from vacuum master diagrams (dual to <i>f</i> -graphs): by cutting 2 ≤ <i>n</i> &lt; <i>N</i> propagators, one can also extract the planar integrand of <i>n</i>-point FF squared at (<i>N</i> – <i>n</i>) loops, thus our results automatically include integrands of 2-point (Sudakov) FF up to four loops (where the squaring is trivial), 3-point FF squared up to three loops, and so on. Our ansatz is completely fixed using soft limits of (tree and loop) FF squared and the multi-collinear limit which reduces it to the splitting function, without any other inputs such as unitarity cuts. This method opens up the exciting possibility of a <i>graphical bootstrap</i> for FF squared for higher <i>N</i> (which contains e.g. planar Sudakov FF to <i>N</i> – 2 loops) similar to that for the amplitude squared via <i>f</i> -graphs. We also comment on applications to the computation of leading order energy correlators where new structures are expected after performing phase-space integrations.</p>