AtCoder Regular Contest 182 A~D

A.Chmax Rush!（枚举）

题意：

有一个长度为$N$的整数序列$S$。最初，$S$的所有元素都是$0$。

同时给你两个长度为$Q$的整数序列：$P=(P_1,P_2,\dots ,P_Q)$和$V=(V_1,V_2,\dots ,V_Q)$。

斯努克希望依次对序列$S$执行$Q$次操作。第$i$个操作如下：

• 执行以下操作之一：
  • 将每个元素$S_1,S_2,\dots ,S_{P_i}$替换为$V_i$。但是，在执行此操作前，如果$S_1,S_2,\dots ,S_{P_i}$中存在严格大于$V_i$的元素，斯努克就会开始哭泣。
  • 将每个元素$S_{P_i},S_{P_i+1},\dots ,S_N$替换为$V_i$。但是，在此操作之前，如果$S_{P_i},S_{P_i+1},\dots ,S_N$中有元素严格大于$V_i$，斯努克就会开始哭泣。

求在$Q$的操作序列中，斯努克可以在不哭泣的情况下完成所有操作的次数，答案对$998244353$取模。

当且仅当有$1\le i\le Q$使得第$i$次操作的选择不同时，两个操作序列是不同的。

分析：

考虑对于任意一个$i<j$，如果$V_i>V_j$，那么$i$与$j$就会互相干扰。

如果$P_i=P_j$，那么$j$不管怎么选操作最后在$P_i$这个点都会出现$S_{P_i}>V_j$的现象，此时答案为$0$。

如果$P_i\ne P_j$，那么让$P$小的那个执行操作$1\to P$，大的那个$P\to n$。

记录一个方向数组，如果$s{d}_i=1$代表必须往$n$操作，$s{d}_i=-1$代表必须往$1$操作，$s{d}_i=0$代表无约束，枚举$i,j$即可，如果出现既要一个操作向左又要向右也是无解的。

答案统计就遍历$1\to Q$，$ans$初始设为$1$，如果$s{d}_i=0$就把答案乘$2$。

代码：

#include<bits/stdc++.h>

using namespace std;
const int mod = 998244353;
const int N = 5005;
int P[N], V[N];
int n, q, sd[N], ans = 1;

int main() {
    cin >> n >> q;
    for (int i = 1; i <= q; i++) {
        cin >> P[i] >> V[i];
    }
    for (int j = 1; j <= q; j++) {
        for (int i = j + 1; i <= q; i++) {
            if (V[i] >= V[j])
                continue;
            if (P[i] == P[j]) {
                cout << 0;
                return 0;
            }
            if (P[j] < P[i]) {
                if (sd[i] && sd[i] != 1) {
                    cout << 0;
                    return 0;
                }
                if (sd[j] && sd[j] != -1) {
                    cout << 0;
                    return 0;
                }
                sd[j] = -1, sd[i] = 1;
            } else {
                if (sd[i] && sd[i] != -1) {
                    cout << 0;
                    return 0;
                }
                if (sd[j] && sd[j] != 1) {
                    cout << 0;
                    return 0;
                }
                sd[j] = 1, sd[i] = -1;
            }
        }
    }
    for (int i = 1; i <= q; i++)
        ans = ans * ((1 + (sd[i] == 0))) % mod;
    cout << ans << endl;
    return 0;
}

B.|{floor(A_i/2^k)}|（构造）

题意：

给你正整数$N$和$K$。

长度为$N$的整数序列，其中所有元素都在$1$和$2^K-1$之间（含），称为良好序列。

良好序列$A=(A_1,A_2,\dots ,A_N)$的得分定义如下：

• 使用介于$1$和$N$（含）之间的整数$i$和非负整数$k$可以表示为$\lfloor \frac{A_i}{2^k}\rfloor$的不同整数的个数。

例如，对于$A=(3,5)$，五个整数可以表示为$\lfloor \frac{A_i}{2^k}\rfloor$：$0$、$1$、$2$、$3$和$5$，因此得数为$5$。

找出一个得分最高的好序列。

每个输入文件都有$T$个测试用例需要解决。

分析：

看到除以$2^k$，联想到与位运算有关，想要得到更多的分数那么右移后得到的数就要越多，所以需要构造出更多不同的二进制下位数不超过$k$的数。

从最高位开始，使用类似线段树建树的方式，每次为左边的区间按位或上$1$，右边的区间不进行处理。$k$为$0$时直接结束，若递归到叶子节点时未满$k$位，后面全部填$1$即可。这样二分着构造出来的所有数均不同且不断右移后能得到更多的结果。

