Пересечение треугольников

#include <iostream>
#include <cmath>
 
using namespace std;
 
struct D_Point{
    double x;
    double y;
};
 
double s_n(D_Point a[], int n){
    //ïëîùàäü n - óãîëüíèêà
    double sum = 0;
    for (int i = 0; i < n; ++i){
        double det = a[i % n].x * a[(i + 1) % n].y 
                      - a[(i + 1) % n].x * a[i % n].y;
        sum += det;
    }
    return fabs(sum / 2);
}
 
int main(){
 
    const double EPS = 1e-10;
    
    D_Point a[3], b[3];
    
    a[0].x = 0; a[0].y = 5;
    a[1].x = 1; a[1].y = 0;
    a[2].x = 3; a[2].y = 2;
    
    b[0].x = 4; b[0].y = 4;
    b[1].x = 5; b[1].y = 1;
    b[2].x = 1; b[2].y = 0;
    
    double sa = s_n(a, 3);//ïëîùàäü ïåðâîãî òðåóãîëüíèêà
    
    int flag = 0;
    for (int i = 0; i < 3; ++i){//òî÷êà âòîðîãî êîòîðóþ áåðåì
        
        double sum = 0;
        for (int j = 0; j < 3; ++j){//òî÷êà ïåðâîãî òð-êà êîòîðóþ íå áåðåì
            
            D_Point c[3];
            for (int k = 0; k < 3; ++k)
                if (k == j)
                    c[k] = b[i];
                else
                    c[k] = a[k];
            
            double sc = s_n(c, 3);//ïëîùàäü íîâîãî òðåóãîëüíèêà
            
            sum += sc;
        }
        
        if (fabs(sa - sum) < EPS){//çíà÷èò òî÷êà âíóòðè òðåóãîëüíèêà
            flag = 1;
            break;
        }
    }
    
    cout << flag << endl;
}