[Codeforces] number theory (R1600) Part.3

最新推荐文章于 2024-06-12 10:00:23 发布

Hytidel

最新推荐文章于 2024-06-12 10:00:23 发布

阅读量184

点赞数

分类专栏： Codeforces数论专精文章标签：算法

本文链接：https://blog.csdn.net/Hytidel/article/details/127008875

版权

Codeforces数论专精专栏收录该内容

24 篇文章 1 订阅

订阅专栏

[Codeforces] number theory (R1600) Part.3

题单:https://codeforces.com/problemset/page/1?tags=number+theory%2C1201-1600

346A. Alice and Bob

原题指路:https://codeforces.com/problemset/problem/346/A

题意 $(2\ \mathrm{s})$

Alice和Bob两人玩游戏,Alice先手.初始时给定一个整数集合 $S$ .每轮玩家选两个相异的整数 $x,y\in S$ ,使得 $|x-y|\not\in S$ ,并将 $∣ x - y ∣$ 加入 $S$ 中,不能操作者负.问最后谁赢.

第一行输入一个整数 $n\ \ (2\leq n\leq 100)$ .第二行输入 $n$ 个相异的整数 $a_1,\cdots,a_n\ \ (1\leq a_i\leq 1\mathrm{e}9)$ .

思路

类似于闭包,无论中间过程,最后得到的集合都相同.设 $\displaystyle g=\gcd_{1\leq i\leq n}a_i$ ,则最终集合为 $\{d,2d,3d,\cdots,\max\{x_i\}\}$ ,故轮数为 $\dfrac{\max\{x_i\}}{d}$ ,判断其奇偶性即可.

代码

void solve() {
	int n; cin >> n;
	valarray<int> a(n);
	int g = 0;  // gcd
	for (auto& ai : a) {
		cin >> ai;
		g = gcd(g, ai);
	}

	int ans = a.max() / g - n;
	cout << (ans & 1 ? "Alice" : "Bob");
}

int main() {
	solve();
}

353C. Find Maximum

原题指路:https://codeforces.com/problemset/problem/353/C

题意

给定一个长度为 $n$ 的序列 $a_1,\cdots,a_n$ .定义函数 $\displaystyle f(x)=\sum_{i=0}^{n-1} a_i\cdot bit(i)$ ,它以 $0,2^n-1]$ 中的一个整数为参数 $x$ ,其中 $bi t (i) = 1$ 当且仅当参数 $x$ 的二进制表示第 $i$ 位为 $1$ ;否则 $bi t (i) = 0$ .如 $n=4,x=11=2^0+2^1+2^3$ 时, $f(x)=a_0+a_1+a_3$ .给定一个整数 $m$ ,求 $f (x)$ 在 $x\in[0,m]$ 上的最大值.

第一行输入一个整数 $n\ \ (1\leq n\leq 1\mathrm{e}5)$ .第二行输入 $n$ 个整数 $a_1,\cdots,a_n\ \ (0\leq a_i\leq 1\mathrm{e}4)$ .第三行输入一个长度为 $n$ 的 $0 - 1$ 串,表示整数 $m$ 的二进制表示.

思路I

①若 $m$ 的MSB为 $0$ ,则 $\forall x\in[0,m]$ 的二进制表示的第 $(n - 1)$ 位为 $0$ ,则 $a_{n-1}$ 不会出现在 $f (x)$ 中.

②若 $m$ 的MSB为 $1$ ,则对某些 $x$ , $a_{n-1}$ 会出现在 $f (x)$ 中.考察不包含 $a_{n-1}$ 的 $f (x)$ ,显然此时 $x\in[0,2^{n-1}-1]$ .

故 $x=2^{n-1}-1$ 时 $f (x)$ 取得最大值.

从高位到地位枚举 $m$ 的二进制表示的数位,数位为 $1$ 时更新答案即可,过程用前缀和、后缀和加速.

代码I

const int MAXN = 1e5 + 5;
int n;
int a[MAXN];
int pre[MAXN];  // a[]的前缀和
string s;

void solve() {
	cin >> n;
	for (int i = 0; i < n; i++) {
		cin >> a[i];
		pre[i + 1] = pre[i] + a[i];
	}
	cin >> s;

	int ans = 0;
	for (int i = n - 1, suf = 0; i >= 0; i--)  // suf为a[]的后缀和
		if (s[i] == '1') ans = max({ ans,pre[i] + suf,suf += a[i] });
	cout << ans;
}

