bzoj2564 集合的面積

Description

　　對於一個平面上點的集合P={(x _i,y _i )}，定義集合P的面積F(P)為點集P的凸包的面積。
　　對於兩個點集A和B，定義集合的和為：
　　A+B={(x _i ^A+x _j ^B,y _i ^A+y _j ^B ):(x _i ^A,y _i ^A )∈A,(x _j ^B,y _j ^B )∈B}
　　現在給定一個N個點的集合A和一個M個點的集合B，求2F(A+B)。

Input

　第一行包含用空格隔開的兩個整數，分別為N和M；
　　第二行包含N個不同的數對，表示A集合中的N個點的坐標；
　　第三行包含M個不同的數對，表示B集合中的M個點的坐標。

Output

　一共輸出一行一個整數，2F(A+B)。

Sample Input

4 5
0 0 2 1 0 1 2 0
0 0 1 0 0 2 1 2 0 1

Sample Output

18
數據規模和約定
對於30%的數據滿足N ≤ 200，M ≤ 200；
對於100%的數據滿足N ≤ 10^5，M ≤ 10^5，|xi|, |yi| ≤ 10^8。

 

正解：$Minkowski$和。

$Minkowski$和，就是題面的這個東西。。

分別求出兩個點集的凸包，然后貪心地加點就行。

首先，$A$凸包和$B$凸包的第一個點的和肯定會在最終的凸包里。

然后我們設$A$凸包到了$i$點，$B$凸包到了$j$點。

如果$a[i+1]+b[j]$比$a[i]+b[j+1]$更凸，那么就用$A$凸包更新，否則用$B$凸包更新。

最后求出的這個就是新凸包，直接用叉積算面積就行了。

 

#include <bits/stdc++.h>
#define il inline
#define RG register
#define ll long long
#define N (500010)
 
using namespace std;
 
struct point{
  ll x,y;
  il point operator + (const point &a) const{
    return (point){x+a.x,y+a.y};
  }
  il point operator - (const point &a) const{
    return (point){x-a.x,y-a.y};
  }
}p[N],t1[N],t2[N],st[N];
 
int n,m,top;
ll S;
 
il int gi(){
  RG int x=0,q=1; RG char ch=getchar();
  while ((ch<'0' || ch>'9') && ch!='-') ch=getchar();
  if (ch=='-') q=-1,ch=getchar();
  while (ch>='0' && ch<='9') x=x*10+ch-48,ch=getchar();
  return q*x;
}
 
il int cmp(const point &a,const point &b){
  if (a.x==b.x) return a.y<b.y; return a.x<b.x;
}
 
il ll cross(RG point a,RG point b){ return a.x*b.y-a.y*b.x; }
 
il void graham(point *p,point *t,RG int n){
  sort(p+1,p+n+1,cmp);
  for (RG int i=1;i<=n;++i){
    while (top>=2 && cross(p[i]-t[top-1],t[top]-t[top-1])>=0) --top;
    t[++top]=p[i];
  }
  for (RG int i=n-1,la=top;i>=1;--i){
    while (top>la && cross(p[i]-t[top-1],t[top]-t[top-1])>=0) --top;
    t[++top]=p[i];
  }
  --top; return;
}
 
int main(){
#ifndef ONLINE_JUDGE
  freopen("area.in","r",stdin);
  freopen("area.out","w",stdout);
#endif
  n=gi(),m=gi();
  for (RG int i=1;i<=n;++i) p[i].x=gi(),p[i].y=gi(); graham(p,t1,n),n=top,top=0;
  for (RG int i=1;i<=m;++i) p[i].x=gi(),p[i].y=gi(); graham(p,t2,m),m=top,top=0;
  st[top=1]=t1[1]+t2[1];
  for (RG int i=1,j=1;i<=n || j<=m;){
    RG point x=t1[(i-1)%n+1]+t2[j%m+1],y=t1[i%n+1]+t2[(j-1)%m+1];
    if (cross(x-st[top],y-st[top])>=0) st[++top]=x,++j; else st[++top]=y,++i;
  }
  for (RG int i=2;i<top;++i) S+=cross(st[i]-st[1],st[i+1]-st[1]); cout<<S; return 0;
}