代码：

#include<bits/stdc++.h>

using namespace std;
const int N = 1e5 + 5;
int ans[N << 1];

void solve(int l, int r, int k) {
    if (!k)
        return;
    if (l == r) {
        while (k--)
            ans[l] = (ans[l] << 1) + 1;
        return;
    }
    int mid = (l + r) >> 1;
    for (int i = l; i <= mid; i++)
        ans[i] = (ans[i] << 1) + 1;
    for (int i = mid + 1; i <= r; i++)
        ans[i] <<= 1;
    solve(l, mid, k - 1);
    solve(mid + 1, r, k - 1);
}

int n, k;

int main() {
    int T;
    cin >> T;
    while (T--) {
        cin >> n >> k;
        for (int i = 1; i <= n; i++)
            ans[i] = 1;
        solve(1, n, k - 1);
        for (int i = 1; i <= n; i++)
            cout << ans[i] << " ";
        cout << endl;
    }
    return 0;
}

C.Sum of Number of Divisors of Product（数学）

题意：

长度在$1$和$N$之间（含）的整数序列，其中每个元素的长度在$1$和$M$之间（含），称为好序列。

好序列的分数定义为$X$的正除数，其中$X$是序列中各元素的乘积。

好序列有$\sum _{k=1}^N{M}^k$个。求所有这些序列的分数之和对$998244353$取模的结果。

分析：

$m\le 16$，先找出所有质数$2,3,5,7,11,13$，一共有$6$个，考虑状压也只有$64$。

进一步发现$n$巨大，但是状压状态很少，可以想到矩阵快速幂。考虑怎么构造原$DP$序列和矩阵进行转移：

首先要求对于每个长度求和所以肯定要设一维来统计前缀和。

其次我们写出转移，并构造转移矩阵：发现如果设“某一状态为某一质数有没有”无法转移，所以我们先令第$i(0\le i<6)$个质数的编号为$i$，个数为$a_i$（$a_i=0$表示没有这个质数），则显然一个序列的答案为$(a_0+1)(a_1+1)(a_2+1)(a_3+1)(a_4+1)(a_5+1)$而如果我们在序列末尾加入一个$6$则答案就变成$(a_0+2)(a_1+2)(a_2+1)(a_3+1)(a_4+1)(a_5+1)$。

根据$c(a+x)(b+y)=c(ab+ay+bx+xy)$，设$c=(a_2+1)(a_3+1)(a_4+1)(a_5+1)$原式子即为：$(a_0+2)(a_1+2)c=(a_0+1)(a_1+1)c+(a_0+1)c+(a_1+1)c+c$

设某一状态$f_S$表示$(a_i+1)(a_j+1)(a_k+1)\cdots (a_l+1),(i,j,k,\dots ,l\in S)$的和。

那么转移过程就变成我们枚举每一位和$[1,m]$间的每一位数字，来看看对于原序列的系数是多少，然后给转移矩阵增加即可，因为一个数字最多有两个不同的质数，找两个变量记录一下即可。

式子形如：$f_S=f_S+f_{S-i}+f_{S-j}+f_{S-i-j}$。

代码：

#include<bits/stdc++.h>

typedef long long LL;
using namespace std;
const LL N = 66;
const LL mod = 998244353;

struct Matrix {
    LL n, m, num[N][N];

    Matrix() {
        n = m = 0;
        memset(num, 0, sizeof(num));
    }
};

Matrix operator*(const Matrix &x, const Matrix &y) {
    Matrix c;
    c.n = x.n;
    c.m = y.m;
    for (LL i = 0; i < x.n; i++)
        for (LL j = 0; j < y.m; j++)
            for (LL k = 0; k < x.m; k++)
                (c.num[i][j] += x.num[i][k] * y.num[k][j] % mod) %= mod;
    return c;
}

Matrix ksm(Matrix a, LL b) {
    Matrix t;
    t.n = t.m = a.n;
    for (LL i = 0; i < a.n; i++)t.num[i][i] = 1;
    for (; b; b >>= 1, a = a * a)if (b & 1)t = t * a;
    return t;
}

LL n, m;
LL id[17];
Matrix a, b;

