OO’s Sequence

5
1 2 3 4 5

23

记录l[i]，r[i]，i的约数存在的位置，要求最接近i的。可以用数组记录下出现的最接近当前位置的数。得到l[i],r[i]后，那么i可以被计数的区间也就有了，左边界为l[i]到i，右边界为i到r[i]，那么i一共会被记录(i-l[i])*(r[i]-i)次，累加所有的计数。

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std ;
#define LL __int64
const int MOD = 1e9+7 ;
int a[100010] , l[100010] , r[100010] ;
int Map[10010] ;
int main() {
    int i , j , n ;
    LL ans ;
    while( scanf("%d", &n) != EOF ) {
        for(i = 1 ; i <= n ; i++)
            scanf("%d", &a[i]) ;
        memset(Map,0,sizeof(Map)) ;
        for(i = 1 ; i <= n ; i++) {
            for(j = 1 , l[i] = 0 ; j*j <= a[i] ; j++) {
                if( a[i]%j ) continue ;
                if( Map[j] ) l[i] = max(l[i],Map[j]) ;
                if( Map[a[i]/j] ) l[i] = max(l[i],Map[a[i]/j]) ;
            }
            Map[a[i]] = i ;
        }
        memset(Map,0,sizeof(Map)) ;
        for(i = n ; i > 0 ; i--) {
            for(j = 1 , r[i] = n+1 ; j*j <= a[i] ; j++) {
                if( a[i]%j ) continue ;
                if( Map[j] ) r[i] = min(r[i],Map[j]) ;
                if( Map[a[i]/j] ) r[i] = min(r[i],Map[a[i]/j]) ;
            }
            Map[a[i]] = i ;
        }
        for(i = 1 , ans = 0 ; i <= n ; i++) {
            ans = (ans + (i-l[i])*(r[i]-i) ) % MOD ;
        }
        printf("%I64d\n", ans) ;
    }
    return 0 ;
}

hdu5288(2015多校1)OO’s Sequence

原文：http://blog.csdn.net/winddreams/article/details/46991181