HDU 4961 Boring Sum 因数分解

题意：给出序列，对于每个数ai，求出分别求出他的左边和右边离他最近的且是他倍数的bi,ci，如果没有令bi = ai,或ci = ai. 求出所有bi * ai的和。

思路：对于每个数，他只能成为他的因子的倍数。所以，我们从前往后更新时，对于每个数，将它的所有的因子的坐标跟新为它的坐标，查询离它最近的倍数，也就是查询它作为因子，在它左面的数中，是否有数成为它的倍数，即：在记录数组中，对应的位置是否有合法的坐标。

代码如下：

#include<cstdio>
#include<cstring>
#include<vector>
#include<cmath>
using namespace std;
const int maxn=100005;
int a[maxn];
int b[maxn];
int c[maxn];
int num[maxn];
int n;
vector<int> f[maxn];
void init()
{
int i,j;
for(i=1;i<maxn;i++)
for(j=1;j<=sqrt((double)i);j++)
{
if(i%j==0)
{
f[i].push_back(j);
f[i].push_back(i/j);
}

}
}
int main()
{
init();
while(scanf("%d",&n)&&n)
{
int i,j;
memset(num,0,sizeof(num));
for(i=1;i<=n;i++)
{
scanf("%d",&a[i]);
}
for(i=1;i<=n;i++)
{
int sz=f[a[i]].size();
if(num[a[i]]==0) b[i]=a[i];
else b[i]=a[num[a[i]]];
for(j=0;j<sz;j++)
{
num[f[a[i]][j]]=i;
}
}
memset(num,0,sizeof(num));
for(i=n;i>0;i--)
{
int sz=f[a[i]].size();
if(num[a[i]]==0) c[i]=a[i];
else c[i]=a[num[a[i]]];
for(j=0;j<sz;j++)
{
num[f[a[i]][j]]=i;
}
}
long long ans=0;
for(i=1;i<=n;i++)
{
ans=(long long)ans+((long long)b[i]*(long long)c[i]);
}
printf("%I64d\n",ans);
}
return 0;
}