題目來源：http://poj.org/problem?id=2318

一個(gè) 矩陣 被分成多個(gè) 區(qū)域， 然后輸入 多個(gè)點(diǎn)， 輸出 每個(gè)區(qū)域點(diǎn)的 個(gè)數(shù)。 當(dāng)尋找點(diǎn) 落在某個(gè)區(qū)域時(shí)， 用二分法。

?

#include <cstdlib> #include <cstring> #include <algorithm> #include <cstdio> #include <cmath> #include <iostream> #include <vector> using namespace std; typedef long long ll; const int N =6000; const double PI = 3.1415927;struct Point{int x,y;Point(){}Point(int x,int y):x(x),y(y){} // 構(gòu)造函數(shù)，方便代碼編寫Point(const Point & p):x(p.x),y(p.y){}Point operator +(Point p){return Point(x+p.x,y+p.y);}Point operator-(Point p){return Point(x-p.x,y-p.y);}Point operator*(double d){return Point(x*d,y*d);}int operator*(Point p){ // 內(nèi)積 點(diǎn)乘return x*p.x+ y*p.y;}int operator^(Point p){// 外積 叉乘return x*p.y-y*p.x;}friend ostream& operator<<(ostream& os,const Point& p ){os<<p.x<<" "<<p.y<<endl;return os;}friend istream& operator>>(istream& is, Point& p) {// Point 不能是常量，是變量is>>p.x>>p.y;return is;} }; Point up[N]; Point low[N]; int num[N]; int n; void bin_search(Point p) // 二分法，尋找點(diǎn)所在的區(qū)域，用叉乘，注意結(jié)束條件是 left+1== right {int left=0,mid;int right = n+1;while(left<=right){mid=(left+right) /2;int t=(up[mid]-low[mid])^(p-low[mid]);if( t >0 )right=mid;else if( t <0 )left=mid;if(left+1 == right ){num[left]++;break;}} } int main() {int m;int a,b,count=0;Point p;while(scanf("%d",&n)!=EOF && n){if(count) printf("\n");count=1;scanf("%d",&m);scanf("%d%d%d%d",&up[0].x,&up[0].y,&low[n+1].x,&low[n+1].y);low[0].x=up[0].x;up[n+1].x=low[n+1].x;for(int i=0;i<=n+1;i++){up[i].y=up[0].y;low[i].y=low[n+1].y;}for(int i=1;i<=n;i++){scanf("%d%d",&a,&b);up[i].x=a;low[i].x=b;}memset(num,0,sizeof(num));for(int i=1;i<=m;i++){scanf("%d%d",&p.x,&p.y);bin_search(p);}for(int i=0;i<=n;i++)printf("%d: %d\n",i,num[i]);}return 0; }

?

轉(zhuǎn)載于:https://www.cnblogs.com/zn505119020/p/3623613.html

總結(jié)

以上是生活随笔為你收集整理的寻找区域中有几个点 叉乘+二分 poj 2318的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。