BZOJ 2127: happiness

標籤：自己   接下來   分配   while   false   sample   i++   color   center   

2127: happinessTime Limit: 51 Sec  Memory Limit: 259 MB
Submit: 2084  Solved: 1028
Description

高一一班的座位表是個n*m的矩陣，經過一個學期的相處，每個同學和前後左右相鄰的同學互相成為了好朋友。這學期要分文理科了，每個同學對於選擇文科與理科有著自己的喜悅值，而一對好朋友如果能同時選文科或者理科，那麼他們又將收穫一些喜悅值。作為電腦競賽教練的scp大老闆，想知道如何分配可以使得全班的喜悅值總和最大。

Input

第一行兩個正整數n，m。接下來是六個矩陣第一個矩陣為n行m列 此矩陣的第i行第j列的數字表示座位在第i行第j列的同學選擇文科獲得的喜悅值。第二個矩陣為n行m列 此矩陣的第i行第j列的數字表示座位在第i行第j列的同學選擇理科獲得的喜悅值。第三個矩陣為n-1行m列 此矩陣的第i行第j列的數字表示座位在第i行第j列的同學與第i+1行第j列的同學同時選擇文科獲得的額外喜悅值。第四個矩陣為n-1行m列 此矩陣的第i行第j列的數字表示座位在第i行第j列的同學與第i+1行第j列的同學同時選擇理科獲得的額外喜悅值。第五個矩陣為n行m-1列 此矩陣的第i行第j列的數字表示座位在第i行第j列的同學與第i行第j+1列的同學同時選擇文科獲得的額外喜悅值。第六個矩陣為n行m-1列 此矩陣的第i行第j列的數字表示座位在第i行第j列的同學與第i行第j+1列的同學同時選擇理科獲得的額外喜悅值。

Output

輸出一個整數，表示喜悅值總和的最大值

Sample Input1 2
1 1
100 110
1
1000
Sample Output1210
【範例說明】
兩人都選理，則獲得100+110+1000的喜悅值。
【資料規模】
對於100%以內的資料，n,m<=100 所有喜悅值均為小於等於5000的非負整數
 首先我們來吐槽一下

什麼狗屁BZOJ，錢收的越來越多，事卻辦不利索，TLE了一晚上，一次比一次慢，最後一次是STD；

 

#include<iostream>#include<cstring>#include<cstdio>#include<queue>#include<algorithm>using namespace std;const int S=0;#define T (m*n+1)#define maxm 10100#define INF 30000000struct Edge{ int to,value,next; }e[1001001];int tot=1,head[maxm],ans,n,m;int a[110][110],b[110][110];void Add_Edge(int u,int v,int w){    e[++tot].to=v;e[tot].value=w;    e[tot].next=head[u];head[u]=tot;}queue<int> q;int pre(int x,int y){ return ((x*n)-n+y); }namespace Max_Flow{    int dep[maxm];    bool BFS(){        memset(dep,-1,sizeof dep );        while(!q.empty()) q.pop();        dep[S]=0;q.push(S);        while(!q.empty()){            int u=q.front();q.pop();            for(int i=head[u];i;i=e[i].next){                int v=e[i].to;                if(dep[v]==-1&&e[i].value){                    dep[v]=dep[u]+1; q.push(v);                    if(v==T) return true;                }            }        }        return false;    }    int Dinic(int u,int flow){        if(u==T) return flow;        int ret=flow;        for(int i=head[u];i;i=e[i].next){            int v=e[i].to;            if(dep[v]==dep[u]+1&&e[i].value){                int x=Dinic(v,min(e[i].value,ret));                ret-=x;e[i].value-=x;e[i^1].value+=x;            }        }        if(ret) dep[u]=-1;        return flow-ret;    }}int main(){    scanf("%d%d",&m,&n);    for(int i=1;i<=m;i++)        for(int j=1;j<=n;j++){            scanf("%d",&a[i][j]);            ans+=a[i][j]; a[i][j] <<= 1;        }    for(int i=1;i<=m;i++)        for(int j=1;j<=n;j++){            scanf("%d",&b[i][j]);            ans+=b[i][j]; b[i][j] <<= 1;        }    for(int i=1;i<m;i++)        for(int x,j=1;j<=n;j++){            scanf("%d",&x);ans+=x;            a[i][j]+=x; a[i+1][j]+=x;            int u=pre(i,j),v=pre(i+1,j);            Add_Edge(u,v,x);Add_Edge(v,u,x);        }    for(int i=1;i<m;i++)        for(int x,j=1;j<=n;j++){            scanf("%d",&x);ans+=x;            b[i][j]+=x; b[i+1][j]+=x;            int u=pre(i,j),v=pre(i+1,j);            Add_Edge(u,v,x);Add_Edge(v,u,x);        }    for(int i=1;i<=m;i++)        for(int x,j=1;j<n;j++){            scanf("%d",&x);ans+=x;            a[i][j]+=x; a[i][j+1]+=x;            int u=pre(i,j),v=pre(i,j+1);            Add_Edge(u,v,x);Add_Edge(v,u,x);        }    for(int i=1;i<=m;i++)        for(int x,j=1;j<n;j++){            scanf("%d",&x);ans+=x;            b[i][j]+=x; b[i][j+1]+=x;            int u=pre(i,j),v=pre(i,j+1);            Add_Edge(u,v,x);Add_Edge(v,u,x);        }    for(int i=1;i<=m;i++)        for(int j=1;j<=n;j++){            int P=pre(i,j);            Add_Edge(S,P,a[i][j]);            Add_Edge(P,S,a[i][j]);            Add_Edge(P,T,b[i][j]);            Add_Edge(T,P,b[i][j]);        }    using namespace Max_Flow;    ans <<= 1;    while(BFS()) ans-=Dinic(S,INF);    printf("%d\n",(ans>>1));    return 0;}

 

BZOJ 2127: happiness