2018南京Gym - 101981J - Prime Game（计数）

第一个元素的素因子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;
}