首页 > 代码库 > 【BZOJ 1185】 凸包+旋转卡壳
【BZOJ 1185】 凸包+旋转卡壳
Description
【分析】
打计算几何真的可以哭出来。。。
跟那个求线段最远点差不多,这题弄三个东西转一转,一个表示左端最远点,一个表示右端最远点,一个表示上面最远点。
左右两边的最远点用点积判断,上方最远点用差积判断。
【向量是最好的解决平面几何问题的,不要画图!!!哭!!!
输出那里,用的是向量加减法,乘一个比例系数,然后还用到了垂直单位向量。
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cstring> 4 #include<iostream> 5 #include<algorithm> 6 #include<cmath> 7 using namespace std; 8 #define Maxn 50010 9 10 const double eps=0.0000001; 11 const double INF=999999999; 12 13 struct P 14 { 15 double x,y; 16 }a[Maxn],t[Maxn]; 17 int len,n; 18 19 P operator - (P x,P y) 20 { 21 P tt; 22 tt.x=x.x-y.x; 23 tt.y=x.y-y.y; 24 return tt; 25 } 26 27 P operator * (P x,double y) 28 { 29 P tt; 30 tt.x=x.x*y; 31 tt.y=x.y*y; 32 return tt; 33 } 34 35 P operator + (P x,P y) 36 { 37 P tt; 38 tt.x=x.x+y.x; 39 tt.y=x.y+y.y; 40 return tt; 41 } 42 43 double myabs(double x) {return x>0?x:-x;} 44 45 int fbs(double x) 46 { 47 if(myabs(x)<=eps) return 0; 48 return x<0?-1:1; 49 } 50 51 double Dot(P x,P y) {return x.x*y.x+x.y*y.y;} 52 double Cross(P x,P y) {return x.x*y.y-x.y*y.x;} 53 bool cmp(P x,P y) {return fbs(x.x-y.x)==0?(x.y<y.y):(x.x<y.x);} 54 55 void chull() 56 { 57 sort(a+1,a+1+n,cmp); 58 len=0; 59 for(int i=1;i<=n;i++) 60 { 61 while(len>1&&Cross(t[len]-t[len-1],a[i]-t[len])<=0) len--; 62 t[++len]=a[i]; 63 } 64 int k=len; 65 for(int i=n-1;i>=1;i--) 66 { 67 while(len>k&&Cross(t[len]-t[len-1],a[i]-t[len])<=0) len--; 68 t[++len]=a[i]; 69 }len--; 70 t[0]=t[len]; 71 } 72 73 void output() 74 { 75 for(int i=0;i<len;i++) 76 { 77 printf("%lf %lf\n",t[i].x,t[i].y); 78 }printf("\n"); 79 } 80 81 double ans; 82 P op[5]; 83 84 void RC() 85 { 86 double L,R,H,D; 87 ans=INF; 88 int l=1,r=1,h=1; 89 for(int i=0;i<len;i++) 90 { 91 int j=(i+1)%len; 92 // l=i;r=j;h=j; 93 // printf("%lf %lf\n",(t[h+1]-t[i]).x,(t[h+1]-t[i]).y); 94 // printf("%lf %lf\n",(t[h+1]-t[j]).x,(t[h+1]-t[j]).y); 95 while(Cross(t[h]-t[i],t[h]-t[j])<=Cross(t[h+1]-t[i],t[h+1]-t[j])) h=(h+1)%len; 96 while(Dot(t[i]-t[j],t[r]-t[j])>=Dot(t[i]-t[j],t[r+1]-t[j])) r=(r+1)%len; 97 if(i==0) l=r; 98 99 // printf("%lf %lf\n",Dot(t[j]-t[i],t[l]-t[i]),Dot(t[j]-t[i],t[l+1]-t[i])); 100 while(Dot(t[j]-t[i],t[l]-t[i])>=Dot(t[j]-t[i],t[l+1]-t[i])) l=(l+1)%len; 101 102 D=sqrt(Dot(t[i]-t[j],t[i]-t[j])); 103 H=Cross(t[h]-t[i],t[h]-t[j])/D; 104 L=-Dot(t[j]-t[i],t[l]-t[i])/D; 105 R=-Dot(t[i]-t[j],t[r]-t[j])/D; 106 107 // printf("**%lf %lf %lf %lf\n",D,H,L,R); 108 109 double now=(R+L+D)*H; 110 // printf("%lf %lf %d %d %d %lf\n",t[i].x,t[i].y,l,r,h,now); 111 if(now<ans) 112 { 113 ans=now; 114 op[0]=t[i]-(t[j]-t[i])*(L/D); 115 op[1]=t[j]-(t[i]-t[j])*(R/D); 116 P tt; 117 tt.x=(t[i]-t[j]).y;tt.y=-(t[i]-t[j]).x; 118 op[2]=op[1]+tt*(H/D); 119 op[3]=op[0]+tt*(H/D); 120 // op[0]=t[i];op[1]=t[j];op[2]=t[r];op[3]=t[l]; 121 } 122 } 123 } 124 125 int main() 126 { 127 scanf("%d",&n); 128 for(int i=1;i<=n;i++) 129 { 130 scanf("%lf%lf",&a[i].x,&a[i].y); 131 } 132 chull(); 133 // output(); 134 RC(); 135 int fr=0; 136 for(int i=1;i<4;i++) if(op[i].y<op[fr].y||fbs(op[fr].y-op[i].y)==0&&op[i].x<op[i].x) fr=i; 137 printf("%.5lf\n",ans); 138 for(int i=0;i<4;i++) 139 { 140 if(fbs(op[i].x)==0) op[i].x=0; 141 if(fbs(op[i].y)==0) op[i].y=0; 142 } 143 // for(int i=0;i<4;i++) op[i].x=myabs(op[i].x),op[i].y=myabs(op[i].y); 144 for(int i=0;i<4;i++) 145 { 146 printf("%.5lf %.5lf\n",op[(fr+i)%4].x,op[(fr+i)%4].y); 147 } 148 return 0; 149 }
调试过程不删了,心酸。。
2016-12-15 16:55:18
【BZOJ 1185】 凸包+旋转卡壳
声明:以上内容来自用户投稿及互联网公开渠道收集整理发布,本网站不拥有所有权,未作人工编辑处理,也不承担相关法律责任,若内容有误或涉及侵权可进行投诉: 投诉/举报 工作人员会在5个工作日内联系你,一经查实,本站将立刻删除涉嫌侵权内容。