[Codeforces] number theory (R1600) Part.8

[Codeforces] number theory (R1600) Part.8

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

1242A. Tile Painting

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

题意

编号$1\sim n$的$n$的格子排成一行.现对这些格子染色,定义一种染色方案是好的,如果对$\forall i,j(i\ne j)s.t.|j-i|$是$n$的一个$>1$的约数,则它们颜色相同.给定$n(1\le n\le 1\mathrm{e}12)$,问至多能用多少种颜色染色,使得染色方案是好的.

思路

①若$n=p^k(p\in primes)$,因$n$模$p$有$p$种余数,故$ans=p$.

②若$n=pq(p,q>1)$且$\gcd (p,q)=1$,则$ans=1$.

​ [证] 对$\forall i,j(i\ne j)$,由CRT:$\exists x\in [1,n]s.t.x\equiv i(\mathrm{mod}p),x\equiv j(\mathrm{mod}p)$,

​ 则下标$i,j$与下标$x$颜色相同,进而所有格子颜色相同.

代码

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

  for (int d = 2; (ll)d * d <= n; d++) {
    if (n % d == 0) {
      while (n % d == 0) n /= d;
      if (n > 1) {
        cout << 1;
        return;
      }
      else {
        cout << d;
        return;
      }
    }
  }
  cout << n;
}

int main() {
  solve();
}

---

---

1263C. Everyone is a Winner!

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

题意

有$t(1\le t\le 10)$组测试数据.每组测试数据输入一个整数$n(1\le n\le 1\mathrm{e}9)$.对所有正整数$k$,问$\lfloor \frac{n}{k}\rfloor$有多少种不同的取值,第一行输出不同取值的个数,第二行升序输出所有取值.

思路

所有$0\le x<\lfloor \sqrt{n}\rfloor$都是解.

[证] 方程$\lfloor \frac{n}{k}\rfloor =x$等价于$x\le \frac{n}{k}<x+1$,解得$k\in \left(\frac{n}{x+1},\frac{n}{x}\right]$,区间长度为$\frac{n}{x^2+x}$.

注意到$x<\lfloor \sqrt{n}\rfloor$时,$\frac{n}{x^2+x}>1$,则区间内至少存在一个整数$k=\lfloor \frac{n}{x}\rfloor$,故所有$x\in \left[0,\lfloor \sqrt{n}\rfloor \right)$都是解.

整除分块求所有解即可.

代码

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

  vi ans;
  for (int l = 1, r, t; l <= n; l = r + 1) {
    r = (t = n / l) ? min(n / t, n) : n;
    ans.push_back(n / l);
  }
  ans.push_back(0);

  reverse(all(ans));  // 注意翻转
  cout << ans.size() << endl;
  for (int i : ans) cout << i << ' ';
  cout << endl;
}

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

---

---

1266C. Diverse Matrix

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

题意

设$A$是一个$r\times c$的整型矩阵,其中元素未必不同.构造一个长度为$(r+c)$的序列$b[]$,其中前$r$个元素是$A$每行的$\gcd$,后$c$个元素是$A$每列的$\gcd$.称序列$b[]$是好的,如果其中的元素两两相异.定义好的序列的权值为其中元素的最大值.

给定两整数$r$和$c(1\le r,c\le 500)$,构造一个$r\times c$矩阵$A$使得其产生的序列$b[]$是好的,且其权值最小.无解输出$0$.若有多组解,输出任一组解.

思路

①$r=c=1$时无解.

②$r=1$时,最优解为$A=[23\cdots c+1]$,不取$1\sim c$是因为$1$会重复.

③$c=1$时,最优解为$A=[23\cdots c+1{]}^T$.

④$r,c\ge 2$时,令$a_{ij}=i(j+r)$,则 b i = gcd ⁡ { i ( r + 1 ) , ⋯   , i ( r + c ) } = i ⋅ gcd ⁡ { r + 1 , ⋯   , r + c } = i b_i=\gcd\{i(r+1),\cdots,i(r+c)\}=i\cdot \gcd\{r+1,\cdots,r+c\}=i bi​=gcd{i(r+1),⋯,i(r+c)}=i⋅gcd{r+1,⋯,r+c}=i.

