bzoj3630 镜面通道 网络流

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

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

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

AC代码如下：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#define N 200005
#define inf 1000000000
using 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