#include <bits/stdc++.h>

using namespace std;

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

vector <int> v[100100], vT[100100],sol[100100];

int M,t,vectorT[100100],viz[100100];

void DFS_T(int nod)
{
    int vecin;
    viz[nod]=0;

    for(int i=0;i<vT[nod].size();i++)
    {
        vecin=vT[nod][i];
        if(viz[vecin]==1)
        {
            viz[vecin]=0;
            DFS_T(vecin);
        }
    }

    sol[M].push_back(nod);
}

void DFS(int nod)
{
    int vecin;

    viz[nod]=1;

    for(int i=0;i<v[nod].size();i++)
    {
        vecin=v[nod][i];
        if(viz[vecin]==0)
        {
            viz[vecin]=1;
            DFS(vecin);
        }
    }

    t++;
    vectorT[t]=nod;
}

int n,m,x,y,j,i;
int main()
{
    f>>n>>m;
    for(i=1;i<=m;i++)
    {
        f>>x>>y;
        v[x].push_back(y);
        vT[y].push_back(x);
    }

    for(i=1;i<=n;i++)
    {
        if(viz[i]==0)DFS(i);
    }



    for(i=n;i>=1;i--)
    {
        if(viz[vectorT[i]]==1)
        {
            M++;
            DFS_T(vectorT[i]);
        }
    }

    g<<M<<'\n';
    for(i=1;i<=M;i++)
    {
        for(j=0;j<sol[i].size();j++)
            g<<sol[i][j]<<" ";
        g<<'\n';
    }
    return 0;
}