HDU 1565 1569 方格取數 (最小割)
阿新 • • 發佈:2019-01-26
連結:http://acm.hdu.edu.cn/showproblem.php?pid=1565
這兩題其實是同一題,一個數據小,一個數據大,據說資料小的時候能用狀壓DP。哇!好神奇呀,可是我不會。就用dinic搞了一下。
思路:假設選了(x,y)點,那麼(x,y)周圍的點都不能選了。那麼將圖分為兩部分,一部分是(x+y)%2==0,另一部分是(x+y)%2==1。所以在一個點和他周圍的點中,會做一個選擇,要麼選這個點,要麼選他周圍的某些點。建圖就可以按照這個思路建了,源點跟(x+y)是偶數的點連邊,流量是該點的權值。(x+y)是奇數的點跟匯點連邊,流量是權值。然後每個(x+y)是偶數的點,跟他周圍的點連邊,權值是INF。這樣就是一個最小割模型,可以用最大流求解。但是這個解出來的答案是最小割,是權值最小的選擇方案,怎麼辦呢?
選擇不是最小割中的點。如果選擇最小權值時沒有非法選擇的話,那麼選擇其餘所有點也不會非法,於是將所有權值加起來減去最小割的權值就是答案了。
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<cctype>
#include<string>
#include<iostream>
#include<vector>
#include<map>
#include<queue>
using namespace std;
#define PB push_back
typedef long long ll;
typedef pair<int,int> pii;
const int INF=1e9+7;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
struct EE{
int to,cap,next;
EE(){}
EE(int to,int cap,int next):to(to),cap(cap),next(next){}
}edge[100000];
int n;
int st,ed,Ecnt,head[410],d[410];
int go[4][2]={{1,0 },{-1,0},{0,1},{0,-1}};
inline int id(int x,int y){
return x*n+y+1;
}
inline bool judge(int x,int y){
return x>=0&&x<n&&y>=0&&y<n;
}
inline void add(int s,int t,int v){
edge[Ecnt]=EE(t,v,head[s]);
head[s]=Ecnt++;
edge[Ecnt]=EE(s,0,head[t]);
head[t]=Ecnt++;
}
bool BFS(){
memset(d,-1,sizeof d);
d[st]=0;
queue<int> Q;Q.push(st);
while(!Q.empty()){
int s=Q.front();Q.pop();
for(int i=head[s];~i;i=edge[i].next){
EE e=edge[i];
if(e.cap>0&&d[e.to]<0){
d[e.to]=d[s]+1;Q.push(e.to);
}
}
}
return d[ed]!=-1;
}
int DFS(int s,int t,int flow){
if(s==ed||!flow)return flow;
int ans=0;
for(int i=head[s];~i;i=edge[i].next){
EE &e=edge[i];
if(e.cap>0&&d[e.to]==d[s]+1){
int ff=DFS(e.to,t,min(flow,e.cap));
if(ff>0){
e.cap-=ff;
edge[i^1].cap+=ff;
flow-=ff;
ans+=ff;
if(!flow)break;
}
}
}
if(!ans)d[s]=-1;
return ans;
}
int dinic(){
int ans=0;
while(BFS()){
ans+=DFS(st,ed,INF);
}
return ans;
}
int main(){
while(scanf("%d",&n)!=EOF){
st=0;ed=n*n+1;
int sum=0;
Ecnt=0;memset(head,-1,sizeof head);
for(int i=0;i<n;i++){
for(int j=0;j<n;j++){
int x;scanf("%d",&x);
sum+=x;
if((i+j)&1){
add(id(i,j),ed,x);
}
else{
add(st,id(i,j),x);
for(int k=0;k<4;k++){
int x=i+go[k][0],y=j+go[k][1];
if(judge(x,y)){
add(id(i,j),id(x,y),INF);
}
}
}
}
}
printf("%d\n",sum-dinic());
}
return 0;
}