#include <bits/stdc++.h>
using namespace std;
const long double pi = acos((long double)-1);
const int inf = 0x3f3f3f3f;
const double eps = 1e-8;
const int maxn = 1e6 + 5;
const int maxm = 1e6 + 5;
const long long mod = 998244353;
void in(initializer_list<int *> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) scanf("%d", *ptr);
}
void in(initializer_list<long long *> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) scanf("%lld", *ptr);
}
void in(initializer_list<double *> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) scanf("%lf", *ptr);
}
void out(initializer_list<int> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) printf("%d ", *ptr);
}
void outln(initializer_list<int> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) printf("%d\n", *ptr);
}
void out(initializer_list<long long> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) printf("%lld ", *ptr);
}
void outln(initializer_list<long long> li) {
  for (auto ptr = li.begin(); ptr != li.end(); ptr++) printf("%lld\n", *ptr);
}
int in(int &a, int &b, int &c, int &d) {
  return scanf("%d%d%d%d", &a, &b, &c, &d);
}
int in(int &a, int &b, int &c) { return scanf("%d%d%d", &a, &b, &c); }
int in(int &a, int &b) { return scanf("%d%d", &a, &b); }
int in(long long &a) { return scanf("%lld", &a); }
int in(int &a) { return scanf("%d", &a); }
int in(char *s) { return scanf("%s", s); }
int in(char &c) { return scanf("%c", &c); }
void out(int a, bool ln) { printf("%d%c", a, " \n"[ln]); }
void out(long long a, bool ln) { printf("%lld%c", a, " \n"[ln]); }
void out(double a, int digit, bool ln) {
  printf("%.*f%c", digit, a, " \n"[ln]);
}
long long gcd(long long a, long long b) { return b ? gcd(b, a % b) : a; }
long long lcm(long long a, long long b) { return a / gcd(a, b) * b; }
map<int, long long> ha, hb;
map<pair<int, int>, long long> hab;
map<int, long long> temp;
long long ans = 0, suma = 1, sumb = 1, sumab = 1;
struct node {
  int a, b;
} no[maxn];
bool cmp1(node a, node b) {
  if (a.a == b.a) return a.b < b.b;
  return a.a < b.a;
}
long long qpl[maxn];
void init(int n) {
  qpl[1] = 1;
  long long ans = 1;
  for (int i = 2; i <= n; ++i) {
    ans *= i;
    ans %= mod;
    qpl[i] = ans;
  }
}
int main() {
  int n;
  in(n);
  init(n);
  for (int i = 0; i < n; ++i) {
    in(no[i].a, no[i].b);
    ha[no[i].a]++, hb[no[i].b]++;
    hab[make_pair(no[i].a, no[i].b)]++;
  }
  map<int, long long>::iterator it = ha.begin();
  for (; it != ha.end(); it++) {
    suma *= qpl[it->second];
    suma %= mod;
  }
  it = hb.begin();
  for (; it != hb.end(); it++) {
    sumb *= qpl[it->second];
    sumb %= mod;
  }
  ans = (qpl[n]) % mod;
  sort(no, no + n, cmp1);
  bool flag = 1;
  for (int i = 1; i < n; ++i) {
    if (no[i].b < no[i - 1].b) {
      flag = 0;
      break;
    }
  }
  if (!flag) {
    ans = (ans - suma - sumb + mod * 3) % mod;
    cout << ans << "\n";
  } else {
    map<pair<int, int>, long long>::iterator it = hab.begin();
    for (; it != hab.end(); it++) {
      sumab *= qpl[it->second];
      sumab %= mod;
    }
    ans = (ans + sumab) % mod;
    ans = (ans - suma - sumb + mod * 3) % mod;
    cout << ans << "\n";
  }
  return 0;
}
