【数学】【求每个值对结果的贡献】Prime Game

题意：

　　　定义了mul(l,r) 为a[l]*a[l+1]*……*a[r]， fac(l,r) 为mul(l,r)中素因子的个数

           求  ，(其中i从1-n,j从i-n)

---

题解:

       首先要考虑的是每个素因子能够贡献多少个区间

　　比如说对第p个位置的数a[p]的素因子，它所贡献的区间共有p*(n-p+1)　　(应该好理解吧。。）

　　但是如果第p个之前，比如说第q个存在了一个一样的素因子，那么就会有重复计算的区间，重复计算的区间的个数为q*(n-p+1)

       所以总的因子贡献为(p-q)*(n-p+1)

       

       所以可以线性筛出1e6内的素数，然后分解质因数即可

       其中用了vector数组存储，vec[i][j]=p表示，对素因子 i 的vector数组第 j 个数的在原数组的位置为 p　,（ 原数组是输入的数组）

---

Code :

//#include<bits/stdc++.h>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<stdio.h>
#include<string.h>
#include<vector>
#include<map>
#include<set>
#include<stack>
#include<unordered_map>
#pragma GCC optimize(2)
#define debug cout<<"*"<<endl;
#define in(x) scanf("%d",&x)
#define out(x) printf("%d\n",x)
#define llin(x) scanf("%lld",&x)
#define llout(x) printf("%lld\n",x)
using namespace std;
typedef long long ll;
typedef __int128 Biginteger;
typedef unsigned long long ull;
typedef pair<int,int> pii;
typedef pair<ll,ll> pll;
typedef pair<string,int> psi;
typedef pair<char,char> pcc;
const ll mod=14937;
const ll INF=0x3f3f3f3f3f;
const int maxn=1e6+17;
const int base=233;
const double PI=acos(-1);
const double euler_cons=0.57721566490153286060651209;
const int N=4;



int a[maxn];
vector<int> vec[maxn];    
int prime[maxn];
int cot=0;
bool vis[maxn];

void pre()    //埃氏筛筛出素数
{

    for(int i=2;i<=maxn;++i)
    {
        if(!vis[i]) prime[++cot]=i;
        for(int j=1;j<=cot&&i*prime[j]<=maxn;++j)
        {
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)   break;
        }
    }
}

void dec(int pos)    //分解质因子
{
    int m=a[pos];
    int i;
    for(i=1;prime[i]*prime[i]<=m&&i<=cot;i++)
    {
        if(m%prime[i]==0)
        {
            vec[prime[i]].push_back(pos);
            while(m%prime[i]==0)    //每个质因子只用记录一次
            {
                m/=prime[i];
            }
        }
    }
    if(m>1)
        vec[m].push_back(pos);
}
int main()
{
    pre();
    for(int i=1;i<=cot;i++)    //防止越界
    {
        vec[prime[i]].push_back(0);
    }
    int n;
    in(n);
    for(int i=1;i<=n;i++)
    {
        in(a[i]);
        dec(i);
    }

    ll ans=0;
    for(int i=1;i<=cot;i++)
    {
        int Size=vec[prime[i]].size();
        for(int j=1;j<Size;j++)
        {
          　ans+=(ll)(n-vec[prime[i]][j]+1)*(vec[prime[i]][j]-vec[prime[i]][j-1]);
        }
    }
    llout(ans);
    return 0;
}