[Codeforces] number theory (R1600) Part.7

[Codeforces] number theory (R1600) Part.7

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

1113B. Sasha and Magnetic Machines

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

题意

给定一个长度为$n(2\le n\le 5\mathrm{e}4)$的序列$a_1,\cdots ,a_n(1\le a_i\le 100)$.现有操作:选择一对$a_i$和$a_j$,选择$a_i$的一个约数$x$,令$a_i/=x,a_j\ast =x$.问经过至多一次操作后得到的序列的元素之和的最小值.

思路

为使得序列元素之和最小,约数$x$应乘在序列的最小元素上.因至多进行一次操作,枚举每个$a_i$进行操作即可.

代码

void solve() {
  int n; cin >> n;
  vi a(n);
  for (auto& ai : a) cin >> ai;

  int minnum = *min_element(all(a));
  int sum = accumulate(all(a), 0);
  int ans = sum;

  for (auto ai : a) {
    for (int d = 1; d <= ai; d++) {
      if (ai % d == 0) {
        int res = sum - minnum - ai;
        res += minnum * d + ai / d;
        ans = min(ans, res);
      }
    }
  }
  cout << ans;
}

int main() {
  solve();
}

---

---

1133D. Zero Quantity Maximization

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

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

给定两长度为$n(1\le n\le 2\mathrm{e}5)$的序列$a_1,\cdots ,a_n$和$b_1,\cdots ,b_n(-1\mathrm{e}9\le a_i,b_i\le 1\mathrm{e}9)$.现按如下方式构造一个长度为$n$的序列$c[]$:选定一个实数$d$,令$c_i=d\cdot a_i+b_i(1\le i\le n)$.求$c[]$中$0$的个数的最大值.

思路

①若$a_i=0$,则$d$的选择对答案无影响,只需检查$b_i$是否为$0$即可.

②若$a_i\ne 0$,应取$d=-\frac{b_i}{a_i}$.将分数化为最简分数后,记录其出现的次数,最后$d$取出现次数最多的分数即可.

代码

map<pii, int> cnt;

void norm(pii& x) {
  if (!x.first) {
    x.first = 0, x.second = 1;
    return;
  }

  if (x.first / abs(x.first) * x.second / abs(x.second) < 0)
    x.first = -abs(x.first), x.second = abs(x.second);
  else x.first = abs(x.first), x.second = abs(x.second);

  int g = abs(gcd(x.first, x.second));
  x.first /= g, x.second /= g;
}

void solve() {
  int n; cin >> n;
  vi a(n), b(n);
  for (auto& ai : a) cin >> ai;
  for (auto& bi : b) cin >> bi;

  int cnt0 = 0;  // 0的个数
  int ans = 0;
  for (int i = 0; i < n; i++) {
    if (a[i] == 0) cnt0 += !b[i];
    else{
      pii d = { -b[i],a[i] };
      norm(d);
      ans = max(ans, ++cnt[d]);
    }
  }
  cout << cnt0 + ans;
}

int main() {
  solve();
}

---

---

1155C. Alarm Clocks Everywhere

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

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

给定一长度为$n(2\le n\le 3\mathrm{e}5)$的序列$x_1,\cdots ,x_n(1\le x_i\le 1\mathrm{e}18)$和一长度为$m(1\le m\le 3\mathrm{e}5)$的序列$p_1,\cdots ,p_m(1\le p_i\le 1\mathrm{e}18)$.构造一个等差数列,使得其首项$y\in \mathbb{Z}$,公差$\in p[]$,且$x[]$中的所有数都是等差数列中的项.若有解,第一行输出"YES",第二行输出任一组首项和公差的下标;否则输出"NO".

思路

显然可取$y=x_i$.考察$x[]$的差分序列$diff[]$,则公差$p$需为所有$dif{f}_i$的约数.

代码

