2017-04-04 水題信心 03選美比賽
阿新 • • 發佈:2019-02-12
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;
}