#include <cstdio>
#include <vector>
#include <bitset>

#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]);
		puts("");
	}
}