#include <iostream>
#include <cstdio>
#include <vector>
#include <algorithm>
#include <stack>
#define NMAX 100005

using namespace std;

int n, m, viz[NMAX], vizt[NMAX], nrSol;
vector <int> g[NMAX], gt[NMAX], sol[NMAX];
stack <int> s;

void read()
{
    scanf("%d %d", &n, &m);
    for(int i=1; i<=m; ++i)
    {
        int x, y;
        scanf("%d %d", &x, &y);
        g[x].push_back(y);
        gt[y].push_back(x); /// graf transpus
    }
}

void dfs(int node)
{
    viz[node] = 1;
    vector <int>::iterator it;
    for(it = g[node].begin(); it != g[node].end(); ++it)
        if(viz[*it] == 0)
        {
            dfs(*it);
            s.push(*it);
        }
}

void dfsTrans(int node)
{
    vizt[node] = 1;
    vector <int>::iterator it;
    for(it = gt[node].begin(); it != gt[node].end(); ++it)
        if(vizt[*it] == 0)
        {
            sol[nrSol].push_back(*it);
            dfsTrans(*it);
        }
}

void solve()
{
    for(int i=1; i<=n; ++i)
        if(viz[i] == 0)
        {
            dfs(i);
            while(!s.empty())
            {
                while(!s.empty() && vizt[s.top()] == 1)
                    s.pop();
                if(!s.empty())
                {
                    int x = s.top();
                    s.pop();
                    sol[nrSol].push_back(x);
                    dfsTrans(x);
                    nrSol++;
                }
            }
        }
}

void print()
{
    printf("%d\n", nrSol);
    for(int i=0; i<nrSol; ++i)
    {
        vector <int>::iterator it;
        for(it = sol[i].begin(); it != sol[i].end(); ++it)
            printf("%d ", *it);
        printf("\n");
    }
}

int main()
{
    freopen("ctc.in", "r", stdin);
    freopen("ctc.out", "w", stdout);
    read();
    solve();
    print();
    return 0;
}