#include <bits/stdc++.h>

using namespace std;

int n,m,nr;
vector<int>G[100005],GT[100005],stiva,elem[100005];
bool fost[100005];

void dfs(int nod)
{
    fost[nod]=1;
    for(int i=0; i<G[nod].size(); i++)
    {
        int x=G[nod][i];
        if(!fost[x]) dfs(x);
    }

    stiva.push_back(nod);
}

void dfsT(int nod)
{
    fost[nod]=0;
    elem[nr].push_back(nod);


    for(int i=0; i<GT[nod].size(); i++)
    {
        int x=GT[nod][i];
        if(fost[x]) dfsT(x);
    }
}

int main()
{
    ifstream f("ctc.in");
    ofstream g("ctc.out");

    f>>n>>m;

    for(int i=1; i<=m; i++)
    {
        int x,y;
        f>>x>>y;
        G[x].push_back(y);
        GT[y].push_back(x);
    }

    for(int i=1; i<=n; i++)
        if(!fost[i]) dfs(i);

    for(int i=n-1; i>=0; i--)
    {
        int nod=stiva[i];
        if(fost[nod])
        {
            nr++;
            dfsT(nod);
        }
    }

    g<<nr<<'\n';
    for(int i=1; i<=nr; i++)
    {
        for(int j=0; j<elem[i].size(); j++)
            g<<elem[i][j]<<' ';
        g<<'\n';
    }

    f.close();
    g.close();

    return 0;
}