题意
现在给你一个含有
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;
}