1. 程式人生 > >HDU 1565 1569 方格取數 (最小割)

HDU 1565 1569 方格取數 (最小割)

連結: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; }