题目大意:给定一个座位图,相邻两人之间是朋友,每个人选择学文或学理会有相应的喜悦值,一对朋友同时选择学文/学理也会有相应的喜悦值,求喜悦值之和最大的方案
这个题的模型显然是最小割- -
看到矩阵上相邻点之间的关系 很容易想到黑白染色 随后就能想到将某种颜色的点源汇对调
但是很可惜我没建出图来- - 自己YY了一种做法 但是喜闻乐见地发现建图方法不对- -
还是OTZ HZWER一下吧- - http://hzwer.com/2422.html
我不赘述了- -
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define M 10100
#define S 0
#define T (m*n+1)
#define INF 0x3f3f3f3f
#define P(i,j) ((i)*n-n+(j))
using namespace std;
struct abcd{
int to,f,next;
}table[1001001];
int head[M],tot=1;
int m,n,ans;
int a[110][110],b[110][100];
void Add(int x,int y,int z)
{
table[++tot].to=y;
table[tot].f=z;
table[tot].next=head[x];
head[x]=tot;
}
void Link(int x,int y,int z)
{
Add(x,y,z);
Add(y,x,z);
}
namespace Max_Flow{
int dpt[M];
bool BFS()
{
static int q[M];
int i,r=0,h=0;
memset(dpt,-1,sizeof dpt);
dpt[S]=1;q[++r]=S;
while(r!=h)
{
int x=q[++h];
for(i=head[x];i;i=table[i].next)
if(table[i].f&&!~dpt[table[i].to])
{
dpt[table[i].to]=dpt[x]+1;
q[++r]=table[i].to;
if(table[i].to==T)
return true;
}
}
return false;
}
int Dinic(int x,int flow)
{
int i,left=flow;
if(x==T) return flow;
for(i=head[x];i&&left;i=table[i].next)
if(table[i].f&&dpt[table[i].to]==dpt[x]+1)
{
int temp=Dinic(table[i].to,min(left,table[i].f) );
left-=temp;
table[i].f-=temp;
table[i^1].f+=temp;
}
if(left) dpt[x]=-1;
return flow-left;
}
}
int main()
{
int i,j,x;
cin>>m>>n;
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
{
scanf("%d",&a[i][j]);
ans+=a[i][j];
a[i][j]<<=1;
}
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
{
scanf("%d",&b[i][j]);
ans+=b[i][j];
b[i][j]<<=1;
}
for(i=1;i<m;i++)
for(j=1;j<=n;j++)
{
scanf("%d",&x);ans+=x;
a[i][j]+=x;a[i+1][j]+=x;
Link(P(i,j),P(i+1,j),x);
}
for(i=1;i<m;i++)
for(j=1;j<=n;j++)
{
scanf("%d",&x);ans+=x;
b[i][j]+=x;b[i+1][j]+=x;
Link(P(i,j),P(i+1,j),x);
}
for(i=1;i<=m;i++)
for(j=1;j<n;j++)
{
scanf("%d",&x);ans+=x;
a[i][j]+=x;a[i][j+1]+=x;
Link(P(i,j),P(i,j+1),x);
}
for(i=1;i<=m;i++)
for(j=1;j<n;j++)
{
scanf("%d",&x);ans+=x;
b[i][j]+=x;b[i][j+1]+=x;
Link(P(i,j),P(i,j+1),x);
}
for(i=1;i<=m;i++)
for(j=1;j<=n;j++)
Link(S,P(i,j),a[i][j]),Link(P(i,j),T,b[i][j]);
using namespace Max_Flow;
ans<<=1;
while( BFS() )
ans-=Dinic(S,INF);
cout<<(ans>>1)<<endl;
return 0;
}