P-variation calculation and statistical inference for the discretely observed differential equation driven by fractional Brownian motion