1. 程式人生 > >2017-04-04 水題信心 03選美比賽

2017-04-04 水題信心 03選美比賽

Description

約翰農夫的母牛貝西,剛剛在牛選美比賽贏得第一名,贏得了“世界牛小姐”稱號。為了傳播善意,貝西將參觀世界各地的農場 N(2 < = N < = 50000)。為簡單起見,世界將被表示成一個二維平面,其中每個農場位於一對整數座標(x,y),其值在-2000000~2000000 範圍內。沒有兩個農場共享相同的座標。
儘管貝西旅行中走農場之間的直線,但一些農場之間的距離可能很大,所以她想帶著一個手提箱,足夠她裝食物,即使最遠的距離。她想確定旅行中最大可能的距離,這樣她才能決定箱子的大小。幫助貝西計算所有成對的農場間最大距離。(最遠點對)

Input

第一行一個整數 N
以下 N 行,每行兩個整數,表示一個農場的座標。

Output

一行,一個整數表示最遠一對農場距離的平方。

Sample Input

4
0 0
0 1
1 1
1 0

Sample Output

2

題解

裸的最遠點對。
旋xuan2轉zhuan3卡qia3殼ke2即可。

#include<cmath>
#include<cstdio>
#include<iostream>
#include<algorithm>
using namespace std;

typedef long long ll;

inline int read(){
    int x  = 0, f = 1
; char c = getchar(); while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); } while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); } return x * f; } const int N = 50000 + 10; struct P{ ll x, y; P(){}; P(ll a, ll b){x = a, y = b;} }p[N], s[N]; int n, top; ll ans; ll sqr(ll x){return
x*x;} ll dis(P a, P b){return sqr(a.x - b.x) + sqr(a.y - b.y);} P operator - (P a, P b){ return P(a.x - b.x, a.y - b.y); } ll operator * (P a, P b){ return a.x * b.y - a.y * b.x; } bool operator < (P a, P b){ ll t = (a - p[1]) * (b - p[1]); return t > 0 || (t == 0 && dis(p[1], a) < dis(p[1], b)); } void init(){ n = read(); int k = 1; for(int i = 1; i <= n; i++){ p[i].x = read(), p[i].y = read(); if(p[i].x < p[k].x || (p[i].x == p[k].x && p[i].y < p[k].y)) k = i; } swap(p[1], p[k]); sort(p + 2, p + n + 1); } void graham(){ s[++top] = p[1]; s[++top] = p[2]; for(int i = 3; i <= n; i++){ while(top > 1 && (s[top] - s[top-1]) * (p[i] - s[top-1]) <= 0) top--; s[++top] = p[i]; } } void RC(){ s[top+1] = s[1]; int now = 2; for(int i = 1; i <= top; i++){ while((s[i+1] - s[i]) * (s[now] - s[i]) < (s[i+1] - s[i]) * (s[now+1] - s[i])){ now++; if(now == top + 1) now = 1; } ans = max(ans, dis(s[now], s[i])); } } void work(){ graham(); RC(); cout<<ans<<endl; } int main(){ freopen("contest.in", "r", stdin); freopen("contest.out", "w", stdout); init(); work(); return 0; }