bzoj1145

题目：http://www.lydsy.com/JudgeOnline/problem.php?id=1145

神题。。。。。。

定义f(abcd)为高度排名为abcd的个数，例如闪电的个数为f(1324)。

我们化简一下：

 f(1324)-f(1243)-f(1432)

=f(1x2x)-f(1423)-f(12xx)+f(1234)-f(14xx)+f(1423)

=f(1x2x)+f(1234)-f(12xx)-f(14xx)

=f(1x2x)+f(1234)-[f(1xxx)-f(13xx)]

=f(1x2x)+f(13xx)+f(1234)-f(1xxx)

实验变成求f(1x2x)，f(13xx)，f(1234)和f(1xxx)。

其中f(1234)和f(1xxx)都比较好求，我们来重点讨论f(1x2x)和f(13xx)。

首先预处理出两个比较重要的数组l[i]和r[i]。l[i]表示1到i-1小于a[i]的个数；r[i]表示i+1到N中小于a[i]的个数。

--------------------------------------------------------------------------------------------------------------------------------------------

f(1x2x)：

我们在2的位置进行统计，不妨记2的位置为i。为方便识别，我们在位置i那里涂成红色，即f(1xx)。

容易发现，如果我们知道了f(13)，那么f(1xx)=f(13)*(N-i-r[i])。

于是变成了求f(13)。

13

=(13+31+12+21)-(31+12+21)

观察第1个括号：132 312 123 21

我们发现他们左边两个黑色的数中，其中一个小于红色的数，另一个可以大于或小于红色的数。

所以第1个括号为：l[i]*(l[i]-1)/2+l[i]*(i-1-l[i])

观察第2个括号：312 123 21

我们发现他们第2个黑色的数小于红色的数，第1个黑色的数没有限制。

所以第2个括号为：∑(j-1)(j<i且a[j]<a[i])

这个可以用树状数组实现。

--------------------------------------------------------------------------------------------------------------------------------------------

f(13xx)：

我们在3的位置进行统计，记3的位置为i，把位置i那里涂成红色，即f(1xx)。

同样容易发现，如果我们知道了f(12)，那么f(1xx)=f(12)*(N-i-r[i])。

注意这里的f(12)和上面的f(13)的不同之处在于统计的位置（红色）。

于是变成求f(12)。

12

=(12+12+21)-(12+21)

观察第1个括号：12 12 21

我们发现3的右边一定有一个2，然后1的位置随便放。

所以第1个括号为：∑(a[j]-1)(i<j且a[i]>a[j])

观察第2个括号：12 21

其实就是 xx

所以第2个括号为：r[i]*(r[i]-1)/2

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<fstream>
#include<algorithm>
#include<cstring>
#include<string>
#include<cmath>
#include<queue>
#include<stack>
#include<map>
#include<utility>
#include<set>
#include<bitset>
#include<vector>
#include<functional>
#include<deque>
#include<cctype>
#include<climits>
#include<complex>
//#include<bits/stdc++.h>适用于CF,UOJ，但不适用于poj

using namespace std;

typedef long long LL;
typedef double DB;
typedef pair<int,int> PII;
typedef complex<DB> CP;

#define mmst(a,v) memset(a,v,sizeof(a))
#define mmcy(a,b) memcpy(a,b,sizeof(a))
#define re(i,a,b)  for(i=a;i<=b;i++)
#define red(i,a,b) for(i=a;i>=b;i--)
#define fi first
#define se second
#define m_p(a,b) make_pair(a,b)
#define SF scanf
#define PF printf
#define two(k) (1<<(k))

template<class T>inline T sqr(T x){return x*x;}
template<class T>inline void upmin(T &t,T tmp){if(t>tmp)t=tmp;}
template<class T>inline void upmax(T &t,T tmp){if(t<tmp)t=tmp;}

const DB EPS=1e-;
inline int sgn(DB x){if(abs(x)<EPS)return ;return(x>)?:-;}
const DB Pi=acos(-1.0);

inline int gint()
  {
        int res=;bool neg=;char z;
        for(z=getchar();z!=EOF && z!='-' && !isdigit(z);z=getchar());
        if(z==EOF)return ;
        if(z=='-'){neg=;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*+z-'',z=getchar());
        return (neg)?-res:res;
    }
inline LL gll()
  {
      LL res=;bool neg=;char z;
        for(z=getchar();z!=EOF && z!='-' && !isdigit(z);z=getchar());
        if(z==EOF)return ;
        if(z=='-'){neg=;z=getchar();}
        for(;z!=EOF && isdigit(z);res=res*+z-'',z=getchar());
        return (neg)?-res:res;
    }

const int maxN=;
const LL Mod=;

int N;
int a[maxN+];
LL l[maxN+],r[maxN+];
LL ans;

LL tree[maxN+];
#define lowbit(a) (a&(-a))
inline void update(int a,LL v){for(;a<=N;a+=lowbit(a))tree[a]=(tree[a]+v)%Mod;}
inline LL ask(int a){LL res=;for(;a>=;a-=lowbit(a))res=(res+tree[a])%Mod;return res;}

LL g[maxN+];

/*
 f(1324)-f(1243)-f(1432)
=f(1x2x)-f(1423)-f(12xx)+f(1234)-f(14xx)+f(1423)
=f(1x2x)+f(1234)-f(12xx)-f(14xx)
=f(1x2x)+f(1234)-[f(1xxx)-f(13xx)]
=f(1x2x)+f(13xx)+f(1234)-f(1xxx)
*/

int main()
  {
      freopen("bzoj1145.in","r",stdin);
      freopen("bzoj1145.out","w",stdout);
      int i;
      N=gint();
      re(i,,N)a[i]=gint();
      mmst(tree,);re(i,,N)l[i]=ask(a[i]-),update(a[i],);
      mmst(tree,);red(i,N,)r[i]=ask(a[i]-),update(a[i],);
      ans=;

      //+f(1x2x)
      re(i,,N)g[i]=(LL(l[i])*LL(l[i]-)/+LL(l[i])*LL(i--l[i]))%Mod;
      mmst(tree,);
      re(i,,N)
        {
            g[i]=(g[i]-ask(a[i]-))%Mod;
            update(a[i],i-);
        }
      re(i,,N)ans=(ans+g[i]*LL(N-i-r[i]))%Mod;

      //+f(13xx)
      mmst(tree,);
      red(i,N,)
        {
            LL t=LL(N-i-r[i]);
            LL p=(ask(a[i]-)-r[i]*(r[i]-)/)%Mod;
            ans=(ans+p*t)%Mod;
            update(a[i],a[i]-);
        }

      //+f(1234)
      re(i,,N)g[i]=l[i];
      mmst(tree,);re(i,,N)update(a[i],g[i]),g[i]=ask(a[i]-);
        mmst(tree,);re(i,,N)update(a[i],g[i]),g[i]=ask(a[i]-);
        re(i,,N)ans=(ans+g[i])%Mod;

        //-f(1xxx)
        re(i,,N)
          {
              LL t=LL(N-i-r[i]);
              ans=(ans-t*(t-)*(t-)/)%Mod;
            }

        ans=(ans%Mod+Mod)%Mod;
        cout<<ans<<endl;
        return ;
    }