HDU 5211 筛法求约数

给出n个数a1,a2...an，定义函数 f[i]=j，(i<j)，表示aj mod ai=0 的最小j，其中j大于i，如果不存在这样的数，则f[i]=0

求n个数所有f[]值的和

先用筛法o(nlogn)求出每个数的约数

然后每读入一个数x，先找出所有的约数，再看看之前有没有出现过这些约数。

这个回看的过程可以用一个数组维护。

维护watch[]，watch[aj]=j 表示aj还没有找到函数值，如果aj是x的约数，那么说明aj的函数值为x的位置。

#include<cstdio>
#include<cstring>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<set>
#include<map>
#include<stack>
#include<vector>
#include<queue>
#include<string>
#include<sstream>
#define eps 1e-9
#define ALL(x) x.begin(),x.end()
#define INS(x) inserter(x,x.begin())
#define FOR(i,j,k) for(int i=j;i<=k;i++)
#define MAXN 1005
#define MAXM 40005
#define INF 0x3fffffff
#define PB push_back
#define MP make_pair
#define X first
#define Y second
#define lc (k<<1)
#define rc ((k<<1)1)
using namespace std;
typedef long long LL;
int i,j,k,n,m,x,y,T,ans[],big,cas,num,len;
bool flag;
int sum,wat[];

vector <int> fac[];
void GetFact(int size)
{
    for (int i=;i<=size;i++)
    {
        for (int j=i;j<=size;j+=i)
        {
            fac[j].push_back(i);
        }
    }
}

int main()
{
    GetFact();
    while (~scanf("%d",&n))
    {
        memset(wat,,sizeof(wat));
        memset(ans,,sizeof(ans));
        for (i=;i<=n;i++)
        {
            scanf("%d",&x);
            for (j=;j<fac[x].size();j++)
            {
                int u=fac[x][j];
                if (wat[u])
                {
                    k=wat[u];
                    ans[k]=i;
                    wat[u]=;
                }
            }
            wat[x]=i;
        }
        sum=;
        for (i=;i<=n;i++)
        {
            sum+=ans[i];
        }
        printf("%d\n",sum);
    }

    return ;
}