A bstract It has been shown recently that the anomalies observed in B ¯ → D * τ ν ¯ τ and B ¯ → K ¯ ℓ + ℓ − decays could be resolved with just one scalar leptoquark. Fitting to the current data on R ( D (∗)) along with acceptable q 2 distributions in B ¯ → D ∗ τ ν ¯ τ decays, four best-fit solutions for the operator coefficients have been found. In this paper, we explore the possibilities of how to discriminate these four solutions. Firstly, we find that two of them are already excluded by the decay B c − → τ − ν ¯ τ , because the predicted decay widths have already overshot the total width Γ B c . It is then found that the remaining two solutions result in two effective Hamiltonians governing b → c τ ν ¯ τ transition, which differ by a sign and enhance the absolute value of the coefficient of c ¯ L γ μ b L τ ¯ L γ μ ν τ L operator by about 12%. However, they give nearly the same predictions as in the SM for the D ∗ and τ longitudinal polarizations as well as the lepton forward-backward asymmetries in B ¯ → D ∗ τ ν ¯ τ decays. For the other observables like ℬ B c − → τ − ν ¯ τ , ℬ B c − → γ τ − ν ¯ τ , R D ∗ q 2 , d ℬ B ¯ → D ∗ τ ν ¯ τ / d q 2 and ℬ B c − → X 3 τ ν ¯ τ , on the other hand, the two solutions give sizable enhancements relative to the SM predictions. With measurement of B c − → τ − ν ¯ τ at LHCb and refined measurements of observables in B ¯ → D ∗ τ ν ¯ τ at both LHCb and Belle-II, such a specific NP scenario could be further deciphered.