hdu 5288 OO’s Sequence（2015多校第一场第1题）枚举因子

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=5288

 

题意：在闭区间[l,r]内有一个数a[i]，a[i]不能整除 除去自身以外的其他的数，f(l,r)表示在这区间内a[i]这样的数的个数，，现给你n个数，求所有区间的f(l,r)的和。

 

思路：对于每个数a[i]求出他的左右侧最靠近他的且是他的因子的位置L、R，并记录，那么对于每个数a[i]都有了他的L，R，而对于每个a[i]在f(l,r)有价值的次数之和就是(i-L+1)*(R-i+1)

 

代码：

 

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <iostream>
#include <queue>
#include <algorithm>
#include <vector>
using namespace std;
#define LL __int64
#define INF 0x3f3f3f3f
const int MAXN=100005;
#define mod 1000000007
int l[MAXN];
int r[MAXN];
int a[MAXN];
int visl[MAXN];
int visr[MAXN];
LL ans;
int main()
{
    int n,i,j;
    while(~scanf("%d",&n))
    {
        memset(visl,0,sizeof(visl));
        memset(visr,0,sizeof(visr));
        for(i=1;i<=n;i++)
        {
            scanf("%d",&a[i]);
            l[i]=1;
            r[i]=n;
        }
        for(i=1;i<=n;i++)
        {
            for(j=a[i];j<=10000;j+=a[i])
            {
                if(visr[j]&&r[visr[j]]==n)
                    r[visr[j]]=i-1;
            }
            visr[a[i]]=i;
        }
        for(i=n;i>0;i--)
        {
            for(j=a[i];j<=10000;j+=a[i])
            {
                if(visl[j]&&l[visl[j]]==1)
                    l[visl[j]]=i+1;
            }
            visl[a[i]]=i;
        }
        ans=0;
        for(i=1;i<=n;i++)
        {
            //printf("%d %d\n",l[i],r[i]);
            ans+=(LL)(i+1-l[i])*(r[i]-i+1);
            ans%=mod;
        }
        printf("%I64d\n",ans%mod);
    }
    return 0;
}