<< NormalOrder.m

a = Op["a"]

Format[Op[n_][x_, y_]] = (Subscript[n, x]^y)

a[3, 4]