AcWing Round #21

A. 整数幂

题目链接

---

给定两个整数 $k$ 和 $l$，请判断是否存在一个正整数 $n$，满足 $kn=l$。

输入格式

第一行包含整数 $T$，表示共有 $T$ 组测试数据。

每组数据占两行，第一行包含整数 $k$，第二行包含整数 $l$。

输出格式

每组数据输出一行结果，如果存在 $n$，则输出 YES，否则输出 NO。

数据范围

前三个测试点满足，$2\le k,l\le 100$。
所有测试点满足，$1\le T\le 10$，$2\le k,l\le 2^{31}-1$。

输入样例：

2
5
25
3
8

输出样例：

YES
NO

题目分析：

手速题，暴力枚举。

Code

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;

LL get_c(int a, int k)
{
    LL res = 1;
    for (int i = 0; i < k; i ++ )
        res *= a;
    return res;
}

int main()
{
    int T;
    cin >> T;
    while (T -- )
    {
        int k, l;
        cin >> k >> l;
        bool flag = false;
        for (int i = 0; i <= 50; i ++ )
            if (get_c(k, i) == l) 
            {
                flag = true;
                break;
            }
        if (flag) puts("YES");
        else puts("NO");
    }
    return 0;
}

---

B. 变成1

题目链接

---

给定一个二进制数 $x$，在它变为 $1$ 之前，不断对它进行如下操作：

• 如果 $x$ 为奇数，则将 $x$ 加 $1$。

• 如果 $x$ 为偶数，则将 x 除以 $2$。

请问，多少次操作后，$x$ 会变为 $1$。

输入格式

共一行，一个 01 字符串，表示二进制数 $x$。

输出格式

一个整数，表示所需操作次数。

数据范围

前六个测试点满足，$x$ 的位数不超过 $11$。
所有测试点满足，$x$ 的首位不为 $0$，且位数不超过 $10^6$。

输入样例1：

1

输出样例1：

0

输入样例2：

1001001

输出样例2：

12

输入样例3：

101110

输出样例3：

8

题目分析：

由于 $x$ 最多可以是 $10^6$ 位，因此需要使用高精度根据题意进行模拟。

但实际上可以不用高精度的方式，转而对二进制数串进行操作，利用双指针的思想处理。

例如对于一个给定的二进制数串 101110，从后往前操作：

手动模拟一下可以发现规律：若当前这一位是 $0$ 时，直接操作数 + 1，然后右移一位；若当前这一位是 $1$，需要加 $1$，那就需要找到所有连续的 $1$ 的个数，因为加上 $1$ 之后这些连续的 $1$ 都会由于进位的原因全部变成 $0$，那么操作数就等于连续 $1$ 的个数加上 $1$。

这样，我们就可以设置两个指针 $i$ 和 $j$，然后累加操作步数即可，具体步骤见代码。

Code

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

int main()
{
    string str;
    cin >> str;
    reverse(str.begin(), str.end()); // 将 01 串倒过来操作
    int cnt = 0;
      
    for (int i = 0; i < str.size() - 1; i ++ )
    {
        if (str[i] == '0') ++ cnt; 
        else
        {
            int j = i;
            while (j < str.size() && str[j] == '1') j ++ ;
            cnt += j - i + 1;
            str[j] = '1';
            i = j - 1;
        }
    }
    cout << cnt << endl;
    return 0;
}

---

C. 最大公约数

题目链接

---

给定两个正整数 $a,m$，其中 $a<m$。

请你计算，有多少个小于 $m$ 的非负整数 $x$ 满足：

$$\gcd (a,m)=\gcd (a+x,m)$$

输入格式

第一行包含整数 $T$，表示共有 $T$ 组测试数据。

每组数据占一行，包含两个整数 $a,m$。

输出格式

每组数据输出一行结果，一个整数，表示满足条件的非负整数 $x$ 的个数。

数据范围

前三个测试点满足，$1\le T\le 10$。
所有测试点满足，$1\le T\le 50$，$1\le a<m\le 10^{10}$。

输入样例：

3
4 9
5 10
42 9999999967

输出样例：

6
1
9999999966

题目分析：

$a<m$。令 $d=\gcd (a,m)=\gcd (a+x,m)$，则 $d\mid a,d\mid m$，而又 $\because$ $d\mid a,d\mid m$，且 $d\mid (a+x)$，$\therefore d\mid x$。

令 $a^{\prime }=\frac{a}{d},m^{\prime }=\frac{m}{d},x^{\prime }=\frac{x}{d}$，则 $\gcd (a^{\prime }+x^{\prime },m^{\prime })=1$，问题转化为：有多少个满足要求的 $x^{\prime }$ 使得 $(a^{\prime }+x^{\prime })$ 与 $m^{\prime }$ 互质。

$\because 0\le x<m$，$\therefore 0\le x^{\prime }<m$，两边同加上 $a^{\prime }$，得 $a^{\prime }\le x^{\prime }<a^{\prime }+m^{\prime }$，即求 $[a^{\prime },a^{\prime }+m^{\prime })$ 中共有多少个数与 $m^{\prime }$ 互质。

图中红色的两段 $[0,a^{\prime })$ 和 $[m^{\prime },a^{\prime }+m^{\prime })$ 在 $\mathrm{m}\mathrm{o}\mathrm{d}{m}^{\prime }$ 的意义下数值是一样的，那么与之互质的数的情况也就是一样的，故我们求 $[a^{\prime },a^{\prime }+m^{\prime })$ 中有多少个与 $m^{\prime }$ 互质的数，就可以转化为求 $[0,m^{\prime })$ 中与 $m^{\prime }$ 互质的数的个数，求一下 $m^{\prime }$ 的欧拉函数 $\phi (m^{\prime })$ 即可。

关于欧拉函数的求解和证明过程，可以参考：AcWing 873. 欧拉函数。

注意，本题不可使用 筛法求欧拉函数 进行求解，筛法求欧拉函数时间复杂度是线性的，可以求出 $1\sim n$ 中每个数的欧拉函数，但本题的部分数据可能达到 $10^{10}$，线性筛不能够在时限内求出所需要数的欧拉函数。

时间复杂度分析：

每组数据只需要求出其欧拉函数即可，复杂度为 $\sqrt{n}$，一共有 $50$ 组数据，值最大可达 $10^{10}$，因此总的计算量约为 $50\times \sqrt{10^{10}}\approx 5\times 10^6$。

Code

#include <iostream>
#include <cstring>
#include <algorithm>

using namespace std;

typedef long long LL;

LL gcd(LL a, LL b)
{
    return b ? gcd(b, a % b) : a;
}

LL phi(LL m)
{
    LL res = m;
    for (int i = 2; i <= m / i; i ++ )
        if (m % i == 0)
        {
            res = res / i * (i - 1);
            while (m % i == 0) m /= i;
        }
    if (m > 1) res = res / m * (m - 1);
    return res;
}

int main()
{
    int T;
    cin >> T;
    while (T -- )
    {
        LL a, m;
        cin >> a >> m;
        m = m / gcd(a, m);
        cout << phi(m) << endl;
    }
    return 0;
}