bzoj1340: [Baltic2007]Escape逃跑问题

题意：http://www.lydsy.com/JudgeOnline/problem.php?id=1340

挺好的一道题...

有曲线不能从左下到右上等价于，不存在从左上到右下的曲线，满足曲线一直在圆内或矩形边界上

所以删除数量最小的圆，使得上下不连通即可

 

#include<bits/stdc++.h>
using namespace std;
#define rep(i,l,r) for(int i=l;i<=r;++i)
const int N=202333,inf=2147483647;
int head[N],n,m,dis[N],tot=1,sum,T,K;
struct node{
    int x,y;
}s[N];
struct zs{
    int to,next,w;
}e[N];
inline void ins(int u,int v,int w){
    e[++tot].to=v; e[tot].next=head[u]; head[u]=tot; e[tot].w=w;
}
inline void insert(int u,int v,int w){
    ins(u,v,w); ins(v,u,0);
}
inline bool bfs(){
     rep(i,0,T) dis[i]=-1; queue<int>q; q.push(0); dis[0]=0;
     while(!q.empty()) {
          int x=q.front(); q.pop();
          for(int k=head[x];k;k=e[k].next) 
             if(dis[e[k].to]<0 && e[k].w>0) {
                   dis[e[k].to]=dis[x]+1; q.push(e[k].to);
             }
     }
     if(dis[T]>0) return 1;else return 0;
}
int find(int x,int low){
     if(x==T) return low;
     int delta=low,now;
     for(int k=head[x];k;k=e[k].next) 
       if(e[k].w>0 && dis[e[k].to]==dis[x]+1){ 
           now=find(e[k].to,min(e[k].w,delta));
           e[k].w-=now; e[k^1].w+=now;   delta-=now;
           if(!delta) return low;
        } 
     dis[x]=-1;
     return low-delta;
}
int main(){
    scanf("%d%d%d",&n,&m,&K);
    T=K*2+1;
    rep(i,1,K) {
        scanf("%d%d",&s[i].x,&s[i].y);
        if(s[i].y<=100) insert(0,i,inf);
        if(abs(s[i].y-m)<=100) insert(i+K,T,inf);
        insert(i,i+K,1);
    }
    rep(i,1,K) rep(j,1,K) if(i!=j&&(s[i].x-s[j].x)*(s[i].x-s[j].x)+(s[i].y-s[j].y)*(s[i].y-s[j].y)<=200*200)
                            insert(i+K,j,inf);
    while(bfs()) sum+=find(0,inf);
    printf("%d\n",sum);
}