Codeforces Round 757 (Div.2) 解题报告

传送门

$\mathrm{C}\mathrm{o}\mathrm{d}\mathrm{e}\mathrm{f}\mathrm{o}\mathrm{r}\mathrm{c}\mathrm{e}\mathrm{s}\mathrm{R}\mathrm{o}\mathrm{u}\mathrm{n}\mathrm{d}757(\mathrm{D}\mathrm{i}\mathrm{v}.2)$

A Divan and a Store

有 $n$ 个物品，第 $i$ 个商品的价格为 $a_i$。

给你 $l$ 和 $r$，表示商人只会选择价格在 $[l,r]$ 内的物品。

商人有 $k$ 元钱，问你商人最多可以买多少个商品。

共 $T$ 组数据。

对于 $100\%$ 的数据，保证 $1\le T\le 100,1\le n\le 100,1\le l\le r\le 10^9,1leqk\le 10^9,1\le a_i\le 10^9$。

sol

贪心，每次买剩下的价格最小的商品。

首先选出满足条件的物品，放到优先队列里（最小值优先）。

然后一件一件弹出来，计算钱数，若钱数 $>k$ 或队列为空，结束程序。

时间复杂度 $\mathcal{O}(Tn\log n)$。

#include <bits/stdc++.h>

using namespace std;

#define int long long

int t;

int n, l, r, k;

signed main()
{
    scanf("%lld", &t);
    while (t--)
    {
        scanf("%lld%lld%lld%lld", &n, &l, &r, &k);
        priority_queue<int, vector<int>, greater<int>> q;
        for (int i = 1, a; i <= n; ++i)
        {
            scanf("%lld", &a);
            if (a < l || a > r)
                continue;
            q.push(a);
        }
        int ans = 0, sum = 0;
        while (!q.empty())
        {
            int now = q.top();
            q.pop();
            sum += now;
            // cout << now << "\n";
            if (sum <= k)
                ans++;
            else
                break;
        }
        printf("%lld\n", ans);
    }
    return 0;
}

B Divan and a New Project

在一条街上，盖 $n+1$ 栋楼，每栋楼坐标为 $x_i$，满足 $\forall j\ne i,x_i\ne x_j$ 且 $x_i$ 是一个整数。

将所有楼从 $0$ 到 $n$ 标号。

有一人在编号为 $0$ 的楼，分别要去编号为 $i$ 的建筑 $a_i$ 次，这个人往返编号为 $i$ 的建筑一趟花费的时间为 $2\times (|x_i-x_0|)$。

请安排每一栋楼的坐标，使得这个人花费的时间最短。

共 $T$ 组数据。

对于 $100\%$ 的数据，保证 $1\le T\le 10^3,1\le n\le 2\ast 10^5,0\le a_i\le 10^6$。

sol

首先，我们把编号为 $0$ 的楼的坐标定为 $0$ 方便去做。

记答案为 $ans$，则 $ans=\sum _{i=1}^{n}2\ast a_i\ast |x_0-x_i|$。

要使得 $ans$ 最少，就要让 $a_i$ 大的楼放到距离 $x_0$ 近的地方。

那么按照 $a_i$ 从大到小排序，在 $x_0$ 按照左右左右左右…的顺序放楼即可。

具体实现见代码。

时间复杂度 $\mathcal{O}(Tn)$。

#include <bits/stdc++.h>

using namespace std;

#define int long long

const int _ = 2e5 + 7;

int t;

int n;

struct abc
{
    int v, id;
} a[_];

int b[_];

int ans;

bool cmp(abc a, abc b)
{
    return a.v > b.v;
}

signed main()
{
    scanf("%lld", &t);
    while (t--)
    {
        ans = 0;
        scanf("%lld", &n);
        for (int i = 1; i <= n; ++i)
            scanf("%lld", &a[i].v), a[i].id = i;
        sort(a + 1, a + n + 1, cmp);
        int k = 0;
        for (int i = 1; i <= n; i += 2)
        {
            b[a[i].id] = ++k;
        }
        k = 0;
        for (int i = 2; i <= n; i += 2)
        {
            b[a[i].id] = --k;
        }
        for (int i = 1; i <= n; ++i)
        {
            ans += 2 * a[i].v * abs(b[a[i].id]);
        }
        cout << ans << "\n";
        cout << "0 ";
        for (int i = 1; i <= n; ++i)
            cout << b[i] << " ";
        cout << "\n";
    }
    return 0;
}

C Divan and bitwise operations

