This Letter addresses the physics-consistent optimization of reconfigurable intelligent surfaces (RISs) with mutual coupling (MC) and 1-bit-programmable RIS elements. This combination of constraints is typical of current prototypes but unexplored in theoretical work. First, we present a simple statistical generator for multiport-network-theory (MNT) parameters of rich-scattering, RIS-parametrized channels. We account for reciprocity, passivity, and coherent backscattering; then, we add a simple hyper-parameter to control the MC strength. Second, we benchmark model-agnostic (dictionary search, coordinate descent, genetic algorithm) and model-based (temperature-annealed back-propagation) strategies under varying MC strength, with and without intelligent initialization. Except when MC is negligible, coordinate descent with random initialization offers the best trade-off in performance, runtime, and memory. Our insights can guide wireless practitioners who optimize RIS prototypes and other reconfigurable wave systems.