#include <iostream>
#include <vector>
#include <deque>
using namespace std;

struct adat
{
    vector<int> be;
    vector<int> ki;
    int be_db;
    int ki_db;
    int szint;
    int lat = 0;
};

deque<int> kov;
vector<adat> x;
vector<int> szint_st, sol1, sol2;

int n, m, a, b, szint_max, db;

void szelessegi(int akt)
{
    szint_max = max(szint_max, x[akt].szint);
    for (auto &e : x[akt].ki)
    {
        if (x[e].szint == 0 || x[e].szint < x[akt].szint + 1)
        {
            x[e].szint = x[akt].szint + 1;
            kov.push_back(e);
        }
    }
}

void keres()
{
    for (int i = 1; i <= n; ++i)
    {
        if (x[i].szint == 1)
        {
            a = i;
            break;
        }
    }
    for (int i = a + 1; i <= n; ++i)
    {
        if (x[i].szint == 1)
        {
            b = i;
            break;
        }
    }

    if (b == 0)
        b = a;
}

void bejar(int akt1, int akt2)
{
    pair<int, int> kov = {-1, -1};
    sol1.push_back(akt1);
    sol2.push_back(akt2);
    for (auto &e : x[akt1].ki)
    {
        if (szint_st[x[e].szint] >= 2)
        {
            if (x[e].lat == 0)
            {
                kov.first = e;
                x[e].lat = 1;
                break;
            }
        }
        else
        {
            kov.first = e;
            x[e].lat = 1;
            break;
        }
    }

    for (auto &e : x[akt2].ki)
    {
        if (szint_st[x[e].szint] >= 2)
        {
            if (x[e].lat == 0)
            {
                kov.second = e;
                x[e].lat = 1;
                break;
            }
        }
        else
        {
            kov.second = e;
            x[e].lat = 1;
            break;
        }
    }

    if(kov.first != -1 && kov.second != -1) bejar(kov.first, kov.second);
}

int main()
{
    cin >> n >> m;
    x.resize(n + 1);
    for (int i = 1; i <= m; ++i)
    {
        cin >> a >> b;

        x[a].ki.push_back(b);
        x[a].ki_db = x[a].ki.size();

        x[b].be.push_back(a);
        x[b].be_db = x[b].be.size();
    }

    for (int i = 1; i <= n; ++i)
    {
        if (x[i].be_db == 0)
        {
            x[i].szint = 1;
            kov.push_back(i);
        }
    }

    while (!kov.empty())
    {
        szelessegi(kov[0]);
        kov.pop_front();
    }

    szint_st.resize(szint_max + 1);
    for (int i = 1; i <= n; ++i)
    {
        ++szint_st[x[i].szint];
    }

    a = 0;
    b = 0;
    keres();
    bejar(a, b);
    for (int i = 1; i <= szint_max; ++i) // feltetelezve, hogy mindig az utolso szinten levo csp egy celpont
    {
        if (szint_st[i] >= 2)
        {
            db += 2;
        }
        else
            db += szint_st[i];
    }

    cout << db << endl;

    for (auto &e : sol1)
        cout << e << " ";
    cout << 0 << endl;

    for (auto &e : sol2)
        cout << e << " ";
    cout << 0;
}