有一个未知的序列 $a$，现知道 $m$ 个信息，每个星信息为 l r x 的形式给出，表示区间 $[l,r]$ 的按位或为 $x$，保证 $a$ 中每一个数都被包含在区间 $[l,r]$ 至少一次。

请输出序列 $a$ 的所有子序列的异或和 $\mathrm{m}\mathrm{o}\mathrm{d}(10^9+7)$。

共 $T$ 组数据。

对于 $100\%$ 的数据，保证 $1\le T\le 10^3,1\le n,m\le 2\ast 10^5,1\le l\le r\le n,0\le x\le 2^{30}-1$。

sol

显然，我们可以得到整个序列的按位或就是所有 $x$ 的按位或，设为 $S$。

如果 $S$ 的第 $i$ 位为 $0$，贡献即为 $0$。

否则总有一个 $1$，当中恰有一个对应贡献为 $2^i$，总贡献为 $2^i\times 2^{n-1}$。

那么 $Ans=S\times 2^{n-1}$。

时间复杂度 $\mathcal{O}(Tn)$。

具体实现见代码。

#include <bits/stdc++.h>

#define int long long

using namespace std;

inline int read()
{
    int s = 0, w = 1;
    char c = getchar();
    for (; !isdigit(c); c = getchar())
        if (c == '-')
            w = -1;
    for (; isdigit(c); c = getchar())
        s = (s << 1) + (s << 3) + (c ^ 48);
    return s * w;
}

const int mod = 1e9 + 7;

inline int qpow(int x, int y)
{
    int res = 1;
    while (y)
    {
        if (y & 1)
            res = res * x % mod;
        x = x * x % mod;
        y >>= 1;
    }
    return res;
}

int T, n, m;

signed main()
{
    T = read();
    while (T--)
    {
        n = read(), m = read();
        int ans = 0;
        while (m--)
        {
            read(), read();
            ans = ans | read();
        }
        ans = ans * qpow(2, n - 1) % mod;
        printf("%lld\n", ans);
    }
    return 0;
}

D Divan and Kostomuksha

给定序列 $a_1,a_2,...,a_n$，要求重排 $a$，使得

$$\sum _{i=1}^n\gcd (a_1,a_2,...,a_i)$$

最大。

输出这个最大值。

sol

显然，一道 dp。

设 $cn{t}_i$ 表示数组 $a$ 中是 $i$ 的倍数的元素个数，$d{p}_i$ 表示时能得到已含有因数 $i$ 结尾排列能获得的最大值。

此时有转移方程
$$d{p}_i=\underset{j=1}{\max }{f}_{i\times pr{i}_j}+i\times (cn{t}_i-cn{t}_{i\times pr{i}_j})$$
最后答案为满足 $cn{t}_i=n$ 的 $d{p}_i$ 的最大值。

求 $cnt$ 可先筛出所有质数，再枚举 $i$ 和质数集 $pri$，计算 $cn{t}_i$。

时间复杂度 $\mathcal{O}(w\log \log w)$，其中 $w$ 为值域。

具体实现见代码。

#include <bits/stdc++.h>

using namespace std;

#define int long long

inline int read()
{
    int x = 0, f = 1;
    char c = getchar();
    while (c < '0' || c > '9')
    {
        if (c == '-')
            f = -1;
        c = getchar();
    }
    while (c >= '0' && c <= '9')
    {
        x = x * 10 + c - '0';
        c = getchar();
    }
    return x * f;
}

const int _ = 2e7 + 1;

int n;

int ans;

int cnt[_ + 7];

int dp[_ + 7];

bool vis[_];

vector<int> primes;

void init()
{
    for (int i = 2; i < _; ++i)
    {
        if (!vis[i])
            primes.push_back(i);
        for (auto p : primes)
        {
            if (p * i > _)
                break;
            vis[p * i] = 1;
            if (i % p == 0)
                break;
        }
    }
}

signed main()
{
    init();
    n = read();
    for (int i = 1, x; i <= n; ++i)
        cnt[read()]++;
    for (auto p : primes)
        for (int j = _ / p; j >= 1; --j)
            cnt[j] += cnt[j * p];
    for (int i = _ - 1; i >= 1; --i)
    {
        dp[i] = cnt[i] * i;
        for (auto p : primes)
        {
            int j = p * i;
            if (j > _)
                break;
            dp[i] = max(dp[i], dp[j] + i * (cnt[i] - cnt[j]));
        }
        if (cnt[i] == n)
            ans = max(ans, dp[i]);
    }
    printf("%lld\n", ans);
    return 0;
}

E Divan and a Cottage

