bzoj3630 镜面通道 网络流

        上(xiao)帝(mao)说：“要有光！”于是只要前后是连通的光就能够通过。。（表问我为什么。这是一个定理）

       所以题目就变成了在相交的两个元件之间连边，S向所有和底部相交的元件连边，T向所有和顶部相交的元件连边。如果光不能通过，相当于存在一条向S->T的路径。然后就变成最小割了。

       判断相交真是太烦了啊啊啊啊！！>_<

AC代码如下：

#include<iostream>#include<cstdio>#include<cstring>#include<cstdlib>#include<cmath>#define N 200005#define inf 1000000000using namespace std;int n,tot=1,fst[N],sta,gol,pnt[N],len[N],nxt[N],d[N],h[N];struct point{ double x,y; };struct obj{ int k; point p,q; double r; }a[N];double dist(point u,point v){return sqrt((u.x-v.x)*(u.x-v.x)+(u.y-v.y)*(u.y-v.y));}void add(int x,int y,int z){pnt[++tot]=y; len[tot]=z; nxt[tot]=fst[x]; fst[x]=tot;}void ins(int x,int y,int z){ add(x,y,z); add(y,x,0); }bool bfs(){memset(d,-1,sizeof(d)); d[sta]=1;int head=0,tail=1; h[1]=sta;while (head<tail){int x=h[++head],p;for (p=fst[x]; p; p=nxt[p]) if (len[p]){int y=pnt[p];if (d[y]==-1){ d[y]=d[x]+1; h[++tail]=y; }}}return d[gol]!=-1;}int dfs(int x,int rst){if (x==gol || !rst) return rst; int p,flow=0;for (p=fst[x]; p; p=nxt[p]) if (len[p]){int y=pnt[p]; if (d[x]+1!=d[y]) continue;int tmp=dfs(y,min(rst,len[p])); if (!tmp) continue;flow+=tmp; len[p]-=tmp;len[p^1]+=tmp; rst-=tmp;if (!rst) break;}if (!flow) d[x]=-1; return flow;}bool ok(obj u,obj v){if (u.k==1 && v.k==1) return dist(u.p,v.p)<=u.r+v.r;if (u.k==2 && v.k==2) return u.p.x<=v.q.x && u.p.y<=v.q.y && u.q.x>=v.p.x && u.q.y>=v.p.y;if (u.k>v.k) swap(u,v);if (dist(u.p,v.p)<=u.r || dist(u.p,v.q)<=u.r) return 1;point t; t.x=v.p.x; t.y=v.q.y; if (dist(u.p,t)<u.r) return 1;t.x=v.q.x; t.y=v.p.y; if (dist(u.p,t)<u.r) return 1;if (v.p.x<=u.p.x && u.p.x<=v.q.x && (fabs(v.p.y-u.p.y)<=u.r || fabs(v.q.y-u.p.y)<=u.r)) return 1;if (v.p.y<=u.p.y && u.p.y<=v.q.y && (fabs(v.p.x-u.p.x)<=u.r || fabs(v.q.x-u.p.x)<=u.r)) return 1;return v.p.x<=u.p.x && u.p.x<=v.q.x && v.p.y<=u.p.y && u.p.y<=v.q.y;}int main(){point t; scanf("%lf%lf%d",&t.x,&t.y,&n); int i,j;for (i=1; i<=n; i++){scanf("%d",&a[i].k);if (a[i].k==1) scanf("%lf%lf%lf",&a[i].p.x,&a[i].p.y,&a[i].r);else scanf("%lf%lf%lf%lf",&a[i].p.x,&a[i].p.y,&a[i].q.x,&a[i].q.y);}for (i=1; i<=n; i++) ins(i,i+n,1);a[0].k=2; a[0].p.y=-1; a[0].q.x=a[n+1].q.x=t.x;a[n+1].k=2; a[n+1].p.y=t.y; a[n+1].q.y=t.y+1;for (i=1; i<=n; i++)for (j=1; j<=n; j++)if (i!=j && ok(a[i],a[j])) ins(i+n,j,1);sta=0; gol=n<<1|1;for (i=1; i<=n; i++){if (ok(a[i],a[0])) ins(0,i,1);if (ok(a[i],a[n+1])) ins(i+n,gol,1);}int ans=0; while (bfs()) ans+=dfs(sta,inf);printf("%d\n",ans);return 0;}

by lych

2016.3.3