日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

POJ 3608 旋转卡壳

發布時間:2023/11/29 编程问答 31 豆豆
生活随笔 收集整理的這篇文章主要介紹了 POJ 3608 旋转卡壳 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

思路:

旋轉卡殼應用

注意點&邊? 邊&邊? 點&點 三種情況

//By SiriusRen #include <cmath> #include <cstdio> #include <algorithm> using namespace std; const double eps=1e-5; const int N=10050; typedef double db; int n,m; struct P{db x,y;P(){}P(db X,db Y){x=X,y=Y;}}p1[N],p2[N]; P operator-(P a,P b){return P(a.x-b.x,a.y-b.y);} db operator*(P a,P b){return a.x*b.y-a.y*b.x;} db operator^(P a,P b){return a.x*b.x+a.y*b.y;} db dis(P c){return sqrt(c.x*c.x+c.y*c.y);} db dotc(P a,P b,P c){return (b-a)^(c-a);} db difc(P a,P b,P c){return (b-a)*(c-a);} db p2seg(P a,P b,P c){if(dotc(a,b,c)<-eps)return dis(a-c);if(dotc(b,a,c)<-eps)return dis(b-c);return abs(difc(a,b,c)/dis(a-b)); } db seg2seg(P a,P b,P c,P d){return min(min(p2seg(a,b,c),p2seg(a,b,d)),min(p2seg(c,d,a),p2seg(c,d,b))); } double solve(P p[],P q[],int np,int nq){int sp=1,sq=1;for(int i=1;i<=np;i++)if(p[i].y<p[sp].y)sp=i;for(int i=1;i<=nq;i++)if(q[i].y>q[sq].y)sq=i;p[np+1]=p[1],q[nq+1]=q[1];double tmp,ans=1e9;for(int i=1;i<=np;i++){while(tmp=(difc(p[sp+1],q[sq+1],p[sp])-difc(p[sp+1],q[sq],p[sp]))>eps)sq=sq%nq+1;if(tmp<-eps)ans=min(ans,p2seg(p[sp],p[sp+1],q[sq]));else ans=min(ans,seg2seg(p[sp],p[sp+1],q[sq],q[sq+1]));sp=sp%np+1;}return ans; } int main(){while(scanf("%d%d",&n,&m)&&(n||m)){for(int i=1;i<=n;i++)scanf("%lf%lf",&p1[i].x,&p1[i].y);for(int i=1;i<=m;i++)scanf("%lf%lf",&p2[i].x,&p2[i].y);printf("%lf\n",min(solve(p1,p2,n,m),solve(p2,p1,m,n)));} }

?

轉載于:https://www.cnblogs.com/SiriusRen/p/9393188.html

創作挑戰賽新人創作獎勵來咯,堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的POJ 3608 旋转卡壳的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。