【题解】LuoGu5307：[COCI2019] Mobitel

最新推荐文章于 2022-08-09 20:22:23 发布

ModestCoder_

最新推荐文章于 2022-08-09 20:22:23 发布

阅读量339

点赞数 1

分类专栏：题解 LuoGu DP 文章标签：题解 LuoGu DP

本文链接：https://blog.csdn.net/ModestCoder_/article/details/99936669

版权

题解同时被 3 个专栏收录

435 篇文章 2 订阅

订阅专栏

LuoGu

255 篇文章 0 订阅

订阅专栏

125 篇文章 0 订阅

订阅专栏

原题传送门
一个很普通的dp
$dp_{i,j,k}$ 表示第i行第j列乘积为k的路径条数
状态冗余太多，并且不优，需要优化一下
首先可以滚掉 $i$ 这一维，接下来考虑把 $k$ 这一维优化掉

一个定理 $[\frac{[\frac{n}{x}]}{y}]=[\frac{n}{xy}]$
dp状态修改为 $dp_{i,j,k}$ 表示第i行第j列，还需要乘上k才能>=n
这个k的值域是n除以所有数上取整
然后就可以愉快地dp了

Code：

#include <bits/stdc++.h>
#define maxn 310
#define maxm 2010
#define qy 1000000007
using namespace std;
int r, s, n, a[maxn][maxn], b[1000010], tot, val[maxm], pos[1000010], dp[2][maxn][maxm];

inline int read(){
	int s = 0, w = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') w = -1;
	for (; isdigit(c); c = getchar()) s = (s << 1) + (s << 3) + (c ^ 48);
	return s * w;
}

int calc(int x, int y){ return x % y == 0 ? x / y : x / y + 1; }
void upd(int &x, int y){ if ((x += y) >= qy) x -= qy; }

int main(){
	r = read(), s = read(), n = read();
	for (int i = 1; i <= r; ++i)
		for (int j = 1; j <= s; ++j) a[i][j] = read();
	for (int i = 1; i <= n; ++i) b[i] = calc(n, i);
	for (int i = 1; i <= n; ++i)
		if (b[i] != b[i - 1]) pos[val[++tot] = b[i]] = tot;
	dp[1][1][pos[calc(n, a[1][1])]] = 1; 
	for (int i = 1; i <= r; ++i){
		int now = i & 1, nxt = now ^ 1;
		for (int j = 1; j <= s; ++j)
			for (int k = 1; k <= tot; ++k){ //printf("%d\n", dp[now][j][k]);
				if (i != r) upd(dp[nxt][j][pos[calc(val[k], a[i + 1][j])]], dp[now][j][k]);
				if (j != s) upd(dp[now][j + 1][pos[calc(val[k], a[i][j + 1])]], dp[now][j][k]);
				if (i != r || j != s || k != tot) dp[now][j][k] = 0;
			}
	}
	printf("%d\n", dp[r & 1][s][tot]);
	return 0;
}