bzoj2127 happiness

题目

看看数据范围，貌似是网络流，so，我们来建建图。

把源点接向学理科的学生，容量为他学习理科的价值。
把学文科的学生接向汇点，容量为他学习文科的价值。
学生之间接双向边，容量为互相学习文理科价值的一半。

这样后，先把所有价值加起来，再减掉多余的。其实，我们需要把文理科分成两个集合，就是割，求个最小割。

#include<bits/stdc++.h>
#define N 10005
#define inf ~0u>>1
using namespace std;
int n,m,S,T,tot,x;
int first[N+1],dis[N+1],siz,id[105][105];
int A[105][105],B[105][105];
vector <int> nex,to,fl;
inline char nc()
{
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline int read()
{
    int x=0,b=1;
    char c=nc();
    for(;!(c<='9'&&c>='0');c=nc())if(c=='-')b=-1;
    for(;c<='9'&&c>='0';c=nc())x=x*10+c-'0';
    return x*b;
}
inline void add(int x,int y,int z)
{
    nex.push_back(first[x]);
    first[x]=siz;
    to.push_back(y);
    fl.push_back(z);
    siz++;
}
inline void insert(int x,int y,int z)
{
    add(x,y,z),add(y,x,0);
}
inline void insert2(int x,int y,int z)
{
    add(x,y,z),add(y,x,z);
}
inline bool BFS(int s,int t)
{
    queue <int> Q;
    for(int i=S;i<=T;i++)dis[i]=inf;
    Q.push(s),dis[s]=0;
    while(!Q.empty())
    {
        int x=Q.front();Q.pop();
        for(int i=first[x];i!=-1;i=nex[i])
        {
            int u=to[i];
            if(fl[i]&&dis[x]+1<dis[u])
            {
                dis[u]=dis[x]+1;
                if(u==t)return true;
                Q.push(u);
            }
        }
    }
    if(dis[t]==inf)return false;return true;
}
inline int DFS(int s,int f,int t)
{
    if(s==t)return f;
    int rest=f;
    for(int i=first[s];i!=-1;i=nex[i])
    {
        int u=to[i];
        if(rest&&fl[i]&&dis[u]==dis[s]+1)
        {
            int tmp=DFS(u,min(rest,fl[i]),t);
            if(tmp==0)dis[u]=0;
            fl[i]-=tmp;
            fl[i^1]+=tmp;
            rest-=tmp;
        }
    }
    return f-rest;  
}
inline int Dinic(int s,int t)
{
    int ans=0;
    while(BFS(s,t))ans+=DFS(s,inf,t);
    return ans;
}
inline void build()
{
    for(int i=1;i<n;i++)
    for(int j=1;j<=m;j++)
    {
        x=read();tot+=x;
        A[i][j]+=x,A[i+1][j]+=x;
        insert2(id[i+1][j],id[i][j],x);
    }
    for(int i=1;i<n;i++)
    for(int j=1;j<=m;j++)
    {
        x=read();tot+=x;
        B[i][j]+=x,B[i+1][j]+=x;
        insert2(id[i+1][j],id[i][j],x);
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<m;j++)
    {
        x=read(),tot+=x;
        A[i][j]+=x,A[i][j+1]+=x;
        insert2(id[i][j],id[i][j+1],x);
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<m;j++)
    {
        x=read(),tot+=x;
        B[i][j]+=x,B[i][j+1]+=x;
        insert2(id[i][j],id[i][j+1],x);
    }
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
        insert(S,id[i][j],A[i][j]),insert(id[i][j],T,B[i][j]);
}
int main()
{
    //freopen("in.txt","r",stdin);
    memset(first,-1,sizeof(first));
    n=read(),m=read();
    S=0,T=m*n+1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            A[i][j]=read(),tot+=A[i][j],A[i][j]<<=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            B[i][j]=read(),tot+=B[i][j],B[i][j]<<=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)    
            id[i][j]=(i-1)*m+j;
    build();
    cout<<tot-(Dinic(S,T)>>1);
    return 0;
}

日常语文不好。

看建图把，除以2可能会有小数，我们就把所有权值扩大两倍。