void solve() {
  int n, m; cin >> n >> m;
  vl x(n), p(m);
  for (auto& xi : x) cin >> xi;
  for (auto& pi : p) cin >> pi;

  ll g = x[1] - x[0];  // 序列x[]的差分序列的gcd
  for (int i = 2; i < n; i++) g = gcd(g, x[i] - x[i - 1]);

  for (int i = 0; i < m; i++) {
    if (g % p[i] == 0) {
      cout << "YES" << endl;
      cout << x[0] << ' ' << i + 1;
      return;
    }
  }
  cout << "NO";
}

int main() {
  solve();
}

---

---

1165D. Almost All Divisors

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

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

有$t(1\le t\le 25)$组测试数据.每组测试数据给定一个长度为$n$的序列$d_1,\cdots ,d_n(2\le d_i\le 1\mathrm{e}6)$,表示某数$x$的所有非平凡约数,其中$d_i$两两相异.若存在这样的数$x$,输出$x$;否则输出$-1$.

思路

对数$x$,若存在约数$d$,则也存在约数$\frac{n}{d}$,则$d[]$中最小元素与最大元素之积即$x$.因$x$的约数一定,则原序列$d[]$是$x$的非平凡约数序列的某个排列,将它们排序后判断是否相等即可.

代码

void solve() {
  int n; cin >> n;
  vl d(n);
  for (auto& di : d) cin >> di;

  sort(all(d));
  ll x = (ll)d[0] * d[n - 1];

  vl divisors;
  for (int d = 2; (ll)d * d <= x; d++) {
    if (x % d == 0) {
      divisors.push_back(d);
      if (d != x / d) divisors.push_back(x / d);
    }
  }
  sort(all(divisors));
  cout << (d == divisors ? x : -1) << endl;
}

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

---

---

1174C. Ehab and a Special Coloring Problem

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

题意

给定一个整数$n(2\le n\le 1\mathrm{e}5)$.构造一个序列$a_2,\cdots,a_n\ \ (1\leq a_i\leq n)\ s.t.\ $对$\forall i\in[2,n]$,若$\gcd(i,j)=1$,则$a_i\neq a_j$,且$a[]$的最大值最小.

思路

显然可让每个数与其最小素因子取值相等,且每个素因子下标的$a[i]$应取尽量小的不相等的数.类似于筛法的过程染色即可.

代码

const int MAXN = 1e5 + 5;
int ans[MAXN], idx;

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

  for (int i = 2; i <= n; i++) {
    if (!ans[i]) {
      ans[i] = ++idx;
      for (int j = 2 * i; j <= n; j += i) ans[j] = ans[i];
    }
    cout << ans[i] << ' ';
  }
}

int main() {
  solve();
}

---

---

1178D. Prime Graph

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

题意

给定一个整数$n$,构造一个包含$n$个节点的简单无向图(无重边无自环),要求边数为素数,每个节点的度数也为素数.

如下图,图$1$和图$3$是$n=4$时的合法解,而图$2$不是.

第一行输入一个整数$n(3\le n\le 1000)$.

若无解,则输出$-1$;否则第一行输出一个素数$m\left(2\le m\le \frac{n(n-1)}{2}\right)$,表示边数.接下来$m$行每行输出两个整数$u,v(1\le u,v\le n)$,表示节点$u$与$v$间存在无向边.若有多组解,输出任一组.

思路

一定有解,下面给出构造:

①$n=3$时,三角形是唯一解.

②$n\ge 4$时,先构造一个环$1\leftrightarrow 2\leftrightarrow \cdots \leftrightarrow n\leftrightarrow 1$,此时每个节点的度数都为$2$,但边数$n$可能非素数.

​ 若此时$n$非素数,设$n$的下一个素数为$m$,则只需对任意$i\in [1,m-n]$,加一条边$i\leftrightarrow \left(i+\frac{n}{2}\right)$.

代码

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

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

  int m = n;
  while (!check(m)) m++;

  cout << m << endl;
  for (int i = 1; i < n; i++) cout << i << ' ' << i + 1 << endl;
  cout << "1 " << n << endl;

  for (int i = 1; i <= m - n; i++) cout << i << ' ' << i + n / 2 << endl;
}

int main() {
  solve();
}

---

---

