Function bfs_with_parents 

pub fn bfs_with_parents(
    adj: &CsrMatrix,
    sources: &[NodeId],
) -> Result<(Vec<Distance>, Vec<NodeId>)>

Multi-source BFS returning parent pointers for path reconstruction.