#include <cstdio>
#include <vector>
#include <stack>
#define N 100001
using namespace std;
vector <int> g[N],gt[N],sol[N];
stack <int> st;
bool v[N];
int n,m,lg;
void R()
{
    int x,y;
    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);
        g[x].push_back(y);
        gt[y].push_back(x);
    }
}

void DFS(int k)
{
    v[k]=1;
    vector <int>::iterator it;
    for(it=g[k].begin(); it!=g[k].end(); it++)
    if(!v[*it])
    {
        DFS(*it);
    }
    st.push(k);
}

void DFST(int k)
{
    v[k]=1;
     sol[lg].push_back(k);
    vector <int>::iterator it;
    for(it=gt[k].begin(); it!=gt[k].end(); it++)
    if(!v[*it])
        DFST(*it);

}

void Rez()
{
    int nr;
    for(int i=1; i<=n; i++)
    v[i]=0;
    while(!st.empty())
    {
        nr=st.top();
        st.pop();
        if(v[nr]==0)
        {
            lg++;
            //printf("%d ",nr);
            DFST(nr);
        }
    }
}

int main()
{
    R();
    int i,nr;
    for( i=1; i<=n; i++)
    if(!v[i])
    DFS(i);
    Rez();
    printf("%d\n",lg);
    vector <int>::iterator it;
    for(i=1; i<=lg; i++)
    {
        for(it=sol[i].begin(); it!=sol[i].end(); it++)
        {
            nr=*it;
            printf("%d ",nr);
        }
        printf("\n");
    }
    return 0;
}