#include <iostream>
#include <fstream>
#include <map>
#include <list>
#include <algorithm>
#include <vector>
using namespace std;
ifstream f("ctc.in");
ofstream g("ctc.out");
list<int>lext,lint;
map<int,int>ge[110000],gi[110000];
bool ok[110000];
int n,m,x,y, nrsol, kmax;
vector<int> solfin[110000];
void dfse(int x,bool viz[])
{
    viz[x]=1;
    lext.push_back(x);
    for(auto it=ge[x].begin();it!=ge[x].end();it++)
    {
        if(!ok[it->first] && !viz[it->first])
            dfse(it->first,viz);
    }
}

void dfsi(int x,bool viz[])
{
    viz[x]=1;
    for(auto it=gi[x].begin();it!=gi[x].end();it++)
    {
        if(!ok[it->first] && !viz[it->first])
            dfsi(it->first,viz);
    }
}

void ctc(int x)
{
    lext.clear();
    bool vizext[110000]={0};
    dfse(x,vizext);
    while(!lext.empty())
    {
        bool vizint[110000]={0};
        int sol[110000]={0}, k=0;
        dfsi(*lext.begin(),vizint);
        for(auto it=lext.begin(); it!=lext.end(); it++)
            if(vizint[*it])
            {
                ok[*it]=1;
                sol[++k]=*it;
                lext.erase(it);
                it--;
            }
        nrsol++;
        for(int i=1;i<=k;i++)
            solfin[nrsol].push_back(sol[i]);
    }
}
int main() {
    f>>n>>m;
    for(int i=1;i<=m;i++)
    {
        f>>x>>y;
        ge[x][y]=gi[y][x]=1;
    }
    for(int i=1;i<=n;i++)
        if(!ok[i])
            ctc(i);
    g<<nrsol<<'\n';
    for(int i=1;i<=nrsol;i++)
    {
        sort(solfin[i].begin(),solfin[i].end());
        for(auto &j: solfin[i])
        g<<j<<' ';
        g<<'\n';
    }
    return 0;
}