洛谷 1373

      最近要补一补自己的DP了，感觉自己DP好弱...这道题首先我们比较容易设计出转移方程，令f[i][j][k1][k2][opt]表示opt取完了i,j格的魔液，此时小A有k1个，uim有k2个的方案数。那么我们容易写出转移方程，可以用递推的形式来转移，由于每个点都可以是起点，那么初始状态就需要把所有的f[i][j][a[i][j]][0][0]赋值为1，这样做复杂度为n*m*k*k，过不了所有点。怎么办呢？对于一个DP的优化，大概只有2种方式，一个是优化转移，一个是优化状态。那么这题的转移方程没有最值，没有单调性，那么我们只能优化状态。考虑状态f[i][j][k][opt]，k表示小a和uim的收集的差值。这样复杂度就降到n*m*k，就可以AC了

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 805
#define mod 1000000007
int f[maxn][maxn][16][2];
int n,m,k,a[maxn][maxn],ans;

void zy(int &x,int &y) 
{
	x=(x+y)%mod;	
}

int main()
{
	scanf("%d%d%d",&n,&m,&k);k++;
	for (int i=1;i<=n;i++)
		for (int j=1;j<=m;j++) scanf("%d",&a[i][j]);
	for (int i=1;i<=n;i++) 
		for (int j=1;j<=m;j++) f[i][j][a[i][j]][0]=1;
	for (int i=1;i<=n;i++) 
		for (int j=1;j<=m;j++) 
			for (int k1=0;k1<k;k1++)  
			{
				zy(f[i+1][j][(k1+a[i+1][j])%k][0],f[i][j][k1][1]);
				zy(f[i+1][j][(k1-a[i+1][j]+k)%k][1],f[i][j][k1][0]);
				zy(f[i][j+1][(k1+a[i][j+1])%k][0],f[i][j][k1][1]);
				zy(f[i][j+1][(k1-a[i][j+1]+k)%k][1],f[i][j][k1][0]);
			}
	for (int i=1;i<=n;i++) 
		for (int j=1;j<=m;j++) 
			ans=(ans+f[i][j][0][1])%mod;
	printf("%d\n",ans);			
	return 0;	
}