​ 此时$b[]$的权值为$(r+c)$,显然这是最小的.

代码

void solve() {
  int r, c; cin >> r >> c;

  if (r * c == 1) {
    cout << 0;
    return;
  }

  if (r == 1) {
    for (int i = 2; i <= c + 1; i++) cout << i << ' ';
    return;
  }
  else if (c == 1) {
    for (int i = 2; i <= r + 1; i++) cout << i << endl;
    return;
  }

  for (int i = 1; i <= r; i++)
    for (int j = 1; j <= c; j++) cout << i * (j + r) << " \n"[j == c];
}

int main() {
  solve();
}

---

---

1285C. Fadi and LCM

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

题意

给定一个整数$x(1\le x\le 1\mathrm{e}12)$,求两正整数$a,bs.t.\mathrm{lcm}(a,b)=x$且 max ⁡ { a , b } \max\{a,b\} max{a,b}最小.若有多组解,输出任一组.

思路I

将$x$素因数分解后,枚举将素因子分配给$a$和$b$的所有方案,求 max ⁡ { a , b } \max\{a,b\} max{a,b}的最小值.

代码I

vl factors;

void get(ll n) {  // 对x素因数分解
  for (ll d = 2; d * d <= n; d++) {
    if (n % d == 0) {
      ll cur = 1;
      while (n % d == 0) n /= d, cur *= d;
      factors.push_back(cur);
    }
  }
  if (n > 1) factors.push_back(n);
}

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

  get(x);

  ll ansa = INFF, ansb = INFF;
  int n = factors.size();
  for (int i = 0; i < (1 << n); i++) {
    ll a = 1, b = 1;
    for (int j = 0; j < n; j++) {
      if (i >> j & 1) a *= factors[j];
      else b *= factors[j];
    }

    if (max(a, b) < max(ansa, ansb)) ansa = a, ansb = b;
  }
  cout << ansa << ' ' << ansb;
}

int main() {
  solve();
}

---

思路II

枚举$x$的约数$d$,检查是否有$\mathrm{lcm}\left(d,\frac{x}{d}\right)=x$,更新 max ⁡ { a , b } \max\{a,b\} max{a,b}的最小值.

代码II

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

  ll ans = 0;
  for (ll d = 1; d * d <= x; d++) 
    if (x % d == 0 && lcm(d, x / d) == x) ans = d;
  cout << ans << ' ' << x / ans;
}

int main() {
    solve();
}

---

---

1294C. Product of Three Numbers

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

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

有$t(1\le t\le 100)$组测试数据.每组测试数据输入一个整数$n(2\le n\le 1\mathrm{e}9)$.求三个相异的正整数$a,b,cs.t.a,b,c\ge 2$且$n=abc$,若有解则第一行输出"YES",第二行输出任一组解;否则输出"NO".

思路

不妨设$a<b<c$,下面最小化$a$并最大化$c$.设$a$是$n$的$>1$的最小约数,$b$是$\frac{n}{a}$的非$a$的$>1$的最小约数,若$\frac{n}{ab}\ne a,b,1$,则找到一组解.

代码

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

  uset<int> s;
  for (int d = 2; (ll)d * d <= n; d++) {  // 求n的>1的最小约数a
    if (n % d == 0) {
      n /= d;
      s.insert(d);
      break;
    }
  }
  for (int d = 2; (ll)d * d <= n; d++) {  // 求n/a的非a的>1的最小约数b
    if (n % d == 0 && !s.count(d)) {
      n /= d;
      s.insert(d);
      break;
    }
  }

  if (s.size() < 2 || s.count(n) || n == 1) cout << "NO" << endl;
  else {
    cout << "YES" << endl;
    s.insert(n);
    for (auto i : s) cout << i << ' ';
    cout << endl;
  }
}

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

---

---

1305C. Kuroni and Impossible Calculation

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

题意

