bzoj1091 切割多边形 枚举&计算几何

       发现边很少只有8条。。那么我们可以枚举每一次用多边形的那一条边去割外面的矩形。显然每一条边一定且只会割一次。然么就是直线求交点了，然后把直线外侧的点全都删掉。

       小心割线刚好和外面的多边形在定点相交的特殊情况；以及割最后一条边的特殊情况。

       嘴巴AC还是很简单的。。实现就233（其实不算太复杂，细节坑爹啊）

       最后吐槽：woc C++什么破编译器调试的时候加一句cout<<i答案居然会不一样。。。真是让人吐血。。。希望有高手看到向我解释为什么~\(≧▽≦)/~辣！！

AC代码如下：

#include<iostream>#include<cstdio>#include<cstring>#include<cmath>#define N 25using namespace std;int n,cnt; double ans=1e200; bool bo[N];struct point{ double x,y; }a[N],b[N];point operator -(point u,point v){u.x-=v.x; u.y-=v.y; return u;}double crs(point u,point v){ return u.x*v.y-u.y*v.x; }double dist(point u,point v){return sqrt((u.x-v.x)*(u.x-v.x)+(u.y-v.y)*(u.y-v.y));}bool upon(point u,point v,point w){return !crs(u-v,w-v) && (u.x-v.x)*(u.x-w.x)<=0 && (u.y-v.y)*(u.y-w.y)<=0;}bool havitr(point p,point q,point u,point v){return !upon(p,u,v) && crs(u-p,q-p)<=0 && crs(v-p,q-p)>=0;}point itr(point u,point s,point v,point t){double tmp=crs(t,u-v)/crs(s,t);u.x+=s.x*tmp; u.y+=s.y*tmp; return u;}void dfs(int k,double now){if (now>=ans) return; int i,j;if (k>n){ ans=now; return; }point c[N],p1,p2; int len=cnt;for (i=1; i<=cnt+1; i++) c[i]=b[i];for (i=1; i<=n; i++) if (bo[i]){bo[i]=0; cnt=0; bool ok=1;for (j=1; j<=len; j++)if (havitr(a[i+1],a[i],c[j],c[j+1])){p1=itr(a[i+1],a[i]-a[i+1],c[j],c[j+1]-c[j]); break;}for (j=1; j<=len; j++)if (havitr(a[i],a[i+1],c[j],c[j+1])){p2=itr(a[i],a[i+1]-a[i],c[j],c[j+1]-c[j]); break;}for (j=1; j<=len; j++)if (crs(p2-p1,c[j]-p1)<0) b[++cnt]=c[j];else if (ok){ ok=0; b[++cnt]=p1; b[++cnt]=p2; }b[cnt+1]=b[1];dfs(k+1,now+dist(p1,p2));cnt=len; for (j=1; j<=cnt+1; j++) b[j]=c[j];bo[i]=1;}}int main(){scanf("%lf%lf",&b[3].x,&b[3].y); int i;b[4].x=b[3].x; b[2].y=b[3].y; cnt=4; scanf("%d",&n);for (i=1; i<=n; i++) scanf("%lf%lf",&a[i].x,&a[i].y);a[n+1]=a[1]; b[cnt+1]=b[1];memset(bo,1,sizeof(bo));dfs(1,0); printf("%.3f\n",ans);return 0;}

by lych

2016.3.3