int main() {
	solve();
}

371B. Fox Dividing Cheese

原题指路:https://codeforces.com/problemset/problem/371/B

题意

给定两整数 $a, b$ .现有如下三种操作,选定一个数 $x$ :①若 $2\mid x$ ,令 $x-=\dfrac{x}{2}$ ;②若 $3\mid x$ ,令 $x-=\dfrac{2}{3}x$ ;③若 $5\mid x$ ,令 $x-=\dfrac{4}{5}x$ .若经若干次操作后能使得 $a = b$ ,输出最小操作次数;否则输出 $- 1$ .

思路

显然操作①等价于 $x / = 2$ ,操作②等价于 $x / = 3$ ,操作③等价于 $x / = 5$ .

设 $a=x\cdot 2^{\alpha_1}\cdot 3^{\alpha_2}\cdot 5^{\alpha_3},b=y\cdot 2^{\beta_1}\cdot 3^{\beta_2}\cdot 5^{\beta_3}$ .

①若 $x\neq y$ ,显然无解.

②若 $x = y$ ,则 $ans=|\alpha_1-\beta_1|+|\alpha_2-\beta_2|+|\alpha_3-\beta_3|$ .

代码

vi get(int n) {
	vi res(4, 0);  // 素因子2的次数、3的次数、5的次数、余下的因子之积
	while (n % 2 == 0) {
		n /= 2;
		res[0]++;
	}
	while (n % 3 == 0) {
		n /= 3;
		res[1]++;
	}
	while (n % 5 == 0) {
		n /= 5;
		res[2]++;
	}
	res[3] = n;
	return res;
}

void solve() {
	int a, b; cin >> a >> b;

	auto res1 = get(a), res2 = get(b);
	if (res1[3] != res2[3]) cout << -1;
	else cout << abs(res1[0] - res2[0]) + abs(res1[1] - res2[1]) + abs(res1[2] - res2[2]);
}

int main() {
	solve();
}

375A. Divisible by Seven

原题指路:https://codeforces.com/problemset/problem/375/A

题意

给定一个至少包含包含数码 $1, 6, 8, 9$ 的整数 $a$ ,重排其数码使得其能被 $7$ 整除,输出重排后的数字(无前导零),无解输出 $0$ .

思路

考察 $1, 6, 8, 9$ 是否能拼出模 $7$ 后余数为 $0\sim 6$ 的所有数.枚举知: $1869\% 7=0,8961\% 7=1,1689\% 7=2,6198\% 7=3,1698\% 7=4,9861\% 7=5,1896\% 7=6$ .故前面的数码任意排,后面补上对应余数的 $1, 6, 8, 9$ 的排列即可.

因不输出前导零,将所有 $0$ 放在最后输出即可,显然这不影响余数.

代码

const int MAXN = 10;
const int last[] = { 1869,8961,1689,6198,1698,9861,1896,1869 };
int cnt[MAXN];  // 每个数码出现的次数

void solve() {
	string s; cin >> s;
	for (auto ch : s) cnt[ch & 15]++;

	cnt[1]--, cnt[6]--, cnt[8]--, cnt[9]--;

	int r = 0;  // 余数
	for (int i = 1; i < 10; i++) {
		while (cnt[i]--) {
			r = (r * 10 + i) % 7;
			cout << i;
		}
	}
	cout << last[7 - r * 10000 % 7];
	while (cnt[0]--) cout << 0;
}

int main() {
	solve();
}

414A. Mashmokh and Numbers

原题指路:https://codeforces.com/problemset/problem/414/A

题意

给定一个长度为 $n$ 的序列 $a_1,\cdots,a_n$ .现有操作:若当前序列元素个数 $\geq 2$ ,取序列的前两个元素 $x, y$ ,将它们删去,同时获得 $\gcd(x,y)$ 的分数.构造一个元素相异的序列$a_1,\cdots,a_n\ \ (1\leq a_i\leq 1\mathrm{e}9)\ s.t.\ $操作后总得分为$ k $. 无解输出$ -1$.

第一行输入两个整数 $n,k\ \ (1\leq n\leq 1\mathrm{e}5,0\leq k\leq 1\mathrm{e}8)$ .

思路

显然每次操作至少获得 $1$ 分,则 $k<\left\lfloor\dfrac{n}{2}\right\rfloor$ 时无解.

