[Codeforces] number theory (R1600) Part.4

[Codeforces] number theory (R1600) Part.4

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

566F. Clique in the Divisibility Graph

原题指路:https://codeforces.com/problemset/problem/566/F

题意

无向图中的一个团是一些节点构成的集合,其中任意两点都有边相连.特别的,空集和只包含一个节点的集合也视为团.

对一个元素为正整数的集合 A = { a 1 , ⋯   , a n } A=\{a_1,\cdots,a_n\} A={a1​,⋯,an​},定义可分性图为:对图中以$a_i,a_j\in A(i\ne j)$为编号的两节点,它们间有边相连当且仅当$a_i\mid a_j$或$a_j\mid a_i$.

给定一个包含$n(1\le n\le 1\mathrm{e}6)$个正整数的集合 A = { a 1 , ⋯   , a n }    ( 1 ≤ a i ≤ 1 e 6 ) A=\{a_1,\cdots,a_n\}\ \ (1\leq a_i\leq 1\mathrm{e}6) A={a1​,⋯,an​}  (1≤ai​≤1e6),元素升序给出.对$A$,求其可分性图中的最大团中的节点数.

思路

将团中的节点升序排列为$x_1,\cdots ,x_k$,则$x_i\mid x_{i+1}(1\le i\le k-1)$.

$dp[x]$表示以$x$结尾的能构成团的序列长度的最大值.状态转移方程$dp[x]=\max \left\{1,\underset{d\mid x}{\max }dp[d]+1\right\}$.

总时间复杂度$\frac{n}{1}+\frac{n}{2}+\cdots +\frac{n}{n}=O(n\log n)$.

代码

const int MAXN = 1e6 + 5;
int n;
bool cnt[MAXN];
int dp[MAXN];  // dp[x]表示以x结尾的能构成团的序列长度的最大值

void solve() {
  CaseT{
    int x; cin >> x;
    cnt[x] = true;
  }

  int ans = 0;
  for (int i = 1; i <= 1e6; i++) {
    if (cnt[i]) {
      ans = max(ans, ++dp[i]);

      for (int j = i; j <= 1e6; j += i) dp[j] = max(dp[j], dp[i]);
    }
  }
  cout << ans;
}

int main() {
  solve();
}

---

---

568A. Primes or Palindromes?

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

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

对正整数$n$,定义$\pi (n)$为不大于$n$的素数的个数,$rub(n)$为不大于$n$的回文数的个数.给定$A=\frac{p}{q}\left(1\le p,q\le 1\mathrm{e}4,\frac{1}{42}\le \frac{p}{q}\le 42\right)$,求最大的正整数$ns.t.\pi (n)\le A\cdot rub(n)$,无解输出"Palindromic tree is better than splay tree".

思路

$\pi (n)\sim \frac{n}{\ln n}$.长度为$k$的回文数的个数约$10^{\lfloor \frac{k+1}{2}\rfloor }$,是$O(\sqrt{n})$的级别.故$A\sim \frac{\sqrt{n}}{\ln n}$,则$n<1\mathrm{e}7$.

筛出$1\mathrm{e}7$内的素数后求$\pi (n)$的前缀和,再求出$1\mathrm{e}7$内的数的$rub(n)$的前缀和即可,总时间复杂度$O(ans\cdot \log \log ans)$.

事实上,对$A\le 42,ans\le 2\mathrm{e}6$.

代码

const int MAXN = 2e6 + 5;
bool prime[MAXN] = { true };
int prepi[MAXN], prerub[MAXN];  // pi(n)和rub(n)的前缀和

bool check(int n) {  // 判断n是否为回文数
  string s = to_string(n), tmps = s;
  reverse(all(s));
  return s == tmps;
}

void init() {
  for (int i = 2; i < MAXN; i++) prime[i] = true;
  for (int i = 2; i < MAXN; i++) {
    if (!prime[i]) continue;

    for (int j = 2 * i; j < MAXN; j += i) prime[j] = false;
  }

  for (int i = 2; i < MAXN; i++) prepi[i] = prepi[i - 1] + prime[i];
  for (int i = 1; i < MAXN; i++) prerub[i] = prerub[i - 1] + check(i);
}

void solve() {
  init();

  int p, q; cin >> p >> q;
  for (int i = MAXN - 1;; i--) {
    if (q * prepi[i] <= p * prerub[i]) {
      cout << i;
      return;
    }
  }
}

int main() {
  solve();
}

---

---

573A. Bear and Poker

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

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

给定一个长度为$n(2\le n\le 1\mathrm{e}5)$的序列$a_1,\cdots ,a_n(1\le a_i\le 1\mathrm{e}9)$.现有操作:选择一个下标$i(1\le i\le n)$,令$a_i\ast =2$或$a_i\ast =3$.问经过若干次操作后能否使得序列元素相等,若能则输出"Yes";否则输出"No".

