1
  2
  3
  4
  5
  6
  7
  8
  9
 10
 11
 12
 13
 14
 15
 16
 17
 18
 19
 20
 21
 22
 23
 24
 25
 26
 27
 28
 29
 30
 31
 32
 33
 34
 35
 36
 37
 38
 39
 40
 41
 42
 43
 44
 45
 46
 47
 48
 49
 50
 51
 52
 53
 54
 55
 56
 57
 58
 59
 60
 61
 62
 63
 64
 65
 66
 67
 68
 69
 70
 71
 72
 73
 74
 75
 76
 77
 78
 79
 80
 81
 82
 83
 84
 85
 86
 87
 88
 89
 90
 91
 92
 93
 94
 95
 96
 97
 98
 99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
#include <bits/stdc++.h>

using namespace std;

typedef long long LL;
typedef long double LD;
typedef unsigned long long ULL;
typedef vector<int> VI;
typedef vector<LL> VLL;
typedef set<int> SI;
typedef pair<int, int> PII;
typedef pair<LL, LL> PLL;
typedef vector<PII> VPII;
typedef vector<PLL> VPLL;

const int INF = 1000000001;
const LD EPS = 1e-9;
const int MOD = 1000000007;
const LL LLINF = 1000000000000000001;

#define FOR(i, b, e) for(int i = b; i <= e; i++)
#define FORD(i, b, e) for(int i = b; i >= e; i--)
#define ALL(c) (c).begin(), (c).end()
#define SIZE(x) ((int)(x).size())
#define DPRINT(x) cout << #x << ": " << x << endl
#define BOOST ios_base::sync_with_stdio(false); cin.tie(0)

#define MP make_pair
#define PB push_back
#define ST first
#define ND second
#define GGL(x) "Case #" << x << ": "

// ゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴゴ

static const int C = 5;
using arr_t = array<int, C>;
static const arr_t primes = {1000002139, 1000002581, 1000003271, 1000003889, 1000025863};
static const arr_t bases = {3, 5, 2, 2137, 7};

const int maxn = 500042;
arr_t vals[maxn];

LL solve(int n, int X, vector<PII> v)
{
    v.PB({0, X});

    map<arr_t, int> ans;

    arr_t act;
    act.fill(0);

    vector<PII> events;
    events.reserve(2*SIZE(v));
    FOR(i, 0, SIZE(v)-1)
    {
        events.PB({v[i].ST, i+1});
        events.PB({v[i].ND, -i-1});
    }

    sort(ALL(events));

    FOR(i, 0, SIZE(events)-2)
    {
        //cout << events[i].ST << ' ' << events[i].ND << endl;
        int mul = 1;
        if(events[i].ND < 0)
            mul = -1;
        int aa = events[i].ND;
        if(aa < 0)
            aa = -aa;
        FOR(j, 0, SIZE(act)-1)
        {
            act[j] = ((act[j] + (vals[aa][j] * mul)) % primes[j] + primes[j]) % primes[j];
        }
        ans[act] += events[i+1].ST - events[i].ST;
    }

    int r = 0;
    for(auto &i : ans)
        r = max(r, i.ND);

    //cout << "returning " << r << endl;
    return r;
}

int main()
{
    BOOST;

    int n, X, Y;

    cin >> n >> X >> Y;
    /*
    n = 500000;
    X = 1000000000;
    Y = 1000000000;
    */

    vals[0].fill(1);
    FOR(i, 1, n+2)
        FOR(j, 0, SIZE(bases)-1)
        {
            vals[i][j] = (LL(vals[i-1][j]) * bases[j]) % primes[j];
        }

    mt19937 mt(42);

    VPII x, y;
    FOR(i, 1, n)
    {
        int x1, y1, x2, y2;
        cin >> x1 >> y1 >> x2 >> y2;
        
        /*
        x1 = mt() % X;
        x2 = mt() % X;
        y1 = mt() % X;
        y2 = mt() % X;
        */
        
        if(x2 < x1)
            swap(x1, x2);
        if(y2 < y1)
            swap(y1, y2);
        x.PB({x1, x2});
        y.PB({y1, y2});
    }

    cout << solve(n, X, x) * solve(n, Y, y) << '\n';
}