#include <bits/stdc++.h>
using namespace std;
int n, m, l, r;
bool a[2020], b[2020];
double f[2020][2020];
int main() {
  scanf("%d%d", &n, &m);
  l = n;
  r = n;
  int i, j, x, y;
  for (i = 1; i <= m; ++i) {
    scanf("%d%d", &x, &y);
    if (!a[x]) --l;
    if (!b[y]) --r;
    a[x] = true;
    b[y] = true;
  }
  for (i = 1; i <= l; ++i) f[i][0] = (n + i * f[i - 1][0]) / i;
  for (j = 1; j <= r; ++j) f[0][j] = (n + j * f[0][j - 1]) / j;
  for (i = 1; i <= l; ++i)
    for (j = 1; j <= r; ++j)
      f[i][j] = (n * n + i * (n - j) * f[i - 1][j] + (n - i) * j * f[i][j - 1] +
                 i * j * f[i - 1][j - 1]) /
                (n * n - (n - i) * (n - j));
  printf("%0.8lf\n", f[l][r]);
  return 0;
}