$k\geq \dfrac{n}{2}$ 时,考察 $x=k-\left\lfloor\dfrac{n-2}{2}\right\rfloor$ ,取 $a_1=x,a_2=2x$ ,其余元素取一个连续的升序序列即可.

注意特判 $n = 1$ 的情况,此时有解当且仅当 $k = 0$ .

代码

void solve() {
	int n, k; cin >> n >> k;
	
	if (n == 1) cout << (k ? -1 : 1);
	else if (k < n / 2) cout << -1;
	else {
		int x = k - (n - 2) / 2;
		cout << x << ' ' << 2 * x << ' ';
		for (int i = 2 * x + 1; i <= 2 * x  + n - 2; i++) cout << i << ' ';
	}
}

int main() {
	solve();
}

414B. Mashmokh and ACM

原题指路:https://codeforces.com/problemset/problem/414/B

题意

称一个整数序列 $a_1,\cdots,a_n$ 是好的,如果 $a_i\mid a_{i+1}\ \ (1\leq i\leq n-1)$ .给定整数 $n,k\ \ (1\leq n,k\leq 2000)$ ,求长度为 $k$ 的、元素范围为 $[1, n]$ 的好的整数序列的个数,答案对 $1\mathrm{e}9+7$ 取模.

思路

$d p [i] [j]$ 表示长度为 $i$ 的、以 $j$ 结尾的好的序列的个数.枚举 $j$ 的约数 $d$ ,状态转移方程 $\displaystyle dp[i][j]=\sum_{d\mid x}dp[i-1][d]$ .初始条件可取 $d p [0] [1] = 1$ .

暴力转移时间复杂度 $O(kn\sqrt{n})$ ,会TLE.考虑优化,先预处理出 $2000$ 内所有数的约数,总时间复杂度 $O(kn\log n)$ .

代码

const int MAXN = 2005;
const int MOD = 1e9 + 7;
vi divisors[MAXN];  // 存每个数的约数
int dp[MAXN][MAXN];  // dp[i][j]表示长度为i的、以j结尾的好的序列的个数

void init() {
	for (int i = 1; i < MAXN; i++)
		for (int j = i; j < MAXN; j += i) divisors[j].push_back(i);
}

void solve() {
	init();

	int n, k; cin >> n >> k;
	
	dp[0][1] = 1;
	for (int i = 1; i <= k; i++) {
		for (int j = 1; j <= n; j++) {
			for (auto d : divisors[j])
				dp[i][j] = ((ll)dp[i][j] + dp[i - 1][d]) % MOD;
		}
	}

	int ans = 0;
	for (int i = 1; i <= n; i++) ans = ((ll)ans + dp[k][i]) % MOD;
	cout << ans;
}

int main() {
	solve();
}

460B. Little Dima and Equation

原题指路:https://codeforces.com/problemset/problem/460/B

题意

给定系数 $a,b,c\ \ (1\leq a\leq 5,1\leq b\leq 1\mathrm{e}4,-1\mathrm{e}4\leq c\leq 1\mathrm{e}4)$ ,求方程 $x=b\cdot s(x)^a+c$ 在 $x\in(0,1\mathrm{e}9)$ 中的所有整数解,其中 $s (x)$ 表示 $x$ 的十进制表示的数位之和.第一行输出解的个数,第二行升序输出所有解.

思路

注意到对 $\forall x\in(0,1\mathrm{e}9),s(x)\in[1,81]$ ,故可枚举 $s$ ,求出对应的 $x$ 后,判断是否 $s (x) = s$ 且 $x<1\mathrm{e}9$ 即可.

代码

int get_sum(int x) {
  int res = 0;
  while (x) {
    res += x % 10;
    x /= 10;
  }
  return res;
}

ll qpow(int a, int k) {
  ll res = 1;
  while (k) {
    if (k & 1) res = res * a;
    k >>= 1;
    a = a * a;
  }
  return res;
}

void solve() {
  int a, b, c; cin >> a >> b >> c;

  vi ans;
  for (int s = 1; s <= 81; s++) {
    ll x = qpow(s, a) * b + c;
    if (get_sum(x) == s && x < 1e9) ans.push_back(x);
  }

  cout << ans.size() << endl;
  for (auto i : ans) cout << i << ' ';
}

int main() {
  solve();
}

495B. Modular Equations

原题指路:https://codeforces.com/problemset/problem/495/B

