2127: happiness

题目链接

题目大意：给定一个座位图，相邻两人之间是朋友，每个人选择学文或学理会有相应的喜悦值，一对朋友同时选择学文/学理也会有相应的喜悦值，求喜悦值之和最大的方案

题解：Orz神奇的解不定方程建图

我的收获：跪啊

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

const int M=10005;
#define INF 0x7fffffff
#define rep(x,y) for(int i=1;i<=x;i++) for(int j=1;j<=y;j++)

bool Exit;
int n,m,t,st,ed,cnt,sum,x;
int a[105][105],b[105][105],id[105][105];
int num[M],d[M],head[M],last[M]; 

struct edge{int to,nex,c;}e[300005];

void add(int i,int j,int w){e[t].to=j,e[t].c=w,e[t].nex=head[i],last[i]=head[i]=t++;}
void ins(int x,int y,int z){add(x,y,z),add(y,x,0);}
void ins2(int x,int y,int z){add(x,y,z),add(y,x,z);}

int dfs(int x,int in)
{
    if(x==ed) return in;
    int ans=0,t;
    for(int i=last[x];i!=-1;last[x]=i=e[i].nex)
    {
        int v=e[i].to;
        if(e[i].c&&d[v]==d[x]-1){
            t=dfs(v,min(in-ans,e[i].c));
            ans+=t;
            e[i].c-=t; e[i^1].c+=t;
            if(Exit||ans==in) return ans;
        }
    }
    if(--num[d[x]]==0) Exit=1;
    d[x]++,num[d[x]]++,last[x]=head[x];
    return ans;
}

int ISAP()
{
    Exit=0;int flow=0;
    while(!Exit) flow+=dfs(st,INF);
    return flow;
}

void work(){
    cout<<sum-ISAP()/2<<endl;
}

void build()
{
    rep(n-1,m) scanf("%d",&x),sum+=x,a[i][j]+=x,a[i+1][j]+=x,ins2(id[i][j],id[i+1][j],x);
    rep(n-1,m) scanf("%d",&x),sum+=x,b[i][j]+=x,b[i+1][j]+=x,ins2(id[i][j],id[i+1][j],x);
    rep(n,m-1) scanf("%d",&x),sum+=x,a[i][j]+=x,a[i][j+1]+=x,ins2(id[i][j],id[i][j+1],x);
    rep(n,m-1) scanf("%d",&x),sum+=x,b[i][j]+=x,b[i][j+1]+=x,ins2(id[i][j],id[i][j+1],x);
    rep(n,m) ins(st,id[i][j],a[i][j]),ins(id[i][j],ed,b[i][j]);//将起点&终点相同的边合并 
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(last,-1,sizeof(last));
    cin>>n>>m;
    rep(n,m) id[i][j]=++cnt;
    rep(n,m) scanf("%d",&a[i][j]),sum+=a[i][j],a[i][j]<<=1;
    rep(n,m) scanf("%d",&b[i][j]),sum+=b[i][j],b[i][j]<<=1;
    st=0,ed=cnt+1;num[0]=ed+1;
    build();
}

int main()
{
    init();
    work();
    return 0;
}