Scalable Link Prediction in Dynamic Networks via Non-Negative Matrix Factorization