上(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