CodeForces301D 【BIT】

转自：http://www.cnblogs.com/oyking/archive/2013/08/14/3256922.html
题意：
给你n个数，求区间有多少对除数和被除数的对数

思路：
线段树 区间和/ 树状数组 前缀和
2e5的数的除数最多才81个
对于一个数的贡献，也就是1 - Current_Position的贡献。

Sum[pos]: 1到POS的对数。
对于区间[L,R]，如果简单考虑 Sum[1, R] - Sum[1, L] 的话，对于[1,L-1]在[L,R]的有效值很难处理
有种方法：
对区间排序，一个把左端点按升序排序，一个右端点按升序排序，
树状数组维护当前位置1-POS的对数；
单点修改：对一个值，维护以他的倍数的位置。
每次遇到左端点，先减去Sum[1, R] - Sum[1, L]（ [1,L-1]在[L,R]的有效值）。（此时不存在[L,R]的有效对数，这个如果想到就简单了吧）
每次遇到右端点，加上 Sum[1, R] - Sum[1, L] （之前该减的已经减了）。

就是 区间端点 排序/从小到大处理

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;

typedef pair<LL,LL>PII;

const int N=2e5+10;
int temp[N],n;
struct asd{
    int Left,Right;
    int id;
};
asd node1[N],node2[N];
int pos[N];
int ans[N];
bool cmp1(asd x,asd y){
    if(x.Left == y.Left) return x.Right < y.Right;
    return x.Left < y.Left;
}
bool cmp2(asd x,asd y){
    if(x.Right == y.Right) return x.Left < y.Left;
    return x.Right < y.Right;
}

int C[N<<1];
int lowbit(int x) {
    return x & (-x);
}
void add(int x, int d) {
    while(x <= n) {
        C[x] += d;
        x += lowbit(x);
    }
}
int Sum(int x) {
    int s = 0;
    while(x > 0) {
        s += C[x];
        x -= lowbit(x);
    }
    return s;
}

void update(int x){
    for(int i=x;i<=n;i+=x)
        add(pos[i],1);
}

int main(){
    int m;
    scanf("%d%d",&n,&m);
    memset(C,0,sizeof(C));

    for(int i=1;i<=n;i++){
        scanf("%d",&temp[i]);
        pos[temp[i]] = i;
    }

    for(int i=1;i<=m;i++){
        scanf("%d%d",&node1[i].Left,&node1[i].Right);
        node2[i].Left = node1[i].Left;
        node2[i].Right = node1[i].Right;
        node2[i].id = node1[i].id = i;
    }
    sort(node1+1,node1+m+1,cmp1);
    sort(node2+1,node2+m+1,cmp2);

    int j,k;
    j = k =1;
    memset(ans,0,sizeof(ans));
    for(int i=1;i<=n;i++){
        while(j <= m && node1[j].Left == i){
            ans[node1[j].id] -= (Sum(node1[j].Right) - Sum(node1[j].Left - 1));
            j++;
        }
        update(temp[i]);
        while(k <= m && node2[k].Right == i){
            ans[node2[k].id] += (Sum(node2[k].Right) - Sum(node2[k].Left - 1));
            k++;
        }
    }
    for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
    return 0;
}