We investigate the long-wavelength evolution of linear perturbations in a homogeneous and anisotropic background with a scalar field coupled to a vector field. Using the spatial gradient expansion in the uniform-$\mathcal{N}$ gauge in which the number of $e$-folds is unperturbed, we derive the complete set of superhorizon solutions and establish their correspondence with infinitesimal variations of the homogeneous anisotropic background. This extends the separate-universe picture, previously known for isotropic FLRW cosmology, to anisotropic spacetimes despite the mixing of scalar, vector, and tensor perturbations induced by broken rotational symmetry. We show that the long-wavelength equations form a self-consistent system and identify a conserved quantity that generalizes the conserved Wronskian of isotropic cosmology. Unlike the isotropic case, the superhorizon modes sourcing the curvature perturbation are governed by three independent conserved channels associated with the scalar field, the background shear, and the gauge-field tilt, together with an additional dynamical shear contribution originating from the anisotropic geometry. This reveals that the evolution of curvature perturbations around anisotropic background is intrinsically richer than in isotropic multi-field models. Our formulation provides a practical prescription for computing the final curvature perturbation directly from horizon-crossing fluctuations, thereby establishing the anisotropic generalization of the $δN$ formalism. We further derive an explicit relation between curvature perturbations and primordial gravitational waves, demonstrating how anisotropic expansion couples scalar and tensor sectors on superhorizon scales. Our framework provides a practical basis for predicting statistical anisotropies in primordial scalar and tensor perturbations.