显然这是一个最小路径覆盖
转化为二分图做
先拆点,把每个点拆成
x1
和
x2
这样
如果
x
到
然后跑一遍二分图最大匹配
匈牙利或网络流都行
我匈牙利忘记了只好dinic
最后答案就是原图的点减二分图最大匹配
#include<bits/stdc++.h>
#define g getchar()
#define T n*2+1
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
inline ll read(){
ll x=0,f=1;char ch=g;
for(;ch<'0'||ch>'9';ch=g)if(ch=='-')f=-1;
for(;ch>='0'&&ch<='9';ch=g)x=x*10+ch-'0';
return x*f;
}
inline void out(ll x){
int a[25],t=0;
if(x<0)putchar('-'),x=-x;
for(;x;x/=10)a[++t]=x%10;
for(int i=t;i;--i)putchar('0'+a[i]);
if(t==0)putchar('0');
putchar('\n');
}
struct re{int v,fl,next;}ed[100001];
int head[505],dep[505],dui[505],mp[201][201];
int e=1,n,ans,m;
inline void ins(int x,int y){
ed[++e]=(re){y,1,head[x]};head[x]=e;
ed[++e]=(re){x,0,head[y]};head[y]=e;
}
bool bfs(){
int tou=1,wei=1;
memset(dep,-1,sizeof(dep));
dui[1]=0,dep[0]=0;
for(;tou<=wei;++tou){
int u=dui[tou];
for(int i=head[u];i;i=ed[i].next)
if(ed[i].fl&&dep[ed[i].v]==-1)
dep[ed[i].v]=dep[u]+1,dui[++wei]=ed[i].v;
}
return dep[T]!=-1;
}
int dfs(int x,int fl){
if(x==T)return fl;
int used=0;
for(int i=head[x];i;i=ed[i].next)
if(ed[i].fl&&dep[ed[i].v]==dep[x]+1){
int w=dfs(ed[i].v,min(fl-used,ed[i].fl));
used+=w;ed[i].fl-=w;ed[i^1].fl+=w;
if(used==fl)return fl;
}
if(!used)dep[x]=-1;
return used;
}
void dinic(){for(;bfs();)ans+=dfs(0,inf);}
int main(){
n=read();m=read();
for(int i=1;i<=m;++i){
int x=read(),y=read();
mp[x][y]=1;
}
for(int i=1;i<=n;++i)mp[i][i]=1;
for(int i=1;i<=n;++i)
for(int j=1;j<=n;++j){
for(int k=1;k<=n;++k)
if(mp[i][j])continue;else mp[i][j]=mp[i][k]&&mp[k][j];
}
for(int i=1;i<=n;++i){
for(int j=1;j<=n;++j)
if(i!=j&&mp[i][j]){
ins(i,j+n);
}
ins(0,i);
ins(i+n,T);
}
dinic();
out(n-ans);
return 0;
}