bzoj1458士兵占领

传送门

和上一题差不多,每行和每列分别看做一个点,障碍点坐标的行和列就不建边,再按照有源汇上下界建图就好了,唯一的区别就是这个题求的是最小流
这个题的数据好水呢,建错图也能A呢

#include<cstdio>
#include<queue>
#include<cstring>
#include<iostream>
using namespace std;
#define rg register
#define min(a,b) (a<b?a:b)
int dis[2001],ans,n,m,k,cnt=1,s,t,ss,tt,sum,inf=1e9,pre[31000],nxt[31000],h[2001],v[30001],cur[2001],in[2001];
queue<int>q;bool mp[101][101];
void read(int &x) {
    char ch; bool ok;
    for(ok=0,ch=getchar(); !isdigit(ch); ch=getchar()) if(ch=='-') ok=1;
    for(x=0; isdigit(ch); x=x*10+ch-'0',ch=getchar()); if(ok) x=-x;
}
void add(int x,int y,int z)
{
    pre[++cnt]=y,nxt[cnt]=h[x],h[x]=cnt,v[cnt]=z;
    pre[++cnt]=x,nxt[cnt]=h[y],h[y]=cnt;
}
bool bfs(int s,int t)
{
    memset(dis,0,sizeof dis);
    q.push(s),dis[s]=1;
    while(!q.empty())
    {
        int x=q.front();q.pop();
        for(int i=h[x];i;i=nxt[i])if(!dis[pre[i]]&&v[i])dis[pre[i]]=dis[x]+1,q.push(pre[i]);
    }
    return dis[t];
}
int dfs(int x,int flow,int t)
{
    if(x==t||!flow)return flow;
    int f=flow;
    for(int &i=cur[x];i;i=nxt[i])
        if(v[i]&&dis[pre[i]]>dis[x])
        {
            int y=dfs(pre[i],min(v[i],f),t);
            f-=y,v[i]-=y,v[i^1]+=y;
            if(!f)return flow;
        }
    if(f==flow)dis[x]=-1;
    return flow-f;
}
int main()
{
    read(n),read(m),read(k);s=0,t=n+m+1,ss=t+1,tt=ss+1;
    for(rg int i=1,x;i<=n;i++)read(x),add(s,i,inf-x),in[i]+=x,in[s]-=x;
    for(rg int i=1,x;i<=m;i++)read(x),add(i+n,t,inf-x),in[t]+=x,in[i+n]-=x;
    for(rg int i=1,x,y;i<=k;i++)read(x),read(y),mp[x][y]=1;
    for(rg int i=1;i<=n;i++)
    	for(rg int j=1;j<=m;j++)
    		if(!mp[i][j])add(i,j+n,1);
    for(int i=0;i<=t;i++)
    {
    	if(in[i]>0)add(ss,i,in[i]),sum+=in[i];
    	if(in[i]<0)add(i,tt,-in[i]);
    }
    add(t,s,inf);
    for(;bfs(ss,tt);ans+=dfs(ss,inf,tt))memcpy(cur,h,sizeof cur);
    if(ans!=sum){printf("JIONG\n");return 0;}
    sum=0,ans=v[cnt];
    for(int i=h[ss];i;i=nxt[i])v[i]=v[i^1]=0;
    for(int i=h[tt];i;i=nxt[i])v[i]=v[i^1]=0;
    v[cnt]=v[cnt^1]=0;
    for(;bfs(t,s);sum+=dfs(t,inf,s))memcpy(cur,h,sizeof cur);
    printf("%d\n",ans-sum);
}
posted @ 2019-01-05 14:04  蒟蒻--lichenxi  阅读(97)  评论(0编辑  收藏  举报