Wannafly 挑战赛12 D

链接：https://www.nowcoder.com/acm/contest/79/D
来源：牛客网

时间限制：C/C++ 1秒，其他语言2秒
空间限制：C/C++ 262144K，其他语言524288K
64bit IO Format: %lld

题目描述

给一个n x m的01矩阵，其中C个格子是1，其余全是0。求有多少全0的子矩阵。
答案对10⁹ + 7取模。

输入描述:

第一行三个数，分别表示n,m,C。
接下来C行，每行两个数x,y（1 <= x <= n, 1 <= y <= m)，表示为1的格子，保证互不相同。

输出描述:

仅一行一个数，表示答案。

示例1

输入

3 3 1
2 3

输出

24

备注:

1 ≤ n,m ≤ 10

⁹

0 ≤ C ≤ 5000
思路：按坐标排序，统计每个1第一次出现的子矩阵个数即可。

#include <iostream>
#include <fstream>
#include <sstream>
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#include <list>
#include <iomanip>
#include <cctype>
#include <cassert>
#include <bitset>
#include <ctime>

using namespace std;

#define pau system("pause")
#define ll long long
#define pii pair<int, int>
#define pb push_back
#define mp make_pair
#define clr(a, x) memset(a, x, sizeof(a))

const double pi = acos(-1.0);
const int INF = 0x3f3f3f3f;
const int MOD = 1e9 + 7;
const double EPS = 1e-9;

/*
#include <ext/pb_ds/assoc_container.hpp>
#include <ext/pb_ds/tree_policy.hpp>

using namespace __gnu_pbds;
tree<pli, null_type, greater<pli>, rb_tree_tag, tree_order_statistics_node_update> T;
*/


struct gg {
    int x, y;
    gg () {}
    gg (int x, int y) : x(x), y(y) {}
    bool operator < (const gg &g) const {
        return x == g.x ? y < g.y : x < g.x;}
} g[5015];
int n, m, c;
ll ans;
int main() {
    scanf("%d%d%d", &n, &m, &c);
    for (int i = 1; i <= c; ++i) {
        scanf("%d%d", &g[i].x, &g[i].y);
    }
    sort(g + 1, g + c + 1);
    g[0] = gg(0, 0);
    for (int i = 1; i <= c; ++i) {
        ll k = (n - g[i].x + 1);
        for (int j = i - 1, l = 1, r = m; ~j; --j) {
            ans += k * (g[i].y - l + 1) % MOD * (r - g[i].y + 1) % MOD * (g[j + 1].x - g[j].x) % MOD;
            if (g[j].y <= g[i].y) l = max(l, g[j].y + 1);
            if (g[i].y <= g[j].y) r = min(r, g[j].y - 1);
        }
    }
    ans = (ll)n * (n + 1) / 2 % MOD * ((ll)m * (m + 1) / 2 % MOD) - ans;
    printf("%d
", (ans % MOD + MOD) % MOD);
    return 0;
}