hdu5288OO’s Sequence

//给一个序列
//定义函数f(l ,r) 为区间[l ,r] 中
//的数ai不是在这个区间其它随意数aj的倍数
//求全部f(l,r)之和
//对于每个数a[i]找其最左的区间l[i]和最右的区间r[i]
//包括a[i]且使得a[i]满足条件的个数为(i-l[i]+1)*(r[i] - i + 1)
//对于每个r[i]仅仅要从左到右遍历
//pre[i] 表示i最后出现的位置
//对于一个数a[i]。枚举其全部倍数值
//在其左边找到其相应的位置，那么 r[pre[j]] = i - 1;
//同理能够得到l[i]
#include<cstdio>
#include<cstring>
#include<iostream>
using namespace std ;
const int maxn = 100010 ;
const int mod = 1e9+7 ;
int a[maxn];
typedef __int64 ll ;
ll last[maxn] ;
ll pre[maxn] ;
ll l[maxn] ;
ll r[maxn] ;
int  main()
{
    //freopen("in.txt" ,"r" ,stdin) ;
    int n ;
    while(~scanf("%d" ,&n))
    {
        memset(last , 0 ,sizeof(last)) ;
        memset(pre , 0 ,sizeof(pre)) ;
        for(ll i =1;i <= n;i++)
        {
            l[i] = 1;r[i] = n ;
            scanf("%d" ,&a[i]) ;
            for(ll j = a[i] ;j < maxn/10;j+=a[i])
            if(pre[j] != 0 && r[pre[j]] == n)
            r[pre[j]] = i - 1;
            pre[a[i]] = i;
        }
        for(ll i = n ;i >= 1;i--)
        {
            for(ll j = a[i];j < maxn/10;j+=a[i])
            if(last[j] != 0 && l[last[j]] == 1)
            l[last[j]] = i + 1;
            last[a[i]] = i ;
        }
        ll ans = 0 ;
        for(ll i = 1;i <= n;i++)
            ans = (ans + (((i-l[i]+1)%mod)*((r[i] - i + 1)%mod))%mod)%mod ;
        printf("%I64d\n" ,ans) ;
    }
    return 0 ;
}