Mam następującą rzadką matrycę, która zawieraO(N) elementy

boost::numeric::ublas::compressed_matrix<int> adjacency (N, N);

Mógłbym napisać podwójną pętlę brutalnej siły, aby przejrzeć wszystkie wpisyO(N^2) czas jak poniżej, ale to będzie zbyt wolne.

for(int i=0; i<N; ++i)
   for(int j=0; j<N; ++j)
       std::cout << adjacency(i,j) std::endl;

Jak mogę zapętlić tylko niezerowe wpisy wO(N) czas? Dla każdego niezerowego elementu chciałbym mieć dostęp do jego wartości i indeksówi,j.