#include<bits/stdc++.h>
using namespace std;
#define FOR(i,s,t) for(int i=(s),_t=(t); i<=_t; ++i)
typedef double db;
typedef long long ll;
struct Point {
    ll x,y;
    Point operator -(const Point &s) const{
        return {x-s.x,y-s.y};
    }
    db Len() const{
        return hypot(x,y);
    }
};
ll det(Point a,Point b) {
    return a.x*b.y-a.y*b.x;
}
ll dot(Point a,Point b) {
    return a.x*b.x+a.y*b.y;
}
const int N=405;
Point a[N];
int w[N];
vector<db> X;
struct node {
    db ans,lmx,rmx,sum;
} T[N<<2];
#define ls p<<1
#define rs p<<1|1
#define lson L,mid,ls
#define rson mid+1,R,rs
void build(int L,int R,int p) {
    T[p]={0,0,0,-(X[R]-X[L])*2};
    if(L==R) return ;
    int mid=L+R>>1;
    build(lson);
    build(rson);
}
void Insert(int L,int R,int p,int x,int v) {
    if(L==R) {
        T[p].ans+=v;
        T[p].lmx+=v;
        T[p].rmx+=v;
        T[p].sum+=v;
        return ;
    }
    int mid=L+R>>1;
    if(x<=mid) Insert(lson,x,v);
    else Insert(rson,x,v);
    db d=(X[mid+1]-X[mid])*2;
    T[p].ans=max({T[ls].ans,T[rs].ans,T[ls].rmx+T[rs].lmx-d});
    T[p].lmx=max(T[ls].lmx,T[ls].sum+T[rs].lmx-d);
    T[p].rmx=max(T[rs].rmx,T[rs].sum+T[ls].rmx-d);
    T[p].sum=T[ls].sum+T[rs].sum-d;
}
int main() {
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    cout << fixed << setprecision(12);
    int n;
    cin >> n;
    // n=400;
    db res=0;
    FOR(i,1,n) {
        int x,y;
        cin >> x >> y >> w[i];
        // x=rand()%10000,y=rand()%10000,w[i]=rand()%10000;
        a[i]={x,y};
        res=max(res,(db)w[i]);
    }
    FOR(i,1,n) FOR(j,1,n) if(i!=j) {
        vector<tuple<db,db,int>> vec;
        X.clear();
        FOR(k,1,n) if(det(a[j]-a[i],a[k]-a[i])>=0) {
            db y=det(a[j]-a[i],a[k]-a[i])/(a[j]-a[i]).Len();
            db x=dot(a[j]-a[i],a[k]-a[i])/(a[j]-a[i]).Len();
            X.push_back(x);
            vec.push_back({y,x,w[k]});
        }
        sort(X.begin(),X.end());
        X.erase(unique(X.begin(),X.end()),X.end());
        sort(vec.begin(),vec.end());
        build(0,X.size()-1,1);
        for(auto [y,x,v]: vec) {
            int p=lower_bound(X.begin(),X.end(),x)-X.begin();
            // cerr << i << ' ' << j << ' ' << v << endl;
            Insert(0,X.size()-1,1,p,v);
            res=max(res,T[1].ans-y*2);
        }
    }
    cout << res << '\n';
    return 0;
}