Improved parallel QR method for large least squares problems involving Kronecker products