bzoj2132: 圈地计划（最小割）

Posted 2021-01-03 bztminamoto

传送门

 

看来以后见到矩形就要黑白染色冷静一下了……

首先，如果它的要求时候相邻的选择相同，那么就是和这一题一样了->这里

然后考虑不同的要怎么做

那就把矩形黑白染色一下吧

然后令其中一种颜色的A，B反过来，那么就和上面那道题一样了

 1 //minamoto
 2 #include<iostream>
 3 #include<cstdio>
 4 #include<cstring>
 5 #include<queue>
 6 #define inf 0x3f3f3f3f
 7 #define id(i,j) ((i-1)*m+j)
 8 using namespace std;
 9 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
10 char buf[1<<21],*p1=buf,*p2=buf;
11 inline int read(){
12     #define num ch-‘0‘
13     char ch;bool flag=0;int res;
14     while(!isdigit(ch=getc()))
15     (ch==‘-‘)&&(flag=true);
16     for(res=num;isdigit(ch=getc());res=res*10+num);
17     (flag)&&(res=-res);
18     #undef num
19     return res;
20 }
21 const int N=10005,M=100005;
22 int head[N],Next[M],ver[M],edge[M],tot=1;
23 int dep[N],cur[N],c[105][105],S,T,n,m,sum;
24 queue<int> q;int dx[]={1,-1,0,0},dy[]={0,0,1,-1};
25 inline void add(int u,int v,int e){
26     ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e;
27     ver[++tot]=u,Next[tot]=head[v],head[v]=tot,edge[tot]=0;
28 }
29 bool bfs(){
30     while(!q.empty()) q.pop();
31     for(int i=S;i<=T;++i) cur[i]=head[i];
32     memset(dep,-1,sizeof(dep));
33     q.push(S),dep[S]=0;
34     while(!q.empty()){
35         int u=q.front();q.pop();
36         for(int i=head[u];i;i=Next[i]){
37             int v=ver[i];
38             if(dep[v]<0&&edge[i]){
39                 dep[v]=dep[u]+1,q.push(v);
40                 if(v==T) return true;
41             }
42         }
43     }
44     return false;
45 }
46 int dfs(int u,int limit){
47     if(u==T||!limit) return limit;
48     int flow=0,f;
49     for(int i=cur[u];i;cur[u]=i=Next[i]){
50         int v=ver[i];
51         if(dep[v]==dep[u]+1&&(f=dfs(v,min(limit,edge[i])))){
52             flow+=f,limit-=f;
53             edge[i]-=f,edge[i^1]+=f;
54             if(!limit) break;
55         }
56     }
57     if(!flow) dep[u]=-1;
58     return flow;
59 }
60 int dinic(){
61     int flow=0;
62     while(bfs()) flow+=dfs(S,inf);
63     return flow;
64 }
65 int main(){
66     //freopen("testdata.in","r",stdin);
67     n=read(),m=read(),S=0,T=n*m+1;
68     for(int i=1;i<=n;++i)
69     for(int j=1;j<=m;++j){
70         int k=read();sum+=k;
71         (i+j&1)?add(S,id(i,j),k):add(id(i,j),T,k);
72     }
73     for(int i=1;i<=n;++i)
74     for(int j=1;j<=m;++j){
75         int k=read();sum+=k;
76         (i+j&1)?add(id(i,j),T,k):add(S,id(i,j),k);
77     }
78     for(int i=1;i<=n;++i)
79     for(int j=1;j<=m;++j) c[i][j]=read();
80     for(int i=1;i<=n;++i)
81     for(int j=1;j<=m;++j){
82         if(i+j&1) continue;
83         for(int k=0;k<4;++k){
84             int x=i+dx[k],y=j+dy[k];
85             if(x<1||x>n||y<1||y>m) continue;
86             add(id(i,j),id(x,y),c[i][j]+c[x][y]);
87             add(id(x,y),id(i,j),c[i][j]+c[x][y]);
88             sum+=c[i][j]+c[x][y];
89         }
90     }
91     printf("%d
",sum-dinic());
92     return 0;
93 }