#include <cstdio>
#include <vector>

using namespace std;

vector<vector<int> > sol;
vector<int> v[100010],q;
int niv[100010],nivmin[100010],nr;
char vaz[100010],bagat[100010];

void dfs(int nod)
{
    vaz[nod]=1;
    niv[nod]=nivmin[nod]=++nr;
    q.push_back(nod);
    bagat[nod]=1;
    for(vector<int>::iterator it=v[nod].begin();it!=v[nod].end();it++)
        if(!vaz[*it])
        {
            dfs(*it);
            nivmin[nod]=min(nivmin[nod],nivmin[*it]);
        }
        else if(bagat[*it]) nivmin[nod]=min(nivmin[nod],nivmin[*it]);
    if(niv[nod]==nivmin[nod])
    {
        int a;
        vector<int> aux;
        do
        {
            a=q.back();
            q.pop_back();
            aux.push_back(a);
            bagat[a]=0;
        }while(a!=nod);
        sol.push_back(aux);
    }
}

int main()
{
    freopen("ctc.in", "r", stdin);
    freopen("ctc.out", "w", stdout);
    int n,m,x,y;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&x,&y);
        v[x].push_back(y);
    }
    for(int i=1;i<=n;i++)
        if(!vaz[i]) dfs(i);
    printf("%d\n",sol.size());
    for(int i=0;i<sol.size();i++)
    {
        for(int j=0;j<sol[i].size();j++) printf("%d ",sol[i][j]);
        printf("\n");
    }
    return 0;
}