hdu 5288 OO’s Sequence

题面求的是 对于每一组l r 存在的i的总和，显然很不好求，可以换一个角度，转换成每个数作为i时 可以选择的区间的数量，只要先预处理出每个说左右最近的可以被它整除的数的位置，左右长度乘一起就是这个数对应的区间数量。预处理左侧区间长度时可以从中向右扫，对每个数在遍历它的因子，看是否出现过，如果出现过，进行取优就可以了，右侧同理。

#include <iostream>
#include <cstring>
#include <string>
#include <cstdio>
#include <cmath>
#include <algorithm>
#include <vector>
#include <queue>
#include <map>
#include<set>
#include<stack>
#define inf 0x3f3f3f3f
#define LL long long
#define bug puts("bugbugbug");
using namespace std;
const int mod=1000000007;
int a[100005];
LL visit[100005];
LL l[100005];
LL r[100005];
int main()
{
    int n;
    while(~scanf("%d",&n))
    {
        memset(visit,0,sizeof(visit));
        for(int i=1; i<=n; i++)
        {
            l[i]=0;
            r[i]=n+1;
            scanf("%d",&a[i]);
            int ll=sqrt(a[i]+1);
            for(int j=1; j<=ll; j++)
                if(a[i]%j==0){
                    l[i]=max(l[i],visit[j]);
                    l[i]=max(l[i],visit[a[i]/j]);
                }
            visit[a[i]]=i;
        }
        for(int i=0;i<=100004;i++)
            visit[i]=n+1;
        for(int i=n; i>0; i--)
        {
            int ll=sqrt(a[i]+1);
            for(int j=1; j<=ll; j++)
                if(a[i]%j==0){
                    r[i]=min(r[i],visit[j]);
                     r[i]=min(r[i],visit[a[i]/j]);
                }
           visit[a[i]]=i;
        }
        LL ans=0;
        for(LL i=1;i<=n;i++)
        {
            ans+=(i-l[i])*(r[i]-i)%mod;
           ans%=mod;
        }
        printf("%I64d\n",ans);

    }
}