int main() {
    id[2] = 1;
    id[3] = 2;
    id[5] = 3;
    id[7] = 4;
    id[11] = 5;
    id[13] = 6;
    cin >> n >> m;
    a.n = 1, a.m = 65;
    for (LL i = 0; i < 64; i++)
        a.num[0][i] = 1;
    b.n = b.m = 65;
    b.num[64][64] = b.num[63][64] = 1;
    for (LL i = 0; i < 64; i++)
        for (LL k = 1; k <= m; k++) {
            LL t = k, t1 = -1, t2 = -1, s1 = 0, s2 = 0;
            for (LL l = 2; l <= t; l++)
                if (t % l == 0) {
                    if (~t1) {
                        t2 = id[l] - 1;
                        while (t % l == 0)t /= l, s2++;
                    } else {
                        t1 = id[l] - 1;
                        while (t % l == 0)t /= l, s1++;
                    }
                }
            b.num[i][i]++;
            if (t1 >= 0 && (i >> t1 & 1))b.num[i ^ (1 << t1)][i] += s1;
            if (t2 >= 0 && (i >> t2 & 1))b.num[i ^ (1 << t2)][i] += s2;
            if (t1 >= 0 && (i >> t1 & 1) && t2 >= 0 && (i >> t2 & 1))b.num[i ^ (1 << t1) ^ (1 << t2)][i] += s1 * s2;
        }
    a = a * ksm(b, n + 1);
    cout << a.num[0][64] - 1 << endl;
    return 0;
}

D.Increment Decrement Again（思维）

题意：

没有两个相邻元素相同的整数序列称为良好序列。

给你两个长度为$N$的良好序列：$A=(A_1,A_2,\dots ,A_N)$和$B=(B_1,B_2,\dots ,B_N)$。$A$和$B$中的每个元素都介于$0$和$M-1$之间。

您可以对$A$执行以下任意次数的运算，可能为零次：

• 在$1$和$N$之间选择一个整数$i$，并执行以下操作之一：
  • 设$A_i\leftarrow (A_i+1)\mathrm{mod}M$。
  • 设$A_i\leftarrow (A_i-1)\mathrm{mod}M$。此处为$(-1)\mathrm{mod}M=M-1$。

但是，不能进行使$A$不再是一个好序列的操作。

判断是否有可能使$A$等于$B$，如果有可能，求这样做所需的最小运算次数。

分析：

先特判$m=2$。

然后我们看作$A$不对$m$取模，要求变为$\forall i<n$，$|a_i-a_{i+1}|<m\wedge a_i\ne a_{i+1}$。

可以发现$a$的大小关系是不会变的，所以如果固定了$a_1^{\prime }$，则后面的都是固定的。可以先令$a_1^{\prime }=b_1$，根据大小关系求出任意一组合法的$a^{\prime }$，接下来问题变为最小化$\sum _{i=1}^n|a_i-a_i^{\prime }-km|$，可以二分类似找中位数的方法解决。复杂度$O(nlognlogV)$。

代码：

#include<bits/stdc++.h>

typedef long long LL;
const LL inf = 0x3f3f3f3f;
using namespace std;

LL n, m, L, R, mid, p, a[200010], b[200010], c[200010];

int main() {
    cin >> n >> m;
    for (LL i = 1; i <= n; ++i)
        cin >> a[i];
    for (LL i = 1; i <= n; ++i)
        cin >> b[i];
    if (m == 2) {
        if (a[1] == b[1])
            cout << "0" << endl;
        else
            cout << "-1" << endl;
        return 0;
    }
    c[1] = b[1];
    for (LL i = 2; i <= n; ++i) {
        LL val = (b[i] - b[i - 1] + m) % m;
        if (a[i] > a[i - 1])
            c[i] = c[i - 1] + val;
        else
            c[i] = c[i - 1] + val - m;
    }
    for (LL i = 1; i <= n; ++i)
        c[i] = a[i] - c[i];
    sort(c + 1, c + 1 + n);
    L = -1000 * inf;
    R = 1000 * inf;
    while (L < R) {
        mid = L + ((R - L + 1) >> 1);
        p = lower_bound(c + 1, c + 1 + n, mid * m) - c;
        if (p <= (n + 1) / 2)
            L = mid;
        else
            R = mid - 1;
    }
    LL ans = 0, s = 0;
    for (LL i = 1; i <= n; ++i)
        ans += abs(c[i] - L * m);
    for (LL i = 1; i <= n; ++i)
        s += abs(c[i] - (L + 1) * m);
    cout << min(ans, s) << endl;
    return 0;
}

赛后交流

在比赛结束后，会在交流群中给出比赛题解，同学们可以在赛后查看题解进行补题。

群号： 704572101，赛后大家可以一起交流做题思路，分享做题技巧，欢迎大家的加入。