In this paper, the ( \(3+1\) )-dimensional Sawada–Kotera (SK) equation is introduced via the Fokas method. Through direct computation, the Nth-order Wronskian determinant and the Nth-order Pfaffian form are rigorously verified to satisfy the Hirota bilinear equation associated with the Hirota bilinear method. Consequently, we are able to construct both Wronskian and Pfaffian solutions for the ( \(3+1\) )-dimensional SK equation. As illustrative examples, we explicitly discuss the Wronskian solutions for the cases \(N = 1, 2, 3\) and the Pfaffian solutions for \(N = 2, 4\) , along with a graphical analysis of the dynamical behaviours of these solutions. In the final section, starting from another bilinear form of the model, we derive a bilinear Bäcklund transformation (BT) involving three arbitrary constants. By applying a gauge transformation to this BT and selecting appropriate parameters, we obtain a new bilinear BT. Using the perturbation method on the new BT, we iteratively derive multi-soliton solutions, explicitly presenting the one-soliton, two-soliton and three-soliton solutions as illustrative examples.