Higher derivative corrections of f(R) gravity with varying equation of state in the case of variable G and Λ