题意

给定整数 $a,b\ \ (0\leq a,b\leq 1\mathrm{e}9)$ ,求关于 $x$ 的方程 $a\ \mathrm{mod}\ x=b$ 的解数,无解输出 $0$ ,有无穷多组解输出"infinity".

思路

① $a < b$ 时,显然无解.

② $a = b$ 时, $\forall x>a$ 都是方程的解.

③ $a > b$ 时,显然 $x\mid (a-b)$ ,且 $x > b$ .故方程的解数即 $(a - b)$ 的 $> b$ 的约数个数.

代码

void solve() {
  int a, b; cin >> a >> b;

  if (a < b) cout << 0;
  else if (a == b) cout << "infinity";
  else {
    int c = a - b;
    int ans = 0;
    for (int i = 1; i <= c / i; i++) 
      if (c % i == 0) ans += (i > b) + (i != c / i && c / i > b);
    cout << ans;
  }
}

int main() {
  solve();
}

515B. Drazil and His Happy Friends

原题指路:https://codeforces.com/problemset/problem/515/B

题意 ( $2\ \mathrm{s}$ )

有 $n$ 个男生和 $m$ 个女生,分别编号 $0\sim (n-1)$ 和 $0\sim (m-1)$ .初始时他们之中有些人开心,有些人不开心.现第 $i$ 天会邀请编号 $(i\ \mathrm{mod}\ n)$ 的男生和编号 $(i\ \mathrm{mod}\ m)$ 的女生吃饭,若两人中有一人开心,则两人都会永久地开心.问最后能否使得所有人开心,若能则输出"Yes";否则输出"No".

第一行输入两个整数 $n,m\ \ (1\leq n,m\leq 100)$ .第二行先输入一个整数 $b\ \ (0\leq b\leq n)$ .接下来输入 $b$ 个整数 $x_1,\cdots,x_b\ \ (0\leq x_i<n)$ ,表示初始时开心的男生.第三行先输入一个整数 $g\ \ (0\leq g\leq m)$ .接下来输入 $g$ 个整数 $y_1,\cdots,y_b\ \ (0\leq y_i<m)$ ,表示初始时开心的女生.

思路

显然至多 $nm$ 天会出现循环,若 $nm$ 天后还有人不开心,则答案为No.可模拟前 $nm$ 天的情况,时间复杂度 $O (nm)$ .

事实上本题有 $O (n + m)$ 的解法.设 $g=\gcd(n,m)$ .若 $i$ 号人开心,则编号为 $x\equiv i\ \ (\mathrm{mod}\ g)$ 的人都会变开心,则只需考察编号 $\mathrm{mod}\ g$ 的每个余数中是否至少有一个人开心即可.

代码

const int MAXN = 105;
bool happy[MAXN];

void solve() {
	int n, m; cin >> n >> m;
	int g = gcd(n, m);

	int xn; cin >> xn;
	while (xn--) {
		int x; cin >> x;
		happy[x % g] = true;
	}
	int yn; cin >> yn;
	while (yn--) {
		int y; cin >> y;
		happy[y % g] = true;
	}

	for (int i = 0; i < g; i++) {
		if (!happy[i]) {
			cout << "No";
			return;
		}
	}
	cout << "Yes";
}

int main() {
	solve();
}

519C. A and B and Team Training

原题指路:https://codeforces.com/problemset/problem/519/C

题意

一支队伍有如下两种组成方式:① $1$ 个A类人和 $2$ 个B类人;② $2$ 个A类人和 $1$ 个B类人.现有 $n$ 个A类人和 $m\ \ (0\leq n,m\leq 5\mathrm{e}5)$ 个B类人,求最多能组成多少支队伍.

思路

设组成①类队 $i\ \ (0\leq i\leq n)$ 支,则剩余 $(n - i)$ 个A类人和 $(m - 2 i)$ 个B类人,他们能组成②类队 $\min\left\{\dfrac{n-i}{2},m-2i\right\}$ .枚举每个 $i$ ,更新答案即可.

代码

void solve() {
  int n, m; cin >> n >> m;

  int ans = 0;
  for (int i = 0; i <= n; i++) {  // 枚举①类队队伍数
    int res = i;
    if (m - 2 * i >= 0) {
      res += min((n - i) / 2, m - 2 * i);
      ans = max(ans, res);
    }
  }
  cout << ans;
}

int main() {
  solve();
}