1195D1. Submarine in the Rybinsk Sea (easy edition)

原题指路:https://codeforces.com/problemset/problem/1195/D1

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

给定两个整数的十进制表示(无前导零)$a_1,\cdots ,a_p$和$b_1,\cdots ,b_q$.定义函数$f(a_1,\cdots ,a_p,b_1,\cdots ,b_q)=\{\begin{array}{l}{a}_1{a}_2\cdots a_{p-q+1}{b}_1{a}_{p-q+2}{b}_2\cdots a_{p-1}{b}_{q-1}{a}_p{b}_q,p\ge q\\ b_1{b}_2\cdots b_{q-p}{a}_1{b}_{q-p+1}{a}_2\cdots a_{p-1}{b}_{q-1}{a}_p{b}_q\end{array}$.如$f(1111,2222)=12121212,f(7777,888)=7787878,f(111,2222)=(2121212)$.

给定一个长度为$n(1\le n\le 1\mathrm{e}5)$的序列$a_1,\cdots ,a_n(1\le a_i\le 1\mathrm{e}9)$,其中每个$a_i$的十进制表示(无前导零)等长.求$\sum _{i=1}^n\sum _{j=1}^{n}f(a_i,a_j)$,答案对$998244353$取模.

思路

显然有$f(a_i,a_j)$也会有$f(a_j,a_i)$.因每个$a_i$的十进制表示等长,注意到$f(a_1\cdots a_p,b_1\cdots ,b_p)=a_1{b}_1{a}_2{b}_2\cdots b_p{b}_p,f(b_1\cdots ,b_p,a_1\cdots a_p)=b_1{a}_1{b}_2{a}_2\cdots b_p{a}_p$,即$f(a_i,a_j)+f(a_j,a_i)=f(a_i,a_i)+f(a_j,a_j)$,则$\sum _{i=1}^n\sum _{j=1}^{n}f(a_i,a_j)=n\sum _{i=1}^{n}f(a_i,a_i)$.

代码

const int MOD = 998244353;

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

  int ans = 0;
  for (int i = 0; i < n; i++) {
    string s; cin >> s;
    int len = s.length();
    int cur = 0;
    for (int i = 0; i < len; i++) {
      cur = ((ll)cur * 10 + s[i] - '0') % MOD;
      cur = ((ll)cur * 10 + s[i] - '0') % MOD;
    }
    ans = (ans + cur) % MOD;
  }

  ans = (ll)ans * n % MOD;
  cout << ans;
}

int main() {
  solve();
}

---

---

1200C. Round Corridor

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

题意

双层圆形空间分为若干个房间,内层等分为$n$部分,相邻两部分间有墙;外层等分为$m$部分,相邻两部分间有墙.如下图为$n=4,m=6$时的房间分布:

现有$q$个询问,每个询问输入四个整数$s_x,s_y,e_x,e_y(1\le s_x,e_x\le 2)$,若$s_x=1$,则$s_y\in [1,n]$;否则$s_y\in [1,m]$,$e_y$的范围同理,表示询问是否能从房间$(s_x,s_y)$走到$(e_x,e_y)$,若能则输出"YES";否则输出"NO".

第一行输入三个整数$n,m,q(1\le n,m\le 1\mathrm{e}18,1\le q\le 1\mathrm{e}4)$.接下来$q$行每行输入四个整数$s_x,s_y,e_x,e_y(1\le s_x,e_x\le 2)$,若$s_x=1$,则$s_y\in [1,n]$;否则$s_y\in [1,m]$,$e_y$的范围同理,表示一个询问.

思路

墙的位置$\left(1,\frac{1}{n}\right),\left(1,\frac{2}{n}\right),\cdots ,\left(1,\frac{n}{n}\right),\left(2,\frac{1}{m}\right),\left(2,\frac{2}{m}\right),\cdots ,\left(2,\frac{m}{m}\right)$.对某个分数$x$,若同时有墙$(1,x)$和$(2,x)$,则无法从房间$y<x$到房间$z>x$.

