C++ 如何打印Delaunay图的边？

以下是：，我现在有：

#include <iostream>
#include <fstream>
#include <cassert>

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Segment_Delaunay_graph_2.h>
#include <CGAL/Segment_Delaunay_graph_adaptation_policies_2.h>
#include <CGAL/Segment_Delaunay_graph_traits_2.h>

typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Segment_Delaunay_graph_traits_2<K> Gt;
typedef CGAL::Segment_Delaunay_graph_2<Gt> DT;

int main() {

    std::ifstream ifs("data.cin");
    assert( ifs );
    DT         vd;
    DT::Site_2 site;
    // read the sites from the stream and insert them in the diagram
    while ( ifs >> site ) { vd.insert( site ); }
    ifs.close();
    // validate the diagram
    assert( vd.is_valid(true, 1) );
    std::cout << std::endl << std::endl;
    // Iterate over edges
    DT::Finite_edges_iterator eit = vd.finite_edges_begin();
    for (int k = 1; eit != vd.finite_edges_end(); ++eit, ++k) {
        DT::Edge e = *eit;
        //std::cout << e << std::endl;
    }
}

---

可以按如下方式打印形成边缘的站点：

for (int k = 1; eit != vd.finite_edges_end(); ++eit, ++k) {
    DT::Edge e = *eit;
    std::cout << "k = " << k << std::endl;
    if ( vd.is_infinite(e.first->vertex( vd.ccw(e.second) )) )
        std::cout << "infinite\n";
    else
        std::cout << e.first->vertex( vd.ccw(e.second) )->site() << std::endl;
    if ( vd.is_infinite(e.first->vertex( vd.cw(e.second) )) )
        std::cout << "infinite\n";
    else
        std::cout << e.first->vertex( vd.cw(e.second) )->site() << std::endl;
    if ( vd.is_infinite(e.first->vertex( e.second )) )
        std::cout << "infinite\n";
    else
        std::cout << e.first->vertex( e.second )->site() << std::endl;
    if ( vd.is_infinite(vd.tds().mirror_vertex(e.first, e.second) ) )
        std::cout << "infinite\n";
    else
        std::cout << vd.tds().mirror_vertex(e.first, e.second)->site() << std::endl;
}

for（int k=1；eit！=vd.finite_edges_end（）；++eit，++k）{
DT:：边e=*eit；
标准：：cout
for (int k = 1; eit != vd.finite_edges_end(); ++eit, ++k) {
    DT::Edge e = *eit;
    std::cout << "k = " << k << std::endl;
    if ( vd.is_infinite(e.first->vertex( vd.ccw(e.second) )) )
        std::cout << "infinite\n";
    else
        std::cout << e.first->vertex( vd.ccw(e.second) )->site() << std::endl;
    if ( vd.is_infinite(e.first->vertex( vd.cw(e.second) )) )
        std::cout << "infinite\n";
    else
        std::cout << e.first->vertex( vd.cw(e.second) )->site() << std::endl;
    if ( vd.is_infinite(e.first->vertex( e.second )) )
        std::cout << "infinite\n";
    else
        std::cout << e.first->vertex( e.second )->site() << std::endl;
    if ( vd.is_infinite(vd.tds().mirror_vertex(e.first, e.second) ) )
        std::cout << "infinite\n";
    else
        std::cout << vd.tds().mirror_vertex(e.first, e.second)->site() << std::endl;
}