2015 Multi-University Training Contest 1（OO’s Sequence）

题意：求每一个[L, R]区间中   j!=i  a[i]对任意a[j]取余为0 的i的个数的总和

思路：讨论每个数，L表示i左边离i最近的i的因子的位置，R表示i右边离i最近的i的因子，则包括a[i]的有（i-L）*（R-i）个区间。

代码：

#include<set>
#include<map>
#include<cmath>
#include<queue>
#include<stack>
#include<math.h>
#include<string>
#include<vector>
#include<cstdio>
#include<time.h>
#include<cstring>
#include<iomanip>
#include<iostream>
#include<string.h>
#include<algorithm>
#define ll long long
#define PI acos(-1.0)
using namespace std;
const ll mod=1e9+7;
vector<ll>mp[10005];
ll pre[10005];
ll l[100005],r[100005];
ll n;
ll a[100005];
void init()
{
    for(int i=1;i<=10000;i++)
    {
        for(int j=1;j*j<=i;j++)
        {
            if(i%j==0)
            {
                mp[i].push_back(j);
                if(j!=(i/j))mp[i].push_back(i/j);
            }
        }
    }
}
int main()
{
    init();
    while(scanf("%lld",&n)!=EOF)
    {
        for(int i=1;i<=n;i++) scanf("%lld",&a[i]);
        memset(pre,0,sizeof(pre));
        for(int i=1;i<=n;i++)
        {
            ll maxn=0;
            for(int j=0;j<mp[a[i]].size();j++) maxn=max(maxn,pre[mp[a[i]][j]]);
            l[i]=maxn;
            pre[a[i]]=i;
        }
        memset(pre,0x3f3f3f,sizeof(pre));
        for(int i=n;i>=1;i--)
        {
            ll minn=n+1;
            for(int j=0;j<mp[a[i]].size();j++) minn=min(minn,pre[mp[a[i]][j]]);
            r[i]=minn;
            pre[a[i]]=i;
        }
        ll ans=0;
        for(int i=1;i<=n;i++)
        {
            ans=(ans+(i-l[i])*(r[i]-i)%mod)%mod;
        }
        cout<<ans<<endl;
    }
}