思路

判断所有$a_i$除$2,3$外的素因子之积是否相等即可.

代码

void solve() {
  int n; cin >> n;
  vi a(n);
  for (auto& ai : a) {
    cin >> ai;
    while (ai % 2 == 0) ai /= 2;
    while (ai % 3 == 0) ai /= 3;
  }

  cout << (a == vi(n, a[0]) ? "Yes" : "No");
}

int main() {
  solve();
}

---

---

576A. Vasya and Petya’s Game

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

题意

A和B玩游戏.A想一个整数$x\in [1,n](1\le n\le 1000)$,B猜该数.B每次可询问$x$是否整除某整数$y(1\le y\le n)$,A会回答其问题.问B至少问多少个问题能确定$x$,先输出问题数,再输出每个问题的$y$.若有多组问题数最少的解,输出任一组.

思路

对$n$的某个素因子$p$,若B不问$p^k(k\ge 1)$,则他无法分辨$p^{k-1}$和$p^k$,故B需问$\le n$的所有素数及其不超过$n$的幂次.

代码

const int MAXN = 1005;
int primes[MAXN], cnt;
bool state[MAXN];

void init() {
  state[0] = state[1] = true;
  for (int i = 2; i * i < MAXN; i++) {
    if (!state[i]) 
      for (int j = i * i; j < MAXN; j += i) state[j] = true;
  }
}

void solve() {
  init();

  int n; cin >> n;

  vi ans;
  for (int i = 1; i <= n; i++) {
    if (!state[i]) {
      int cur = 1;
      while (cur <= n / i) {
        cur *= i;
        ans.push_back(cur);
      }
    }
  }

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

int main() {
  solve();
}

---

---

626C. Block Towers

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

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

有A和B两种砖,分别由$2$、$3$个方块组成.现分别用两种砖横放堆叠造塔,使得任意两塔包含的方块数不同,共造出了$n$座用A型砖构成的塔和$m$座用B型砖构成的塔.给定两整数$n,m(0\le n,m\le 1\mathrm{e}6,n+m>0)$问最高的塔的高度的最小值.

思路

设答案为$ans$,则至少有$n$个不超过$ans$的$2$的倍数、$m$个不超过$ans$的$3$的倍数、$(n+m)$个$2$或$3$的倍数,即$n\le \lfloor \frac{ans}{2}\rfloor ,m\le \lfloor \frac{ans}{3}\rfloor ,n+m\le \lfloor \frac{ans}{2}\rfloor +\lfloor \frac{ans}{3}\rfloor -\lfloor \frac{ans}{6}\rfloor$,二分即可.

移项知$ans\ge \frac{3}{2}(n+m)$,故二分上界可取$3\mathrm{e}6$.

代码

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

  int l = 1, r = 3e6;
  while (l < r) {
    int mid = l + r >> 1;
    if (n <= mid / 2 && m <= mid / 3 && n + m <= mid / 2 + mid / 3 - mid / 6) r = mid;
    else l = mid + 1;
  }
  cout << l;
}

int main() {
  solve();
}

---

---

630K. Indivisibility

原题指路:https://codeforces.com/problemset/problem/630/K

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

给定整数$n(1\le n\le 1\mathrm{e}18)$,求$[1,n]$中有多少个数不能被$2\sim 10$整除.

思路

容斥,$ans=n-\lfloor \frac{n}{2}\rfloor -\lfloor \frac{n}{3}\rfloor -\lfloor \frac{n}{5}\rfloor -\lfloor \frac{n}{7}\rfloor$

​ $+\lfloor \frac{n}{2\cdot 3}\rfloor +\lfloor \frac{n}{2\cdot 5}\rfloor +\lfloor \frac{n}{2\cdot 7}\rfloor +\lfloor \frac{n}{3\cdot 5}\rfloor +\lfloor \frac{n}{3\cdot 7}\rfloor +\lfloor \frac{n}{5\cdot 7}\rfloor$

​ $-\lfloor \frac{n}{2\cdot 3\cdot 5}\rfloor -\lfloor \frac{n}{2\cdot 3\cdot 7}\rfloor -\lfloor \frac{n}{2\cdot 5\cdot 7}\rfloor -\lfloor \frac{n}{3\cdot 5\cdot 7}\rfloor +\lfloor \frac{n}{2\cdot 3\cdot 5\cdot 7}\rfloor$.

代码

