Problem
- 方格上有一些给定点,每次可以覆盖一行或者覆盖一列,求覆盖所有给定点的最小次数。
Solution
- 将行作为左边的点,列作为右边的点,对于每个给定点行向列连边,二分图最大匹配即为答案。
- 合法性:最大匹配的基础上所有的点一定会被覆盖。反证如果存在这样的点,所求一定不是最大匹配。
- 最优性:最大匹配一定是最优的。等价于在最大匹配合法性的基础上证明比最大匹配小的方案一定不合法:比最大匹配小,必然存在<u,v>没有选择,在原问题中等价于没有覆盖到所有点。
- 还是row左col右的建图。
- 多了个输出方案,在这卡了半天,最后发现一开始想的本质上是没问题的,必然性处理上不好。
- 通过dfs遍历残余网络,最后方案是左边没被遍历到的部分和右边遍历到的部分。
#include<iostream>
#include<queue>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int inf=0x7f7f7f7f;
long long n,m,s,t,cnt=1,last[2005],x,y,z,cur[2005],ceng[2005],k;
bool inq[2005],vis[2005];
struct edge{
long long v,w,next;
}e[2100000];
inline void add(int u,int v,int w)
{
cnt++;
e[cnt].w=w;
e[cnt].v=v;
e[cnt].next=last[u];
last[u]=cnt;
}
bool bfs()
{
queue<int>q;
memset(ceng,0,sizeof(ceng));
ceng[s]=1;
q.push(s);
while(!q.empty())
{
int u=q.front();q.pop();
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v,w=e[i].w;
if(ceng[v]==0&&w>0)
{
ceng[v]=ceng[u]+1;
q.push(v);
}
}
}
return ceng[t];
}
long long dfs(int u,long long dis)
{
if(u==t)
return dis;
for(int i=cur[u];i;i=e[i].next)
{
cur[u]=e[i].next;
long long v=e[i].v,w=e[i].w;
if(ceng[v]==ceng[u]+1&&w!=0)
{
long long di=dfs(v,min(dis,w));
if(di>0)
{
e[i].w-=di;
e[i^1].w+=di;
return di;
}
}
}
return 0;
}
long long dinic()
{
long long ans=0;
while(bfs())
{
for(int i=0;i<=t;i++)
cur[i]=last[i];
while(long long d=dfs(s,inf))
ans+=d;
}
return ans;
}
void init()
{
memset(inq,0,sizeof(inq));
memset(last,0,sizeof(last));
memset(vis,0,sizeof(vis));
cnt=1;
}
void find(int u)
{
vis[u]=1;
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v,w=e[i].w;
if(w&&!vis[v])
find(v);
}
}
int main()
{
while(scanf("%lld%lld%lld",&n,&m,&k)&&n)
{
s=0;t=n+m+1;
for(int i=1,x,y;i<=k;i++)
{
scanf("%d%d",&x,&y);
inq[x]=inq[n+y]=1;
add(x,y+n,1);
add(y+n,x,0);
}
for(int i=1;i<=n;i++)
{
if(!inq[i]) continue;
add(s,i,1);
add(i,s,0);
}
for(int i=1;i<=m;i++)
{
if(!inq[i+n]) continue;
add(i+n,t,1);
add(t,i+n,0);
}
cout<<dinic();
find(s);
for(int i=1;i<=n;i++)
if(!vis[i]&&inq[i])
cout<<" r"<<i;
for(int i=n+1;i<=n+m;i++)
if(vis[i]&&inq[i])
cout<<" c"<<i-n;
cout<<endl;
init();
}
return 0;
}
- 多了障碍物,重新划分行和列即可。
- 有了最小顶点覆盖的模型这题就好做。
#include<iostream>
#include<queue>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
const int inf=0x7f7f7f7f;
long long n,m,s,t,cnt=1,last[20005],x,y,z,cur[20005],ceng[20005],k1,k2,ansn,ansm,tt;
int g[105][105],numr[105][105],numc[105][105];
struct edge{
long long v,w,next;
}e[2100000];
inline void add(int u,int v,int w)
{
cnt++;
e[cnt].w=w;
e[cnt].v=v;
e[cnt].next=last[u];
last[u]=cnt;
}
bool bfs()
{
queue<int>q;
memset(ceng,0,sizeof(ceng));
ceng[s]=1;
q.push(s);
while(!q.empty())
{
int u=q.front();q.pop();
for(int i=last[u];i;i=e[i].next)
{
int v=e[i].v,w=e[i].w;
if(ceng[v]==0&&w>0)
{
ceng[v]=ceng[u]+1;
q.push(v);
}
}
}
return ceng[t];
}
long long dfs(int u,long long dis)
{
if(u==t)
return dis;
for(int i=cur[u];i;i=e[i].next)
{
cur[u]=e[i].next;
long long v=e[i].v,w=e[i].w;
if(ceng[v]==ceng[u]+1&&w!=0)
{
long long di=dfs(v,min(dis,w));
if(di>0)
{
e[i].w-=di;
e[i^1].w+=di;
return di;
}
}
}
return 0;
}
long long dinic()
{
long long ans=0;
while(bfs())
{
for(int i=0;i<=t;i++)
cur[i]=last[i];
while(long long d=dfs(s,inf))
ans+=d;
}
return ans;
}
void init()
{
memset(last,0,sizeof(last));
memset(g,0,sizeof(g));
memset(numr,0,sizeof(numr));
memset(numc,0,sizeof(numc));
cnt=1;
ansn=ansm=0;
}
int main()
{
cin>>tt;
while(tt--)
{
scanf("%lld%lld%lld",&n,&m,&k1);
for(int i=1,x,y;i<=k1;i++)
{
scanf("%d%d",&x,&y);
g[x][y]=1;
}
cin>>k2;
for(int i=1,x,y;i<=k2;i++)
{
scanf("%d%d",&x,&y);
g[x][y]=2;
}
for(int i=1;i<=n;i++)
{
ansn++;
for(int j=1;j<=m;j++)
{
if(g[i][j]==2)
ansn++;
if(g[i][j]==1)
numr[i][j]=ansn;
}
}
for(int i=1;i<=m;i++)
{
ansm++;
for(int j=1;j<=n;j++)
{
if(g[j][i]==2)
ansm++;
if(g[j][i]==1)
numc[j][i]=ansm;
}
}
s=0;t=ansn+ansm+1;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=m;j++)
{
if(g[i][j]==1)
{
add(numr[i][j],numc[i][j]+ansn,1);
add(numc[i][j]+ansn,numr[i][j],0);
}
}
}
for(int i=1;i<=ansn;i++)
{
add(s,i,1);
add(i,s,0);
}
for(int i=1;i<=ansm;i++)
{
add(i+ansn,t,1);
add(t,i+ansn,0);
}
cout<<dinic()<<endl;
init();
}
return 0;
}