bzoj2127 happiness

来源:互联网 发布:淘宝面膜销量排行榜 编辑:程序博客网 时间:2024/05/21 11:06

题目

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

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

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

#include<bits/stdc++.h>#define N 10005#define inf ~0u>>1using 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可能会有小数,我们就把所有权值扩大两倍。