题意
有一个n*n的01矩阵,有m个点为1,每一次可以把一行或一列的1全部变成0,问最少要消多少次。
n<=500,m<=10000
分析
用最小割来做。
若一个点(i,j)为1则若不消第i行则必消第j列,从i到j+n连一条inf的边。从s到1-n,从n+1-n*2连流量为1的边,那么最小割即为答案。
不想打dinic了,所以直接去blog上面复制了一个板子下来……
代码
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<queue>
#define N 1005
#define inf 0x3f3f3f3f
using namespace std;
int cnt,n,m,s,t,last[N],cur[N],dis[N],ans;
struct edge{int to,c,next;}e[N*1000];
queue <int> q;
void addedge(int u,int v,int c)
{
e[++cnt].to=v;e[cnt].c=c;e[cnt].next=last[u];last[u]=cnt;
e[++cnt].to=u;e[cnt].c=0;e[cnt].next=last[v];last[v]=cnt;
}
bool bfs()
{
memset(dis,0,sizeof(dis));
dis[s]=1;
while (!q.empty()) q.pop();
q.push(s);
while (!q.empty())
{
int u=q.front();
q.pop();
for (int i=last[u];i;i=e[i].next)
if (e[i].c&&!dis[e[i].to])
{
dis[e[i].to]=dis[u]+1;
if (e[i].to==t) return 1;
q.push(e[i].to);
}
}
return 0;
}
int dfs(int x,int maxf)
{
if (x==t||!maxf) return maxf;
int ret=0;
for (int &i=cur[x];i;i=e[i].next)
if (e[i].c&&dis[e[i].to]==dis[x]+1)
{
int f=dfs(e[i].to,min(maxf-ret,e[i].c));
e[i].c-=f;
e[i^1].c+=f;
ret+=f;
if (ret==maxf) break;
}
return ret;
}
void dinic()
{
while (bfs())
{
for (int i=s;i<=t;i++)
cur[i]=last[i];
ans+=dfs(s,inf);
}
}
int main()
{
scanf("%d%d",&n,&m);
s=0;t=n*2+1;cnt=1;
for (int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
addedge(x,y+n,inf);
}
for (int i=1;i<=n;i++)
{
addedge(s,i,1);
addedge(i+n,t,1);
}
dinic();
printf("%d",ans);
return 0;
}