题目:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=2414
一个 r*c 的矩阵,有n个方格上有怪物,你可以从边界上发射子弹,子弹会把你发射过去的那一行或者列的全部小怪都清掉。问你最少需要多少发子弹,可以全清掉,然后输出具体方案。来张图片吧:
SampleInput
4 4 3
1 1
1 4
3 2
4 4 2
1 1
2 2
0 0 0
SampleOutput
2 r1 r3
2 r1 r2
最小点覆盖,选出最少的点构成集合S,使得对于图G=(V,E)所有的边都至少有一个端点在集合S里。
本题可以把横纵坐标看做二分图的两边的点,怪物的位置就是联系两点的边信息,这样顺利的转化成最小点覆盖问题,也即最大匹配数。怎样寻找射弹位置?在左侧未覆盖的点拓展匈牙利,树上所有的点都标记,最后结果就是左边(Y轴)未覆盖的点和右边(X轴)覆盖了的点。例如:
数据:
4 4 6
4 1
4 2
4 3
1 3
3 3
2 3
void show(){
cout<<"show: "<<endl;
for(int i=r;i>=1;i--)cout<<S[i]<<endl;
cout<<" "; for(int i=1;i<=c;i++) cout<<T[i]<<" "; cout<<endl;
}
int main()
{
freopen("cin.txt","r",stdin);
while(scanf("%d%d%d",&r,&c,&n)&&(r+c+n))
{
e=0;
memset(head,-1,sizeof(head));
memset(linkl,0,sizeof(linkl));
memset(linkr,0,sizeof(linkr));
for(int i=0;i<n;i++)
{
int u,v;
scanf("%d%d",&u,&v);
addedge(u,v);
}
int ans=0;
for(int i=1;i<=r;i++)
{
show();
//memset(S,0,sizeof(S));
//memset(T,0,sizeof(T));
if(match(i)){
ans++;
}
}
show();
memset(S,0,sizeof(S));
memset(T,0,sizeof(T));
show();
for(int i=1;i<=r;i++)
{
if(!linkr[i])
match(i); //对没有匹配的点再次调用match,排除行边
cout<<"linker: "<<i<<endl;
show();
}
printf("%d",ans);
for(int i=1;i<=r;i++)
if(!S[i]) printf(" r%d",i);
for(int i=1;i<=c;i++)
if(T[i]) printf(" c%d",i);
puts("");
}
return 0;
}
恩,就是这样。
代码:
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std;
const int maxn=2e3+10;
const int maxm=1e6+10;
int e,r,c,n,head[maxn],edgeto[maxm],next[maxm],linkl[maxn],linkr[maxn];
bool S[maxn],T[maxn];
void addedge(int u,int v)
{
edgeto[e]=v;
next[e]=head[u];
head[u]=e++;
}
bool match(int u)
{
S[u]=1;
for(int i=head[u];i!=-1;i=next[i])
if(!T[edgeto[i]])
{
T[edgeto[i]]=1;
if(!linkl[edgeto[i]]||match(linkl[edgeto[i]]))
{
linkl[edgeto[i]]=u;
linkr[u]=edgeto[i];
return true;
}
}
return false;
}
int main()
{
while(scanf("%d%d%d",&r,&c,&n)&&(r+c+n))
{
e=0;
memset(head,-1,sizeof(head));
memset(linkl,0,sizeof(linkl));
memset(linkr,0,sizeof(linkr));
for(int i=0;i<n;i++)
{
int u,v;
scanf("%d%d",&u,&v);
addedge(u,v);
}
int ans=0;
for(int i=1;i<=r;i++)
{
memset(S,0,sizeof(S));
memset(T,0,sizeof(T));
if(match(i))
ans++;
}
memset(S,0,sizeof(S));
memset(T,0,sizeof(T));
for(int i=1;i<=r;i++)
if(!linkr[i]) match(i);
printf("%d",ans);
for(int i=1;i<=r;i++)
if(!S[i]) printf(" r%d",i);
for(int i=1;i<=c;i++)
if(T[i]) printf(" c%d",i);
puts("");
}
return 0;
}