给定一个长度为$n(2\le n\le 2\mathrm{e}5)$的序列$a_1,\cdots ,a_n(0\le a_i\le 1\mathrm{e}9)$和一个整数$m(1\le m\le 1000)$,求$\prod _{1\le i<j\le n}|a_i-a_j|$,答案对$m$取模.

思路

①$n\le m$时,暴力求即可.

②$n>m$时,因模$m$只有$m$种余数,则至少存在两个数$a_i,a_{j}s.t.a_i\equiv a_j(\mathrm{mod}m)$,故$ans=0$.

代码

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

  if (n > m) cout << 0;
  else {
    int ans = 1;
    for (int i = 0; i < n; i++)
      for (int j = i + 1; j < n; j++) ans = (ll)abs(a[i] - a[j]) % m * ans % m;
    cout << ans;
  }
}

int main() {
  solve();
}

---

---

1312C. Adding Powers

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

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

序列$v_1,\cdots ,v_n$初始时都为零.现有操作:在第$i$步(下标从$0$开始)时,可选择跳过此步,也可选择一个下标$pos\in [1,n]$,令$v_{pos}+=k^i$.给定一个序列$a_1,\cdots ,a_n$,问若干次操作后能否将序列$v[]$变为序列$a[]$,若能则输出"YES";否则输出"NO".

有$t(1\le t\le 1000)$组测试数据.每组测试数据第一行输入两个整数$n,k(1\le n\le 30,2\le k\le 100)$.第二行输入$n$个整数$a_1,\cdots ,a_n(0\le a_i\le 1\mathrm{e}16)$.

思路

将每个$a_i$表示成$k$进制数,检查是否有相同的$k$的幂次.

代码

const int MAXN = 35;
int n, k;
ll a[MAXN];

bool check() {
  uset<int> s;  // 记录用过的k的幂次
  for (int i = 0; i < n; i++) {
    ll x = a[i];
    int d = 0;  // k的幂次
    while (x) {
      if (x % k != 0 && x % k != 1) return false;

      if (x % k == 1) {
        if (s.count(d)) return false;
        s.insert(d);
      }
      d++;
      x /= k;
    }
  }
  return true;
}

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

  cout << (check() ? "YES" : "NO") << endl;    
}

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

---

---

1332B. Composite Coloring

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

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

称一个正整数是好的,如果它能表示为两$>1$的正整数之积.称一个序列是好的,如果其所有元素都是好的.

给定一个长度为$n$的好的序列$a_1,\cdots ,a_n$.现选定一个正整数$m\in [1,m]$,用编号$1\sim m$的$m$种颜色给序列染色,使得每种颜色至少用一次,序列中每个元素被染有且只有一种颜色,任意两个相同颜色的元素的$\gcd >1$.注意相等的元素可染不同色.可以证明当所有$a_i\le 1000$时有解.构造任一组染色方案,第一行输出所用颜色的种数,第二行输出一个染色方案.

有$t(1\le t\le 1000)$组测试数据.每组测试数据第一行输入一个整数$n(1\le n\le 1000)$.第二行输入一个长度为$n$的好的序列$a_1,\cdots ,a_n(4\le a_i\le 1000)$.数据保证所有测试数据的$n$之和不超过$1\mathrm{e}4$.

思路

对每个好的数$k$,存在一个素数$p\le \sqrt{k}s.t.p\mid k$.因$\sqrt{1000}$以内共$11$个素数,则每个数与其素因子染相同颜色即可.

代码

const int MAXN = 1005;
vi res[MAXN];  // 含有最小素因子p的元素的下标
int ans[MAXN];

int get(int x) {  // 求x的最小素因子
  for (int i = 2; i * i <= x; i++)
    if (x % i == 0) return i;
  return x;
}

void solve() {
  int n; cin >> n;
  for (int i = 1; i < MAXN; i++) res[i].clear();

  for (int i = 1; i <= n; i++) {
    int a; cin >> a;
    res[get(a)].push_back(i);
  }

  int color = 0;
  for (int i = 1; i < MAXN; i++) {
    if (res[i].size()) {
      color++;
      for (auto j : res[i]) ans[j] = color;
    }
  }

  cout << color << endl;
  for (int i = 1; i <= n; i++) cout << ans[i] << " \n"[i == n];
}

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

