ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

2018-2019 ACM-ICPC, Asia Nanjing Regional Contest J. Prime Game(数论)

2021-09-28 19:59:51  阅读:181  来源: 互联网

标签:Prime Contest int ACM vec fac mul primes mp


题组传送门

题意

现在给你一个含有 n n n个正整数的集合 a [ ] a[ ] a[],现在定义一个函数 m u l ( l , r ) mul(l,r) mul(l,r)代表 a l ∗ a l + 1 ∗ . . . ∗ a r − 1 ∗ a r a_l*a_{l+1}*...*a_{r-1}*a_r al​∗al+1​∗...∗ar−1​∗ar​并定义 f a c ( l , r ) fac(l,r) fac(l,r)代表 m u l ( l , r ) mul(l,r) mul(l,r)中含有不同的素因子的个数。
请你计算出 f a c ( i , j ) fac(i,j) fac(i,j)

思路

这个题其实就是计算每一个素数的贡献,然后进行运算

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N = 1e6 + 7;
int primes[N], cnt; // primes[]存储所有素数
bool st[N];         // st[x]存储x是否被筛掉
int a[2102100];
map<int, int> mp;
void get_primes(int n)
{
    for (int i = 2; i <= n; i++)
    {
        if (!st[i])
            primes[cnt++] = i;
        for (int j = 0; primes[j] <= n / i; j++)
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0)
                break;
        }
    }
}
signed main()
{
    get_primes(N - 7);
    vector<vector<int> >vec;
    int n;
    cin >> n;
    for (int i = 1;i<=n;i++){
        cin >> a[i];
        int t = a[i];
        for (int j = 0; primes[j] <= t / primes[j];j++){
            int ju = 0;
            while(t%primes[j]==0){
                t /= primes[j];
                ju = 1;
            }
            if(ju){
                vec[i].push_back(primes[j]);
            }
        }
        if(t>1){
            vec[i].push_back(t);
        }
    }
    int ans = 0;
    for (int i = 1; i <= n;i++){
        for (int j = 0; j < vec[i].size();j++){
            if(!mp[vec[i][j]]){
                ans += (n - i + 1);
            }else{
                ans += (n - mp[vec[i][j]] + 1);

            }
            mp[vec[i][j]] = i;
        }
    }
    cout << ans << endl;
}

标签:Prime,Contest,int,ACM,vec,fac,mul,primes,mp
来源: https://blog.csdn.net/m0_51841071/article/details/120463524

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有