设$g=\gcd (n,m)$,则内外层同时有墙的位置为$\frac{1}{g},\frac{2}{g},\cdots ,\frac{g}{g}$,这些墙将空间分为$g$个部分,每个部分内部可互相到达,不同部分不能互相到达.

①$x=1$时,$\left(1,1\right),\left(1,2\right),\cdots ,\left(1,\frac{n}{g}\right)$属于第一组;$\left(1,\frac{n}{g}+1\right),\left(1,\frac{n}{g}+2\right),\left(1,\frac{2n}{g}\right)$属于第二组;$\cdots$.

②$x=2$时,$(2,1),(2,1),\cdots ,\left(2,\frac{m}{g}\right)$属于第一组;$\left(2,\frac{m}{g}+1\right),\left(2,\frac{m}{g}+2\right),\cdots ,\left(2,\frac{2m}{g}\right)$属于第二组;$\cdots$.

代码

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

  ll g = gcd(n, m);
  CaseT{
    ll sx, sy, ex, ey; cin >> sx >> sy >> ex >> ey;
    sy--, ey--;

    if (sx == 1) sy /= n / g;
    else sy /= m / g;
    if (ex == 1) ey /= n / g;
    else ey /= m / g;
    cout << (sy == ey ? "YES" : "NO") << endl;
  }
}

int main() {
  solve();
}

---

---

1209B. Koala and Lights

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

题意

编号$1\sim n$的$n$盏灯排成一排,每盏灯有打开和关闭两种状态,其中$i$号灯有两个参数$a_i,b_i$,表示它会在$b_i,b_i+a_i,b_i+2a_i,\cdots$时刻反转状态.经过充分长的时间,求每个时刻打开的灯的数量的最大值.

第一行输入一个整数$n(1\le n\le 100)$.第二行输入一个长度为$n$的$0-1$串,其中第$i$个字符为$0$表示初始时$i$号灯关闭,为$1$表示初始时$i$号灯打开.接下来$n$行每行输入两个整数$a,b(1\le a,b\le 5)$,表示每盏灯的参数.

思路

显然所有灯开始循环切换状态后所有灯有公共周期,且该公共周期即所有等差数列的公差的$\mathrm{lcm}$,再加上最晚进入循环切换状态的灯的$b$值作为预周期,则只需模拟预周期和第一个周期灯的变化情况即可.

$\mathrm{lcm}(1,2,3,4,5)=120$,再加上最长的预周期$5$,故只需模拟前$125$个时刻.

代码

const int MAXN = 125;
int cnt[MAXN];  // 每个时刻打开的灯的数量

void solve() {
  int n; cin >> n;
  string s; cin >> s;
  for (int i = 0; i < n; i++) {
    int a, b; cin >> a >> b;
    int c = s[i] - '0';
    for (int j = 0; j < MAXN; j++) {  // 枚举时刻
      cnt[j] += c;
      if (j == b) {  // 进入循环切换状态
        c ^= 1;
        b += a;
      }
    }
  }
  cout << *max_element(cnt, cnt + MAXN);
}

int main() {
  solve();
}

---

---

1220B. Multiplication Table

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

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

给定一个$n\times n(3\le n\le 1000)$的乘法表$M$,其元素$M_{ij}=a_i{a}_j(1\le M_{ij}\le 1\mathrm{e}9)$,其中$a_1,\cdots ,a_n(1\le a_i\le 1\mathrm{e}9)$是一个正整数序列.现给定删去对角线元素后的乘法表$M$,求原序列$a_1,\cdots ,a_n$.

思路

用$\frac{(xy)\cdot (xz)}{yz}=x^2$计算即可.

代码

const int MAXN = 1005;
int a[MAXN][MAXN];

void solve() {
  int n; cin >> n;
  for (int i = 1; i <= n; i++)
    for (int j = 1; j <= n; j++) cin >> a[i][j];

  for (int i = 1; i <= n; i++) {
    int j = i % n + 1, k = j % n + 1;
    cout << (int)sqrt((ll)a[i][j] * a[i][k] / a[j][k]) << ' ';
  }
}

int main() {
  solve();
}

---

---