#include <iostream>
#include <deque>

#include <boost/geometry.hpp>
#include <boost/geometry/geometries/point_xy.hpp>
#include <boost/geometry/geometries/polygon.hpp>
#include <boost/geometry/geometries/linestring.hpp>

#include <boost/foreach.hpp>


int main()
{
    typedef boost::geometry::model::d2::point_xy<long> point; 
    typedef boost::geometry::model::polygon<point> polygon;
    typedef boost::geometry::model::linestring<point> linestring;

    polygon green;
    linestring blue;

    boost::geometry::read_wkt(
        "POLYGON((0 0, 0 2, 2 2, 2 0, 1 0, 0 0))", green);

    boost::geometry::read_wkt(
        "LINESTRING(1 -1,1 0, 1 1)", blue);

    boost::geometry::correct(green);
    boost::geometry::correct(blue);

    std::deque<linestring> output;
    boost::geometry::intersection(blue, green, output);

    int i = 0;
    std::cout << "green && blue:" << std::endl;
    BOOST_FOREACH(linestring const& p, output)
    {
        std::cout << i++ << ": " << boost::geometry::dsv(p) << std::endl;
    }


    return 0;
}