题目链接:CSU-ACM 2164
2018
Description
Bobo 想统计满足下面条件的矩阵 A 的数量。
矩阵 A 有 n 行 m 列,每个元素都是正整数。第 i 行第 j 列的元素用 Ai, j 表示。
A1, 1 = 2018.
对于所有 2 ≤ i ≤ n, 1 ≤ j ≤ m,Ai, j 是 Ai − 1, j 的约数。
对于所有 1 ≤ i ≤ n, 2 ≤ j ≤ m,Ai, j 是 Ai, j − 1 的约数。
因为满足条件的矩阵 A 数量很多,Bobo 只想统计满足条件的矩阵数量除以 (10e9 + 7) 的余数。
1 ≤ n, m ≤ 2000
数据组数不超过 10e5.
Input
输入文件包含多组数据,请处理到文件结束。
每组数据包含 2 个整数 n 和 m.
Output
对于每组数据输出 1 个整数表示所求的数量除以 (109 + 7) 的余数。
Sample Input
1 1
1 2
2 2
2 3
2000 2000
Sample Output
1
4
25
81
570806941
Hint
对于第二组样例(n = 1, m = 2),满足条件的矩阵 A 有 (2018, 2018),(2018, 1009),(2018, 2),(2018, 1) 共 4 种。
Source
2018湖南省第14届大学生计算机程序设计竞赛
扫描二维码关注公众号,回复:
3382224 查看本文章
Author
ftiasch
题解:不知道别人打表干嘛,显然是组合数学规律,组合数的底数是n+m,上面的数随便取n或者m就行。详见AC代码
// 组合数
#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int maxn = 4005; ///
const int M = 1000000007;
LL fac[maxn + 1] = {1, 1}, inv[maxn + 1] = {1, 1}, f[maxn + 1] = {1, 1};
LL
C(LL a, LL b) {
return (fac[a] * inv[b] % M * inv[a - b] % M) % M;
}
void
init() {
for(int i = 2; i < maxn; i++) {
fac[i] = fac[i - 1] * i % M;
f[i] = (M - M / i) * f[M % i] % M;
inv[i] = inv[i - 1] * f[i] % M;
}
}
int
main() {
int n, m;
init();
while( ~scanf("%d %d", &n, &m) ) {
LL ans = ((C(n + m, n) - 1 + M) % M * (C(n + m, n) - 1 + M ) % M) % M;
printf("%lld\n", ans);
}
return 0;
}