Bilayer transition metal dichalcogenides (TMDs) belong to a class of materials with two unique features, the coupled spin-valley-layer degrees of freedom and the crystal structure that is globally centrosymmetric but locally non-centrosymmetric. In this work, we will show that the combination of these two features can lead to a rich phase diagram for unconventional superconductivity, including intra-layer and inter-layer singlet pairings and inter-layer triplet pairings, in bilayer superconducting TMDs. In particular, we predict that the inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov state can exist in bilayer TMDs under an in-plane magnetic field. We also discuss the experimental relevance of our results and possible experimental signatures. 74.25.Dw Introduction.-Unconventional superconductivity [1][2][3], which is beyond the simple s-wave spin-singlet superconductivity in the Bardeen-Cooper-Schrieffer theory, can emerge in two dimensional (2D) systems, such as surfaces [4][5][6] or interfaces [7], superconducting heterostructures [8] and 2D or quasi-2D superconducting materials [9][10][11][12][13][14]. Recently, it was demonstrated that "Ising" superconductivity can exist in monolayer transition metal dichalcogenides (TMDs), such as MoS 2 [11,13] and NbSe 2 [12], based on experimental observation that in-plane upper critical field H c2, is far beyond the paramagnetic limit. The space symmetry group of the monolayer TMD is the D 3h group without inversion symmetry. Thus, the monolayer superconducting TMDs belong to the so-called non-centrosymmetric superconductors (SCs) [3], for which spin-up and spin-down Fermi surfaces are split by strong spin-orbit coupling (SOC), leading to a mixing of spin singlet and triplet pairings [15,16]. The existence of triplet components can enhance H c2, in non-centrosymmetric SCs [17]. In monolayer TMDs, Ising SOC fixes spin axis along the out-of-plane direction and greatly reduces the Zeeman effect of in-plane magnetic fields, thus explaining the experimental observations of high H c2, . A high H c2, was also observed in bilayer TMDs (e.g. NbSe 2 ) [12]. The crystal structure of bilayer TMDs is described by the symmetry group D 3d with inversion symmetry and the corresponding Fermi surfaces are spin degenerate. This experimental result motivates us to study the difference between bilayer superconducting TMDs and conventional SCs.We first illustrate the difference from symmetry aspect. Although inversion symmetry exists in bilayer TMDs, the inversion center should be chosen at the center between two layers, labeled by "P" in Fig. 1a. As a result, bilayer TMDs belong to a class of materials which are globally centrosymmetric, but locally non-centrosymmetric (for each layer). The absence of local inversion symmetry can lead to the "hidden" spin polarization [18,19], the spinlayer locking [20,21] and other exotic physical phenomena [22]. The superconductivity for these materials has been studied in the CeCoIn 5 /YbCoIn 5 hybrid system [10,23], SrPtAs [23][24][25][26] and other ...