A. Peter and Snow Blower 解析(思維、幾何)

Codeforce 613 A. Peter and Snow Blower 解析(思維、幾何)

今天我們來看看CF613A
題目連結

題目
給你一個點(P)和(n)個點形成的多邊形(照順或逆時針順序給)，求這個多邊形繞著(P)轉最後可以造成的面積。(有關正式的"旋轉"定義請看原題)

前言

儲存點的座標時沒想過要用(pair<long long,long long>)，結果debug超久

想法

首先要注意到:由於題目的旋轉的定義是把每個點都對於點(P)去做旋轉，所以最後的圖形一定是兩個同心圓，而面積就是兩個圓中間的面積，而我們只需要維護最長的半徑和最短的半徑就好。
由於題目是按照順序給多邊形的點，所以我們可以把每條邊單獨拿出來考慮和(P)點的最短和最長距離。


如上圖所示，想要判斷點(P)到線段(overline{SE})的最短距離線段是否在線段(overline{SE})上，我們只需要判斷(overrightarrow{PM})是否被(overrightarrow{PS},overrightarrow{PE})所包住，而其中一種方法就是利用外積(叉積、cross product):
如果(overrightarrow{PM})是被包住的，那麼(sgn(overrightarrow{PM} imesoverrightarrow{PS})=-sgn(overrightarrow{PM} imesoverrightarrow{PE}))
反之如果(sgn(overrightarrow{PM} imesoverrightarrow{PS})=sgn(overrightarrow{PM} imesoverrightarrow{PE}))，那麼代表沒有被包住。以上是利用了外積的性質:(overrightarrow{AB} imesoverrightarrow{CD}=-overrightarrow{CD} imesoverrightarrow{AB})對於任何向量(overrightarrow{AB},overrightarrow{CD})。
而要計算最短距離，我們有兩種方法:

1. 利用內積是投影長度的相乘的性質，我們把線段的法向量和(overrightarrow{PE})作內積，再除以法向量的長度，就是最短距離。
2. 利用外積的絕對值是向量們所展出的四邊形面積，且等於底乘以高，(|overrightarrow{PS} imesoverrightarrow{PE}|/|overrightarrow{SE}|)就是最短距離。

而透過觀察可以發現，(P)點到線段的長度，不是最短距離，那就是端點。有了以上資訊，我們就可以寫了。

程式碼:

const int _n=1e5+10;
int t,n,m;
PII p,prev,ps[_n];
db minn=1e9,maxx=-1e9,pi=acos(-1);
bool sgn(db x){
  return x>=0.0?0:1;
}
db cp(PII u,PII v){
  return (db)(u.fi*v.se-u.se*v.fi);
}
db len(PII u){
  return sqrt(u.fi*u.fi+u.se*u.se);
}
void f(PII x,PII y){
  PII tt2={y.fi-p.fi,y.se-p.se},tt3={x.fi-p.fi,x.se-p.se},tt1={-(tt3.se-tt2.se),tt3.fi-tt2.fi};
  db res1=len(tt2),res2=len(tt3),res3=abs((db)(tt1.fi*tt2.fi+tt1.se*tt2.se))/len(tt1);
  bool z=1;if(sgn(cp(tt1,tt2))==sgn(cp(tt1,tt3)))z=0;
  if(z){
    minn=min(minn,min(res1,min(res2,res3)));
    maxx=max(maxx,max(res1,max(res2,res3)));
  }else{
    minn=min(minn,min(res1,res2));
    maxx=max(maxx,max(res1,res2));
  }
}
main(void) {ios_base::sync_with_stdio(0);cin.tie(0);cout.tie(0);
  //這邊的PII必須是pair<ll,ll>
  cin>>n>>p.fi>>p.se;rep(i,0,n)cin>>ps[i].fi>>ps[i].se; prev=ps[0];
  rep(i,1,n)f(prev,ps[i]),prev=ps[i];
  f(prev,ps[0]);
  cout<<setprecision(20)<<pi*(maxx*maxx-minn*minn)<<'
';
  return 0;
}

標頭、模板請點Submission看
Submission