#include<cstdio>
#include<vector>
using namespace std;
int n,m,k,x,y,a[100004],viz[100004];
vector < int > v[100004],vt[100004],ras[100004];
void dfs(int nod)
{
    viz[nod]=1;
    for(int i=0;i<v[nod].size();i++)
        if(viz[v[nod][i]]==0)
            dfs(v[nod][i]);
    a[++k]=nod;
}
void dfp(int nod)
{
    viz[nod]=0;
    ras[k].push_back(nod);
    for(int i=0;i<vt[nod].size();i++)
        if(viz[vt[nod][i]]==1)
            dfp(vt[nod][i]);
}
int main()
{
    freopen("ctc.in","r",stdin);
    freopen("ctc.out","w",stdout);
    scanf("%d %d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        scanf("%d %d",&x,&y);
        v[x].push_back(y);
        vt[y].push_back(x);
    }
    for(int i=1;i<=n;i++)
        if(viz[i]==0)dfs(i);
    k=0;
    for(int i=n;i>=1;i--)
        if(viz[a[i]]==1)
        {
            k++;
            dfp(a[i]);
        }
    printf("%d\n",k);
    for(int i=1;i<=k;i++&&printf("\n"))
        for(int j=0;j<ras[i].size();j++)
            printf("%d ",ras[i][j]);
    return 0;
}