A protocol is proposed to realize one-step implementation of the N-qubit nonadiabatic holonomic quantum gates with superconducting qubits. The inverse Hamiltonian engineering is applied in designing microwave pulses to drive superconducting qubits. By combining curve fitting, the wave shapes of the designed pulses can be described by simple functions, which are not hard to realize in experiments. To demonstrate the effectiveness of the protocol, a three-qubit holonomic controlled π -phase gate is taken as an example in numerical simulations. The results show that the protocol holds robustness against noise and decoherence. Therefore, the protocol may provide an alternative approach for implementing N-qubit nonadiabatic holonomic quantum gates.The devices for implementing N-qubit nonadiabatic holonomic quantum gates are shown in Figure 1a. As shown in Figure 1a, the devices are composed of N + 2 superconducting qubits (0, 1, 2, . . . , N − 1, N A , N B ) and N + 1 1D coplanar waveguide resonators (CPWR 0 , CPWR 1 , CPWR 2 , . . . ,CPWR N ). The level configuration of superconducting qubit 0 is shown in