Integrable discretization of recursion operators and unified bilinear forms to soliton hierarchies