M_sea

洛谷4313 文理分科
LuoguBZOJ分析看做总收益减去最小割。源点表示文科,汇点表示理科。源点向$(i,j)$连容量为$art[i]...
扫描右侧二维码阅读全文
15
2019/01

洛谷4313 文理分科

Luogu

BZOJ

分析

看做总收益减去最小割。

源点表示文科,汇点表示理科。

源点向$(i,j)$连容量为$art[i][j]$的边,$(i,j)$向汇点连容量为$science[i][j]$的边。

每个点$(i,j)$建一个新点,然后源点向这个新点连容量为$same\_art[i][j]$的边,这个新点向$(i,j)$及与$(i,j)$相邻的点连容量为$\infty$的边。$same\_science[i][j]$的处理类似。

求最小割就行了。

代码

写得略长。

//It is made by M_sea
#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#include <queue>
#define re register
using namespace std;

inline int read() {
    int X=0,w=1; char c=getchar();
    while (c<'0'||c>'9') { if (c=='-') w=-1; c=getchar(); }
    while (c>='0'&&c<='9') X=X*10+c-'0',c=getchar();
    return X*w;
}

const int N=100+10;
const int INF=0x3f3f3f3f;
const int dx[4]={0,1,0,-1};
const int dy[4]={1,0,-1,0};

struct Edge { int v,w,nxt; };
Edge e[5000000];
int head[100000];

inline void addEdge(int u,int v,int w) {
    static int cnt=0;
    e[cnt].v=v,e[cnt].w=w,e[cnt].nxt=head[u],head[u]=cnt++;
    e[cnt].v=u,e[cnt].w=0,e[cnt].nxt=head[v],head[v]=cnt++;
}

int n,m,s,t;
int level[100000];

inline int bfs() {
    memset(level,0,sizeof(level));
    queue<int> Q; Q.push(s); level[s]=1;
    while (!Q.empty()) {
        int u=Q.front(); Q.pop();
        for (re int i=head[u];i!=-1;i=e[i].nxt) {
            int v=e[i].v;
            if (e[i].w&&!level[v])
                level[v]=level[u]+1,Q.push(v);
        }
    }
    return level[t]!=0;
}

inline int dfs(int u,int cpflow) {
    if (u==t) return cpflow;
    int addflow=0;
    for (re int i=head[u];i!=-1;i=e[i].nxt) {
        int v=e[i].v;
        if (e[i].w&&level[v]==level[u]+1) {
            int tmpadd=dfs(v,min(cpflow-addflow,e[i].w));
            e[i].w-=tmpadd,e[i^1].w+=tmpadd;
            addflow+=tmpadd;
        }
    }
    if (!addflow) level[u]=0;
    return addflow;
}

inline int dinic() {
    int maxflow=0;
    while (bfs()) maxflow+=dfs(s,INF);
    return maxflow;
}

int art[N][N],science[N][N];
int same_art[N][N],same_science[N][N];
int tot=0;
inline int g(int i,int j) { return (i-1)*m+j; }

inline void build() {
    memset(head,-1,sizeof(head));
    s=0,t=tot=n*m+1;
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j) {
            addEdge(s,g(i,j),art[i][j]);
            addEdge(g(i,j),t,science[i][j]);
    }
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j) {
            ++tot;
            addEdge(s,tot,same_art[i][j]);
            addEdge(tot,g(i,j),INF);
            for (re int k=0;k<4;++k) {
                int x=i+dx[k],y=j+dy[k];
                if (x<1||x>n||y<1||y>m) continue;   
                addEdge(tot,g(x,y),INF);
            }
        }
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j) {
            ++tot;
            addEdge(tot,t,same_science[i][j]);
            addEdge(g(i,j),tot,INF);
            for (re int k=0;k<4;++k) {
                int x=i+dx[k],y=j+dy[k];
                if (x<1||x>n||y<1||y>m) continue;
                addEdge(g(x,y),tot,INF);
            }
        }
}

int main() {
    n=read(),m=read(); int ans=0;
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j)
            ans+=(art[i][j]=read());
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j)
            ans+=(science[i][j]=read());
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j)
            ans+=(same_art[i][j]=read());
    for (re int i=1;i<=n;++i)
        for (re int j=1;j<=m;++j)
            ans+=(same_science[i][j]=read());
    build(); ans-=dinic();
    printf("%d\n",ans);
    return 0;
}
最后修改:2019 年 05 月 26 日 06 : 34 PM

发表评论