瞬间,地面上出现了一个n*m的巨幅矩阵,矩阵的每个格子上有一坨0~k不等量的魔液。怪物各给了小a和uim一个魔瓶,说道,你们可以从矩阵的任一个格子开始,每次向右或向下走一步,从任一个格子结束。开始时小a用魔瓶吸收地面上的魔液,下一步由uim吸收,如此交替下去,并且要求最后一步必须由uim吸收。魔瓶只有k的容量,也就是说,如果装了k+1那么魔瓶会被清空成零,如果装了k+2就只剩下1,依次类推。怪物还说道,最后谁的魔瓶装的魔液多,谁就能活下来。小a和uim感情深厚,情同手足,怎能忍心让小伙伴离自己而去呢?沉默片刻,小a灵机一动,如果他俩的魔瓶中魔液一样多,不就都能活下来了吗?小a和他的小伙伴都笑呆了!
现在他想知道他们都能活下来有多少种方法。
输入格式:
第一行,三个空格隔开的整数n,m,k
接下来n行,m列,表示矩阵每一个的魔液量。同一行的数字用空格隔开。
输出格式:
一个整数,表示方法数。由于可能很大,输出对1 000 000 007取余后的结果。
输入样例#1:
2 2 3 1 1 1 1
输出样例#1:
4
我们用dp[i][j][l][0/1]表示坐标i,j时魔瓶里有l,0表示a,1表示uims,然后推出状态转移方程即可。
要注意魔瓶的上限为k+1而不是k。
#include<bits/stdc++.h>
#define mod 1000000007
using namespace std;
int Map[810][810],dp[810][810][16][2];
int n,m,k;
int main()
{
cin>>n>>m>>k;
k++;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
cin>>Map[i][j];
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
dp[i][j][Map[i][j]%k][0]=1;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
for(int l=0;l<k;l++)
{
dp[i][j][l][1]=(dp[i][j][l][1]+dp[i][j-1][(l+Map[i][j]+2*k)%k][0]+dp[i-1][j][(l+Map[i][j]+2*k)%k][0])%mod;
dp[i][j][l][0]=(dp[i][j][l][0]+dp[i][j-1][(l-Map[i][j]+2*k)%k][1]+dp[i-1][j][(l-Map[i][j]+2*k)%k][1])%mod;
}
long long Max=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=m;j++)
Max+=(dp[i][j][0][1])%mod;
cout<<Max%mod<<endl;
return 0;
}