bzoj1458: 士兵占领 网络流

Posted 2021-02-12 dsrdsr

链接

https://www.lydsy.com/JudgeOnline/problem.php?id=1458
也可以去luogu

思路

想成倒着删去点，使得依旧满足覆盖！！
左边横，右边列，之间用点链接起来，流量1
S到i为之间总公共 最多去掉m-h[i]-L[i]个，i到T同理。

代码

#include <iostream>
#include <queue>
#include <cstdio>
#include <cstring>
const int N=507,inf=0x3f3f3f3f;
using namespace std;
int read() {
    int x=0,f=1;char s=getchar();
    for(;s>'9'||s<'0';s=getchar()) if(s=='-') f=-1;
    for(;s>='0'&&s<='9';s=getchar()) x=x*10+s-'0';
    return x*f;
}
int n,m,k,S,T,L[N],C[N],dis[N],mp[N][N],h[N],l[N];
struct node {
    int v,nxt,cap;
}e[2000007];
int head[2000007],tot=1;
void add(int u,int v,int cap) {
    e[++tot].v=v;
    e[tot].cap=cap;
    e[tot].nxt=head[u];
    head[u]=tot;
}
queue<int> q;
bool bfs() {
    memset(dis,-1,sizeof(dis));
    q.push(S);
    dis[S]=0;
    while(!q.empty()) {
        int u=q.front();
        q.pop();
        for(int i=head[u];i;i=e[i].nxt) {
            int v=e[i].v;
            if(dis[v]==-1&&e[i].cap) {
                dis[v]=dis[u]+1;
                q.push(v);
            }
        }
    }
    return dis[T]!=-1;
}
int dfs(int u,int f) {
    if(u==T) return f;
    int rest=f;
    for(int i=head[u];i;i=e[i].nxt) {
        int v=e[i].v;
        if(dis[v]==dis[u]+1&&e[i].cap&&rest) {
            int t=dfs(v,min(rest,e[i].cap));
            if(!t) dis[t]=0;
            e[i].cap-=t;
            e[i^1].cap+=t;
            rest-=t;
        }
    }
    return f-rest;
}
int dinic() {
    int ans=0;
    while(bfs()) ans+=dfs(S,inf);
    return ans;
}
int main() {
//  freopen("a.in","r",stdin);
    n=read(),m=read(),k=read();
    for(int i=1;i<=m;++i) L[i]=read();
    for(int i=1;i<=n;++i) C[i]=read();
    for(int i=1;i<=k;++i) {
        int x=read(),y=read();
        h[x]++,l[y]++;
        mp[x][y]=1;
    }
    for(int i=1;i<=n;++i) if(m-h[i]<L[i]) {puts("JIONG!");return 0;}
    for(int i=1;i<=m;++i) if(n-l[i]<C[i]) {puts("JIONG!");return 0;}
    for(int i=1;i<=n;++i)
        for(int j=1;j<=m;++j)
            if(!mp[i][j]) add(i,j+n,1),add(j+n,i,0);
    S=n+m+1,T=n+m+2;
    for(int i=1;i<=n;++i) add(S,i,m-h[i]-L[i]),add(i,S,0);
    for(int i=1;i<=m;++i) add(i+n,T,n-l[i]-C[i]),add(T,i+n,0);
    printf("%d
",n*m-dinic()-k);
    return 0;
}