2127: happiness

Time 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 Input

1 2
1 1
100 110
1
1000

Sample Output

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

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

?

#include<cstdio> #include<iostream> #include<cmath> #include<cstring> #include<algorithm> #define maxn 10005 #define maxm 550000 #define inf 1061109567 using namespace std; char ch; bool ok; void read(int &x){for (ok=0,ch=getchar();!isdigit(ch);ch=getchar()) if (ch=='-') ok=1;for (x=0;isdigit(ch);x=x*10+ch-'0',ch=getchar());if (ok) x=-x; } int n,m,a[105][105],b[105][105],c[105][105],d[105][105],ans; struct flow{int s,t,tot,now[maxn],son[maxm],pre[maxm],val[maxm];int dis[maxn],head,tail,list[maxn];bool bo[maxn];void init(){s=0,t=n*m+1,tot=1,memset(now,0,sizeof(now));}void put(int a,int b,int c){pre[++tot]=now[a],now[a]=tot,son[tot]=b,val[tot]=c;}void add(int a,int b,int c){put(a,b,c),put(b,a,0);}bool bfs(){memset(bo,0,sizeof(bo));head=0,tail=1,list[1]=s,dis[s]=0,bo[s]=1;while (head<tail){int u=list[++head];for (int p=now[u],v=son[p];p;p=pre[p],v=son[p])if (val[p]&&!bo[v]) bo[v]=1,dis[v]=dis[u]+1,list[++tail]=v;}return bo[t];}int dfs(int u,int rest){if (u==t) return rest;int ans=0;for (int p=now[u],v=son[p];p&&rest;p=pre[p],v=son[p])if (val[p]&&dis[v]==dis[u]+1){int d=dfs(v,min(rest,val[p]));val[p]-=d,val[p^1]+=d,ans+=d,rest-=d;}if (!ans) dis[u]=-1;return ans;}int dinic(){int ans=0;while (bfs()) ans+=dfs(s,inf);return ans;} }f; inline int calc(int a,int b){return (a-1)*m+b;} int main(){read(n),read(m),f.init();for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) read(a[i][j]);for (int i=1;i<=n;i++) for (int j=1;j<=m;j++) read(b[i][j]);for (int i=1;i<=n;i++) for (int j=1;j<=m;j++){int u=calc(i,j);f.add(f.s,u,b[i][j]<<1),f.add(u,f.t,a[i][j]<<1);ans+=a[i][j]+b[i][j];}for (int i=1;i<n;i++) for (int j=1;j<=m;j++) read(c[i][j]);for (int i=1;i<n;i++) for (int j=1;j<=m;j++) read(d[i][j]);for (int i=1;i<n;i++) for (int j=1;j<=m;j++){int u=calc(i,j),v=calc(i+1,j);f.add(f.s,u,d[i][j]),f.add(f.s,v,d[i][j]);f.add(u,f.t,c[i][j]),f.add(v,f.t,c[i][j]);f.add(u,v,c[i][j]+d[i][j]),f.add(v,u,c[i][j]+d[i][j]);ans+=c[i][j]+d[i][j];}for (int i=1;i<=n;i++) for (int j=1;j<m;j++) read(c[i][j]);for (int i=1;i<=n;i++) for (int j=1;j<m;j++) read(d[i][j]);for (int i=1;i<=n;i++) for (int j=1;j<m;j++){int u=calc(i,j),v=calc(i,j+1);f.add(f.s,u,d[i][j]),f.add(f.s,v,d[i][j]);f.add(u,f.t,c[i][j]),f.add(v,f.t,c[i][j]);f.add(u,v,c[i][j]+d[i][j]),f.add(v,u,c[i][j]+d[i][j]);ans+=c[i][j]+d[i][j];}printf("%d\n",ans-(f.dinic()>>1));return 0; }

?

?

轉載于:https://www.cnblogs.com/suishiguang/p/6827509.html

總結

以上是生活随笔為你收集整理的BZOJ 2127: happiness的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。