9月10日模拟赛题解

前言

一时爆炸一时爽，一直爆炸一直爽！

因为是教师节，所以中午就放学了，但是我们敬爱的 tc 还给我们准备了一场模拟赛！

$\text{T1}$：签到。直接 $\mathrm{D}\mathrm{P}$。评分 $10$。

$\text{T2}$：签到 $2$。分解一下质因数再取个 $\max$ 就没了（我竟然搞出了 $\mathcal{O}(n+\sum \sqrt{a_i})$ 比题解解快的方法）。评分 $30$。

update：
算法被卡了。

$\text{T3}$：签到 $3$。在 $\mathrm{S}\mathrm{P}\mathrm{F}\mathrm{A}$ 的模板上进行更改即可。评分 $30$.

$\text{T4}$：$\mathrm{D}\mathrm{P}+$ 容斥。评分 $60$。

然后出现了玄学错误，导致从 rk2 掉到 rk5。

$\text{T3}$ 的读入我是这样写的：

for (int i = 1; i <= n; i++)
{
	getchar();
	for (int j = 1; j <= m; j++)
	{
		char c = getchar();

这个在本地测都没问题，但评测机和 $\text{LuoguIDE}$ 都输出了另一个答案。

后来 大桥封印术！ 告诉我千万不要用 $\mathrm{g}\mathrm{e}\mathrm{t}\mathrm{c}\mathrm{h}\mathrm{a}\mathrm{r}$，于是改成了：

for (int i = 1; i <= n; i++)
{
	char c;
	for (int j = 1; j <= m; j++)
	{
		cin >> c;

然后就有 $70$ 分了。

正文

$\text{T1奶酪工厂}$

$\text{Description}$

有一个奶酪工厂，在 $N$ 周里，第 $i$ 周，生产一单位奶酪需要 $C_i$ 便士。有一个货栈，保存一单位奶酪，每周需要 $S$ 便士。货栈可以存无限量的奶酪，而且保证它们不变质。

每周有一个订单，在第 $i$ 周需要交付 $Y_i$ 单位的奶酪，第 $i$ 周生产的奶酪以及之前的存货都可以交付。

请计算完成这 $N$ 周的订单花的最小钱数。

$\text{Solution}$

发现第 $i$ 周的订单可以选择用第 $i$ 周生产的奶酪或之前的存货交付，这是一个递推的过程。

我们设 $d{p}_i$ 表示第 $i$ 周生产一单位奶酪需要的最少钱数，那么第 $i$ 周可以：

1. 用第 $i-1$ 周的存货，存在货仓一周；
2. 直接生产。
   那么方程就是：
   $$d{p}_i=\min (d{p}_{i-1}+s,C_i)$$
   但是注意第一周只能直接生产，所以我们可以初始化 $d{p}_0=INF$。

时间复杂度为 $\mathcal{O}(n)$。

$\text{Code}$

#include <iostream>
#include <cstdio>
#define int long long
using namespace std;

const int MAXN = 10005;

int c[MAXN], y[MAXN], dp[MAXN];

signed main()
{
	int n, s, ans = 0;
	scanf("%lld%lld", &n, &s);
	dp[0] = 0x3f3f3f3f3f3f3f3fll;
	for (int i = 1; i <= n; i++)
	{
		scanf("%lld%lld", c + i, y + i);
		dp[i] = min(dp[i - 1] + s, c[i]);
		ans += dp[i] * y[i];
	}
	printf("%lld\n", ans);
	return 0;
}

$\text{T2阶乘}$

$\text{Description}$

有 $n$ 个数 $a_1,a_2\dots ,a_n\in {\mathbb{Z}}^{+}$，设它们的积为 $p$，你可以找到一个 $q\in {\mathbb{Z}}^{+}$，使 $p\times q$ 刚好为 $m!$（$m\in {\mathbb{Z}}^{+}$），求 $m$ 的最小值。

$\text{Solution}$

题解用了玄学二分，复杂度比我的方法多了个 $\log$，然而因为 $n$ 很小所以显不出优势。

和 大桥封印术！ 交换思路后发现两人凑一凑就能凑出答案。

先在读入时将 $p$ 分解质因数，假设长这样子：$p=2^a\times 3^b\times 5^c\times \cdots$，设 $m_1$ 为满足 $m!$ 中含 $2^a$ 的 $m$ 的最小值，$m_2$ 为满足 $m!$ 中含 $3^b$ 的 $m$ 的最小值，$m_3$ 为满足 $m!$ 中含 $5^c$ 的 $m$ 的最小值……

那么显然当 $m$ 取 max ⁡ { m i } \max\{m_i\} max{mi​} 时，$m$ 是同时含所有的最小值。

于是问题转化为如何求 $m_i$。

现在要让 $m_i!$ 中含 $t^k$，求 $m_i$ 的最小值，我们来推推柿子：

$\begin{array}{rl}k& =\lfloor \frac{m_i}{t}\rfloor +\lfloor \frac{m_i}{t^2}\rfloor +\lfloor \frac{m_i}{t^3}\rfloor +\cdots \\ & \le \frac{m_i}{t}+\frac{m_i}{t^2}+\frac{m_i}{t^3}+\cdots \\ & =\frac{m_i}{t-1}\end{array}$

$\therefore k\le \frac{m_i}{t-1}$

现在是知道 $t,k$ 要求 $m_i$。

移项得 $m_i\ge (t-1)k$

但这个可能有 bug，以 $t=5,k=3$ 为例：此时算出来的 $m_i=(5-1)\times 3=12$。

这个显然是错的，$12$ 本来就不含质因数 $5$，怎么可能是最小的？

可以发现 $12!$ 只含 $5^2$，而 $15!$ 含 $5^3$，才是正确的 $m_i$。

我们设 $x=(t-1)k$，那么 $m_i$ 应该是不小于 $x$ 的数中最小的 $t$ 的倍数，我们直接用 $x+t-x\mathrm{m}\mathrm{o}\mathrm{d}t$ 即可。

分解质因数用 $\mathcal{O}(n)+\sum \sqrt{a_i}$，求 $m_i$ 用 $o(n)$，总时间复杂度为 $\mathcal{O}(n+\sum \sqrt{a_i})$。

$\text{Code}$

#include <iostream>
#include <cstdio>
using namespace std;

const int MAXN = 1e5 + 5;

int tot;
int a[MAXN], p[MAXN], cnt[MAXN], x[MAXN];
bool vis[MAXN];

void pre() //预处理
{
	for (int i = 2; i <= MAXN; i++)
	{
		if (!vis[i])
		{
			p[++tot] = i;
			x[i] = i;
		}
		for (int j = 1; j <= tot && i * p[j] <= MAXN; j++)
		{
			vis[i * p[j]] = true;
			x[i * p[j]] = p[j];
			if (i % p[j] == 0)
			{
				break;
			}
		}
	}
}

void dyy(int n) //分解质因数
{
	while (n != 1)
	{
		cnt[x[n]]++;
		n /= x[n];
	}
}

int main()
{
	int n, maxx = 0, ans = 0;
	scanf("%d", &n);
	pre();
	for (int i = 1; i <= n; i++)
	{
		scanf("%d", a + i);
		maxx = max(maxx, a[i]);
		dyy(a[i]);
	}
	for (int i = 2; i <= maxx; i++)
	{
		if (cnt[i])
		{
			int k = (i - 1) * cnt[i];
			ans = max(ans, k + i - k % i);
		}
	}
	printf("%d\n", ans);
	return 0;
}

update：
做法被 hack 了，取到不小于 $x$ 的数中最小的 $t$ 的倍数还是有可能不够，例如：

Input
5
2 2 2 2 2
Output
6
Answer
8

$\text{Solution}$

所以还是得用二分。

上界为 $\sum a_i$，二分 $m$，$\mathrm{check}(mid)$ 判断每个质因数的个数够不够。

时间复杂度为 $\mathcal{O}(\log [{\sum }_{1\le i\le n}{a}_i\times {\sum }_{1\le 质数p\le mid}(\lfloor \frac{mid}{p}\rfloor +\lfloor \frac{mid}{p^2}\rfloor +\lfloor \frac{mid}{p^3}\rfloor +\cdots )])\approx \mathcal{O}(\log ({\sum }_{1\le i\le n}{a}_i\times {\sum }_{1\le 质数p\le mid}\lfloor \frac{mid}{p-1}\rfloor ))$。

$\text{Code}$

#include <iostream>
#include <cstdio>
#define int long long
using namespace std;

const int MAXN = 1e5 + 5;

int tot;
int a[MAXN], p[MAXN], cnt[MAXN];
bool vis[MAXN];

void pre(int n)
{
	for (int i = 2; i <= n; i++)
	{
		if (!vis[i])
		{
			p[++tot] = i;
		}
		for (int j = 1; j <= tot && i * p[j] <= n; j++)
		{
			vis[i * p[j]] = true;
			if (i % p[j] == 0)
			{
				break;
			}
		}
	}
}

void dyy(int n) 
{
	for (int i = 1; p[i] * p[i] <= n; i++)
	{
		while (n % p[i] == 0)
		{
			n /= p[i];
			cnt[p[i]]++;
		}
	}
	if (n != 1)
	{
		cnt[n]++;
	}
}

bool check(int n)
{
	for (int i = 1; i <= tot; i++)
	{
		int t = n, res = 0;
		while (t > 1)
		{
			t /= p[i];
			res += t;
		}
		if (res < cnt[p[i]])
		{
			return false;
		}
	}
	return true;
}

signed main()
{
	int n, sum = 0;
	scanf("%lld", &n);
	pre(n);
	for (int i = 1; i <= n; i++)
	{
		scanf("%lld", a + i);
		dyy(a[i]);
		sum += a[i];
	}
	int l = 1, r = sum, ans;
	while (l <= r)
	{
		int mid = (l + r) >> 1;
		if (check(mid))
		{
			ans = mid;
			r = mid - 1;
		}
		else
		{
			l = mid + 1;
		}
	}
	printf("%lld\n", ans);
	return 0;
}

$\text{T3小S练跑步}$

$\text{Description}$

有一张地图，每一个位置上都写了到达该位置后不能往哪一个方向移动。有 $5$ 种表示的符号：U代表不能向上移动，D代表不能向下移动，L代表不能向左移动，R代表不能向右移动，S代表该位置有障碍物，到达那里后无法继续训练。整个公园共有 $n$ 行 $m$ 列，要从第 $1$ 行第 $1$ 列出发，到第 $n$ 行第 $m$ 列结束练习。

求从起点到终点途中最少要改变几次方向。

$\text{Solution}$

直接用 $\mathrm{S}\mathrm{P}\mathrm{F}\mathrm{A}$ 搜，用 $lst$ 记录当前方向，注意更改了松弛条件，$dis$ 数组也需要开成三维。

时间复杂度为 $\mathcal{O}(nm)$。

$\text{Code}$

#include <iostream>
#include <cstdio>
#include <queue>
#include <cstring>
using namespace std;

const int MAXN = 505;

int n, m;
int mp[MAXN][MAXN], fx[5][2] = {{0, 0}, {-1, 0}, {1, 0}, {0, -1}, {0, 1}}, dis[MAXN][MAXN][5];
bool vis[MAXN][MAXN][5];

struct que
{
	int x, y, lst;
};

void spfa()
{
	memset(dis, 0x3f, sizeof(dis));
	queue<que> q;
	for (int i = 1; i <= 4; i++)
	{
		dis[1][1][i] = 0;
		q.push(que{1, 1, i});
	}
	while (!q.empty())
	{
		int x = q.front().x, y = q.front().y, lst = q.front().lst;
		q.pop();
		vis[x][y][lst] = false;
		if (mp[x][y] == -1)
		{
			continue;
		}
		for (int i = 1; i <= 4; i++)
		{
			int tx = x + fx[i][0], ty = y + fx[i][1];
			if (tx < 1 || ty < 1 || tx > n || ty > m || mp[x][y] == i)
			{
				continue;
			}
			if (dis[x][y][lst] + (lst != i) < dis[tx][ty][i])
			{
				dis[tx][ty][i] = dis[x][y][lst] + (lst != i);
				if (!vis[tx][ty][i])
				{
					vis[tx][ty][i] = true;
					q.push(que{tx, ty, i});
				}
			}
		}
	}
}

int main()
{
	scanf("%d%d", &n, &m);
	for (int i = 1; i <= n; i++)
	{
		char c;
		for (int j = 1; j <= m; j++)
		{
			cin >> c;
			if (c == 'S')
			{
				mp[i][j] = -1;
			}
			else if (c == 'U')
			{
				mp[i][j] = 1;
			}
			else if (c == 'D')
			{
				mp[i][j] = 2;
			}
			else if (c == 'L')
			{
				mp[i][j] = 3;
			}
			else
			{
				mp[i][j] = 4;
			}
		}
	}
	spfa();
	int ans = 0x3f3f3f3f;
	for (int i = 1; i <= 4; i++)
	{
		ans = min(ans, dis[n][m][i]);
	}
	if (ans == 0x3f3f3f3f)
	{
		puts("No Solution");
	}
	else
	{
		printf("%d\n", ans);
	}
	return 0;
}

$\text{T4轨道}$

$\text{Description}$

对于公式
$$v=\frac{a_1\times a_2\times \cdots \times a_n}{k}$$
给出 $n,m,k$，求这 $n$ 个正整数 $a_1,a_2,\dots ,a_3$ 在都不大于 $m$ 的情况下有多少种选择方式，使得 $v$ 为与 $k$ 互质的正整数，答案对 $10007$ 取模。

$\text{Solution}$

$\text{Code}$