#include <bits/stdc++.h>
using namespace std;
int n, m, a, b, D[3004][3004], tt = -1;
vector<int> e[3005];
vector<pair<int, int> > d1[3005], d2[3005];
struct RES {
  int a, b, c, d;
} res;
inline void bfs(int st) {
  queue<int> q;
  D[st][st] = 0;
  q.push(st);
  while (!q.empty()) {
    int r = q.front();
    q.pop();
    for (vector<int>::iterator it = e[r].begin(); it != e[r].end(); it++)
      if (D[st][*it] == -1) {
        D[st][*it] = D[st][r] + 1;
        q.push(*it);
      }
  }
  return;
}
inline void init() {
  for (int i = 1; i <= n; i++) {
    for (int j = 1; j <= n; j++)
      if (D[i][j] != -1) d1[i].push_back(make_pair(D[i][j], j));
    for (int j = 1; j <= n; j++)
      if (D[j][i] != -1) d2[i].push_back(make_pair(D[j][i], j));
    sort(d1[i].begin(), d1[i].end());
    sort(d2[i].begin(), d2[i].end());
    reverse(d1[i].begin(), d1[i].end());
    reverse(d2[i].begin(), d2[i].end());
  }
}
int main() {
  cin >> n >> m;
  while (m--) {
    cin >> a >> b;
    if (a != b) e[a].push_back(b);
  }
  memset(D, -1, sizeof(D));
  for (int i = 1; i <= n; i++) bfs(i);
  init();
  for (int i = 1; i <= n; i++) {
    for (int j = 1; j <= n; j++) {
      if (i != j && D[i][j] != -1) {
        for (int q = 0; q < d2[i].size() && q < 3; q++) {
          if (d2[i][q].second != j && d2[i][q].second != i)
            for (int p = 0; p < d1[j].size() && p < 4; p++) {
              if (d1[j][p].second == i || d1[j][p].second == j ||
                  d1[j][p].second == d2[i][q].second)
                continue;
              int tmp = D[i][j] + d1[j][p].first + d2[i][q].first;
              if (tmp > tt)
                res.a = d2[i][q].second, res.b = i, res.c = j,
                res.d = d1[j][p].second, tt = tmp;
            }
        }
      }
    }
  }
  cout << res.a << " " << res.b << " " << res.c << " " << res.d << endl;
  return 0;
}
