#include <fstream>
#include <vector>
#include <bitset>
#include <string.h>
#define Nmax 100099
using namespace std;
ifstream f("ctc.in");
ofstream g("ctc.out");

int N,M,sol,aranjat[Nmax];
vector < int > Graph[Nmax],GT[Nmax],Comp[Nmax];
vector < bool > Plus,Minus;


inline void ReadInput()
{
    f>>N>>M;
    for(int i=1;i<=M;++i)
    {
        int x,y;
        f>>x>>y;
        Graph[x].push_back(y);
        GT[y].push_back(x);
    }
}

void DFS_PLUS(int nod)
{
    if(aranjat[nod])return;
    Plus[nod]=1;
    vector < int > :: iterator it;
    for(it=Graph[nod].begin();it!=Graph[nod].end();++it)
        if(!Plus[*it])DFS_PLUS(*it);
}

void DFS_MINUS(int nod)
{
    if(aranjat[nod])return;
    Minus[nod]=1;
    if(Plus[nod]==1)
    {
        Comp[sol].push_back(nod);
        aranjat[nod]=1;
    }
    else return;
    vector < int > :: iterator it;
    for(it=GT[nod].begin();it!=GT[nod].end();++it)
        if(!Minus[*it])DFS_MINUS(*it);
}

int main()
{
    ReadInput();
    for(int i=1;i<=N;++i)
        if(!aranjat[i])
        {

            Plus.resize(N+5);Minus.resize(N+5);
            //Plus.assign(0); Minus.assign(0);
            ++sol;
            DFS_PLUS(i);
            DFS_MINUS(i);
            Plus.clear();Minus.clear();
        }
    g<<sol<<'\n';
    for(int i=1;i<=sol;++i,g<<'\n')
        for(int j=0;j<Comp[i].size();++j)g<<Comp[i][j]<<' ';

    f.close();g.close();
    return 0;
}