#include <iostream>
#include <fstream>
#include <vector>

using namespace std;

vector<int> v, g[1001], _g[1001], sol[1001];
bool viz[1001]={0}, _viz[1001]={0};
int n, nr, nrc;

void dfs(int nod)
{
    viz[nod]=true;
    v.push_back(nod);
    for(int i=0; i<g[nod].size(); i++)
        if(!viz[g[nod][i]])
            dfs(g[nod][i]);
}

void dfs2(int nod)
{
    _viz[nod]=true;
    sol[nrc].push_back(nod);
    int l=_g[nod].size();
    for(int i=0; i<l; i++)
        if(viz[_g[nod][i]] && !_viz[_g[nod][i]])
            dfs2(_g[nod][i]);
}

int main()
{
    ifstream f("ctc.in");
    ofstream gg("ctc.out");
    int m, a, b, lg, i, j;
    f>>n>>m;
    for(i=1; i<=m; i++)
    {
        f>>a>>b;
        g[a].push_back(b);
        _g[b].push_back(a);
    }
    for(i=1; i<=n; i++)
        if(!viz[i])
        {
            dfs(i);
            lg=v.size();
            for(j=0; j<lg; j++)
            {
                if(!_viz[v[j]])
                {
                    nrc++;
                    dfs2(v[j]);
                }
            }
        }
    g<<nrc<<"\n";
    for(i=1; i<=nrc; i++)
    {
        lg=sol[i].size();
        for(j=0; j<lg; j++)
            gg<<sol[i][j]<<" ";
        gg<<"\n";
    }
}