[Codeforces] number theory (R1600) Part.5

[Codeforces] number theory (R1600) Part.5

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

742B. Arpa’s obvious problem and Mehrdad’s terrible solution

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

题意

给定一个长度为$n(1\le n\le 1\mathrm{e}5)$的序列$a_1,\cdots ,a_n(1\le a_i\le 1\mathrm{e}5)$和一个整数$x(0\le x\le 1\mathrm{e}5)$.求使得$1\le i<j\le n$且$a_i\mathrm{xor}{a}_j=x$的数对$(i,j)$的个数.

思路

$cnt[x]$表示数$x$的出现个数,扫一遍即可.

注意$x\mathrm{xor}y$最大可取到 2 max ⁡ { x , y } 2\max\{x,y\} 2max{x,y}.注意答案可能爆int.

代码

const int MAXN = 2e5 + 5;  // 两倍空间
int cnt[MAXN];

void solve() {
  int n, x; cin >> n >> x;
  ll ans = 0;
  while (n--) {
    int a; cin >> a;
    ans += cnt[a ^ x];
    cnt[a]++;
  }
  cout << ans;
}

int main() {
  solve();
}

---

---

743C. Vladik and fractions

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

题意

给定一个整数$n(1\le n\le 1\mathrm{e}4)$,求三个相异的整数$x,y,z(1\le x,y,z\le 1\mathrm{e}9)s.t.\frac{2}{n}=\frac{1}{x}+\frac{1}{y}+\frac{1}{z}$,无解则输出$-1$.

思路

注意到$\frac{2}{n}=\frac{1}{n}+\frac{1}{n+1}+\frac{1}{n(n+1)}$,取$x=n,y=n+1,z=n(n+1)$,显然$x,y,z\le 1\mathrm{e}9$.

注意$n=1$时无解.

代码

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

  if (n == 1) cout << -1;
  else cout << n << ' ' << n + 1 << ' ' << n * (n + 1);
}

int main() {
  solve();
}

---

---

757B. Bash’s Big Day

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

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

给定$n(1\le n\le 1\mathrm{e}5)$个数$a_1,\cdots ,a_n(1\le a_i\le 1\mathrm{e}5)$,从中取最多的数$a_1,\cdots ,a_{k}s.t.\gcd (a_1,\cdots ,a_k)>1$,特别地,可单独取一个$a_i=1$的数.问最多能取多少个数.

思路

对每个$a_i$素因数分解,预处理出每个素因子的出现次数,取最多的即可.

代码

const int MAXN = 1e5 + 5;
umap<int, int> cnt;  // 素因子出现次数

void solve() {
  CaseT{
    int a; cin >> a;
    for (int i = 2; i <= sqrt(a); i++) {
      if (a % i == 0) {
        cnt[i]++;
        while (a % i == 0) a /= i;
      }
    }
    if (a > 1) cnt[a]++;
  }

  int ans = 1;
  for (auto [u, v] : cnt) ans = max(ans, v);
  cout << ans;
}

int main() {
  solve();
}

---

---

762A. k-th divisor

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

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

给定两整数$n,k(1\le n\le 1\mathrm{e}15,1\le k\le 1\mathrm{e}9)$,求$n$的第$k$小的约数,无解输出$-1$.

思路

$O\left(\sqrt{n}\right)$预处理出$n$的约数即可.

代码

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

  vl divisors;  // 注意ll
  for (int i = 1; i <= sqrt(n); i++) {
    if (n % i == 0) {
      divisors.push_back(i);
      if (i != n / i) divisors.push_back(n / i);
    }
  }

  if (divisors.size() < k) cout << -1;
  else {
    sort(all(divisors));
    cout << divisors[k - 1];
  }
}

int main() {
  solve();
}

---

---

840A. Leha and Function

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

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

考察从$n$个元素$[1,\cdots ,n]$中选$k$个元素构成一个子集,定义$f(n,k)$表示所有子集中最小元素的期望.

给定两长度为$n(1\le n\le 2\mathrm{e}5)$的序列$a_1,\cdots ,a_n$和$b_1,\cdots ,b_n(1\le a_i,b_i\le 1\mathrm{e}9)$.将序列$a[]$重排为序列$a^{\prime }[]s.t.\sum _{i=1}^{n}f(a_i,b_i)$最大,输出任一满足条件的序列$a^{\prime }[]$.

思路

从$n$个元素中选$k$个元素的方案数为$C_n^k$.枚举$i$,为使得其是子集中最小的元素,应在比它大的$(n-i)$个元素中选$(k-1)$个元素,故$f(n,k)=\frac{\sum _{i=1}^{n-k+1}i\cdot C_{n-i}^{k-1}}{C_n^k}$.

注意到$\frac{f(n,k)}{f(n,k-1)}=\frac{k}{k+1}$,则$f(n,k)=\frac{n+1}{k+1}$.为使得$f(n,k)$尽量大,应让较大的$n$与较小的$k$配对.

将序列$a[]$降序排列、$b[]$升序排列后配对即可.

代码

void solve() {
  int n; cin >> n;
  vi a(n);
  vii b(n);
  for (auto& ai : a) cin >> ai;
  for (int i = 0; i < n; i++) {
    cin >> b[i].first;
    b[i].second = i;
  }

  sort(all(a), greater<int>()), sort(all(b));

  vi ans(n);
  for (int i = 0; i < n; i++) ans[b[i].second] = a[i];
  for (int i = 0; i < n; i++) cout << ans[i] << " \n"[i == n - 1];
}

