Codeforces Round #182 (Div. 1) D. Yaroslav and Divisors

数学分析很重要啊~~

求l，r之间有多少个序列对：

道题给出的是1-N的一个排列，因此是一个非常特殊的问题，倍数关系的对只有NLOGN(N/1 + N/2 + ..+N/N)
预处理出所有的序列对，可以统计出前缀和sum[l]表示1~l之间的对数

对于每个查询(X,Y)，我们将所有关系对XI>=X的，YI <= Y进行统计，

可以用总数sum[y]-sum[x] - (xI<x但是yI>=x&&yI<=y的对数)

需要对所有对数排序，所有询问也排个序

原文大神链接：http://hi.baidu.com/zyz913614263/item/86338144aba62ab8df2a9f7f

自己写的代码 比较龊

#include<cstdio>
#include<cstring>
#include<climits>
#include<string>
#include<algorithm>
#include<cmath>
#include<iostream>
#include<queue>
#include<vector>
#include<stack>
#include<map>
#include<set>
using namespace std;
typedef long long ll;

template<class T> void show(T *a, int n) {  for(int i = 0; i < n; i++) cout<<a[i]<<"  "; puts("");  }
template<class T> void queue_clear(T *a, int n) {   for(int i = 0; i < n; i++) {while(!a->empty()) a->pop();}   }      
template<class T> void Bit_add(T *a, int x, int p, int MAX) {   while(p < MAX){a[p] += x; p += p&(-p);}   }
template<class T> T Bit_sum(T *a, int p) { T ret = 0; while(p){ret += a[p];p -= p&(-p);}   return ret; }

const int MAXN = 200007;
struct node {
    int l, r, p;
    bool operator < (const node &no) const {
        return l < no.l;
    }
}q[MAXN];
vector<int> pai[MAXN];
int sum[MAXN];
int d[MAXN], p[MAXN];
int del[MAXN];
int ans[MAXN];

int main() {
    int n, m;
    while(~scanf("%d%d", &n, &m)) {
        int i, j, t;
        int l, r;
        for(i = 1; i <= n; i++) {
            scanf("%d", &t);
            p[t] = i;
        }
        for(i = 1; i <= m; i++) {
            scanf("%d%d", &q[i].l, &q[i].r);
            q[i].p = i;
        }
        
        memset(sum, 0, sizeof(sum));
        for(i = 1; i <= n; i++) pai[i].clear();
        for(i = 1; i <= n; i++) {
            for(j = i; j <= n; j += i) {
                l = p[i], r = p[j];
                if(l > r) swap(l, r);
                pai[l].push_back(r);
                sum[r] ++;
            }
        }
        for(i = 1; i <= n; i++) sum[i] += sum[i-1];
        
        sort(q+1, q+1+m);
        memset(d, 0, sizeof(d));
        int nq = 1;
        for(i = 1; i <= n; i++) {
            for(j = 0; j < pai[i-1].size(); j++) {
                Bit_add(d, 1, pai[i-1][j], n+1);
            }
            while(nq <= m && q[nq].l == i) {
                del[q[nq].p] = Bit_sum(d, q[nq].r) - Bit_sum(d, q[nq].l-1);
                nq ++;
            }
        }
        
        for(i = 1; i <= m; i++) ans[q[i].p] = sum[q[i].r] - sum[q[i].l-1] - del[q[i].p];
        for(i = 1; i <= m; i++) printf("%d\n", ans[i]);
    }
    return 0;
}