bzoj1340: [Baltic2007]Escape逃跑问题

传送门
我们大力拆点，在相邻两个距离小于等于200的点连边
然后S向离上边界小于等于100的点连边
离下边界小于等于100的点向T连边
大力最小割一点不虚。

#include<bits/stdc++.h>
#define N 505
using namespace std;
int head[N],q[N],dis[N];
int tot=1,l,w,n,S,T,ans;
struct edge{int to,next,f;}e[N*N];
struct P{int x,y;}p[N];
long long sqr(int x){
    return 1ll*x*x;
}
void add(int x,int y,int z){
    e[++tot]=(edge){y,head[x],z};
    head[x]=tot;
    e[++tot]=(edge){x,head[y],0};
    head[y]=tot;
}
bool spfa(){
    memset(dis,-1,sizeof(dis));
    dis[S]=0; q[1]=S;
    for (int h=0,t=1;h!=t;)
        for (int x=q[++h],i=head[x];i;i=e[i].next)
            if (dis[e[i].to]==-1&&e[i].f)
                dis[e[i].to]=dis[x]+1,q[++t]=e[i].to;
    return dis[T]!=-1;
}
int dfs(int x,int flow){
    if (x==T) return flow;
    int k,rest=flow;
    for (int i=head[x];i&&rest;i=e[i].next)
        if (dis[e[i].to]==dis[x]+1&&e[i].f){
            k=dfs(e[i].to,min(rest,e[i].f));
            rest-=k; e[i].f-=k; e[i^1].f+=k;
        }
    if (rest) dis[x]=-1;
    return flow-rest;
}
int main(){
    scanf("%d%d%d",&l,&w,&n);
    for (int i=1;i<=n;i++)
        scanf("%d%d",&p[i].x,&p[i].y);
    S=n*2+1; T=S+1;
    for (int i=1;i<=n;i++){
        add(i,i+n,1);
        if (p[i].y<=100) add(S,i,1e9);
        if (w-p[i].y<=100) add(i+n,T,1e9);
    }
    for (int i=1;i<=n;i++)
        for (int j=1;j<=n;j++)
            if (sqr(p[i].x-p[j].x)+sqr(p[i].y-p[j].y)<=40000&&i!=j)
                add(i+n,j,1e9);
    while (spfa()) ans+=dfs(S,1e9);
    printf("%d",ans);
}