[JZOJ3462] 休息

　　休息的时候，可以放松放松浑身的肌肉，打扫打扫卫生，感觉很舒服。在某一天，某LMZ 开始整理他那书架。已知他的书有n 本，从左到右按顺序排列。他想把书从矮到高排好序，而每一本书都有一个独一无二的高度Hi。他排序的方法是：每一次将所有的书划分为尽量少的连续部分，使得每一部分的书的高度都是单调下降，然后将其中所有不少于2 本书的区间全部翻转。重复执行以上操作，最后使得书的高度全部单调上升。可是毕竟是休息时间，LMZ 不想花太多时间在给书排序这种事上面。因此他划分并翻转完第一次书之后，他想计算，他一共执行了多少次翻转操作才能把所有的书排好序。LMZ 惊奇地发现，第一次排序之前，他第一次划分出来的所有区间的长度都是偶数。

Input

　　第一行一个正整数n, 为书的总数。

　　接下来一行n个数，第i个正整数Hi，为第i 本书的高度。

Output

　　仅一个整数，为LMZ 需要做的翻转操作的次数。

Sample Input

6

5 3 2 1 6 4

Sample Output

3

【样例解释】

第一次划分之后，翻转(5,3,2,1)，（6,4）。之后，书的高度为1 2 3 5 4 6，然后便是翻转（5,4）即可。

Data Constraint

　　对于10%的数据：n<=50

　　对于40%的数据：n<=3000

　　对于100%的数据：1<=n<=100000, 1<=Hi<=n

Summary

　　先模拟第一次翻转，之后只剩下区间连接处两个数字可能没有递增，变成了求逆序对的个数，将数列中依次放入线段树中，记录每个区间的个数，放入数字a[i]，逆序对的贡献就是已经放入且在区间(a[i]+1,H_max)里的数字个数。

 1 #include<cstdio>
 2 using namespace std;
 3 int n,a[1000000],x,b[1000000],c[1000000],h;
 4 long long ans,s[1000000];
 5 int bb(int w,int x,int y)
 6 {
 7     b[w]=x;
 8     c[w]=y;
 9     if (x==y) return 0;
10     bb(w*2,x,(x+y)/2);
11     bb(w*2+1,(x+y)/2+1,y);
12 }
13 int cc(int w,int x)
14 {
15     if (b[w]==c[w])
16     {
17         s[w]++;
18         return 0;
19     }
20     int mid=(b[w]+c[w])/2;
21     if (x>mid) cc(w*2+1,x);
22     if (x<=mid) cc(w*2,x);
23     s[w]=s[w*2]+s[w*2+1];
24 }
25 long long find(int w,int x,int y)
26 {
27     if (b[w]==x&&c[w]==y)
28         return s[w];
29     int mid=(b[w]+c[w])/2;
30     if (x>mid)
31         return find(w*2+1,x,y);
32     else
33     if (y<=mid)
34         return find(w*2,x,y);
35     else
36     return (find(w*2,x,mid)+find(w*2+1,mid+1,y));
37 }
38 int main()
39 {
40     scanf("%d",&n);
41     for (int i=1;i<=n;i++)
42     {
43         scanf("%d",&a[i]);
44         if (a[i]>h) h=a[i];
45     }
46     int r=0,l[10000],y[10000];
47     bool f=true;a[n+1]=0xffffff;
48     for (int i=2;i<=n+1;i++)
49         if (a[i]<a[i-1]) {
50             if (f) {
51                 x=i-1;f=false;
52             }
53         }
54         else
55         if (not f&&i-x>1) {
56             f=true;
57             l[++r]=x;
58             y[r]=i-1;
59         }
60     for (int j=1;j<=r;j++)
61         for (int i=l[j];i<=(l[j]+y[j])/2;i++) {
62             int c=y[j]-i+l[j];
63             int t=a[i];a[i]=a[c];a[c]=t;
64         }
65     ans=ans+r;
66     bb(1,1,h);
67     for (int i=1;i<=n;i++)
68     {
69         cc(1,a[i]);
70         if (a[i]+1<=n)
71             ans=ans+find(1,a[i]+1,h);
72     }
73     printf("%lld",ans);
74 }