2021ICPC江西省赛G.Magic Number Group莫队

博客介绍了如何通过预处理每个数的质因子来解决区间[L,R]内最多有多少数能被质数p整除的问题。采用莫队算法离线维护众数板子,对区间进行操作,实现高效求解。代码中展示了预处理、维护区间质因子个数以及解答查询的详细过程。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

题目

G.Magic Number Group

给定一个正整数序列,每次询问区间 [ L , R ] [L,R] [L,R],任意选择一个大于1的正整数p,该区间内最多有多少数能被p整除。

解题思路

对于p,我们一定是选择一个质数,问题就转化成了区间 [ L , R ] [L,R] [L,R],要使尽可能多的数包含质因子p。

考虑对于所有数都分解出所有质因子。每个数的质因子数量很少,最后就是莫队离线维护众数板子。

代码

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

const int M = 1e6 + 6;
vector<int> fac[M];
bool vis[M];
void init() //预处理每个数的质因子
{
    for (int i = 2; i < M; i++)
    {
        if (!vis[i])
        {
            for (int j = i; j < M; j += i)
            {
                vis[j] = 1;
                fac[j].emplace_back(i);
            }
        }
    }
}

const int N = 5e4 + 5;
vector<int> a(N);
int block;
struct node
{
    int l, r, id;
    bool operator<(const node &o) const
    {
        if (l / block != o.l / block)
            return l < o.l;
        if (l / block & 1)
            return r > o.r;
        return r < o.r;
    }
} q[N];
int cnt[M];
int num[M];
vector<int> ans(N);
int maxx;
void del(int x)
{
    for (auto p : fac[a[x]])
    {
        num[cnt[p]]--;
        if (maxx == cnt[p] && num[cnt[p]] == 0)
            maxx--;
        cnt[p]--;
        num[cnt[p]]++;
    }
}
void add(int x)
{
    for (auto p : fac[a[x]])
    {
        num[cnt[p]]--;
        cnt[p]++;
        num[cnt[p]]++;
        maxx = max(maxx, cnt[p]);
    }
}

int main()
{
    init();
    int t;
    scanf("%d", &t);
    while (t--)
    {
        int n, m;
        maxx = 0;
        scanf("%d %d", &n, &m);
        for (int i = 1; i <= n; ++i)
            scanf("%d", &a[i]);
        block = sqrt(n);
        for (int i = 1; i <= m; ++i)
        {
            scanf("%d %d", &q[i].l, &q[i].r);
            q[i].id = i;
        }
        sort(q + 1, q + 1 + m);
        int l = 1, r = 0;
        for (int i = 1; i <= m; ++i)
        {
            while (l < q[i].l)
                del(l++);
            while (l > q[i].l)
                add(--l);
            while (r < q[i].r)
                add(++r);
            while (r > q[i].r)
                del(r--);
            ans[q[i].id] = maxx;
        }
        for (int i = 1; i <= m; ++i)
            printf("%d\n", ans[i]);
        while (l <= r) //清空
            del(l++);
    }
    return 0;
}

### 关于2023年ICPC江西省题目解析 目前针对2023年ICPC江西省的比,具体的官方题解尚未广泛传播或公开发布。然而,在准备此类竞时,通常可以通过多种途径获取帮助和资源: - **在线评测平台(OJ)**:已有50多道来自各大公司的最新真题被录入至特定的OJ平台上供免费练习[^1]。虽然这些题目并非直接来源于2023年的ICPC江西省事,但对于熟悉相似类型的编程挑战非常有益。 对于希望深入理解往年类似事的情况,可参考如下实例: - 2022年江西理工大学软件工程学院举办的程序设计竞提供了详细的C++题解案例[^2]。 另外,有关区间选择优化问题的一个例子可以在2020年ICPC江西省大学生程序设计竞中找到,其中涉及到了如何挑选若干区间以实现`min(tot,x)`的最大化策略[^3]。 尽管上述资料不是专门针对2023年度的具体情况,但它们能够提供有价值的背景信息和技术指导,有助于参者更好地应对可能遇到的问题类型。 ```cpp // 示例代码展示了一个简单的字符串处理逻辑,这可能是解决某些ACM/ICPC风格问题的基础部分之一 #include <iostream> using namespace std; int main() { string s1 = "acmer ICPCAK"; string s2 = "CCPCAK"; int length = min(s1.size(), s2.size()); string result; for (int i = 0; i < length; ++i){ if (s1[i] == s2[i]){ result += s1[i]; } else break; } cout << result << endl; } ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

hesorchen

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值