「GMOI」Round 1

P8924 「GMOI R1-T1」Perfect Math Class

思路

1. 因为 $n$ 和 $m$ 均不超过 $100$，所以我们用 $vis[i][j]$ 表示第 $i$ 行第 $j$ 列是否函数上。当然，用 map 去维护点的信息也可以。
2. 用因为 $k\le 7$，故计算函数值 $y$ 的时候要开 long long。

代码

#include <bits/stdc++.h>
using namespace std;
const int maxk = 15;
map <pair <int, int>, int> vis; 
int n, m, k, a[maxk];
int main ()
{
	scanf ("%d %d %d", &n, &m, &k);
	for (int i = 1; i <= k + 1; i++)
		scanf ("%d", &a[i - 1]);
	for (int x = 0; x <= n - 1; x ++) 
	{
		long long y = 0, p = 1;
		for (int i = 1; i <= k + 1; i++) 
		{
			y += p * a[i - 1];
			p *= x;
		}
		if (y < 0 || y > m - 1) continue;
		vis.insert (make_pair (make_pair (x, y), 1)); 
	}
	for (int i = m - 1; i >= 0; i--) 
	{
		for (int j = 0; j <= n - 1; j++) 
			if (vis.find (make_pair (j, i)) == vis.end ())
			printf (".");
			else printf ("*");
		printf ("\n");
	}
	return 0;
}

P8925 「GMOI R1-T2」Light

思路

1. 我们设左边的第 $n$ 项是 $a_n$，右边第 $n$ 项是 $b_n$。易得，$a_1=-2L$，$b_1=2R$。
2. 根据题意，$\frac{a_n+b_{n-1}}{2}=-L$ 且 $\frac{b_n+a_{n-1}}{2}=R$。故易得，$a_n-a_{n-2}=-2R-2L$。分奇偶讨论，当 $n$ 为奇数时，$a_n=-(n+1)L-(n-1)R$；当 $n$ 为偶数时，$a_n=-n(R+L)$。$b_n$ 同理可得。

代码

#include <bits/stdc++.h>
using namespace std;
long long l, r, n, ans;
int T; 
int main ()
{
	scanf ("%d", &T);
	scanf ("%lld %lld", &l, &r);
	while (T --) 
	{
		char turn; 
		cin >> turn, scanf ("%lld", &n);
		if (turn == 'L') 
		{
			if (n & 1ll) ans = -(n + 1) * l - (n - 1) * r;
			else ans = -n * l - n * r;
		} 
		else 
		{
			if (n & 1ll) ans = (n + 1) * r + (n - 1) * l;
			else ans = n * r + n * l;
		}
		printf ("%lld\n", ans);
	}
	return 0;
}

P8926 「GMOI R1-T3」Number Pair

思路

1. 设 $x=\gcd (x,y)x_1$ 与 $y=\gcd (x,y)y_1$，则易得 $\gcd (x_1,y_1)=1$。由 $k\gcd (x,y)=\mathrm{lcm}(x,y)$ 得，$k{\gcd (x,y)}^2=xy$，进而 $k=x_1{y}_1$，故欲求满足条件的 $(x,y)$ 的对数，我们需要找到 $k$ 分解成两个互质的数相乘的方案数。
2. 根据唯一分解定理把 $k$ 分解成 $p_1^{r_1}{p}_2^{r_2}\cdots p_n^{r^n}$。若 $p_i\mid x_1$，则 $p_i\nmid y_1$。故可知，$k$ 分解成两个互质的数相乘的方案数是 $2^n$，可由快速幂得到。
3. 如何得到 $k$ 的标准分解式。因为 p i , i ∈ { 1 , 2 … n } p_i,i\in \{1,2\dots n\} pi​,i∈{1,2…n} 中至多有一个大于 $\sqrt{k}$ 的素数，所以我们可以先通过线性筛得到 $1$ 到 $10^8$ 内所有的素数。因为 $1$ 到 $10^8$ 中只有 $28724216$ 个素数，故我们要卡着这个值开数组，避免MLE。我们将 $n$ 初始化为 $0$，然后枚举这些素数，如果某个素数可以整除 $k$，那么 $n+1$，并不断用这个素数除 $k$ 直到这个素数 $\nmid k$。若退出循环后 $k\ne 1$，则说明 $k$ 的标准分解式中有一个大于 $\sqrt{k}$ 的素数，那么 $n+1$。
4. 因为 $P\le \gcd (x,y)\le Q$，所以最后的答案就是 $(Q-P+1)2^n$。
5. 但是这个题卡常了，于是我们要优化。我们发现在枚举素数进行试除的时候，若某个素数 $p$ 满足了条件 $p^2>k$，则没有必要在枚举了，这时候退出循环就可以了。
6. 最坏时间复杂度为 $\mathcal{O}(10^8+28724216T)$。

代码

