we discuss branching ratios for the LFV decays of heavier Higgs bosons, $H/A \to \tau^\pm \ell_i^\mp$, using Eqs. (\[eq:BrhtmH\]) and (\[eq:BrhtmA\]) under the current data of LFV rare tau decays. In the THDM, there are many possible decay modes for $H$ depending on the mass spectrum. In the numerical analysis, we included contributions from one-loop induced $Z\gamma$, $\gamma\gamma$ and $gg$ modes in addition to all the tree level modes. Those for $A$ are one-loop induced modes of $hh$, $hH$, $HH$, $h\gamma$, $H\gamma$, $W^\pm W^\mp$, $ZZ$, $Z\gamma$, $\gamma\gamma$ and $gg$ in addition to all the tree level modes. The branching ratios for $H/A \to \tau^\pm \ell_i^\mp$ are sensitive to the masses of all the Higgs bosons. Here we consider the case of $\sin(\alpha-\beta)=-1$ and $m_H^{}=m_A^{}=m_{H^\pm}^{}$ $(\equiv m_\Phi^{})$. As discussed in Sec. II, the $\rho$ parameter constraint is satisfied for this choice. From the $b \to s \gamma$ results, $m_\Phi^{}$ is taken to be greater than $350$ GeV. As also discussed in Sec. II, $M$ determines the decoupling property of heavier Higgs bosons. Although the branching ratios ${\rm Br}(H/A \to \tau^\pm \ell_i^\mp)$ are insensitive to $M$ in the present parameter set, its value strongly affects the allowed parameter region under the theoretical 