#include <fstream>
#include <vector>
#include <bitset>
#include <stdio.h>

#define maxN 100100

using namespace std;

vector <int> G[maxN], GT[maxN], S[maxN], A;
int N, M, T[maxN], Sol, Nr;

bitset <maxN> viz;


void df (int x) 
{
    viz[x] = 1;
    for (int i = 0; i < (int) G[x].size(); ++ i)
        if (!viz[G[x][i]])
            df(G[x][i]);
    T[x] = ++ Nr;
}


void df2 (int x) 
{
    S[Sol].push_back(x);
    viz[x] = 1;
    for (int i = 0; i < (int) GT[x].size(); ++ i)
        if (!viz[GT[x][i]])
            df2(GT[x][i]);
}

inline bool cmp (int a, int b) 
{
    return T[a] > T[b];
}

int main () 
{
        int i, a, b, 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])
              df (i);

    for (i = 1; i <= N; ++ i)
        A.push_back(i);
      
    sort(A.begin(), A.end(), cmp);

    viz.reset ();

    for (i = 0; i < N; ++ i)
        if (!viz[A[i]]) 
        {
            Sol ++;
            df2 (A[i]);      
        }
    printf("%d\n",Sol);

    for (i = 1; i <= Sol; ++ i) 
    {
        
        for (j = 0; j < (int) S[i].size(); ++ j)
        
            printf("%d ",S[i][j]);
            printf("\n");
    }
 
 return 0;
}