第一个元素的素因子2:
它能贡献的区间有[1,1],[1,2],……,[1,10] 10个区间
第一个元素的素因子3:
它能贡献的区间有[1,1],[1,2],……,[1,10] 10个区间
当前sum = 10+10
第二个元素的素因子7:
它能贡献的区间有[1,2],[1,3],……,[1,10] 9个区间
它能贡献的区间有[2,2],[2,3],……,[2,10] 9个区间
当前sum = 10+10 +9*2
同理第三个元素的素因子5算好后
sum = 10+10+9*2+8*3
当考虑第四元素的素因子5时,发现i = 3时的素因子5
在区间:
[1,3],[1,4],……,[1,10]
[2,3],[2,4],……,[2,10]
[3,3],[3,4],……,[3,10]
这3*8个区间中已经贡献过,所以我们从当前i = 4位置向后考虑 这个位置的素因子5对区间的贡献为 7
sum = 10+10+9*2+8*3+7
最终到n = 10,sum = 10+10+9*2+8*3+7+6*4+5*5+4+0+2*4+1+3=134
这种计数题算贡献,怎么算贡献呢?直接算有多少区间包含某一个素数的倍数。
例如,记录一下所有2的倍数的位置,3的倍数的位置,之后算有多少个区间包含2倍数的位置,这个就是2这个素数的倍数的贡献,然后算有多少个区间包含3这个素数的倍数的位置,这个就是3这个素数的贡献。
例如样例二:
2这个素数的倍数的位置有:1 5 9 10
3:1 6 7 10
5:3 4
7:2
那么2贡献的区间就有:(5-1)*1 + (9-5)*5 + (10-9)*9 + (10+1-10)*10
可以理解为:包含1不包含5的区间有多少,包含1和5不包含9的区间有多少。依次类推。
AC代码:
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<vector>
#include<stdlib.h>
#include<queue>
#include<set>
#include<iomanip>
#include<math.h>
using namespace std;
typedef long long ll;
typedef double ld;
const ll mod=1e9+7;
int a[1000004];
vector<int>v[1000004];
int mp[1000004];
int i,j,k;
int main()
{
int n;
cin>>n;
for(i=1; i<=n; ++i)
{
scanf("%d",&a[i]);
for(j=2; j*j<=a[i]; ++j)
{
if(a[i]%j==0)
{
v[i].push_back(j);
while(a[i]%j==0)
a[i]/=j;
}
}
if(a[i]!=1)
v[i].push_back(a[i]);
}
for(i=1; i<=1000000; ++i)
{
mp[i]=n+1;
}
ll ans=0;
ll sum=0;
for(i=n; i>=1; --i)
{
for(j=0; j<v[i].size(); ++j)
{
sum+=mp[v[i][j]]-i;
mp[v[i][j]]=i;
}
ans+=sum;
}
cout<<ans<<endl;
}