把入侵者看作边,每一行每一列都是点,选取某一行某一列都有费用,这样问题就是选总权最小的点集覆盖所有边,就是最小点权覆盖。
此外,题目的总花费是所有费用的乘积,这时有个技巧,就是取对数,把乘法变为加法运算,最后再还原。
另外还可以从最小割的思路去这么理解:
每一行与源点相连,容量为该行的花费;每一列与汇点相连,容量为该列的花费;对于每个入侵者的坐标,该行该列连接一条容量INF的边。
要让源点汇点不连通,割边集必然与所有入侵者的行或列相关,而这样建模后的最小割就是最小的花费(容量INF的边必然不是最小割的一部分,其余的必然会选择某行或某列)。
1 #include2 #include 3 #include 4 #include 5 #include 6 using namespace std; 7 #define INF (1<<30) 8 #define MAXN 111 9 #define MAXM 1111 10 11 struct Edge{ 12 int v,cap,flow,next; 13 }edge[MAXM]; 14 int vs,vt,NE,NV; 15 int head[MAXN]; 16 17 void addEdge(int u,int v,int cap){ 18 edge[NE].v=v; edge[NE].cap=cap; edge[NE].flow=0; 19 edge[NE].next=head[u]; head[u]=NE++; 20 edge[NE].v=u; edge[NE].cap=0; edge[NE].flow=0; 21 edge[NE].next=head[v]; head[v]=NE++; 22 } 23 24 int level[MAXN]; 25 int gap[MAXN]; 26 void bfs(){ 27 memset(level,-1,sizeof(level)); 28 memset(gap,0,sizeof(gap)); 29 level[vt]=0; 30 gap[level[vt]]++; 31 queue que; 32 que.push(vt); 33 while(!que.empty()){ 34 int u=que.front(); que.pop(); 35 for(int i=head[u]; i!=-1; i=edge[i].next){ 36 int v=edge[i].v; 37 if(level[v]!=-1) continue; 38 level[v]=level[u]+1; 39 gap[level[v]]++; 40 que.push(v); 41 } 42 } 43 } 44 45 int pre[MAXN]; 46 int cur[MAXN]; 47 int ISAP(){ 48 bfs(); 49 memset(pre,-1,sizeof(pre)); 50 memcpy(cur,head,sizeof(head)); 51 int u=pre[vs]=vs,flow=0,aug=INF; 52 gap[0]=NV; 53 while(level[vs]