Codeforces Round #286 (Div. 2) D. Mr. Kitayuta Technology （拓扑序，并查集，思维）

题目链接

题面：

题意：
给定n个点和m条单向边,边有传递性,即 a->b,b->c 可实现 a->c, 问满足 输入的m条单向边的实现 , 最少需要多少条单向边。
即最少需要多少条单向边可以满足给定有向图中点的可到达关系。

弱连通图：有向图的基图联通。
基图：有向图的所有有向边全部替换成无向边所形成的图。

题解：
若有向图的某个弱连通分量没有环，那么通过拓扑序可用 x-1 条边把这 x 个点连接起来使之符合要求。
若有向图的某个弱连通分量有环，那么至少需要 x 条边（使这x个点形成一个强连通分量，就是形成一个环）才能符合要求。不可能使用更少的边符合要求。

先使用并查集找出所有的弱联通分量，然后通过拓扑序找出强连通分量。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<string>
#include<queue>
#include<bitset>
#include<map>
#include<set>
#define ll long long
#define pr make_pair
#define pb push_back
#define ui unsigned int
#define lc (cnt<<1)
#define rc (cnt<<1|1)
#define tlen(x)  (t[(x)].r-t[(x)].l+1)
#define tmid (l+r)>>1
#define forhead(x) for(int i=head[(x)];i;i=nt[i])
#define max(x,y) ((x)>(y)?(x):(y))
#define min(x,y) ((x)>(y)?(y):(x))
using namespace std;
const int inf=0x3f3f3f3f;
const ll lnf=0x3f3f3f3f3f3f3f3f;
const double dnf=1e18;
const int mod=1000000007;
const double eps=1e-8;
const double pi=acos(-1.0);
const int maxn=100100;
const int maxm=100100;
const int up=100000;

int head[maxn],ver[maxn],edge[maxn],nt[maxn],d[maxn],tot=1;
int f[maxn],q[maxn];
bool ha[maxn];

void add(int x,int y,int z=0)
{
    ver[++tot]=y,edge[tot]=z;
    nt[tot]=head[x],head[x]=tot;
}

int fi(int x)
{
    if(f[x]!=x)
        f[x]=fi(f[x]);
    return f[x];
}

void mer(int x,int y)
{
    int xx=fi(x);
    int yy=fi(y);
    if(xx==yy) return;
    f[xx]=yy;
}

int main(void)
{
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        f[i]=i;
    int x,y;
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        add(x,y);
        mer(x,y);
        d[y]++;
    }
    int ans=n;
    int l=1,r=0;
    for(int i=1;i<=n;i++)
    {
        if(ha[fi(i)]==false) ans--;
        ha[fi(i)]=true;
        if(d[i]==0) q[++r]=i;
    }
    while(l<=r)
    {
        int x=q[l++];
        forhead(x)
        {
            int y=ver[i];
            if(--d[y]==0) q[++r]=y;
        }
    }
    for(int i=1;i<=n;i++)
    {
        if(d[i]==0) continue;
        if(ha[fi(i)]==true) ans++;
        ha[fi(i)]=false;
    }
    printf("%d\n",ans);
    return 0;
}