ABC324题解

发布时间 2023-10-15 10:47:28作者: adam01

A/B

赛时没打。

C

暴力判断是相等 s[i] == t 还是替换了一个字符,或者是添加/删除了一个字符。

最后两个判断只需要交换一下 \(s\)\(t\) 的顺序就可以共用一个函数了。

D

注意到 \(N\le 13\),所以平方数不会超过 \(v=10^{13}\),很容易想到暴力枚举 \(\sqrt v\) 以内的数,判断是否能够被表示出来。

时间复杂度 \(O(N\sqrt v)\),或者说 \(O(N10^{\frac{N}{2}})\)

E

注意到如果两个字符串 \(S1,S2\) 拼起来可以包含 \(T\) 那么一定可以被前后分割成 \(T=T1+T2\) 使得 \(T1\subset S1, T2\subset S2\)

预处理出每个 \(S_i\) 包含 \(T\) 的最长前后缀的长度,使用双指针解决。

复杂度 \(O(\sum_i |S_i|)\)

查询时计算前缀后缀和大于 \(|T|\)\((i,j)\) 对数。

最后对计算出的前后缀长度做一个后缀和来优化查询。

时间复杂度 \(O(N+\sum_i |S_i|)\)

F

分数规划裸题。

二分答案,然后在 DAG 上跑个最长路 dp 就完事了。

> 什么是分数规划?<

G

将排列上的数放在二维平面上,每个数可以表示为 \((i,a_i)\) 的形式,对于每个序列可以写作 \((l,r,u,d)\),表示包含了 \(l\le i\le r, u\le a_i\le d\) 的所有数。注意到这相当于包含了一个矩形内的所有点,可以使用主席树来维护。

具体的,把点按照第二维排序,依次加入主席树中。

对于两个操作,相当于在矩形 \(s\) 里划一条平行于坐标轴的分界线(自己模拟一下)。

怎么找这条分界线呢?

对于每个 1 操作,在主席树上找到 \(l_s-1\) 的排名 \(rk\),然后主席树上二分出 \(rk+x\) 的值 \(v\) 就得到了分界线。

\((l,r,u,d)\) 变成了 \((l,v,u,d)(v+1,r,u,d)\)

对于每个 2 操作,分界线就是 \(x\),直接分割就行了。

\((l,r,u,d)\) 变成了 \((l,r,u,v)(l,r,v+1,d)\)

查询是容易的。

附 C,D,E,F,G 代码

C

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

const int N = 5e5 + 5;
string s[N], t;
int n;

bool chk1(string s1, string s2)
{
    int cnt = 0;
    for(int i = 0; i < s1.size(); i ++)
        cnt += s1[i] != s2[i];
    return cnt <= 1;
}

bool chk2(string s1, string s2)
{
    int dt = 0;
    for(int i = 0; i < s1.size(); i ++)
    {
        while(!dt && s1[i] != s2[i + dt]) dt ++;
        if(s1[i] != s2[i + dt]) return false;
    }
    return true;
}

signed main()
{
    ios::sync_with_stdio(0);cin.tie(0);
    cin >> n >> t;
    vector<int> ans;
    for(int i = 1; i <= n; i ++)
    {
        cin >> s[i];
        if(s[i].size() < t.size() - 1 || s[i].size() > t.size() + 1) continue;
        if(s[i] == t)
        {
            ans.push_back(i);
            continue;
        }
        if(s[i].size() == t.size() && chk1(s[i], t))
        {
            ans.push_back(i);
            continue;
        }
        if(s[i].size() == t.size() - 1 && chk2(s[i], t))
        {
            ans.push_back(i);
            continue;
        }
        if(s[i].size() == t.size() + 1 && chk2(t, s[i]))
        {
            ans.push_back(i);
            continue;
        }
    }
    cout << ans.size() << "\n";
    for(int i : ans) cout << i << " ";

    return 0;
}

D

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

#define int ll

const int N = 15;
int a[N], n, p[N];
string s;

inline bool chk(int x)
{
    int b[N];
    memcpy(b, a, sizeof a);
    while(x)
    {
        b[x % 10] --;
        if(b[x % 10] < 0) return false;
        x /= 10;
    }
    for(int i = 1; i < 10; i ++)
        if(b[i] > 0) return false;
    return true;
}

signed main()
{
    ios::sync_with_stdio(0);cin.tie(0);
    cin >> n >> s;
    p[0] = 1;
    for(int i = 1; i <= n; i ++) p[i] = p[i - 1] * 10ll;
    for(int i = 0; i < n; i ++) a[s[i] - '0'] ++;
    int ans = 0;
    for(int i = 0; i * i < p[n]; i ++)
        ans += chk(i * i);
    cout << ans;

    return 0;
}

E

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;

#define int ll

const int N = 5e5 + 5;
string s[N], t;
int n, l[N], r[N], sl[N], sr[N];

