题目链接:点击打开链接
题意很好理解,实际上转换一下,b[i]代表的就是在a[i]前面的第一个倍数,如果没有就等于a[i],,c[i]就是a[i]后面的第一个倍数。
这样先预处理出1到100000所有数的约数存入p[ ],用k[ ]表示1到100000中每个数紧挨着的倍数。先从前往后扫扫到每个数就根据这个数的约数更新k[ ]得到b[ ],再从后往前得到c[ ]。
代码:
#include <cstdio>
#include <iostream>
#include <vector>
#define MAX 100010
#include <cstring>
using namespace std;
int a[MAX];
int b[MAX];
int c[MAX];
int k[MAX];
bool vis[MAX];
vector <int > p[MAX];
int main(){
int n;
for(int i=1;i<=100000;i++){
for(int j=i;j<=100000;j+=i){
p[j].push_back(i);
}
}
while(~scanf("%d",&n),n){
memset(vis,0,sizeof(vis));
for(int i=1;i<=n;i++){
scanf("%d",&a[i]);
}
for(int i=1;i<=n;i++){
if(!vis[a[i]]) b[i]=a[i];
else b[i]=k[a[i]];
for(int j=0;j<p[a[i]].size();j++){
k[p[a[i]][j]]=a[i];
vis[p[a[i]][j]]=1;
}
}
memset(vis,0,sizeof(vis));
for(int i=n;i>=1;i--){
if(!vis[a[i]]) c[i]=a[i];
else c[i]=k[a[i]];
for(int j=0;j<p[a[i]].size();j++){
k[p[a[i]][j]]=a[i];
vis[p[a[i]][j]]=1;
}
}
long long res=0;
for(int i=1;i<=n;i++){
res+=(long long)b[i]*(long long)c[i];
}
printf("%I64d\n",res);
}
return 0;
}