func_code_string: #include <bits/stdc++.h>   using namespace std;   template <class T>   void _R(T &x) {     cin >> x;   }   void _R(int &x) { scanf( %d , &x); }   void _R(int64_t &x) { scanf( %lld , &x); }   void _R(double &x) { scanf( %lf , &x); }   void _R(char &x) { scanf(  %c , &x); }   void _R(char *x) { scanf( %s , x); }   void R() {}   template <class T, class... U>   void R(T &head, U &...tail) {     _R(head);     R(tail...);   }   template <class T>   void _W(const T &x) {     cout << x;   }   void _W(const int &x) { printf( %d , x); }   void _W(const int64_t &x) { printf( %lld , x); }   void _W(const double &x) { printf( %.16f , x); }   void _W(const char &x) { putchar(x); }   void _W(const char *x) { printf( %s , x); }   template <class T, class U>   void _W(const pair<T, U> &x) {     _W(x.first);     putchar(   );     _W(x.second);   }   template <class T>   void _W(const vector<T> &x) {     for (auto i = x.begin(); i != x.end(); _W(*i++))       if (i != x.cbegin()) putchar(   );   }   void W() {}   template <class T, class... U>   void W(const T &head, const U &...tail) {     _W(head);     putchar(sizeof...(tail) ?     :   n );     W(tail...);   }   template <class T, class... U>   void DEBUG(const T &head, const U &...tail) {}   int MOD = 1e9 + 7;   void ADD(long long &x, long long v) {     x = (x + v) % MOD;     if (x < 0) x += MOD;   }   const int SIZE = 5e5 + 10;   long long mypow(long long x, long long y) {     x %= MOD;     long long res = 1 % MOD;     while (y) {       if (y & 1) res = res * x % MOD;       y >>= 1;       x = x * x % MOD;     }     return res;   }   int ll[2][SIZE];   int d[SIZE], dn;   vector<pair<int, int> > seg[2];   template <class T>   void maa(T &x, T y) {     if (x < y) x = y;   }   template <class T>   void mii(T &x, T y) {     if (x > y) x = y;   }   int get_v(int x) { return lower_bound(d, d + dn, x) - d; }   int get_len(int x) { return d[x] - d[x - 1]; }   struct BIT {     int _n;     long long _d[SIZE];     void init(int __n) {       _n = __n;       memset(_d, 0, sizeof(long long) * (_n + 1));     }     void ins(int x, long long v) {       x++;       for (; x <= _n; x += x & -x) ADD(_d[x], v);     }     long long qq(int x) {       x++;       long long res = 0;       for (; x; x -= x & -x) ADD(res, _d[x]);       return res;     }     long long qq_range(int x, int y) {       long long ret = qq(y);       if (x) ADD(ret, -qq(x - 1));       return ret;     }   } bit[2];   int main() {     int K, N, M;     R(K, N, M);     if (N + M == 0) {       W(mypow(2, K));       return 0;     }     d[dn++] = 0;     for (int i = 0; i < (N); ++i) {       int x, y;       R(x, y);       d[dn++] = x - 1;       d[dn++] = y;       seg[0].push_back(make_pair(x - 1, y));     }     for (int i = 0; i < (M); ++i) {       int x, y;       R(x, y);       d[dn++] = x - 1;       d[dn++] = y;       seg[1].push_back(make_pair(x - 1, y));     }     sort(d, d + dn);     dn = unique(d, d + dn) - d;     for (int j = 0; j < (2); ++j) {       for (int i = 0; i < (((int)(seg[j]).size())); ++i) {         maa(ll[j][get_v(seg[j][i].second)], get_v(seg[j][i].first) + 1);       }     }     for (int i = 0; i < (2); ++i) bit[i].init(dn + 2);     bit[0].ins(1, mypow(2, get_len(1)) - 2);     if (ll[0][1] == 0) {       bit[1].ins(0, 1);     }     if (ll[1][1] == 0) {       bit[0].ins(0, 1);     }     int bb[2] = {};     for (int i = (2); i < (dn); ++i) {       long long both[2] = {};       int len = get_len(i);       for (int j = 0; j < (2); ++j) ADD(both[j], bit[j].qq(i - 1));       while (bb[0] < ll[1][i]) {         long long v = bit[0].qq_range(bb[0], bb[0]);         bit[0].ins(bb[0], -v);         bb[0]++;       }       if (ll[1][i] != i) bit[0].ins(i - 1, both[1]);       while (bb[1] < ll[0][i]) {         long long v = bit[1].qq_range(bb[1], bb[1]);         bit[1].ins(bb[1], -v);         bb[1]++;       }       if (ll[0][i] != i) bit[1].ins(i - 1, both[0]);       bit[0].ins(i, (both[0] + both[1]) * (mypow(2, len) - 2) % MOD);     }     long long ans = 0;     for (int i = 0; i < (2); ++i) ADD(ans, bit[i].qq(dn));     W(ans * mypow(2, (K - d[dn - 1])) % MOD);     return 0;   }   