//Algoritmul plus-minus cu 2 DFS


#include <bits/stdc++.h>
#define NMAX 100005
using namespace std;
ifstream fin("ctc.in");
ofstream fout("ctc.out");
int sol;
int n, m;
bool del[NMAX];
vector<int> G[NMAX];
vector<int> GT[NMAX];

vector<int> s[NMAX];
int crt;
vector<int> aux;
struct plus_minus
{
    bool PLUS, MINUS;
};
plus_minus fr[NMAX];
void read();
void edge()
{
    for(int i=1; i<=n; ++i)
        fr[i].MINUS=fr[i].PLUS=0;
    aux.clear();
}
void DFS(int nod, int nr);
int counter()
{
    int temp=0;
    for(int i=1; i<=n; ++i)
        if(fr[i].PLUS && fr[i].MINUS)
        {
            ++temp;
            del[i]=1;
            aux.push_back(i);
        }
    return temp;
}
int main()
{
    read();
    for(int i=1; i<=n; ++i)
    {
        if(del[i])
            continue;
        edge();
        DFS(i, 1);
        DFS(i, 2);
        if(counter())
        {
            ++sol;
            s[crt].insert(s[crt].begin(), aux.begin(), aux.end());
            ++crt;
        }
    }
    fout<<sol<<'\n';
    for(int i=0; i<crt; ++i)
    {
        for(int j=0; j<s[i].size(); ++j)
            fout<<s[i][j]<<' ';
        fout<<'\n';
    }
    return 0;
}
void DFS(int nod, int nr)
{
    if(nr==1)
    {
        if(fr[nod].PLUS)
            return;
        fr[nod].PLUS=1;
    }
    else
    {
        if(fr[nod].MINUS)
            return;
        fr[nod].MINUS=1;
    }
    if(nr==1)
    {
        for(int i=0; i<G[nod].size(); ++i)
            DFS(G[nod][i], nr);
    }
    else
    {
        for(int i=0; i<GT[nod].size(); ++i)
            DFS(GT[nod][i], nr);
    }
}
void read()
{
    fin>>n>>m;
    for(int i=1; i<=m; ++i)
    {
        int x, y;
        fin>>x>>y;
        G[x].push_back(y);
        GT[y].push_back(x);
    }
}