bzoj 2127

3894: 文理分科

Time Limit: 10 Sec   Memory Limit: 512 MB
Submit: 663   Solved: 385
[ Submit][ Status][ Discuss]

Description

 文理分科是一件很纠结的事情！（虽然看到这个题目的人肯定都没有纠

结过）

 小P所在的班级要进行文理分科。他的班级可以用一个n*m的矩阵进行

描述，每个格子代表一个同学的座位。每位同学必须从文科和理科中选择

一科。同学们在选择科目的时候会获得一个满意值。满意值按如下的方式

得到：

1．如果第i行第秒J的同学选择了文科，则他将获得art[i][j]的满意值，如

  果选择理科，将得到science[i][j]的满意值。

2．如果第i行第J列的同学选择了文科，并且他相邻（两个格子相邻当且

  仅当它们拥有一条相同的边）的同学全部选择了文科，则他会更开

  心，所以会增加same_art[i][j]的满意值。

3．如果第i行第j列的同学选择了理科，并且他相邻的同学全部选择了理

  科，则增加same_science[i]j[]的满意值。

  小P想知道，大家应该如何选择，才能使所有人的满意值之和最大。请

告诉他这个最大值。

Input

第一行为两个正整数：n，m

接下来n术m个整数，表示art[i][j]；

接下来n术m个整数．表示science[i][j]；

接下来n术m个整数，表示same_art[i][j]；

Output

输出为一个整数，表示最大的满意值之和

Sample Input

3 4
13 2 4 13
7 13 8 12
18 17 0 5

8 13 15 4
11 3 8 11
11 18 6 5

1 2 3 4
4 2 3 2
3 1 0 4

3 2 3 2
0 2 2 1
0 2 4 4

Sample Output

152

一如既往的网络流神题，我一开始按照bzoj2127来写，那真是写得我。。

后来感觉这不对啊，于是去看了写题解，ORZ，有更简单建图的方法，在2127的基础下，一下就看懂了，于是调了半天的点的编号。。 MDZZ

看来网络流还是要多刷

#include<cstdio>
#include<cstring>
const int N=105*105*10;
const int MAX=1<<30;
int n,m;

int ans=0;
struct qq
{
    int x,y,z;
    int last;
    int other;
}s[600001];
int num,last[N];
int init1 (int x,int y,int z)
{
    num++;
    s[num].x=x;s[num].y=y;s[num].z=z;
    s[num].last=last[x];
    last[x]=num;
    return num;
}
void init (int x,int y,int z1)
{
    int num1=init1(x,y,z1),num2=init1(y,x,0);
    s[num1].other=num2;s[num2].other=num1;
    return ;
}
int st1,ed1;
int a[105][105],b[105][105];
int Get (int x,int y)
{
	return (x-1)*m+y;
}
int dx[5]={1,-1,0,0,0};
int dy[5]={0,0,1,-1,0};
void Init ()
{
    num=0;memset(last,-1,sizeof(last));
    scanf("%d%d",&n,&m);st1=3*n*m+1;ed1=st1+1;
    for (int u=1;u<=n;u++)
    	for (int i=1;i<=m;i++)
    	{
    		int a;
    		scanf("%d",&a);
    		init(st1,Get(u,i),a);
    		ans=ans+a;
    	}
    for (int u=1;u<=n;u++)
    	for (int i=1;i<=m;i++)
    	{
    		int a;
    		scanf("%d",&a);
    		init(Get(u,i),ed1,a);
    		ans=ans+a;
    	}
    for (int u=1;u<=n;u++)
    	for (int i=1;i<=m;i++)
    	{
    		int a;
    		scanf("%d",&a);
    		ans+=a;
    		init(st1,Get(u,i)+n*m,a);
			for(int k=0;k<5;k++){
                int x=dx[k]+u;
                int y=dy[k]+i;
                if(x<1||x>n||y<1||y>m)continue;
                init(Get(u,i)+m*n,Get(x,y),MAX);
            }
		}
    for (int u=1;u<=n;u++)
    	for (int i=1;i<=m;i++)
    	{
    		int a;
    		scanf("%d",&a);
    		ans+=a;
    		init(Get(u,i)+n*m*2,ed1,a);
			for(int k=0;k<5;k++){
                int x=dx[k]+u;
                int y=dy[k]+i;
                if(x<1||x>n||y<1||y>m)continue;
                init(Get(x,y),Get(u,i)+m*n*2,MAX);
            }
		}
	/*printf("%d\n",ans);
	for (int u=1;u<=num;u++)
		printf("%d %d %d\n",s[u].x,s[u].y,s[u].z);*/
}
int q[N];
int h[N];
bool bt ()
{
    memset(h,-1,sizeof(h));
    int st=1,ed=2;
    h[st1]=1;q[st]=st1;
    while (st!=ed)
    {
        int x=q[st];
        for (int u=last[x];u!=-1;u=s[u].last)
        {
            int y=s[u].y;
            if (s[u].z>0&&h[y]==-1)
            {
                h[y]=h[x]+1;
                q[ed]=y;
                ed++;
            }
        }
        st++;
    }
    if (h[ed1]==-1) return false;
    return true;
}
int mymin (int x,int y)
{
    return x<y?x:y;
}
int find_max (int x,int f)
{
    if (x==ed1) return f;
    int s1=0;
    for (int u=last[x];u!=-1;u=s[u].last)
    {
        int y=s[u].y;
        if (h[y]==(h[x]+1)&&s[u].z>0&&s1<f)
        {
            int t=find_max(y,mymin(s[u].z,f-s1));
            s1+=t;
            s[u].z-=t;
            s[s[u].other].z+=t;
        }
    }
    if (s1==0) h[x]=0;
    return s1;
}
void solve ()
{
    while (bt()==true)
        ans=ans-find_max(st1,MAX);
    printf("%d\n",ans);       
}
int main()
{
    Init();
    solve();
    return 0;
}