int main() {
  solve();
}

---

---

870C. Maximum splitting

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

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

有$t(1\le t\le 1\mathrm{e}5)$组测试数据.每组测试数据输入一个整数$n(1\le n\le 1\mathrm{e}9)$,将$n$分解为尽量多个合数之和,输出合数的个数;无解输出$-1$.

思路

先求出$[1,16]$内的数的方案数,对$>16$的数,显然最优策略是拆出若干个$4$,将其减小为$\le 16$的数.

代码

const int MAXN = 20;
int ans[MAXN] = { 0,-1,-1,-1,1,-1,1,-1,2,1,2,-1,3,2,3,2,4 };

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

    if (n <= 16) cout << ans[n] << endl;
    else {
        int tmp = (n - 16) / 4 + 1;
        cout << tmp + ans[n - 4 * tmp] << endl;
    }
}

int main() {
	CaseT  // 单测时注释掉该行
	solve();
}

---

---

876B. Divisiblity of Differences

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

题意

给定$n$个数$a_1,\cdots ,a_n(0\le a_i\le 1\mathrm{e}9)$,从中选出$m(1\le m\le 1\mathrm{e}5)$个元素,使得其中任意两元素之差是$k(2\le k\le n\le 1\mathrm{e}5)$的倍数,若有解,第一行输出"Yes",第二行输出任意方案;否则输出"No".

思路

将所有数按模$k$的余数分类,输出任一元素个数$\ge m$的类中的$m$个元素即可.

代码

const int MAXN = 1e5 + 5;
vi nums[MAXN];

void solve() {
  int n, m, k; cin >> n >> m >> k;
  while (n--) {
    int a; cin >> a;
    nums[a % k].push_back(a);
  }

  for (int i = 0; i < MAXN; i++) {
    if (nums[i].size() >= m) {
      cout << "Yes" << endl;
      for (int j = 0; j < m; j++) cout << nums[i][j] << ' ';
      return;
    }
  }
  cout << "No" << endl;
}

int main() {
  solve();
}

---

---

891A. Pride

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

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

给定一个长度为$n(1\le n\le 2000)$的序列$a_1,\cdots ,a_n(1\le a_i\le 1\mathrm{e}9)$.现有操作:选定$a[]$中两相邻的元素$x,y$,将它们之一替换为$\gcd (x,y)$.若经若干次操作可将$a[]$的元素都变为$1$,输出最小步数;否则输出$-1$.

思路

①若初始时序列中有$1$,设$1$的个数为$cnt$,显然$ans=n-cnt$.

②若初始时序列中无$1$,$O(n^2)$暴力找是否存在一段区间$[l,r]\ s.t.\ $其中元素的$\gcd=1$,若存在,则$ans=(r-l)+(n-1)$.枚举起点$l$即可.

代码

void solve() {
  int n; cin >> n;
  vi a(n + 1);
  int cnt = 0;  // 1的个数
  for (int i = 1; i <= n; i++) {
    cin >> a[i];
    cnt += a[i] == 1;
  }

  if (cnt) {
    cout << n - cnt;
    return;
  }

  int ans = INF;
  for (int i = 1; i <= n; i++) {
    int g = a[i];  // gcd
    for (int j = i + 1; j <= n; j++) {
      g = gcd(g, a[j]);

      if (g == 1) {
        ans = min(ans, j - i + n - 1);
        break;
      }
    }
  }
  cout << (ans == INF ? -1 : ans);
}

int main() {
  solve();
}

---

---

900B. Position in Fraction

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

题意

给定三个整数$a,b,c(1\le a<b\le 1\mathrm{e}5,0\le c\le 9)$,求分数$\frac{a}{b}$的十进制表示中数码$c$在小数点后第一次出现的位数,下标从$1$开始,无解输出$-1$.

思路

分数是有限或无限循环小数,显然循环节不超过$b$,模拟前$b$次除法即可.

代码

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

  int cur = a;
  for (int i = 1; i <= b; i++) {
    cur *= 10;
    int tmp = cur / b;
    if (tmp == c) {
      cout << i;
      return;
    }

    cur %= b;
  }
  cout << -1;
}

int main() {
  solve();
}

---

---

912B. New Year’s Eve

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

题意

给定两整数$n,k(1\le k\le n\le 1\mathrm{e}18)$,从$[1,n]$中选至多$k$个数使得它们异或之和最大,输出最大异或和.

思路

①$k=1$时,显然$ans=n$.

②$k\ge 2$时,设$n$的MSB为第$p$位(下标从$0$开始),则$ans\le 2^p-1$.

​ 注意到$2^p\mathrm{xor}{2}^p-1=2^{p+1}-1$,故取不大于$n$的$2^p$和$2^p-1$时异或和最大.

​ 可遍历$n$的二进制数位求$p$再求答案,也可用若干次$\times 2+1$实现.

代码

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

  if (k == 1) cout << n;
  else {
    ll ans = 1;
    while (ans < n) ans = ans * 2 + 1;
    cout << ans;
  }
}

int main() {
  solve();
}

---

---