hdu 4961 Boring Sum n*logn

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
using namespace std;

int a[100002],b[100002],c[100002];
int vis[100002];

int main()
{
int n;
while(~scanf("%d",&n)&&n)
{
for(int i=1;i<=n;i++)
scanf("%d",&a[i]);
memset(vis,0,sizeof(vis));  //值i倍数且是在左边的中对应的最大下标
for(int i=1;i<=n;i++)
{
if(vis[a[i]])
b[i]=a[vis[a[i]]];
else b[i]=a[i];
for(int j=1;j<(int)sqrt((double)a[i])+1;j++)
{
if(a[i]%j==0)
{
vis[a[i]/j]=i;
vis[j]=i;
}
}//更新的都是i右边的约数，因为i和i前面的都已经得出来了 。由于i是小到大，所有得到的值也是最大的。
}

memset(vis,0,sizeof(vis)); //值i倍数且是在右边的对应的最小下标
for(int i=n;i>=1;i--)
{
if(vis[a[i]])
c[i]=a[vis[a[i]]];
else c[i]=a[i];

for(int j=1;j<(int)sqrt((double)a[i])+1;j++)
{
if(a[i]%j==0)
{
vis[a[i]/j]=i;
vis[j]=i;
}
}//更新的都是i左边的约数的vis值，因为i和i后面的都已经得出来了。由于i是小到大，所以得到的值也是最小的
}
__int64 ans=0;
for(int i=1;i<=n;i++)
ans+=(__int64 )b[i]*c[i];
printf("%I64d\n",ans);
}
return 0;
}