CodeForce 538-F A Heap of Heaps(主席树)

题目: 传送门
思路:
       咱用 离线+树状数组 逃了好多次课 ，终究还是 被抓回来补课了；
       对于每个k，每个点我们都去查询 [k*(I-1)+2,min(n,k*I+1)] 区间中小于 a[I]的个数，我们写出这种区间的总数，用调和级数可知范围在 nlog n 以内，因为查询的区间太大了，而且离线查询我们还有一个排序的操作，所以对时间的损耗太大，使用主席树更优.
       刚刚学习主席树，感觉这个东西很奇妙，给我的感觉有点像动态开点,不过二者目的不同，前者是为了记录更多的历史版本，后者是为了节省更多的空间。

AC_Code

const int maxn = 2e5+10;

int n,len,m,cnt;
int a[maxn];
int v[maxn];
int tree[maxn*40+10];
int ls[maxn*40+10],rs[maxn*40+10];
int root[maxn];
int ans[maxn];

void modify(int l,int r,int& x,int y,int pos) {
    x = ++cnt;
    tree[x] = tree[y];
    ls[x] = ls[y];
    rs[x] = rs[y];
    tree[x]++;
    if(l==r) return ;
    int mid =(l+r) >>1;
    if(pos<=mid) modify(l,mid,ls[x],ls[y],pos);
    if(pos> mid) modify(mid+1,r,rs[x],rs[y],pos);
    return ;
}

void quriy(int l,int r,int x,int y,int k,int& ans) {
    if(k<l) return ;
    if(r <=k) {
        ans += tree[x]-tree[y];
        return ;
    }
    int mid = (l+r) >>1;
    quriy(l,mid,ls[x],ls[y],k,ans);
    if(mid<k) quriy(mid+1,r,rs[x],rs[y],k,ans);
    return ;
}

int main() {
    scanf("%d",&n);
    for(int i=1;i<=n;i++) {
        scanf("%d",&a[i]);
        v[i-1]=a[i];
    }
    sort(v,v+n);
    len = unique(v,v+n)-v;
    for(int i=1;i<=n;i++) modify(1,len,root[i],root[i-1],lower_bound(v,v+len,a[i])-v+1);

    for(int k=1;k<n;k++) {
        for(int i=1;i<=n;i++) {
            long long l = 1ll*k*(i-1)+2;
            if(l>n) break;
            long long r = min(n,k*i+1);
            //cout<<i<<' '<<l<<' '<<r<<endl;
            quriy(1,len,root[r],root[l-1],lower_bound(v,v+len,a[i])-v,ans[k]);
        }
    }
    for (int i = 1; i < n; ++i){
        printf("%d",ans[i]);
        if(i != n-1) printf(" ");
        else printf("\n");
    }
    return 0;
}