#include <fstream>
#include <vector>
#include <algorithm>

using namespace std;

ifstream cin ("ctc.in");
ofstream cout ("ctc.out");

const int N = 1e5;

bool viz[N + 1];

vector <int> g[N + 1], gt[N + 1], v, sol[N + 1];

int n, m, comp, x, y;

void dfs_plus (int node)
{
    viz[node] = 1;
    for (auto it : g[node])
        if (!viz[it])
            dfs_plus (it);
    v.push_back(node);
}

void dfs_minus (int node)
{
    viz[node] = 1;
    sol[comp].push_back(node);
    for (auto it : gt[node])
        if (!viz[it])
            dfs_minus (it);
}

int main()
{
    cin >> n >> m;
    for (int i = 1; i <= m; ++i)
    {
        cin >> x >> y;
        g[x].push_back(y);
        gt[y].push_back(x);
    }
    for (int i = 1; i <= n; ++i)
        if (!viz[i])
            dfs_plus (i);
    reverse (v.begin(), v.end());
    for (int i = 1; i <= n; ++i)
        viz[i] = 0;
    for (auto it : v)
        if (!viz[it])
            ++comp, dfs_minus(it);
    cout << comp << '\n';
    for (int i = 1; i <= comp; ++i, cout << '\n')
        for (auto it : sol[i])
            cout << it << ' ';
    return 0;
}