#include<stdio.h>
#include<vector>
using namespace std;

vector <int> g[100005], gt[100005],v[100005];
int ctc[100005],viz[100005],nrt,t[100005],n,m,nr_ctc;

void dfs(int nod,int c)
{
    int i,x,nr=gt[nod].size();
    ctc[nod]=c;
    v[c].push_back(nod);
    for(i=0;i<nr;i++)
        if(!ctc[x=gt[nod][i]])
            dfs(x,c);
}

void dfs_t(int nod)
{
    int i,x,nr=g[nod].size();
    viz[nod]=1;
    for(i=0;i<nr;i++)
        if(!viz[x=g[nod][i]])
            dfs_t(x);
    t[++nrt]=nod;
}

int main ()
{
    int i,a,b,x,j;
    freopen("ctc.in","r",stdin);
    freopen("ctc.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(i=1;i<=m;i++)
    {
        scanf("%d%d",&a,&b);
        g[a].push_back(b);
        gt[b].push_back(a);
    }
    for(i=1;i<=n;i++)
        if(!viz[i])
            dfs_t(i);
    for(i=nrt;i>=1;i--)
        if(!ctc[t[i]])
            dfs(t[i],++nr_ctc);
    printf("%d\n",nr_ctc);
    for(i=nr_ctc;i>=1;i--)
    {
        for(j=0,x=v[i].size(); j<x; j++)
            printf("%d ",v[i][j]);
        printf("\n");
    }
    return 0;
}