пжста найдите ошибку в задаче на выпуклую оболочку

#include<iostream> 
#include<vector> 
#include<algorithm> 
#include<string> 
#include<set> 
#include<map> 
#include<iterator> 
#include<iomanip> 
#include<ctime> 
#include<deque>
#include<stack>
#include<queue>
#include<math.h> 
#include<cmath> 
using namespace std;
#pragma warning (disable:4996); 
#define lint long long 
#define endl '\n'
const long long INF = 1e17;
const long double eps = 0.00000000000001;
 
struct vec {
    lint x, y;
    vec(lint a, lint b) : x(a), y(b) {}
    long double len() {
        return sqrt(x * x + y * y);
    }
};
vector<vec>g;
vector<vec>st;
 
vec operator + (vec &a, vec &b) {
    return vec(a.x + b.x, a.y + b.y);
}
vec operator - (vec &a, vec &b) {
    return vec(a.x - b.x, a.y - b.y);
}
vec operator * (vec &a, long double k) {
    return vec(a.x * k, a.y * k);
}
lint operator % (vec &a, vec &b) {
    return a.x * b.y - b.x * a.y;
}
lint operator ^ (vec &a, vec &b) {
    return a.x * b.x + a.y * b.y;
}
vec vect(vec &a, vec &b) {
    return a - b;
}
lint vect1(vec &a, vec &b, vec&c, vec&d) {
    vec v1 = vect(a, b);
    vec v2 = vect(c, d);
    return v1 % v2;
}
bool comp(vec a, vec b) {
    vec a1(a.x - g[0].x, a.y - g[0].y);
    vec b1(b.x - g[0].x, b.y - g[0].y);
    if (a1%b1 > 0)return true;
    if (a1%b1 == 0)return a1.len() < b1.len();
 
    return 0;
}
 
int main() {
    ios_base::sync_with_stdio(false);
    cin.tie(NULL);
    cout.tie(NULL);
#ifdef _DEBUG
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
#endif
    long double a1, b1, c1, n, m, k, a2, b2, c2, a, b, c, d;
    cin >> n;
    lint xx, yy;
    map<pair<lint, lint>, lint>mp;
    for (int i = 0; i < n; ++i) {
        cin >> xx >> yy;
        vec vec1(xx, yy);
        g.push_back(vec1);
        mp[make_pair(xx, yy)] = i + 1;
    }
    lint miny = INF, minx = INF, index = 0;
    for (int i = 0; i < n; ++i) {
        if (g[i].x <= minx) {
            if (g[i].y < miny) miny = g[i].y;
            index = i;
            minx = g[i].x;
        }
    }
    swap(g[0], g[index]);
    sort(g.begin() + 1, g.end(), comp);
    st.push_back(g[0]);
    lint x11 = st[st.size() - 1].x, y11 = st[st.size() - 1].y;
    st.push_back(g[1]);
    st.push_back(g[2]);
    for (int i = 3; i < n; ++i) {
        if (g[i].x == st[st.size() - 1].x&&g[i].y == st[st.size() - 1].y)continue;
        while (st.size() > 1 && vect1(st[st.size() - 1], st[st.size() - 2], g[i], st[st.size() - 1]) <= 0ll)
        st.pop_back();
        st.push_back(g[i]);
    }
    cout << st.size() << endl;
    for (int i = 0; i < st.size(); ++i)
        cout << mp[make_pair(st[i].x, st[i].y)] << " ";
    cout << endl;
    long double sum = 0;
    vec vec1(x11, y11);
    st.push_back(vec1);
    vector<vec>st2 = st;
    while (!st.empty()) {
        lint xx = st[st.size() - 1].x;
        lint yy = st[st.size() - 1].y;
        st.pop_back();
        if (!st.empty())
            sum += sqrt((xx - st[st.size() - 1].x)*(xx - st[st.size() - 1].x) + (yy - st[st.size() - 1].y)*(yy - st[st.size() - 1].y));
    }
    cout << fixed << setprecision(10) << sum << endl;
 
    lint sum1 = 0;
    for (int i = 0; i < st2.size() - 1; ++i) 
        sum1 += st2[i] % st2[i + 1];
    sum1 += st2[st2.size() - 1] % st2[0];
    if (sum1 % 2 == 0)cout << sum1 / 2;
    else cout << sum1 / 2 << ".5";
    return 0;
}