新生赛模拟骗

zstu新生赛 p6
前言：
如果忘记前缀和（esp二维）与差分 可以去
https://blog.csdn.net/k_r_forever/article/details/81775899 学习
1414: 曼哈顿距离
时间限制: 1 Sec 内存限制: 128 MB
提交: 9 解决: 4
[提交] [状态] [讨论版] [命题人:admin]
题目描述
在平面上，坐标（x1, y1）的点P1与坐标（x2, y2）的点P2的曼哈顿距离为：|x1-x2|+|y1-y2|。 现在有一个nXn (1<=n<=1000)的矩阵，一开始矩阵中每个元素的值都为零。 对这个矩阵进行m (1<=m<=100000)次操作，每次操作定义如下：
(x1, y1), (x2, y2): 对于矩阵中满足x1<=x<=x2, y1<=y<=y2的所有元素(x, y)加上(x, y)与(x1, y1)的曼哈顿距离|x1-x|+|y1-y|
例如，当n=5，m=2，两组操作分别为(1, 1), (4, 4)和(2, 2), (5, 5):

现在想让你输出m次操作后的矩阵，但是由于本OJ输出文件不能超过1M，所以输出矩阵所有非零元素的乘积即可 （结果对1000000007取余）,没有非零项则输出0；

输入
第一行两个整数n， m，分别表示正方形网格的边长和操作次数；
接下来m行，每行4个整数x1 y1 x2 y2；

输出
输出一个整数代表答案。

样例输入
5 2
1 1 4 4
2 2 5 5

样例输出
853543927
显然需要用到二维平面的差分思想来完成
(x1, y1), (x2, y2)对于矩阵中满足x1<=x<=x2, y1<=y<=y2的所有元素加一个数的操作
但这个数并不是常值，给了这道题特殊性
一种巧妙的处理方法：
|x1-x|+|y1-y|=(x1+y1)-(x+y)
元素终值=|x1-x|+|y1-y|+|x1-x.1|+|y1-y.1|+…+|x1-x.n|+|y1-y.n|=(x1+y1)*n-（求和）(x.i+y.j)
（1<=i<=n;1<=j<=n)
我们设置两个差分数组，一个用于确定n的值，一个用于确定(x.i+y.i)这部分的值
因此有dm[?][?]+=1 dm_spec[?][?]+=x1+y1
Code：
#include
using namespace std;
int dm[1010][1010],dm_spec[1010][1010];
int main()
{
int n,m,i,j,x1,y1,x2,y2;
cin>>n>>m;
while(m–)
{
cin>>x1>>y1>>x2>>y2;
dm[x1][y1]+=1;
dm[x2+1][y1]-=1;
dm[x1][y2+1]-=1;
dm[x2+1][y2+1]+=1;
dm_spec[x1][y1]+=x1+y1;
dm_spec[x2+1][y1]-=x1+y1;
dm_spec[x1][y2+1]-=x1+y1;
dm_spec[x2+1][y2+1]+=x1+y1;
}
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
{
dm[i][j]+=dm[i][j-1]+dm[i-1][j]-dm[i-1][j-1];
dm_spec[i][j]+=dm_spec[i][j-1]+dm_spec[i-1][j]-dm_spec[i-1][j-1];
}
long long ans=1,flag=0;
for(i=1;i<=n;i++)
for(j=1;j<=n;j++)
{
int temp=-dm_spec[i][j]+(i+j)dm[i][j];
if(temp>0)
{
flag=1;
ans=anstemp%1000000007;
}
}
if(flag)cout<<ans<<endl;else cout<<0<<endl;
}