hdu5288 区间求个数和

对每个区间这种数的个数求和，这种数满足在这个区间内没有数能被它整除，不包括它本身。

#include <iostream>
#include<stdio.h>
#include<string.h>
#include<vector>
#include<algorithm>
#include<cmath>
#include<map>
#include<queue>
using namespace std;
const int N=100005;
const int MAXN=(1<<31)-1;
int INF=0x7f7f7f7f;
int T,n,m,k,tot;
const int MOD=1e9+7;
int cas=1;
vector<int> a[N];
int b[N];
int main()
{
#ifndef  ONLINE_JUDGE
   freopen("aaa","r",stdin);
#endif
    int T;
    while(~scanf("%d",&n)){
        for(int i =0;i<N;i++) a[i].clear();
        for(int i=1;i<=n;i++) {
            scanf("%d",&b[i]);
            a[b[i]].push_back(i);
        }

        long long res=0;
        for(int i=1;i<=n;i++){
             int l=0,r=n+1;
             for(int j=1;j*j<=b[i];j++){
                if(b[i]%j==0){
                  if(a[j].size()>0){
                        int rr=upper_bound(a[j].begin(),a[j].end(),i)-a[j].begin();//第一个大于i的数字
                        if(rr<a[j].size()) r=min(a[j][rr],r);
                        int ll=lower_bound(a[j].begin(),a[j].end(),i)-a[j].begin()-1;//第一个小于i的数字
                        if(ll>=0) l=max(a[j][ll],l);
                    }
                    int jj=b[i]/j;
                    if(a[jj].size()>0){
                        int rr=upper_bound(a[jj].begin(),a[jj].end(),i)-a[jj].begin();
                        if(rr<a[jj].size()) r=min(a[jj][rr],r);
                        int ll=lower_bound(a[jj].begin(),a[jj].end(),i)-a[jj].begin()-1;
                        if(ll>=0) l=max(a[jj][ll],l);
                    }
                }
             }
             res+=(long long)(i-l)*(r-i);
             res%=MOD;
        }
        printf("%d\n",res);
    }
    return 0;
}