void solve() {
  ll n; cin >> n;
  cout << n - n / 2 - n / 3 - n / 5 - n / 7
    + n / (2 * 3) + n / (2 * 5) + n / (2 * 7) + n / (3 * 5) + n / (3 * 7) + n / (5 * 7)
    - n / (2 * 3 * 5) - n / (2 * 3 * 7) - n / (2 * 5 * 7) - n / (3 * 5 * 7)
    + n / (2 * 3 * 5 * 7);
}

int main() {
  solve();
}

---

---

633B. A Trivial Problem

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

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

给定一个整数$m(1\le m\le 1\mathrm{e}5)$,求所有使得$n!$以$m$个零结尾的正整数$n$,第一行输出个数,第二行输出所有$n$.

思路

$n!$后缀零的个数等于$n$的素因子$5$的个数.

注意到每相邻十个数至多会为$n!$贡献一个后缀零,故只需枚举$[1,10m]$即可.

代码

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

  vi ans;
  int cnt = 0;  // 素因子5的个数
  for (int i = 1; i < 10 * n; i++) {
    int cur = i;
    while (cur % 5 == 0) {
      cur /= 5;
      cnt++;
    }
    if (cnt == n) ans.push_back(i);
    else if (cnt > n) break;
  }

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

int main() {
	solve();
}

---

---

678C. Joty and Chocolate

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

题意

有编号$1\sim n$的$n$个物品.现以某种顺序对它们染色,若该物品的编号能被$a$整除,则染红色;若该物品的编号能被$b$整除,则染蓝色;若该物品的编号同时被$a$和$b$整除,则只能染一种颜色.染色结束后,每个染红色的物品可获得$p$的收益,每个染蓝色的物品可获得$q$的收益.求最大收益.

第一行输入五个整数$n,a,b,p,q(1\le n,a,b,p,q\le 1\mathrm{e}9)$.

思路

一个物品可选择染红色或蓝色当且仅当其编号是$\mathrm{lcm}(a,b)$的倍数,这样的物品应染收益更大的颜色.

显然$ans=p\cdot \lfloor \frac{n}{a}\rfloor +q\cdot \lfloor \frac{n}{b}\rfloor -\min (p,q)\cdot \lfloor \frac{n}{\mathrm{lcm}(a,b)}\rfloor$.

代码

void solve() {
  ll n, a, b, p, q; cin >> n >> a >> b >> p >> q;
  cout << n / a * p + n / b * q - n / lcm(a, b) * min(p, q);
}

int main() {
  solve();
}

---

---

707C. Pythagorean Triples

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

题意

给定一个整数$n(1\le n\le 1\mathrm{e}9)$,问是否存在一个边长都为整数的直角三角形,使得$n$为其中一条边的边长.若有解则输出任一组解;否则输出$-1$.

思路

设$n$是直角三角形的短直角边的边长,$n^2+m^2=k^2$,则$n^2=(k+m)(k-m)$.

①$n$为奇数时,取$k+m=n^2,k-m=1$,则$m=\frac{n^2-1}{2},k=\frac{n^2+1}{2}$.

②$n$为偶数时,取$k+m=\frac{n^2}{2},k-m=2$,则$m=\frac{n^2}{4}-1,k=\frac{n^2}{4}+1$.

显然$n=1,2$时无解.

代码

void solve() {
  ll n; cin >> n;

  if (n <= 2) cout << -1;
  else if (n & 1) cout << (n * n - 1) / 2 << ' ' << (n * n + 1) / 2;
  else cout << n * n / 4 - 1 << ' ' << n * n / 4 + 1;
}

int main() {
  solve();
}

---

---

735D. Taxes

原题指路:https://codeforces.com/problemset/problem/735/D

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

若收益为$n(2\le n\le 2\mathrm{e}9)$元,则需交税最大的非$n$的约数元.现将$n$拆成$k$份,即$n_1+\cdots +n_k=n(k\ge 1)$,每份分别交税,问最少要交多少税.

思路

(1)$n$为素数时,$ans=1$.

(2)$n$非素数时,$ans\ge 2$.

​ ①若$n$是$\ge 4$的偶数,由Goldbach猜想:$n$可表示为两素数之和,此时$ans=2$.

​ ②若$n$是奇数,因$n=2+(n-2)$,

​ i)若$(n-2)$是素数,则$ans=2$.

​ ii)若$(n-2)$是合数,由弱Goldbach猜想:$n$可表示为三个奇素数之和,此时$ans=3$.

代码

bool check(int n) {  // 判断是否是素数
  if (n <= 2) return n == 2;
  for (int i = 2; i <= sqrt(n); i++)
    if (n % i == 0) return false;
  return true;
}

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

  if (check(n)) cout << 1;
  else if (n % 2 == 0) cout << 2;
  else cout << (check(n - 2) ? 2 : 3);
}

int main() {
  solve();
}

---

---