给定 $n$ 天的气温 $T_i$，设第 $i$ 天温度为 $P$，则第 $i+1$ 天的温度为：

• $P+1(P<T_i)$

• $P-1(P>T_i)$

• $P(P=T_i)$

对第 $i$ 天有 $k_i$ 个询问，问若第 $0$ 天的温度为 $x$ ，那么第 $i$ 天的温度是多少。

强制在线。

sol

显然，所有初始的室温变化后的结果满足单调。

那用一棵动态开点值域线段树维护所有初始的室温变化后的结果即可。

具体地，每次插入一个气温 $T_i$ 后，在线段树上二分以得到变化后的室温为 $T_i$ 的区间 $\left[l,r\right]$，之后 $[-1,l)$ 区间 $+1$，$\left(r,MaxT+1\right]$ 区间 $-1$ 即可。

时间复杂度 $\mathcal{O}(n\log w)$，其中 $w$ 为值域。

#include <bits/stdc++.h>

using namespace std;

typedef long long ll;

inline ll read()
{
    ll x = 0, f = 1;
    char c = getchar();
    while ((c < '0' || c > '9') && (c != '-'))
        c = getchar();
    if (c == '-')
        f = -1, c = getchar();
    while (c >= '0' && c <= '9')
        x = x * 10 + c - '0', c = getchar();
    return x * f;
}

const int N = 2e5 + 10, mod = 1e9 + 1;

int n, m;

const int M = 4e7 + 10;

int tot, mn[M], mx[M], tag[M], ls[M], rs[M];

inline int New(int l, int r)
{
    ++tot, mn[tot] = l, mx[tot] = r;
    return tot;
}

inline void upd(int u, int x)
{
    mn[u] += x, mx[u] += x, tag[u] += x;
}

inline void push_down(int u, int l, int r, int mid)
{
    if (!tag[u])
        return;
    if (!ls[u])
        ls[u] = New(l, mid);
    if (!rs[u])
        rs[u] = New(mid + 1, r);
    upd(ls[u], tag[u]), upd(rs[u], tag[u]);
    tag[u] = 0;
}

inline void update(int &u, int l, int r, int ql, int qr, int x)
{
    if (ql > qr)
        return;
    if (!u)
        u = New(l, r);
    if (l >= ql && r <= qr)
    {
        upd(u, x);
        return;
    }
    int mid = l + r >> 1;
    push_down(u, l, r, mid);
    if (ql <= mid)
        update(ls[u], l, mid, ql, qr, x);
    if (qr > mid)
        update(rs[u], mid + 1, r, ql, qr, x);
    mn[u] = min(!ls[u] ? l : mn[ls[u]], !rs[u] ? mid + 1 : mn[rs[u]]);
    mx[u] = max(!ls[u] ? mid : mx[ls[u]], !rs[u] ? r : mx[rs[u]]);
}

inline int Findl(int u, int l, int r, int x)
{
    if ((!u ? l : mn[u]) >= x)
        return -1;
    if (l == r)
        return l;
    int mid = l + r >> 1;
    push_down(u, l, r, mid);
    int ret = Findl(rs[u], mid + 1, r, x);
    if (ret != -1)
        return ret;
    return Findl(ls[u], l, mid, x);
}

inline int Findr(int u, int l, int r, int x)
{
    if ((!u ? r : mx[u]) <= x)
        return mod;
    if (l == r)
        return l;
    int mid = l + r >> 1;
    push_down(u, l, r, mid);
    int ret = Findr(ls[u], l, mid, x);
    if (ret != mod)
        return ret;
    return Findr(rs[u], mid + 1, r, x);
}

inline int Query(int u, int l, int r, int ql)
{
    if (l == r)
        return !u ? l : mn[u];
    int mid = l + r >> 1;
    push_down(u, l, r, mid);
    if (ql <= mid)
        return Query(ls[u], l, mid, ql);
    return Query(rs[u], mid + 1, r, ql);
}

signed main()
{
    n = read();
    ll lasans = 0;
    int rt = 0;
    for(int i = 1; i <= n; ++i)
    {
        int x = read();
        int l = Findl(rt, 0, 1e9, x), r = Findr(rt, 0, 1e9, x);
        update(rt, 0, 1e9, 0, l, +1);
        update(rt, 0, 1e9, r, 1e9, -1);
        m = read();
        for(int j = 1; j <= m; ++j)
        {
            int x = (read() + lasans) % mod;
            printf("%lld\n", lasans = Query(rt, 0, 1e9, x));
        }
    }
}