Codeforces 785E. Anton and Permutation

题目链接：http://codeforces.com/problemset/problem/785/E

---

其实可以CDQ分治...

我们只要用一个数据结构支持单点修改，区间查询比一个数大(小)的数字有多少个就可以了。

考虑分块，每段区间之内有排序或者二分查询比一个数大的树的个数的操作。

复杂度${O(qn \sqrt n log_{2}^{\sqrt n})}$

 #include<iostream>
 #include<cstdio>
 #include<algorithm>
 #include<vector>
 #include<cstdlib>
 #include<cmath>
 #include<cstring>
 using namespace std;
 #define maxn 1000100
 #define llg int
 #define LL long long
 #define yyj(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
 llg n,m,SIZE,belong[maxn],cnt,ne[maxn],a[maxn],T;
 LL ans;
 vector<llg>c[maxn];

 inline int getint()
 {
     int w=,q=; char c=getchar();
     while((c<'' || c>'') && c!='-') c=getchar(); if(c=='-') q=,c=getchar();
     while (c>='' && c<='') w=w*+c-'', c=getchar(); return q ? -w : w;
 }

 void init()
 {
     SIZE=sqrt(n);
     for (llg i=;i<=n;i++)
     {
         ne[i]=i+SIZE-; ne[i]=min(ne[i],n); cnt++;
         c[cnt].push_back();
         for (llg j=i;j<=ne[i];j++)
         {
             ne[j]=ne[i];
             belong[j]=cnt;
             c[cnt].push_back(j);
         }
         i=ne[i];
     }
     for (llg i=;i<=n;i++) a[i]=i;
 }

 llg erfen(llg x,llg v)
 {
     llg l=,r=c[x].size()-,mid,wz;
     while (l<=r)
     {
         mid=(l+r)>>;
         if (c[x][mid]<=v) {l=mid+; wz=mid;}else r=mid-;
     }
     return c[x].size()-wz;
 }

 llg more(llg l,llg r,llg v)
 {
     if (r<l) return ;
     llg tot=;
     llg stk=belong[l],endk=belong[r];
     for (llg i=l;i<=r;i++)
     {
         if (belong[i]==stk || belong[i]==endk)
         {
             if (a[i]>v) tot++;
             continue;
         }
         tot+=erfen(belong[i],v);
         i=ne[i];
     }
     return tot;
 }

 void change(llg x,llg v1,llg v2)
 {
     llg w=c[x].size();
     for (llg i=;i<w;i++)
         if (c[x][i]==v1)
         {
             c[x][i]=v2;
             sort(c[x].begin(),c[x].end());
             return ;
         }
 }

 int main()
 {
 //    yyj("a");
     cin>>n>>T;
     init();
     while (T--)
     {
         llg l,r;
         l=getint(),r=getint();
         if (l>r) swap(l,r);
         if (l!=r)
         {
             llg morel=more(l+,r-,a[l]);
             llg lessl=r--l-morel;
             llg morer=more(l+,r-,a[r]);
             llg lessr=r--l-morer;
             ans+=morel-lessl+lessr-morer;
             if (a[l]>a[r]) ans--;else ans++;
             change(belong[l],a[l],a[r]); change(belong[r],a[r],a[l]);
             swap(a[l],a[r]);
         }
         printf("%I64d\n",ans);
     }
     return ;
 }

---