Wannafly挑战12 D矩阵计数 - 贪心

题目链接： 点击打开链接

:

题解思路：我们不妨先求出不满足条件的矩阵有多少，最后用全部的减去一下就可以了。如图平面上有四个是1的点，假设先枚举到绿色那块的右上方那个点，首先我们是以x轴坐标大小排好序的，那么我们发现绿色1号那块与4号里面所有的点一对一形成的矩形都是不满足的，其次是2号应该是与3号那块，为什么呢?因为黄色边那个点之前已经把3号上面那块算过了，所以这次只要算3号那块就可以了，之后5和6都是算3号那块面积。此时不难发现我们只需要每次用一个最高y和最小y限制住右上那块面积和左下面积就可以算出不重复的总数量了，实现其实很简单，具体看代码就懂了。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mx = 5e3 + 10,mod = 1e9+7;
int n,m,c;
struct node
{
	ll x, y;
	bool operator < (node A)const
	{
		if(x==A.x) return y < A.y;
		return x < A.x;
	}	
}s[mx];
int main()
{
	scanf("%d%d%d",&n,&m,&c);
	for(int i=1;i<=c;i++) scanf("%lld%lld",&s[i].x,&s[i].y);
	sort(s+1,s+c+1);
	ll cnt = (1ll*n*m)%(2*mod), ans = 0;
	for(int i=1;i<=c;i++){
		ll l = 0,r = m+1;
		for(int j=i-1;~j;j--){
			ans += (r-s[i].y)*(n-s[i].x+1)%mod*(s[j+1].x-s[j].x)%mod*(s[i].y-l)%mod;
			ans %= mod;
			if(s[j].y>s[i].y) r = min(r,s[j].y);
			if(s[j].y<s[i].y) l = max(l,s[j].y);
		}
	} 
	//n*m子矩阵个数就是下面这个公式
	//由n*m+(n-1)*m+(m-1)*n+(n-1)(m-1).....可得 
	cnt = (1ll*n*(n+1)/2%mod)*(1ll*m*(m+1)/2%mod)%mod;
	ans = (cnt-ans+mod)%mod;
	printf("%lld\n",ans);
	return 0;
}