signed main()
{
    ios::sync_with_stdio(0);cin.tie(0);
    cin >> n >> t;
    for(int i = 1; i <= n; i ++)
    {
        cin >> s[i];
        int k = 0;
        for(int j = 0; k < t.size() && j < s[i].size(); j ++)
            if(s[i][j] == t[k]) k ++;
        l[i] = k;
        k = t.size() - 1;
        for(int j = s[i].size() - 1; k >= 0 && j >= 0; j --)
            if(s[i][j] == t[k]) k --;
        r[i] = t.size() - 1 - k;
        sl[l[i]] ++;
        sr[r[i]] ++;
    }
    for(int i = N - 5; i >= 0; i --) sl[i] += sl[i + 1];
    for(int i = N - 5; i >= 0; i --) sr[i] += sr[i + 1];
    int ans = 0;
    for(int i = 1; i <= n; i ++)
    {
        ans += sr[t.size() - l[i]] + sl[t.size() - r[i]];
    }
    cout << ans / 2;
    

    return 0;
}

F

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double D;

#define int ll

const int N = 2e5 + 5;

struct edge{int v, a, b;};
vector<edge> e[N];
int deg[N], n, m;
D f[N];

bool check(D x)
{
    queue<int> q;
    int d[N];
    memcpy(d, deg, sizeof deg);
    memset(f, -0x10, sizeof f);
    f[1] = 0;
    for(int i = 1; i <= n; i ++)
        if(!d[i]) q.push(i);
    while(!q.empty())
    {
        int t = q.front();q.pop();
        for(auto i : e[t])
        {
            if(--d[i.v] == 0) q.push(i.v);
            f[i.v] = max(f[i.v], f[t] + (i.a - x * i.b));
        }
    }
    return f[n] >= 0.L;
}

signed main()
{
    ios::sync_with_stdio(0);cin.tie(0);
    cin >> n >> m;
    for(int i = 1; i <= m; i ++)
    {
        int x, y, a, b;
        cin >> x >> y >> a >> b;
        e[x].push_back({y, a, b});
        deg[y] ++;
    }

    D l = 0, r = 1e4;
    while(l + (D)1e-12 < r)
    {
        D mid = (l + r) / 2.L;
        if(check(mid)) l = mid;
        else r = mid;
    }
    printf("%.17Lf", r);

    return 0;
}

G

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 2e5 + 5;
struct node{int l, r, sz;};
int n, m, a[N];
struct sgtv
{
    node tr[N << 5];
    int root[N], idx;
    int insert(int lstt, int l, int r, int x)
    {
        int nowt = ++idx;
        tr[nowt] = tr[lstt];
        if(l == r)
        {
            tr[nowt].sz++;
            return nowt;
        }
        int mid = (l + r) >> 1;

        if(mid >= x) tr[nowt].l = insert(tr[lstt].l, l, mid, x);
        else tr[nowt].r = insert(tr[lstt].r, mid + 1, r, x);
        tr[nowt].sz = tr[tr[nowt].l].sz + tr[tr[nowt].r].sz;
        return nowt;
    }
    int query(int lt, int rt, int l, int r, int k)
    {
        if(l > k || r < k) return 0;
        if(l == r) return tr[rt].sz - tr[lt].sz;
        int mid = (l + r) >> 1;
        if(mid >= k) return query(tr[lt].l, tr[rt].l, l, mid, k);
        else return tr[tr[rt].l].sz - tr[tr[lt].l].sz + query(tr[lt].r, tr[rt].r, mid + 1, r, k);
    }
    int find(int lt, int rt, int l, int r, int k)
    {
        if(l == r) return r;
        int mid = (l + r) >> 1;
        int cnt = tr[tr[rt].l].sz - tr[tr[lt].l].sz;
        if(cnt >= k) return find(tr[lt].l, tr[rt].l, l, mid, k);
        else return find(tr[lt].r, tr[rt].r, mid + 1, r, k - cnt);
    }
} tr;

struct M{int l, r, u, d;}b[N];
int h = 0;

int main()
{
    ios::sync_with_stdio(0);cin.tie(0);
    cin >> n;
    vector<pair<int, int>> v;
    for (int i = 1; i <= n; i++)
    {
        cin >> a[i];
        v.push_back({a[i], i});
    }

    sort(v.begin(), v.end());

    for(int i = 1; i <= n; i ++)
        tr.root[i] = tr.insert(tr.root[i - 1], 0, n + 1, v[i - 1].second);

    cin >> m;
    b[0] = {1, n, 1, n};
    while(m --)
    {
        int t, s, x;
        cin >> t >> s >> x;
        if(t == 1)
        {
            int L = tr.query(tr.root[b[s].u - 1], tr.root[b[s].d], 0, n + 1, b[s].l - 1);
            int r = tr.find(tr.root[b[s].u - 1], tr.root[b[s].d], 0, n + 1, L + x);
            r = min(r, b[s].r);
            r = max(r, b[s].l - 1);
            b[++h] = {r + 1, b[s].r, b[s].u, b[s].d};
            b[s].r = r;
        }
        else
        {
            x = min(x, b[s].d);
            x = max(x, b[s].u - 1);
            b[++h] = {b[s].l, b[s].r, x + 1, b[s].d};
            b[s].d = x;   
        }
        cout << tr.query(tr.root[b[h].u - 1], tr.root[b[h].d], 0, n + 1, b[h].r) - 
        tr.query(tr.root[b[h].u - 1], tr.root[b[h].d], 0, n + 1, b[h].l - 1) << "\n";
    }

    return 0;
}