hdu 4961 Boring Sum

http://acm.hdu.edu.cn/showproblem.php?pid=4961

先打个表，把每个数的约数存起来，然后从前往后扫一遍，结果存在f[i],然后从后往前扫一遍，结果存在c[i],最后算f[i]*c[i]的和。

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#define maxn 100001
using namespace std;

int f[maxn],c[maxn];
vector<int>g[maxn];
int vis1[maxn],vis2[maxn];
int a[maxn];
void inti()
{
    for(int i=2; i<=maxn; i++)
    {
        for(int j=i; j<=maxn; j+=i)
        {
            g[j].push_back(i);
        }
    }
}
int n;
int main()
{
    inti();
    while(scanf("%d",&n)!=EOF)
    {
        memset(vis1,0,sizeof(vis1));
        memset(vis2,0,sizeof(vis2));
        memset(f,0,sizeof(f));
        memset(c,0,sizeof(c));
        memset(a,0,sizeof(a));
        if(n==0) break;
        for(int i=1; i<=n; i++)
        {
            scanf("%d",&a[i]);
        }
        for(int i=1; i<=n; i++)
        {
            if(vis1[a[i]]==0) f[i]=a[i];
            else f[i]=vis1[a[i]];
            vis1[1]=a[i];
            for(int j=0; j<(int)g[a[i]].size(); j++)
            {
                vis1[g[a[i]][j]]=a[i];
            }
        }
        for(int i=n; i>=1; i--)
        {
            if(vis2[a[i]]==0) c[i]=a[i];
            else c[i]=vis2[a[i]];
            vis2[1]=a[i];
            for(int j=0; j<(int)g[a[i]].size(); j++)
            {
                vis2[g[a[i]][j]]=a[i];
            }
        }
        __int64 sum=0;
        for(int i=1; i<=n; i++)
        {
            sum+=(__int64)f[i]*(__int64)c[i];
        }
        printf("%I64d
",sum);
    }
    return 0;
}