生活随笔
收集整理的這篇文章主要介紹了
[清华集训2017]无限之环（网络流）
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
很妙的一道題
對于每個格子，它合法與否，只跟它上下左右的相鄰格子有關(guān)，所以可以想到黑白染色
（用 (i,j) 表示 i 行 j 列的格子，我把 (i+j) %2 == 0 的格子染成白色，把(i+j)%2 == 1 的格子染成黑色）
關(guān)鍵是怎么描述旋轉(zhuǎn)操作
我樸實的想法是給每個格子建四個點，每個點代表格子可通過旋轉(zhuǎn)達到的一種狀態(tài)，從格子的初始狀態(tài)向格子其它狀態(tài)連邊，然后把格子間可以匹配的狀態(tài)連起來，再然后……就沒有然后了
看了題解，真的被驚到了
還是把每個格子拆成4個節(jié)點，對應四個方向上的接口
從源點向白格的接口連流量上界1，費用0的邊
從黑格的接口向匯點連流量上界1，費用0的邊
從黑格的接口向可以匹配的白格的接口連流量上界1，費用0的邊
然后旋轉(zhuǎn)操作，就可以通過一種神奇方式描述出來：
AD O BC
把這個看成一個白格（黑格類似，只是連邊方向相反）
邊（u->v，f，w）代表從 u 到 v，流量上界 f，費用 w 的邊
下文描述的旋轉(zhuǎn)方向均為順時針
1.此格有1個接口：
A
|
D O BC
（A->B ，1，1)） 對應轉(zhuǎn)90度
（A->C， 1，2)） 對應轉(zhuǎn)180度
（A->D ，1，1)） 對應轉(zhuǎn)270度
2.此格有2個接口：
情況1：
A
|
D O —— BC
（A->C， 1，1)） 對應轉(zhuǎn)90度
（B->D ，1，1)） 對應轉(zhuǎn)270度
（A->C， 1，1)）+（B->D ，1，1)） 對應轉(zhuǎn)180度
情況2：
A
|
D O B
|C
不能轉(zhuǎn)，忽略
3.此格有3個接口：
A
|
D O —— B
|C
（A->D ，1，1)） 對應轉(zhuǎn)270度
（B->D， 1，2)） 對應轉(zhuǎn)180度
（C->D ，1，1)） 對應轉(zhuǎn)90度
4.此格有4個接口：
A
|
D —— O —— B
|C
轉(zhuǎn)了也沒區(qū)別，忽略
建完圖后跑最小費用最大流即可
如果我表達不清的話，這是樣例1的建圖，很丑，湊合著看吧
圖上標的數(shù)代表費用，沒標的邊默認費用為0，所有邊的流量上界默認為1
#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std
;
const int inf
=0x7fffffff;
const int N
=8100;
const int M
=200010;
struct Edge
{int u
,v
,f
,w
,nxt
;
}edge
[M
<<1];
int s
,t
,head
[N
],cnt
,maxflow
,mincost
,dis
[N
],inque
[N
],pre
[N
];
queue
<int> q
;
int n
,m
,id
[2005][2005][4],num
,gr
[2005][2005],tot
;
void add(int u
,int v
,int f
,int w
){edge
[cnt
].u
=u
;edge
[cnt
].v
=v
;edge
[cnt
].f
=f
;edge
[cnt
].w
=w
;edge
[cnt
].nxt
=head
[u
];head
[u
]=cnt
++;edge
[cnt
].u
=v
;edge
[cnt
].v
=u
;edge
[cnt
].f
=0;edge
[cnt
].w
=-w
;edge
[cnt
].nxt
=head
[v
];head
[v
]=cnt
++;
}
bool spfa(){memset(dis
,0x7f,sizeof(dis
));memset(inque
,0,sizeof(inque
));memset(pre
,-1,sizeof(pre
));dis
[s
]=0;q
.push(s
);inque
[s
]=1;while(!q
.empty()){int u
=q
.front();q
.pop();inque
[u
]=0;for(int i
=head
[u
];i
!=-1;i
=edge
[i
].nxt
){int v
=edge
[i
].v
;if(edge
[i
].f
>0&&dis
[v
]>dis
[u
]+edge
[i
].w
){dis
[v
]=dis
[u
]+edge
[i
].w
;pre
[v
]=i
;if(!inque
[v
]){q
.push(v
);inque
[v
]=1;}}}} if(pre
[t
]!=-1) return 1;return 0;
}
void EK(){int flow
;while(spfa()){flow
=inf
;int x
=pre
[t
];while(x
!=-1){flow
=min(edge
[x
].f
,flow
);x
=pre
[edge
[x
].u
];}x
=pre
[t
];while(x
!=-1){edge
[x
].f
-=flow
;edge
[x
^1].f
+=flow
;mincost
+=flow
*edge
[x
].w
;x
=pre
[edge
[x
].u
];}maxflow
+=flow
;}
}
int main(){memset(head
,-1,sizeof(head
));scanf("%d%d",&n
,&m
);for(int i
=1;i
<=n
;i
++)for(int j
=1;j
<=m
;j
++)for(int k
=0;k
<4;k
++)id
[i
][j
][k
]=++num
;s
=++num
,t
=++num
; for(int i
=1;i
<=n
;i
++){for(int j
=1;j
<=m
;j
++){scanf("%d",&gr
[i
][j
]);int tt
=0,pos
[5]={-1,-1,-1,-1,-1},npos
=-1;for(int k
=0;k
<4;k
++){if(gr
[i
][j
]&(1<<k
)){tot
++;pos
[++tt
]=k
;if((i
+j
)%2==0) add(s
,id
[i
][j
][k
],1,0);else add(id
[i
][j
][k
],t
,1,0);}else npos
=k
;}if(tt
==1){if((i
+j
)%2==0){add(id
[i
][j
][pos
[1]],id
[i
][j
][(pos
[1]+1)%4],1,1);add(id
[i
][j
][pos
[1]],id
[i
][j
][(pos
[1]+2)%4],1,2);add(id
[i
][j
][pos
[1]],id
[i
][j
][(pos
[1]+3)%4],1,1);}else{add(id
[i
][j
][(pos
[1]+1)%4],id
[i
][j
][pos
[1]],1,1);add(id
[i
][j
][(pos
[1]+2)%4],id
[i
][j
][pos
[1]],1,2);add(id
[i
][j
][(pos
[1]+3)%4],id
[i
][j
][pos
[1]],1,1);}}else if(tt
==2){if(pos
[2]-pos
[1]==2) continue;if((i
+j
)%2==0){add(id
[i
][j
][pos
[1]],id
[i
][j
][(pos
[1]+2)%4],1,1);add(id
[i
][j
][pos
[2]],id
[i
][j
][(pos
[2]+2)%4],1,1);}else{add(id
[i
][j
][(pos
[1]+2)%4],id
[i
][j
][pos
[1]],1,1);add(id
[i
][j
][(pos
[2]+2)%4],id
[i
][j
][pos
[2]],1,1);}}else if(tt
==3){if((i
+j
)%2==0){add(id
[i
][j
][(npos
+1)%4],id
[i
][j
][npos
],1,1);add(id
[i
][j
][(npos
+2)%4],id
[i
][j
][npos
],1,2);add(id
[i
][j
][(npos
+3)%4],id
[i
][j
][npos
],1,1);}else{add(id
[i
][j
][npos
],id
[i
][j
][(npos
+1)%4],1,1);add(id
[i
][j
][npos
],id
[i
][j
][(npos
+2)%4],1,2);add(id
[i
][j
][npos
],id
[i
][j
][(npos
+3)%4],1,1);}}else continue;}}int dx
[]={-1,0,1,0};int dy
[]={0,1,0,-1};for(int i
=1;i
<=n
;i
++){for(int j
=1;j
<=m
;j
++){if((i
+j
)%2==0){for(int k
=0;k
<=3;k
++){int x
=i
+dx
[k
],y
=j
+dy
[k
];if(x
<1||x
>n
||y
<1||y
>m
) continue;add(id
[i
][j
][k
],id
[x
][y
][(k
+2)%4],1,0);}}}}EK();if(tot
%2!=0||maxflow
!=tot
/2) printf("-1\n"); else printf("%d\n",mincost
);return 0;
}
總結(jié)
以上是生活随笔為你收集整理的[清华集训2017]无限之环（网络流）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。
