传送门
思路
吐槽
吐槽:话说不应该是他俩都得死嘛qwq
咋做
一道
D
P
DP
DP不好题.
我们用 f [ i ] [ j ] [ q ] [ p ] f[i][j][q][p] f[i][j][q][p]来表示:走到第 i i i行第 j j j列魔液差距值为 q q q且当前为 p p p走的方案数( p ∈ { 0 , 1 } p = 0 p\in \{0, 1\} p=0 p∈{0,1}p=0表示目前为小 a a a走, p = 1 p=1 p=1表示目前为 u i m uim uim走)
初始条件: f [ i ] [ j ] [ a [ i ] [ j ] ] [ 0 ] = 1 f[i][j][a[i][j]][0]=1 f[i][j][a[i][j]][0]=1表示小 a a a从每个点开始取,差距值为 a [ i ] [ j ] a[i][j] a[i][j]的方案数为 1 1 1
那么我们可以想出转移方程:
f [ i ] [ j ] [ p ] [ 0 ] + = f [ i − 1 ] [ j ] [ p − a [ i ] [ j ] ] [ 1 ] + f [ i ] [ j − 1 ] [ p − a [ i ] [ j ] ] [ 1 ] f[i][j][p][0]+=f[i-1][j][p-a[i][j]][1]+f[i][j-1][p-a[i][j]][1] f[i][j][p][0]+=f[i−1][j][p−a[i][j]][1]+f[i][j−1][p−a[i][j]][1]
这个式子表示目前在第 i i i行第 j j j列,差距值为 p p p,当前小 a a a走的方案数,因为只能往右走或往下走,且上一步一定是 u i m uim uim走的,所以可以从第 i − 1 i-1 i−1行第 j j j列、第 i i i行第 j − 1 j-1 j−1列转移过来,差距值增大
f [ i ] [ j ] [ p ] [ 1 ] + = f [ i − 1 ] [ j ] [ p + a [ i ] [ j ] ] [ 0 ] + f [ i ] [ j − 1 ] [ p − a [ i ] [ j ] ] [ 0 ] f[i][j][p][1]+=f[i-1][j][p+a[i][j]][0]+f[i][j-1][p-a[i][j]][0] f[i][j][p][1]+=f[i−1][j][p+a[i][j]][0]+f[i][j−1][p−a[i][j]][0]
同理,这个式子表示目前在第 i i i行第 j j j列,差距值为 p p p,当前 u i m uim uim走的方案数,上一步一定是小 a a a走的,所以可以从第 i − 1 i-1 i−1行第 j j j列、第 i i i行第 j − 1 j-1 j−1列转移过来,差距值减小
我们容易想出,最后的答案就是
∑
i
=
1
n
∑
j
=
1
m
f
[
i
]
[
j
]
[
0
]
[
1
]
\sum_{i = 1}^{n}\sum_{j = 1}^{m}f[i][j][0][1]
∑i=1n∑j=1mf[i][j][0][1]最后一维是
1
1
1是因为最后一步只能由
u
i
m
uim
uim走
空间复杂度 n ∗ m ∗ k ∗ 2 n*m*k*2 n∗m∗k∗2刚好可以,时间复杂度 O ( n ∗ m ∗ k ) O(n*m*k) O(n∗m∗k)可以过
注意!
- 因为魔液值到达 k + 1 k+1 k+1之后就会清 0 0 0,所以在过程中要随时对 k + 1 k+1 k+1取模
- 空间开好!!不要开太大!因为最后一维只用 0 、 1 0、1 0、1,所以开 2 2 2的数组就好了就行了,开 3 3 3会 M L E MLE MLE
- 码风丑陋,谨慎阅读 q w q qwq qwq
代码
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
const int B = 20;
const int A = 800 + 7;
const int mod = 1e9 + 7;
const int inf = 0x3f3f3f3f;
inline int read() {
char c = getchar(); int x = 0, f = 1;
for( ; !isdigit(c); c = getchar()) if(c == '-') f = -1;
for( ; isdigit(c); c = getchar()) x = x * 10 + (c ^ 48);
return x * f;
}
int n, m, k, a[A][A], f[A][A][B][2];
int main() {
n = read(), m = read(), k = read() + 1;
for(int i = 1; i <= n; i++) for(int j = 1; j <= m; j++) a[i][j] = read(), f[i][j][a[i][j] % k][0] = 1;
for(int i = 1; i <= n; i++)
for(int j = 1; j <= m; j++)
for(int p = 0; p <= k; p++) {
f[i][j][p][0] = (f[i][j][p][0] + f[i - 1][j][(p - a[i][j] + k) % k][1] + f[i][j - 1][(p - a[i][j] + k) % k][1]) % mod;
f[i][j][p][1] = (f[i][j][p][1] + f[i - 1][j][(p + a[i][j] + k) % k][0] + f[i][j - 1][(p + a[i][j] + k) % k][0]) % mod;
}
int ans = 0;
for(int i = 1; i <= n; i++) for(int j = 1; j <= m; j++) ans += f[i][j][0][1], ans %= mod;
cout << ans << '\n';
return 0;
}