单元格（NHOI2012）

单元格

题目描述：

在一个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

输出：

6

 

样例2：

输入：

3 3 4 7

输出：

0

 

样例3：

输入：

4 6 9 12

输出：

264

 

样例4：

输入：

7 5 13 18

输出：

1212

 

样例5：

输入：

4000 4000 4000 14000

输出：

859690013

 

题解：

       乍一看这道题，我整个人都蒙了。什么？选三个单元格，有几种选法？除了搜索外我想不到其他方法了。但经过大神指点后，我豁然开朗，真是“山重水复疑无路，柳暗花明又一村”哪！大神果然是大神，一语惊醒梦中人哪！

 

       假如我们将每个单元格看做一个点，再将这些点用线连接起来，就会发现：A、B、C三点刚好在一个矩形的边上，并且它的“费用”就是矩形的周长！于是我们就可以这样枚举矩形，计算出它们的边长，将符合条件的矩形中三个的的所有可能的位置数累加，取模，就是最后的答案了。

       问题来了，如何算出每个不同的矩形中三个点位置的方案总数？搜索？暴力搜索肯定会超时，我们试试能不能用数学方法求得不同矩形中三点位置的总数。

       首先我们分类讨论一下三点的位置。

       大致可以分为两大类：

1.两个点分别在矩形对角线所连接的两个顶点上，一个如图：

2.两个点在矩形的边上，另一个点在矩形的顶点上（只有一个点）

 

        最后将两种情况相加。可以得到公式，一个R行C列矩形中三个点的放法总数=（R-2）*（C-2）*6

        这个问题到现在就已经基本上解决了，接下来只需要枚举矩形的长和宽，再用每个矩形中三个点的总放法乘该矩形的总个数，累加即可得到答案。

        这里直接给出求一个 i 行 j 列的矩形 在R行C列的矩形 中数量的公式：数量=(R-i+1)*(C-j+1)

        至于为什么，想必DALAO的您肯定能很快想出来......

 

代码：

 1 //注意！！！
 2 //为了排版美观、增强程序可读性
 3 //此代码省略的丑陋的%mod
 4 //所以，复制代码下来测试或出数据
 5 //一定不要忘记加上%mod
 6 //温馨提示：
 7 //%mod要乘一次就%一次，否则会炸！ 
 8 #include<bits/stdc++.h>
 9 using namespace std;
10 long long r,c,minT,maxT;
11 const long long mod=1000000007;
12 long long ans;
13 int main()
14 {
15     cin>>r>>c>>minT>>maxT;
16     for(long long i=3;i<=r;i++)
17         for(long long j=3;j<=c;j++)
18         {
19             int s=(i-1+j-1)*2;
20             if(s>=minT&&s<=maxT)
21                 ans+=(i-2)*(j-2)*6*(r-i+1)*(c-j+1);
22             else if(s>maxT) break;
23         }
24     cout<<ans<<endl;
25     return 0;
26 }

 

转载于:https://www.cnblogs.com/GDOI2018/p/9029388.html