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 ;

}