---

---

1342C. Yet Another Counting Problem

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

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

给定两整数$a,b$.现给定两整数$l,r$,表示询问有多少个$x\in [l,r]s.t.(x\mathrm{mod}a)\mathrm{mod}b\ne (x\mathrm{mod}b)\mathrm{mod}a$.

有$t(1\le t\le 100)$组测试数据.每组测试数据第一行输入三个整数$a,b,q(1\le a,b\le 200,1\le q\le 500)$,其中$q$表示询问书.接下来$q$行每行输入两个整数$l,r(1\le l\le r\le 1\mathrm{e}18)$,表示一个询问.

思路

注意到$[(x+ab)\mathrm{mod}a]\mathrm{mod}b=(x\mathrm{mod}a)\mathrm{mod}b$,则$x$有上述性质当且仅当$x\mathrm{mod}ab$有上述性质.

预处理出$len=ab$内满足所需性质的数的个数的前缀和$pre[]$,根据$\lfloor \frac{x}{len}\rfloor$和$x\%len$求$[1,x]$中有多少个数满足所需性质.

代码

const int MAXN = 4e4 + 5;
int pre[MAXN];  // 满足所需性质的数的个数的前缀和
int len;  // a*b

void init(int a,int b) {
  len = a * b;
  for (int i = 1; i <= len; i++) 
    pre[i] = pre[i - 1] + ((i % a) % b != (i % b) % a);
}

ll query(ll x) {
  ll cnt = x / len, rest = x % len;
  return cnt * pre[len] + pre[rest];
}

ll query(ll l, ll r) {
  return query(r) - query(l - 1);
}

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

  CaseT{
    ll l, r; cin >> l >> r;
    cout << query(l, r) << ' ';
  }
  cout << endl;
}

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

---

---

1344A. Hilbert’s Hotel

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

题意

有无限个房间,其中编号$0\sim (n-1)$的$n$个房间分别住着客人$a_0,\cdots ,a_{n-1}$.现有操作:对每个整数$k\in [0,k-1]$,将$k$号房间的客人移动到$(k+a_{k\mathrm{mod}n})$号房间.问操作后$0\sim (n-1)$号房间是否有且只有一个客人.

有$t(1\le t\le 1\mathrm{e}4)$组测试数据.每组测试数据第一行输入一个整数$n(1\le n\le 2\mathrm{e}5)$.第二行输入一个长度为$n$的序列$a_0\cdots ,a_{n-1}(-1\mathrm{e}9\le a_i\le 1\mathrm{e}9)$.数据保证所有测试数据的$n$之和不超过$2\mathrm{e}5$.

思路

$0\sim (n-1)$号房间是否有且只有一个客人等价于没有两个客人在同一房间且没有空房间.

设操作后存在两个客人在同一房间,则$x+a_{x\mathrm{mod}n}=y+a_{y\mathrm{mod}n}$至少有一整数解.

设$x=k_{1}n+i,y=k_{2}n+j$,则原方程$\iff x-y=a_{y\mathrm{mod}n}-a_{x\mathrm{mod}n}\iff k_{1}n+i-k_{2}n-j=a_j-a_i$,

则$(k_1-k_2)n=a_j-a_i+j-i=(a_j+j)-(a_i+i)$.

设$b_i=a_i+i$,则等价于求$(k_1-k_2)n=c_j-c_i$的正整数解$c_i,c_j$,显然这等价于$c_i\equiv c_j(\mathrm{mod}n)$.

预处理出所有$(i+a_i)\mathrm{mod}n$,检查余数是否为$0\sim (n-1)$即可.

代码

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

  sort(all(a));
  bool ok = true;
  for (int i = 0; i < n; i++) ok &= a[i] == i;
  cout << (ok ? "YES" : "NO") << endl;
}

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

---

---