#include <bits/stdc++.h>
using namespace std;
template <typename T>
void chkmax(T &x, T y) {
  x = x > y ? x : y;
}
template <typename T>
void chkmin(T &x, T y) {
  x = x > y ? y : x;
}
const int INF = 2139062143;
template <typename T>
void read(T &x) {
  x = 0;
  bool f = 1;
  char ch;
  do {
    ch = getchar();
    if (ch == '-') f = 0;
  } while (ch > '9' || ch < '0');
  do {
    x = x * 10 + ch - '0';
    ch = getchar();
  } while (ch >= '0' && ch <= '9');
  x = f ? x : -x;
}
template <typename T>
void write(T x) {
  if (x < 0) x = ~x + 1, putchar('-');
  if (x > 9) write(x / 10);
  putchar(x % 10 + '0');
}
const int mod = 1e9 + 7;
const int N = 4e5 + 7;
struct Node {
  int l, r;
} a[N], b[N];
struct Segment_Tree {
  struct Tree {
    int l, r;
    long long dat, mul;
  } t[N * 4];
  inline void build(int p, int l, int r) {
    t[p].l = l;
    t[p].r = r;
    t[p].mul = -1;
    if (l == r) return;
    int mid = (t[p].l + t[p].r) >> 1;
    build((p << 1), l, mid);
    build((p << 1 | 1), mid + 1, r);
    pushup(p);
  }
  inline void pushup(int p) {
    t[p].dat = (t[(p << 1)].dat + t[(p << 1 | 1)].dat) % mod;
  }
  inline void pushdown(int p) {
    if (t[p].mul == -1) return;
    t[(p << 1)].dat =
        1ll * (t[(p << 1)].r - t[(p << 1)].l + 1) % mod * t[p].mul % mod;
    t[(p << 1 | 1)].dat = 1ll * (t[(p << 1 | 1)].r - t[(p << 1 | 1)].l + 1) %
                          mod * t[p].mul % mod;
    t[(p << 1)].mul = t[p].mul % mod;
    t[(p << 1 | 1)].mul = t[p].mul % mod;
    t[p].mul = -1;
  }
  inline void modify(int p, int l, int r, int k) {
    if (l > r) return;
    if (l <= t[p].l && t[p].r <= r) {
      t[p].mul = k % mod;
      t[p].dat = 1ll * (t[p].r - t[p].l + 1) % mod * k % mod;
      return;
    }
    pushdown(p);
    int mid = (t[p].l + t[p].r) >> 1;
    if (l <= mid) modify((p << 1), l, r, k);
    if (r > mid) modify((p << 1 | 1), l, r, k);
    pushup(p);
  }
  inline int query(int p, int l, int r) {
    if (l <= t[p].l && t[p].r <= r) return t[p].dat;
    pushdown(p);
    int mid = (t[p].l + t[p].r) >> 1, ans = 0;
    if (l <= mid) ans = (ans + query((p << 1), l, r)) % mod;
    if (r > mid) ans = (ans + query((p << 1 | 1), l, r)) % mod;
    return ans;
  }
} st[2];
int k, n, m, cnt, num, tmp[N * 4], mx[2][N];
inline int qpow(int a, int b) {
  int ret = 1;
  while (b) {
    if (b & 1) ret = 1ll * ret * a % mod;
    a = 1ll * a * a % mod;
    b >>= 1;
  }
  return ret;
}
signed main() {
  read(k);
  read(n);
  read(m);
  for (int i = 1; i <= n; i++) {
    read(a[i].l);
    read(a[i].r);
    a[i].l--;
    tmp[++cnt] = a[i].l;
    tmp[++cnt] = a[i].r;
  }
  for (int i = 1; i <= m; i++) {
    read(b[i].l);
    read(b[i].r);
    b[i].l--;
    tmp[++cnt] = b[i].l;
    tmp[++cnt] = b[i].r;
  }
  tmp[++cnt] = 0;
  tmp[++cnt] = k;
  sort(tmp + 1, tmp + cnt + 1);
  num = unique(tmp + 1, tmp + cnt + 1) - tmp - 1;
  for (int i = 1; i <= n; i++) {
    a[i].l = lower_bound(tmp + 1, tmp + num + 1, a[i].l) - tmp;
    a[i].r = lower_bound(tmp + 1, tmp + num + 1, a[i].r) - tmp;
    chkmax(mx[0][a[i].r], a[i].l);
  }
  for (int i = 1; i <= m; i++) {
    b[i].l = lower_bound(tmp + 1, tmp + num + 1, b[i].l) - tmp;
    b[i].r = lower_bound(tmp + 1, tmp + num + 1, b[i].r) - tmp;
    chkmax(mx[1][b[i].r], b[i].l);
  }
  st[0].build(1, 1, num);
  st[1].build(1, 1, num);
  st[0].modify(1, 1, 1, 1);
  st[1].modify(1, 1, 1, 1);
  for (int i = 2; i <= num; i++) {
    int x = st[0].query(1, 1, i - 1);
    int y = st[1].query(1, 1, i - 1);
    int q = qpow(2, tmp[i] - tmp[i - 1] - 1);
    q = (q - 1 + mod) % mod;
    st[0].modify(1, i - 1, i - 1, (y + st[0].query(1, i - 1, i - 1)) % mod);
    st[1].modify(1, i - 1, i - 1, (x + st[1].query(1, i - 1, i - 1)) % mod);
    q = 1ll * q * (x + y) % mod;
    st[0].modify(1, i, i, q);
    st[1].modify(1, i, i, q);
    st[0].modify(1, 1, mx[1][i], 0);
    st[1].modify(1, 1, mx[0][i], 0);
  }
  printf("%lld\n", 1ll *
                       (st[0].query(1, 1, num) + st[1].query(1, 1, num) % mod) *
                       qpow(2, mod - 2) % mod);
  return 0;
}
