#include <bits/stdc++.h>
#pragma GCC optimize("Ofast")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,avx,avx2")
#pragma GCC optimization("O3")
#pragma GCC optimization("unroll-loops")
using namespace std;
const long long int mod = 998244353;
long long int power(long long int x, long long int y) {
  long long int res = 1;
  x = x % mod;
  if (x == 0) return 0;
  while (y > 0) {
    if (y & 1) res = (res * x) % mod;
    y = y >> 1;
    x = (x * x) % mod;
  }
  return res;
}
long long int inv(long long int a) { return power(a, mod - 2); }
long long int add(long long int a, long long int b) { return (a + b) % mod; }
long long int sub(long long int a, long long int b) {
  return (a - b + mod) % mod;
}
long long int mul(long long int a, long long int b) { return (a * b) % mod; }
long long int divi(long long int a, long long int b) {
  return (a * inv(b)) % mod;
}
long long int fact[1000001];
long long int ifact[1000001];
void pre() {
  fact[0] = 1;
  ifact[0] = 1;
  for (int i = 1; i <= 1000000; i++) {
    fact[i] = mul(fact[i - 1], i);
    ifact[i] = power(fact[i], mod - 2);
  }
}
long long int ncr(int n, int r) {
  if (r > n) {
    return 0;
  }
  return divi(fact[n], mul(fact[r], fact[n - r]));
}
void solve() {
  long long int n, d;
  cin >> n;
  vector<long long int> a;
  for (int i = 0; i < 2 * n; i++) {
    cin >> d;
    a.push_back(d);
  }
  sort(a.begin(), a.end());
  long long int ans = 0;
  for (int i = 0; i < n; i++) {
    ans = sub(ans, mul(ncr(2 * n, n), a[i]));
    ans %= mod;
  }
  for (int i = n; i < 2 * n; i++) {
    ans = add(ans, mul(ncr(2 * n, n), a[i]));
    ans %= mod;
  }
  cout << ans << "\n";
}
int main() {
  ios_base::sync_with_stdio(false);
  cin.tie(NULL);
  pre();
  long long int num = 1;
  for (long long int i = 0; i < num; i++) {
    solve();
  }
}
