#include <bits/stdc++.h>
using namespace std;
ifstream f ("ctc.in");
ofstream g ("ctc.out");
const int N_MAX = 100005;
int n,m,x,y,p;
vector <int> G[N_MAX],GT[N_MAX];
vector <int> sol[N_MAX];
vector <bool> viz;
vector <int> S;
void dfs (int nodStart){
    viz[nodStart] =  true;
    for(auto vecin : G[nodStart]){
        if(!viz[vecin])
            dfs(vecin);


    }
     S.push_back(nodStart);
}
void dfsGT (int nodStart){
    viz[nodStart] =  true;
    for(auto vecin : GT[nodStart]){
        if(!viz[vecin]){
            dfsGT(vecin);
        }
    }

}
int main()
{
    f >> n >> m;
    viz.resize(N_MAX,false);
    for(int i = 1; i <= m; i++){
        f >> x >> y;
        G[x].push_back(y);
        GT[y].push_back(x);
    }
    for(int i = 1; i <= n; i++){
        if (!viz[i]){
            dfs(i);
        }
    }
    viz.assign(N_MAX,false);
    for(vector <int> :: reverse_iterator it = S.rbegin();it != S.rend(); it++){
        if(!viz[*it]){
            p++;
            sol[p].push_back(*it);
            dfsGT(*it);
        }
    }
    g << p << '\n';
    for (int i = 1; i <= p; i++ , g << '\n')
        for(auto j : sol[i])
                g << j << ' ';
    return 0;
}