How do we discretely represent this integral in order to implement it?

Also, is there an efficient analytical method of rendering that computes this integral or something similar instead of sampling a bunch of random directions? If such a method existed, it could be much less noisy.