10.24小练习

题目大意

给定一个长度为n的序列,对于序列中的每个正整数数$a_i$以及n满足:($1\le a_i,n\le 10^6$).
求:${\Sigma }_{i=1}^n{\Sigma }_{j=i}^{n}prime(a_i\ast a_{i+1}\ast ...\ast a_{j-1}\ast a_j)$.
其中$prime(i)$为$i$含有的不同质因数的个数.

输入样例
10
99 62 10 47 53 9 83 33 15 24
输出样例
248

解题思路

先预处理出$prime[i]$(即数字i的最小质因数).
通过找规律可以得出:
$value[i]={\Sigma }_j^{num}=1(n-i+1)\ast (i-pre[t_j])$
($value[i]$为第i个点的总贡献,$t_j$为第i个点的权值的第j小的质因数,$num$为其不同质因数个数，$pre[i]$为i最后出现的位置)
然后维护$pre[i]$更新$ans$即可.

参考代码

#include<bits/stdc++.h>
#define N 1001005
#define ll long long
using namespace std;

int n,prime[N];
int pre[N],x;
ll ans=0;

inline void Init()
{
	int maxx=1000005;
	for(int i=2;i<=maxx;i++){
		if(!prime[i])
		for(int j=i;j<=maxx;j+=i){
			if(!prime[j])
			prime[j]=i;
		}
	}
}
int main()
{
	n=in;
	Init();
	for(int i=1;i<=n;i++)
	{
		x=in;
		while(x>1)
		{
			int t=prime[x];
			while(x%t==0)x/=t;
			ans=ans+(ll)(n-i+1)*(i-pre[t]);//核心部分
			pre[t]=i;
		}
	}
	printf("%lld\n",ans);
	return 0;
}

~咕