#include<fstream>
using namespace std;
ifstream cin("ctc.in");
ofstream cout("ctc.out");
int n,m,p,sol,x,y,Plus[4005],Minus[4005],a[4005][4005],comp[4005];
void dfsplus(int k)
{
    Plus[k]=1;
    for(int i=1;i<=n;i++)
        if(a[k][i] && !Plus[i]) dfsplus(i);
}
void dfsminus(int k)
{
    Minus[k]=1;
    for(int i=1;i<=n;i++)
        if(a[i][k] && !Minus[i]) dfsminus(i);
}
int main()
{
    cin>>n>>m;
    while(m--)
    {
        cin>>x>>y;
        a[x][y]=1;
    }
    for(int i=1;i<=n;i++)
        if(!comp[i])
        {
            ++sol;
            dfsplus(i);
            dfsminus(i);
            for(int j=1;j<=n;j++)
                comp[j]+=Plus[j]*Minus[j]*sol;
            for(int j=1;j<=n;j++)
                if(comp[j]) Plus[j]=Minus[j]=0;
        }
    cout<<sol<<'\n';
    for(int i=1;i<=sol;i++,cout<<'\n')
        for(int j=1;j<=n;j++)
            if(comp[j]==i) cout<<j<<' ';
}