#include<fstream>
using namespace std;
ifstream f("ctc.in");
ofstream g("ctc.out");
struct nod
{
    int inf;
    nod* leg;
};
typedef nod* PNod;
PNod L1[100001],sol[100001];
int n,m,viz1[100001],viz2[100001],nrc;
/*int exista1(int j,int i)
{
    PNod p;
    for( p=L1[i];p;p=p->leg)
        if(p->inf==j)
            return 1;
    return 0;
}*/
/*int exista2(int j,int i)
{
    PNod p;
    for(p=L2[i];p;p=p->leg)
        if(p->inf==j)
            return 1;
    return 0;
}*/
void add1(int x,int y)
{
    PNod p;
    p=new nod;
    p->inf=y;
    p->leg=L1[x];
    L1[x]=p;          
}
void add2(int x,int y)
{
    PNod p;
    p=new nod;
    p->inf=y;
    p->leg=sol[x];
    sol[x]=p;          
}
void df1(int i,int k)
{
     
    viz1[i]=k;
    for(PNod p=L1[i];p;p=p->leg)
        if(viz1[p->inf]==0)
            df1(p->inf,k);
}
void df2(int i,int k)
{
     
    viz2[i]=k;
    for(PNod p=sol[i];p;p=p->leg)
		if(viz2[p->inf]==0)
			df2(p->inf,k);
}
void citire()
{
    int x,y,i;
    f>>n>>m;
    for(i=1;i<=m;i++)
    {
        f>>x>>y;add1(x,y);add2(y,x);
    }
}
void compcon()
{
 int i,k=0;
 for(i=1;i<=n;i++)
 {
	 if(viz1[i]==0)
	 {
		 k++;
		 df1(i,k);
		 df2(i,k);
		 for(int j=1;j<=n;j++)
			 if(viz1[j]==k)
				 if(viz1[j]!=viz2[j])viz1[j]=0;
			 
	 }
 }
  g<<k<<'\n'; 
for (i=1;i<=k;i++)
{
	for(int j=1;j<=n;j++)
		if(viz1[j]==i)	g<<j<<" ";
	g<<'\n';
}
}
int main()
{
    citire();
    compcon();
    g<<nrc;
    return 0;
}