单元格（模拟）

题目描述
在一个R行C列的表格里，我们要选出3个不同的单元格。但要满足如下的两个条件：
（1）选中的任意两个单元格都不在同一行。
（2）选中的任意两个单元格都不在同一列。
假设我们选中的单元格分别是：A，B，C，那么我们定义这种选择的“费用”= f[A][B] + f[B][C] + f[C][A]。 其中f[A][B]是指单元格A到单元格B的距离，即两个单元格所在行编号的差的绝对值 + 两个单元格所在列编号的差的绝对值。例如：单元格A在第3行第2列，单元格B在第5行第1列，那么f[A][B] = |3-5| + |2-1| = 2 + 1 = 3。至于f[B][C], f[C][A]的意义也是同样的道理。现在你的任务是：有多少种不同的选择方案，使得“费用”不小于给定的数minT，而且不大于给定的数maxT，即“费用”在【minT, maxT】范围内有多少种不同的选择方案。答案模1000000007。所谓的两种不同方案是指：只要它们选中的单元格有一个不同，就认为是不同的方案。

输入
一行，4个整数，R、C、minT、maxT。3≤R,C≤4000, 1≤minT≤maxT≤20000。
对于30%的数据, 3 ≤ R ，C ≤ 70。

输出
一个整数，表示不同的选择方案数量模1000000007后的结果。

样例输入
【样例1】
3 3 1 20000
【样例2】
3 3 4 7
【样例3】
4 6 9 12
【样例4】
7 5 13 18
【样例5】
4000 4000 4000 14000

样例输出
【样例1】
6
【样例2】
0
【样例3】
264
【样例4】
1212
【样例5】
859690013

思路
对于一种选择，其实际答案为2*(maxx-minx+maxy-miny)，即a,b,c三点组成的矩形的周长，a,b,c三点一共有六种方式组成一个相同的矩形，因此枚举每种矩形组成方式求解即可

代码实现

#pragma GCC optimize(3,"Ofast","inline")
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
const int N=1e5+5;
const int M=10005;
const int INF=0x3f3f3f;
const ull sed=31;
const ll mod=1e9+7;
const double eps=1e-8;
const double PI=acos(-1.0);
typedef pair<int,int>P;
 
int r,c,mint,maxt;
 
int main()
{
    scanf("%d%d%d%d",&r,&c,&mint,&maxt);
    ll ans=0;
    for(int i=2;i<=r;i++)
    {
        for(int j=2;j<=c;j++)
        {
            int t=(i+j-2)*2;
            if(t>=mint && t<=maxt) ans=(ans+6*((i-2)*(j-2)%mod)*((r-i+1)*(c-j+1)%mod)%mod)%mod;
        }
    }
    printf("%lld\n",ans);
    return 0;
}