Abstract: The discussion of integrability of soliton equation is important for soliton theory. The exact solution of a nonlinear soliton equation is another essential aspect in soliton problem. Among them, the Hirota method and Bäcklund transformation are proved to be the effective approaches to find the exact solutions for the nonlinear soliton equation. In this paper, the supersymmetric form of the variable coefficient Korteweg-de Vries (VCKdV) equation is given. The soliton solutions for the VCKdV equation are derived by Hirota method and Bäcklund transformation. First, how to seek the supersymmetry VCKdV equation was discussed by the direct method. Through variable transformation and bilinear method, the supersymmetry variable coefficient KdV equation can be written in a bilinear form. Soliton solutions for the supersymmetry VCKdV equation are obtained by supersymmetry bilinear derivatives. Then starting from the bilinear form of the Supersymmetry VCKdV equation, the bilinear Bäcklund transformation was obtained. By the commutability of the bilinear Bäcklund transformation, the one soliton solution, two soliton solution, and three soliton solution for the supersymmetry VCKdV equation were given respectively.