#define _CRT_SECURE_NO_WARNINGS
#define _USE_MATH_DEFINES
#include <iostream>
using namespace std;
#include <vector>
#include <set>
#include <map>
#include <cstdio>
#include <cstring>
#include <queue>
#include <cstdlib>
#include <algorithm>
#include <list>
#include <string>
#include <cmath>
#include <bitset>
#include <stack>
#include <math.h>
#include <functional>
#define ull unsigned long long
#define ll long long
#define pll pair<ll ,ll>
#define ft first
#define sd second
#define cy cout<<"YES"<<endl
#define cn cout<<"NO"<<endl
#define forn(i,n) for(ll (i)=0;(i)<(n);(i)++)
#define forne(i,n) for(ll (i)=1;(i)<=(n);(i)++)
#define rforn(i,n) for(ll (i)=(n-1);i>=0;i--)
#define rforne(i,n) for(ll (i)=n;i>=1;i--)
#define vv vector
#define inf32 INT32_MAX
#define inf64 INT64_MAX
#define LD long double
const int N = 4e5 + 100;
const int mod = 998244353;
ll T;
struct Point {
LD x, y;
Point(LD x, LD y) { this->x = x; this->y = y; }
Point() {}
#define line Point
Point operator+(Point b) { return Point(this->x + b.x, this->y + b.y); }
Point operator-(Point b) { return Point(this->x - b.x, this->y - b.y); }
Point operator*(LD t) { return Point(this->x * t, this->y * t); }
LD operator*(Point b) { return this->x * b.x + this->y * b.y; }
bool operator==(Point b) { return this->x == b.x && this->y == b.y; }
};
class Geomtry_tool {
public:
static LD len(line a) { return sqrt(a.x * a.x + a.y * a.y); }
LD dis(Point a, Point b) { return sqrt((a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y)); }
static LD angle(line a, line b) { return acos(a * b / len(a) / len(b)); }
static LD cross(line a, line b) { return a.x * b.y - a.y * b.x; }
static int sgn(LD x) { LD eps = 1e-10; return abs(x) <= eps ? 0 : x > 0 ? 1 : -1; }
static LD cross(Point a, Point b, Point c) { return cross(b - a, c - a); }
static bool l_S_inter(Point a, Point b, Point c, Point d) {
return sgn(cross(a, b, c) * cross(a, b, d)) <= 0;
}
static bool S_S_inter(Point a, Point b, Point c, Point d) {
return l_S_inter(a, b, c, d) && l_S_inter(c, d, a, b);
}
static Point getNode(Point a, line u, Point b, line v) {
LD t = cross((a - b), v) / cross(v, u);
return a + u * t;
}
static Point rotate(Point a, double b) {
return { a.x * cos(b) - a.y * sin(b), a.x * sin(b) + a.y * cos(b) };
}
};
class Andrew {
Geomtry_tool fn;
public:
int top = 0;
struct cmp {
bool operator()(const Point& a, const Point& b) {
return a.x == b.x ? a.y < b.y : a.x < b.x;
}
};
LD Andrew_lenth(vv<Point>& p, vv<Point>& s) {
int n = p.size();
sort(p.begin(), p.end(), cmp());
s.resize(p.size() + 2);
forn(i, n) {
while (top > 1 && fn.sgn(fn.cross(s[top - 1], s[top], p[i]) <= 0))
top--;
s[++top] = p[i];
}
int t = top;
for (int i = n - 2; i >= 0; i--) {
while (top > t && fn.sgn(fn.cross(s[top - 1], s[top], p[i])) <= 0)
top--;
s[++top] = p[i];
}
LD res = 0;
for (int i = 1; i < top; i++)
res += fn.dis(s[i], s[i + 1]);
return res;
}
LD rotating_calipers(vv<Point>& s) {
int n = top - 1;
LD res = 0;
for (int i = 1, j = 2; i <= n; i++) {
while (fn.cross(s[i], s[i + 1], s[j])
< fn.cross(s[i], s[i + 1], s[j % n + 1]))
j = j % n + 1;
res = max(res, max(fn.dis(s[i], s[j])
, fn.dis(s[i + 1], s[j])));
}
return res;
}
};
void solve() {
int n; cin >> n;
vv<Point>nums(n);
forn(i, n)
cin >> nums[i].x >> nums[i].y;
Andrew aw;
vv<Point >rs(n);
aw.Andrew_lenth(nums, rs);
LD ans = aw.rotating_calipers(rs);
cout<<ans*ans<< endl;
}
int main() {
T = 1;
while (T--)
solve();
}