#include<stdio.h>
#include<vector>
#include<stack>

#define maxn 100005
#define pb push_back

using namespace std;

FILE*f=fopen("ctc.in","r");
FILE*g=fopen("ctc.out","w");

int n,m,i,x,y,nivel;
int viz[maxn],niv[maxn],lowest[maxn],inSt[maxn];
vector<int>G[maxn],v;
vector< vector<int> >sol;
stack<int>St;

inline void citire () {
	
	fscanf(f,"%d %d",&n,&m);
	
	for ( i = 1 ; i <= m ; ++i ){
		fscanf(f,"%d %d",&x,&y);
		G[x].pb(y);
	}
}

inline void componenta ( int x ){
	int nod;
	v.clear();
	
	do{
		nod = St.top(); inSt[nod] = 0;
		v.pb(nod);
		St.pop();
	}while(nod != x );
	
	sol.pb(v);
}

void dfs ( int nod , int t ){
	viz[nod] = 1; niv[nod] = lowest[nod] = nivel; St.push(nod); inSt[nod] = 1;
	
	for ( vector<int>::iterator itt = G[nod].begin() ; itt != G[nod].end() ; ++itt ){
		int nodvcn = (*itt);
		
		if ( !viz[nodvcn] ){
			++nivel;
			dfs(nodvcn,nod);
			lowest[nod] = lowest[nod] > lowest[nodvcn] ? lowest[nodvcn] : lowest[nod];
		}
		else{
			if ( inSt[nodvcn] ){
				lowest[nod] = lowest[nod] > lowest[nodvcn] ? lowest[nodvcn] : lowest[nod];
			}
		}
	}
	
	if ( lowest[nod] == niv[nod] ){
		componenta(nod);
	}
}

inline void afisare () {
	
	fprintf(g,"%d\n",sol.size());
	
	for ( i = 0 ; i < sol.size() ; ++i ){
		for ( vector<int>::iterator itt = sol[i].begin() ; itt != sol[i].end() ; ++itt ){
			fprintf(g,"%d ",(*itt));
		}
		fprintf(g,"\n");
	}
}

int main () {
	
	citire();
	
	for ( i = 1 ; i <= n ; ++i ){
		if ( !viz[i] ){
			dfs(i,0);
		}
	}
	
	afisare();
	
	
	fclose(f);
	fclose(g);
	
	return 0;
}