//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];


struct plus_minus
{
    bool PLUS=false, MINUS=false;
};
plus_minus fr[NMAX];
void read();

void DFS(int nod, int nr);
void edge()
{
    for(int i=1; i<=n; ++i)
    {
        if(fr[i].PLUS && fr[i].MINUS)
        {
            del[i]=1;
            s[sol].push_back(i);
        }
        fr[i].PLUS=fr[i].MINUS=0;
    }
}
int main()
{
    read();
    for(int i=1; i<=n; ++i)
    {
        if(del[i])
            continue;
        DFS(i, 1);
        DFS(i, 2);
        edge();
        ++sol;
    }
    fout<<sol<<'\n';
    for(int i=0; i<sol; ++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(del[nod])
        return;
    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);
    }
}