#include <bits/stdc++.h>
using namespace std;
const int MOD = 1e9 + 7;
const int maxn = 1e8 + 10;
int t, p, q, pi, zs[28724221], ans;
bool vis[maxn];
long long k; 
void Init () 
{
	for (int i = 2; i <= 1e8; i++) 
	{
		if (!vis[i]) zs[++ pi] = i;
		for (int j = 1; j <= pi; j++) 
		{
			if (1ll * zs[j] * i > 1e8) break;
			vis[zs[j] * i] = 1;
			if (i % zs[j] == 0) break;
		}
	}
}
int quickpower (int b) 
{
	if (b == 0) return 1;
	if (b == 1) return 2;
	int mid = quickpower (b / 2);
	if (b & 1) return 2ll * mid * mid % MOD;
	else return 1ll * mid * mid % MOD; 
}
int main ()
{
	Init ();
	scanf ("%d", &t);
	while (t --) 
	{
		scanf ("%lld %d %d", &k, &p, &q);
		int n = 0;
		for (int i = 1; i <= pi && k != 1; i++) 
		{
			if (1ll * zs[i] * zs[i] > k) break;
			if (k % zs[i] == 0) 
			{
				n ++;
				while (k % zs[i] == 0 && k != 1) k /= zs[i];
			}
		}
		if (k != 1) n ++;
		ans = 1ll * quickpower (n) * (q - p + 1) % MOD;
		printf ("%d\n", ans);
	}
	return 0;
}

收获

1. 知道了自己以前写的线性筛都是假的。

2. 第一次碰见卡空间复杂度的，$256MB$ 是 $268435456$ 字节。下图是常见数据类型的所占字节数，bool 型是 $1$ 字节。
   

3. 意想不到的卡常优化（思路 5）。

P8927 「GMOI R1-T4」Rain

思路

经过绝对值运算后，将会有 $n$ 项符号为正，$n$ 项符号为负。为了使答案尽可能地大，我们希望，$p{a}_1,p{a}_2\dots p{a}_n,q{a}_1,q{a}_2\dots q{a}_n$ 中前 $n$ 大项的符号全为正，其余较小的项符号全为负。

以 $n=2k$ 和 $p\ge q$ 的情况为例，我们对 $a_1,a_2\dots a_{2k}$ 进行由大到小的排序，得到新的 $a_1,a_2\dots a_{2k}$。令 $b_{1,2\dots 2k}$ 为 $a_{1,2k,2,2k-1,\dots k,k+1}$。下面证明序列 $b$ 满足期望。

先不考虑从 $2k$ 到 $1$ 产生的距离，对其他距离进行讨论。当 $i$ 为奇数时，$|bi-b_{i+1}|=p{b}_i-q{b}_{i+1}$；当 $i$ 为偶数时，由于 $p{b}_i$ 递增且 $q_{i+1}$ 递减，设 $i_0$（如果存在的话），使得 $p{b}_{i_0}\ge q{b}_{i_0+1}$，则当 $i<i_0$ 时，$|p{b}_i-q{b}_{i+1}|=q{b}_{i+1}-p{b}_i$；当 $i\ge i_0$ 时，$|p{b}_i-q{b}_{i+1}|=p{b}_i-q{b}_{i+1}$。

故符号为正的 $2k-1$ 项为，$p{b}_{1,3,\dots ,2k-1},q{b}_{3,5,\dots ,i_0-1},p{b}_{i_0,\dots ,2k-2}$。
即，$p{a}_{1,2,\dots ,k,k+2,2k-\frac{i_0}{2}+1},q{a}_{2,3,\dots ,\frac{i_0}{2}}$。

符号为负的 $2k-1$ 项为，$p{a}_{2k-\frac{i_0}{2}+2,\dots ,2k},q{a}_{\frac{i_0}{2}+1,\dots ,2k}$。

根据 $a$ 的排序方式和 $i_0$ 的意义可知，正项中的每一个大于负项中的每一个。又 $q{a}_1$ 与 $p{a}_{k+1}$ 无论哪项为正，肯定大于上述的负项。故序列 $b$ 符合期望。当 $n$ 为奇数时，同理可得；当 $p<q$ 时，同理可得，但对 $a$ 的排序方式改为“小大小大……”。

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

代码

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e6 + 100;
int n, p, q, a[maxn], cnt, b[maxn];
long long ans;
int read () 
{
	int s = 0, w = 1;
	char ch = getchar ();
	while (ch < '0' || ch > '9') 
	{
		if (ch == '-') w = -1;
		ch = getchar ();
	}
	while (ch >= '0' && ch <= '9') 
	{
		s = s * 10 + ch - '0';
		ch = getchar ();
	}
	return s * w;
}
int main ()
{
	n = read (), p = read (), q = read ();
	for (int i = 1; i <= n; i++)
		a[i] = read ();
	sort (a + 1, a + n + 1);
	for (int i = 1, j = n; i <= j; i++, j--) 
	{
		if (i == j) b[++ cnt] = a[i];
		if (p >= q) 
			b[++ cnt] = a[i], b[++ cnt] = a[j];
		else 
			b[++ cnt] = a[j], b[++ cnt] = a[i];
	}
	for (int i = 1; i <= n; i++) 
	{
		int to = (i == n ? 1 : i + 1);
		ans += abs (1ll * p * b[i] - 1ll * q * b[to]);
	}
	printf ("%lld\n", ans);
	for (int i = 1; i <= n; i++)
		printf ("%d ", b[i]);
	return 0;
}

收获

1. 敢于推测，敢于证明。