#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define db double
#define mp make_pair
#define pii pair<int, int>
#define pb push_back
#define mp make_pair
#define FR first
#define SE second
using namespace std;
inline int read() {
  int x = 0; bool op = false;
  char c = getchar();
  while(!isdigit(c))op |= (c == '-'), c = getchar();
  while(isdigit(c))x = (x << 1) + (x << 3) + (c ^ 48), c = getchar();
  return op ? -x : x;
}
const int N = 5e3 + 10;
const int P = 998244353;
void add(int &a, int b) {a += b; a >= P ? a -= P : 0;}
void sub(int &a, int b) {a -= b; a < 0 ? a += P : 0;}
int n, m, k;
int tot1, tot2;
int ls1[N], ls2[N];
pii a[N];
void disc() {
  for(int i = 1; i <= k; i++)ls1[++tot1] = a[i].FR;
  for(int i = 1; i <= k; i++)ls2[++tot2] = a[i].SE;
  sort(ls1 + 1, ls1 + 1 + tot1);
  tot1 = unique(ls1 + 1, ls1 + 1 + tot1) - (ls1 + 1);
  sort(ls2 + 1, ls2 + 1 + tot2);
  tot2 = unique(ls2 + 1, ls2 + 1 + tot2) - (ls2 + 1);
  for(int i = 1; i <= k; i++) {
    a[i].FR = lower_bound(ls1 + 1, ls1 + 1 + tot1, a[i].FR) - ls1;
    a[i].SE = lower_bound(ls2 + 1, ls2 + 1 + tot2, a[i].SE) - ls2;
  }
  return ;
}
vector<int> A[N];
struct UFS {
  int L[N], R[N];
  void init() {
    for(int i = 1; i <= tot2 + 1; i++)L[i] = R[i] = i;
    return ;
  }
  int ql(int x) {return L[x] == x ? x : L[x] = ql(L[x]);}
  int qr(int x) {return R[x] == x ? x : R[x] = qr(R[x]);}
  void del(int x) {L[x] = x - 1; R[x] = x + 1;}
}ufs;
int sum1[N][N], sum2[N][N], vis[N][N];
int calc2(int l, int r) {
  int i6 = 166374059, res = 0;
  auto _calc = [&](int r) {
    int v1 = 1ll * r * (r + 1) % P, v2 = (2ll * r + 1) % P;
    return 1ll * v1 * v2 % P * i6 % P;
  };
  add(res, _calc(r)); sub(res, _calc(l - 1));
  return res;
}
int calc3(int l, int r) {
  int i2 = 499122177, res = 0;
  auto _calc = [&](int r) {
    int v = 1ll * r * (r + 1) % P * i2 % P;
    return 1ll * v * v % P;
  };
  add(res, _calc(r)); sub(res, _calc(l - 1));
  return res;
}
int calc4(int l, int r) {
  int i30 = 432572553, res = 0;
  auto _calc = [&](int r) {
    int v1 = 1ll * r * (r + 1) % P;
    int v2 = (2ll * r + 1) % P;
    int v3 = (3ll * r % P * r % P + 3ll * r - 1) % P;
    return 1ll * v1 * v2 % P * v3 % P * i30 % P;
  };
  add(res, _calc(r)); sub(res, _calc(l - 1));
  return res;
}
int Calc(int tv1, int td1, int tv2, int td2, int l, int r) {
  int res = 0;
  add(res, 1ll * tv1 * tv2 % P * calc2(l, r) % P);
  add(res, 1ll * (1ll * tv1 * td2 % P + 1ll * tv2 * td1 % P) % P * calc3(l, r) % P);
  add(res, 1ll * td1 * td2 % P * calc4(l, r) % P);
  return res;
}
int calc_coef(int x, int y, int z, int w) {
  int coef = 0;
  if(x == y && z == w)coef = P - 1;
  else if(x == y)coef = (sum1[x][w - 1] - sum1[x][z] ? 0 : 1);
  else if(z == w)coef = (sum2[y - 1][z] - sum2[x][z] ? 0 : 1);
  else {
    vector<pii> p;
    if(vis[x][z])p.pb(mp(0, 0));
    if(vis[x][w])p.pb(mp(0, 2));
    if(vis[y][z])p.pb(mp(2, 0));
    if(vis[y][w])p.pb(mp(2, 2));
    if(sum1[x][w - 1] - sum1[x][z])p.pb(mp(0, 1));
    if(sum1[y][w - 1] - sum1[y][z])p.pb(mp(2, 1));
    if(sum2[y - 1][z] - sum2[x][z])p.pb(mp(1, 0));
    if(sum2[y - 1][w] - sum2[x][w])p.pb(mp(1, 2));
    int f[3][3][3][3] = {0};
    for(auto t : p) {
      int u = t.FR, v = t.SE;
      // printf("%d %d\n", u, v);
      for(int i = 0; i < 3; i++) {
        for(int j = 2; j >= i; j--) {
          for(int k = 0; k < 3; k++) {
            for(int l = 2; l >= k; l--) {
              // printf("qwq:%d\n", f[i][j][k][l]);
              int ni = min(i, u), nj = max(j, u);
              int nk = min(k, v), nl = max(l, v);
              sub(f[ni][nj][nk][nl], f[i][j][k][l]);
            }
          }
        }
      }
      sub(f[u][u][v][v], 1);
    }
    coef = f[0][2][0][2];
  }
  return coef;
}
int calc(int x, int y, int z, int w) {
  int coef = calc_coef(x, y, z, w), res = 0;
  x = ls1[x]; y = ls1[y]; z = ls2[z]; w = ls2[w];
  int st = max(y - x + 2, w - z + 2);
  vector<int> p = {st, y + 1, n - x + 1, w + 1, m - z + 1, min(n, m) + 1};
  sort(p.begin(), p.end());
  p.erase(unique(p.begin(), p.end()), p.end());
  // printf("st:%d\n", st);
  for(int i = 0; i < p.size(); i++) {
    if(p[i] < st)continue;
    if(p[i] > min(n, m))continue;
    int tv1 = 1, td1 = 0, tv2 = 1, td2 = 0;
    int l = p[i], r = p[i + 1] - 1;
    if(l > r)continue;
    (x - 1 < n - p[i] ? add(tv1, x - 1) : (add(tv1, n), sub(td1, 1)));
    (y - p[i] + 1 > 0 ? (sub(tv1, y + 1), add(td1, 1)) : void());
    (z - 1 < m - p[i] ? add(tv2, z - 1) : (add(tv2, m), sub(td2, 1)));
    (w - p[i] + 1 > 0 ? (sub(tv2, w + 1), add(td2, 1)) : void());
    add(res, Calc(tv1, td1, tv2, td2, l, r));
    // printf("range:%d %d %d\n", l, r, Calc(tv1, td1, tv2, td2, 5, 5));
  }
  // printf("calc:%d %d %d %d %d %d\n", x, y, z, w, coef, res);
  return 1ll * res * coef % P;
}
int main() {
  n = read(); m = read(); k = read();
  for(int i = 1; i <= k; i++)a[i].FR = read(), a[i].SE = read();
  disc();
  for(int i = 1; i <= k; i++) {
    vis[a[i].FR][a[i].SE] = 1;
    sum1[a[i].FR][a[i].SE] = 1;
    sum2[a[i].FR][a[i].SE] = 1;
  }
  for(int i = 1; i <= tot1; i++) {
    for(int j = 1; j <= tot2; j++) {
      sum1[i][j] += sum1[i][j - 1];
      sum2[i][j] += sum2[i - 1][j];
    }
  }
  int ans = Calc(n + 1, P - 1, m + 1, P - 1, 1, min(n, m));
  // printf("ans:%d\n", ans);
  sort(a + 1, a + 1 + k);
  vector<array<int, 4> > rec;
  for(int i = 1; i <= k; i++)rec.pb({a[i].FR, a[i].FR, a[i].SE, a[i].SE});
  // for(int i = 1; i <= k; i++)add(ans, calc(a[i].FR, a[i].FR, a[i].SE, a[i].SE));
  for(int i = 1; i <= k; i++) {
    ufs.init();
    vector<int> cnt(tot2 + 1);
    for(int j = i + 1; j <= k; j++)cnt[a[j].SE]++;
    for(int j = 1; j <= tot2; j++)if(!cnt[j])ufs.del(j);
    for(int t = k; t > i; t--) {
      if(!(--cnt[a[t].SE]))ufs.del(a[t].SE);
      if(ufs.ql(a[t].SE - 1) > a[i].SE)continue;
      if(ufs.ql(a[i].SE - 1) > a[t].SE)continue;
      int nl = min(a[i].SE, a[t].SE), nr = max(a[i].SE, a[t].SE);
      int tl = ufs.ql(nl), tr = ufs.qr(nr);
      // printf("id:%d %d %d %d\n", i, t, tl, tr);
      rec.pb({a[i].FR, a[t].FR, nl, nr});
      if(tl)rec.pb({a[i].FR, a[t].FR, tl, nr});
      if(tr <= tot2)rec.pb({a[i].FR, a[t].FR, nl, tr});
      if(tl && tr <= tot2)rec.pb({a[i].FR, a[t].FR, tl, tr});
      // add(ans, calc(a[i].FR, a[t].FR, nl, nr));
      // if(tl)add(ans, calc(a[i].FR, a[t].FR, tl, nr));
      // if(tr <= tot2)add(ans, calc(a[i].FR, a[t].FR, nl, tr));
      // if(tl && tr <= tot2)add(ans, calc(a[i].FR, a[t].FR, tl, tr));
    }
  }
  sort(rec.begin(), rec.end());
  rec.erase(unique(rec.begin(), rec.end()), rec.end());
  for(auto t : rec)add(ans, calc(t[0], t[1], t[2], t[3]));
  printf("%d\n", ans);
  return 0;
}