hdu 5288 数学 ****

给一个序列
定义函数f(l ,r) 为区间[l ,r] 中
的数ai不是在这个区间其他任意数aj的倍数
求所有f(l,r)之和

通过预处理，记录 a[i] 的左右边界（所谓的左右边界时 在从 a[i] 当前位置往左往右找，找到左边第一个和右边第一个能够整除 a[i] 的数，这两个数就是a[i]的左右边界）然后记录到 l[] & r[] 中， 这样 a[i] 对 ans 的贡献是 (i - l[i]) * (r[i] - i);在预处理 l[] 数组时，用pre[j]标记一下 j （表示 j 最后一次出现的位置），如果 j 在之前已经遇到过且右边界没有被更新过，就将 pre[j] 的边界更新到当前 i 所在的位置。

前面做标记，后面遇到时再将记录的值更新到对应的l，r里

 #include<cstdio>
 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cmath>
 #include<queue>
 #include<map>
 using namespace std;
 #define MOD 1000000007
 const int INF=0x3f3f3f3f;
 const double eps=1e-;
 typedef long long ll;
 #define cl(a) memset(a,0,sizeof(a))
 #define ts printf("*****\n");
 const int MAXN=;
 int n,m,tt;
 int a[MAXN];
 ll l[MAXN],r[MAXN],pre[MAXN];
 int main()
 {
     int i,j,k,ca=;
     #ifndef ONLINE_JUDGE
     freopen("1.in","r",stdin);
     #endif
     while(~scanf("%d",&n))
     {
         for(i=;i<=n;i++)
         {
             scanf("%d",a+i);
             l[i]=,r[i]=n+;
         }
         cl(pre);
         for(i=;i<=n;i++)
         {
             for(j=a[i];j<=;j+=a[i])
             {
                 if(pre[j]&&r[pre[j]]==n+)      //保证是第一个遇到的,因此更新过的就不能再更新了
                 {
                     r[pre[j]]=i;
                 }
             }
             pre[a[i]]=i;
         }
         cl(pre);
         for(i=n;i>;i--)
         {
             for(j=a[i];j<=;j+=a[i])
             {
                 if(pre[j]&&l[pre[j]]==)
                 {
                     l[pre[j]]=i;
                 }
             }
             pre[a[i]]=i;
         }
         ll ans=;
         for(i=;i<=n;i++)
         {
             ans+=(ll)(((i-l[i])*(r[i]-i))%MOD);
             ans=(ans%MOD+MOD)%MOD;
         }
         printf("%I64d\n",ans);
     }
 }