bzoj1340: [Baltic2007]Escape逃跑问题

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

#include<bits/stdc++.h>#define N 505using 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);}