P2049 魔术棋子 题解

CSDN同步

原题链接

简要题意：

在一个 (n imes m) 的魔术棋盘中，每个格子中均有一个整数，当棋子走进这个格子中，则此棋子上的数会被乘以此格子中的数。一个棋子从左上角走到右下角，只能向右或向下行动，请问此棋子走到右下角后，模 (\% k) 可以为几？

（原题题意足够简要了吧）

(n,m,k leq 100).

考虑一个很朴素的做法，(mathcal{O}(nmk)) 的那种。

很显然我们不可能算出所有路径的答案（(C_{n+m}^m) 种的级别大家应该都清楚），所以说我们可以从答案入手。

枚举一个答案，看它能不能是合法的答案。

对于 ((i,j) (i > 1 , j > 1)) 点，其数为 (a_{i,j})，如何验证 (a_{i,j} ightarrow l)（表示走到 (a_{i,j}) 的时候数字为 $l）的正确？

考虑前一步的走法，很显然，考虑是否有 (t imes a_{i,j} \% k = l) 且 (t) 是 (a_{i-1,j}) 或 (a_{i,j-1}) 的合法答案。

好，下一步我们考虑，(t) 是否需要用逆元计算？倒推需要逆元，我们可以考虑正推。用当前的答案去更新之后的答案，也是所谓 动态规划 的另一种（不太常见的）方式。

所以说 动态规划 不止可以是 f[i] = ... 的形式，可以是 ... = f[i] ... 的形式，采用“我从哪里来”“我到哪里去”两种均可。

下面就很简单了。初步的方案，用 (h_{i,j,t}) 表示走到 (a_{i,j}) 数字能否为 (t). 一开始只有 (h_{1,1,a_{1,1}} = 1)，考虑如何转移？很简单。

[h_{i,j,t imes a_{i,j} \% k} = max(h_{i,j-1,t , h_{i-1,j,t}}) ]

当前想要表示出 (t imes a_{i,j} \% k) 就去验证 (k)，避免了逆元和转移中记忆化的技巧 ，直奔主题。

时间复杂度：(mathcal{O}(nmk)).

#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;

const int N=1e2+1;

inline int read(){char ch=getchar(); int f=1; while(ch<'0' || ch>'9') {if(ch=='-') f=-f; ch=getchar();}
	int x=0; while(ch>='0' && ch<='9') x=(x<<3)+(x<<1)+ch-'0',ch=getchar(); return x*f;}

inline void write(int x) {
	if(x<0) {putchar('-');write(-x);return;}
	if(x<10) {putchar(char(x%10+'0'));return;}
	write(x/10);putchar(char(x%10+'0'));
}

int n,m,k,a[N][N];
bool h[N][N][N];

int main() {
	n=read(); m=read(); k=read();
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++) a[i][j]=read()%k;
	h[1][1][a[1][1]]=1;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++) {
		if(i==1 && j==1) continue;
		for(int t=0;t<k;t++) h[i][j][t*a[i][j]%k]|=h[i-1][j][t];
		for(int t=0;t<k;t++) h[i][j][t*a[i][j]%k]|=h[i][j-1][t];
	} 
/*	for(int i=1;i<=n;i++)
	for(int j=1;j<=m;j++) {
		printf("%d %d : ",i,j);
		for(int t=0;t<k;t++) if(h[i][j][t]) printf("%d ",t);
		putchar('
');
	}*/
	vector<int> v;
	for(int t=0;t<k;t++) if(h[n][m][t]) v.push_back(t);
	printf("%d
",v.size());
	for(int i=0;i<v.size();i++) printf("%d ",v[i]); putchar('
');
	return 0;
}